Thermodynamic Bethe Ansatz Equations Research Articles

Nonlinear integral equations for the sausage model**Dedicated to the memory of Petr Petrovich Kulish.

The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the sigma model preserving integrability. The target space is deformed from the sphere to ‘sausage’ shape by a deformation parameter ν. This model is defined by a factorizable S-matrix which is obtained by deforming that of the sigma model by a parameter λ. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonlinear integral equations (NLIEs), which are applicable to generic value of λ. Our derivation is based on relations extracted from the truncated TBA equations. For a consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between ν and λ.

On the mass-coupling relation of multi-scale quantum integrable models

We determine exactly the mass-coupling relation for the simplest multi-scale quantum integrable model, the homogenous sine-Gordon model with two independent mass-scales. We first reformulate its perturbed coset CFT description in terms of the perturbation of a projected product of minimal models. This representation enables us to identify conserved tensor currents on the UV side. These UV operators are then mapped via form factor perturbation theory to operators on the IR side, which are characterized by their form factors. The relation between the UV and IR operators is given in terms of the sought-for mass-coupling relation. By generalizing the $\Theta$ sum rule Ward identity we are able to derive differential equations for the mass-coupling relation, which we solve in terms of hypergeometric functions. We check these results against the data obtained by numerically solving the thermodynamic Bethe Ansatz equations, and find a complete agreement.

Yang–Yang equilibrium statistical mechanics: A brilliant method

Yang and Yang in 1969 [J. Math. Phys. 10, 1115 (1969)] for the first time proposed a rigorous approach to the thermodynamics of the one-dimensional system of bosons with a delta-function interaction. This paper was a breakthrough in exact statistical mechanics, after Yang [Phys. Rev. Lett. 19, 1312 (1967)] published his seminal work on the discovery of the Yang–Baxter equation in 1967. Yang and Yang’s brilliant method yields significant applications in a wide range of fields of physics. In this paper, we briefly introduce the method of the Yang–Yang equilibrium statistical mechanics and demonstrate a fundamental application of the Yang–Yang method for the study of thermodynamics of the Lieb–Liniger model with strong and weak interactions in a whole temperature regime. We also consider the equivalence between the Yang–Yang’s thermodynamic Bethe ansatz equation and the thermodynamics of the ideal gas with the Haldane’s generalized exclusion statistics.

Thermodynamics of Inozemtsev's elliptic spin chain

We study the thermodynamic behaviour of Inozemtsev's long-range elliptic spin chain using the Bethe ansatz equations describing the spectrum of the model in the infinite-length limit. We classify all solutions of these equations in that limit and argue which of these solutions determine the spectrum in the thermodynamic limit. Interestingly, some of the solutions are not selfconjugate, which puts the model in sharp contrast to one of the model's limiting cases, the Heisenberg xxx spin chain. Invoking the string hypothesis we derive the thermodynamic Bethe ansatz equations (TBA-equations) from which we determine the Helmholtz free energy in thermodynamic equilibrium and derive the associated Y-system. We corroborate our results by comparing numerical solutions of the TBA-equations to a direct computation of the free energy for the finite-length hamiltonian. In addition we confirm numerically the interesting conjecture put forward by Finkel and González-López that the original and supersymmetric versions of Inozemtsev's elliptic spin chain are equivalent in the thermodynamic limit.

A Riemann-Hilbert formulation for the finite temperature Hubbard model

Inspired by recent results in the context of AdS/CFT integrability, we reconsider the Thermodynamic Bethe Ansatz equations describing the 1D fermionic Hubbard model at finite temperature. We prove that the infinite set of TBA equations are equivalent to a simple nonlinear Riemann-Hilbert problem for a finite number of unknown functions. The latter can be transformed into a set of three coupled nonlinear integral equations defined over a finite support, which can be easily solved numerically. We discuss the emergence of an exact Bethe Ansatz and the link between the TBA approach and the results by Juttner, Klumper and Suzuki based on the Quantum Transfer Matrix method. We also comment on the analytic continuation mechanism leading to excited states and on the mirror equations describing the finite-size Hubbard model with twisted boundary conditions.

Crystal fields and Kondo effect: New results for the magnetic susceptibility

The thermodynamic Bethe ansatz equations for the Coqblin–Schrieffer model have been solved for the first time to obtain the magnetic susceptibility in the presence of crystal fields for non-zero temperatures. For the case of N=4 effective ionic states an analytic expression for the limiting values of the pseudo-energies has been found facilitating the numerical solution for various crystal and magnetic field configurations. The single-impurity model applies to a wide range of dense Kondo systems and has been used before to explain apparent non-Fermi-liquid behavior. The flattening off of the susceptibility curves at a substantially higher temperature than the specific heat is shown to be a general feature of the Coqblin–Schrieffer thermodynamics.

Nontrivial thermodynamics in ’t Hooft’s large-Nlimit

We study the finite volume/temperature correlation functions of the (1+1)-dimensional ${\rm SU}(N)$ principal chiral sigma model in the planar limit. The exact S-matrix of the sigma model is known to simplify drastically at large $N$, and this leads to trivial thermodynamic Bethe ansatz (TBA) equations. The partition function, if derived using the TBA, can be shown to be that of free particles. We show that the correlation functions and expectation values of operators at finite volume/temperature are not those of the free theory, and that the TBA does not give enough information to calculate them. Our analysis is done using the Leclair-Mussardo formula for finite-volume correlators, and knowledge of the exact infinite-volume form factors. We present analytical results for the one-point function of the energy-momentum tensor, and the two-point function of the renormalized field operator. The results for the energy-momentum tensor can be used to define a nontrivial partition function.

Elementary functions in thermodynamic Bethe ansatz**Based on talks given at ‘Infinite Analysis 2014’ (Tokyo, 2014) and at ‘Integrable lattice models and quantum field theories’ (Bad Honnef, 2014).

Some years ago, Fendley found an explicit solution to the thermodynamic Bethe ansatz (TBA) equation for an supersymmetric theory in 2D with a specific F-term. Motivated by this, we seek explicit solutions for other super-potential cases utilizing the idea from the ODE/IM correspondence. We find that the TBA equations, corresponding to a wider class of super-potentials, admit solutions in terms of elementary functions such as modified Bessel functions and confluent hyper-geometric series.

Finite-volume spectra of the Lee-Yang model

We consider the non-unitary Lee-Yang minimal model ${\cal M}(2,5)$ in three different finite geometries: (i) on the interval with integrable boundary conditions labelled by the Kac labels $(r,s)=(1,1),(1,2)$, (ii) on the circle with periodic boundary conditions and (iii) on the periodic circle including an integrable purely transmitting defect. We apply $\varphi_{1,3}$ integrable perturbations on the boundary and on the defect and describe the flow of the spectrum. Adding a $\Phi_{1,3}$ integrable perturbation to move off-criticality in the bulk, we determine the finite size spectrum of the massive scattering theory in the three geometries via Thermodynamic Bethe Ansatz (TBA) equations. We derive these integral equations for all excitations by solving, in the continuum scaling limit, the TBA functional equations satisfied by the transfer matrices of the associated $A_{4}$ RSOS lattice model of Forrester and Baxter in Regime III. The excitations are classified in terms of $(m,n)$ systems. The excited state TBA equations agree with the previously conjectured equations in the boundary and periodic cases. In the defect case, new TBA equations confirm previously conjectured transmission factors.

Thermodynamic properties of one-dimensional Gaudin-Yang model at finite temperature

The one-dimensional system interacting via a delta-function interparticle interaction is a very important one in cold atomic systems and has fundamental importance in many-body physics. In one dimension, due to the geometric confinement induced quantum correlations and quantum fluctuations, there may exist a number of unusual phenomena, such as spin-charge separation, effective fermionization and quantum criticality. This paper studies the basic properties of a uniform one-dimensional Gaudin-Yang model for fermions by solving the thermodynamic Bethe-ansatz equations by a numerical method. Numerically, we use the many-variable Newton’s method to solve the coupled equations. We analyze the physical properties, including density, interaction, temperature and entropy at a given temperature and a given interaction, separately. We know that a lot of researches are limited to zero temperature. However, we cannot reach the absolute zero temperature in the real cold atomic experiment. So it is important to deal with the finite temperature problems. We study the density and entropy as a function of the chemical potential, temperature and interaction and, then give the phase diagrams, respectively. We found that there is a quantum critical zone in the phase diagram of entropy, including the high temperature zone with thermal fluctuations and the Luttinger liquid zone with quantum fluctuations. For a given temperature and low chemical potential, the thermal fluctuations are the main factor in the entropy. With the increase of chemical potential, the system enters the quantum critical zone where the competitive effect between the thermal fluctuations and the quantum fluctuations exists. When the chemical potential is large enough, the quantum fluctuations become the main factor in the system’s entropy, and we get the Luttinger liquid phase. Our results can be further used in the finite temperature density-functional theory and to analyze the collective phenomena at a finite temperature.

Quantum criticality of one-dimensional multicomponent Fermi gas with strongly attractive interaction

Quantum criticality of strongly attractive Fermi gas with symmetry in one dimension is studied via the thermodynamic Bethe ansatz (TBA) equations. The phase transitions driven by the chemical potential , effective magnetic field H1, H2 (chemical potential biases) are analyzed at the quantum criticality. The phase diagram and critical fields are analytically determined by the TBA equations in the zero temperature limit. High accurate equations of state, scaling functions are also obtained analytically for the strong interacting gases. The dynamic exponent and correlation length exponent read off the universal scaling form. It turns out that the quantum criticality of the three-component gases involves a sudden change of density of states of one cluster state, two or three cluster states. In general, this method can be adapted to deal with the quantum criticality of multicomponent Fermi gases with symmetry.

Winding vacuum energies in a deformed [formula omitted] sigma model

We consider the problem of calculating the Casimir energies in the winding sectors of Fateev's SS-model, which is an integrable two-parameter deformation of the O(4) non-linear sigma model in two dimensions. This problem lies beyond the scope of all traditional methods of integrable quantum field theory including the thermodynamic Bethe ansatz and non-linear integral equations. Here we propose a solution based on a remarkable correspondence between classical and quantum integrable systems and express the winding energies in terms of certain solutions of the classical sinh-Gordon equation.

Crystal fields and Kondo effect: Specific heat for Cerium compounds

The thermodynamic Bethe ansatz equations for the N=6 Coqblin-Schrieffer model with crystal fields have been solved numerically. The realistic case of three Kramers doublets with arbitrary splittings has been studied for the first time. The specific heat has been calculated for representative combinations of the ionic energy splittings providing ample material for comparison with experimental results for Cerium impurities and compounds.

Quantum Wronskian approach to six-point gluon scattering amplitudes at strong coupling

We study the six-point gluon scattering amplitudes in N=4 super Yang-Mills theory at strong coupling based on the twisted Z_4-symmetric integrable model. The lattice regularization allows us to derive the associated thermodynamic Bethe ansatz (TBA) equations as well as the functional relations among the Q-/T-/Y-functions. The quantum Wronskian relation for the Q-/T-functions plays an important role in determining a series of the expansion coefficients of the T-/Y-functions around the UV limit, including the dependence on the twist parameter. Studying the CFT limit of the TBA equations, we derive the leading analytic expansion of the remainder function for the general kinematics around the limit where the dual Wilson loops become regular-polygonal. We also compare the rescaled remainder functions at strong coupling with those at two, three and four loops, and find that they are close to each other along the trajectories parameterized by the scale parameter of the integrable model.

Mayer-cluster expansion of instanton partition functions and thermodynamic bethe ansatz

In arXiv:0908.4052, Nekrasov and Shatashvili pointed out that the N=2 instanton partition function in a special limit of the Omega-deformation parameters is characterized by certain thermodynamic Bethe ansatz (TBA) like equations. In this work we present an explicit derivation of this fact as well as generalizations to quiver gauge theories. To do so we combine various techniques like the iterated Mayer expansion, the method of expansion by regions, and the path integral tricks for non-perturbative summation. The TBA equations derived entirely within gauge theory have been proposed to encode the spectrum of a large class of quantum integrable systems. We hope that the derivation presented in this paper elucidates further this completely new point of view on the origin, as well as on the structure, of TBA equations in integrable models.

Reformulating the TBA equations for the quark anti-quark potential and their two loop expansion

The boundary thermodynamic Bethe Ansatz (BTBA) equations introduced in arXiv:1203.1913 and arXiv:1203.1617 to describe the cusp anomalous dimension contain imaginary chemical potentials and singular boundary fugacities, which make its systematic expansion problematic. We propose an alternative formulation based on real chemical potentials and additional source terms. We expand our equations to double wrapping order and find complete agreement with the direct two-loop gauge theory computation of the cusp anomalous dimension.

Thermodynamic Bethe ansatz for non-equilibrium steady states: exact energy current and fluctuations in integrable QFT

We evaluate the exact energy current and scaled cumulant generating function (related to the large-deviation function) in non-equilibrium steady states with energy flow, in any integrable model of relativistic quantum field theory (IQFT) with diagonal scattering. Our derivations are based on various recent results of Bernard and Doyon. The steady states are built by connecting homogeneously two infinite halves of the system thermalized at different temperatures Tl, Tr, and waiting for a long time. We evaluate the current J(Tl, Tr) using the exact QFT density matrix describing these non-equilibrium steady states and using Zamolodchikov’s method of the thermodynamic Bethe ansatz (TBA). The scaled cumulant generating function is obtained from the extended fluctuation relations which hold in integrable models. We verify our formula in particular by showing that the conformal field theory (CFT) result is obtained in the high-temperature limit. We analyze numerically our non-equilibrium steady-state TBA equations for three models: the sinh-Gordon model, the roaming trajectories model, and the sine-Gordon model at a particular reflectionless point. Based on the numerics, we conjecture that an infinite family of non-equilibrium c-functions, associated with the scaled cumulants, can be defined, which we interpret physically. We study the full scaled distribution function and find that it can be described by a set of independent Poisson processes. Finally, we show that the ‘additivity’ property of the current, which is known to hold in CFT and was proposed to hold more generally, does not hold in general IQFT—that is, J(Tl, Tr) is not of the form f(Tl) − f(Tr).

Exact results for the low energy AdS 4 × $ \mathbb{C}\mathbb{P} $ 3 string theory

We derive the Thermodynamic Bethe Ansatz equations for the relativistic sigma model describing the AdS(4)XCP(3) string II A theory at strong coupling (i.e. in the Alday-Maldacena decoupling limit). The corresponding Y-system involves an infinite number of Y functions and is of a new type, although it shares a peculiar feature with the Y-system for AdS(4)XCP(3). A truncation of the equations at level p and a further generalisation to generic rank N allow us an alternative description of the theory as the N=4, p= \infty representative in an infinite family of models corresponding to the conformal cosets CP(N-1)_p X U(1), perturbed by a relevant composite field \phi(N,p) =\phi_[CP(N-1)_p] X \phi[U(1)] that couples the two independent conformal field theories. The calculation of the ultraviolet central charge confirms the conjecture by Basso and Rej and the conformal dimension of the perturbing operator, at every N and p, is obtained using the Y-system periodicity. The conformal dimension of \phi[CP(N-1)_p] matches that of the field identified by Fendley while discussing integrability issues for the purely bosonic CP(N-1) sigma model.

Discontinuity relations for the AdS4/CFT3 correspondence

We study in detail the analytic properties of the Thermodynamic Bethe Ansatz (TBA) equations for the anomalous dimensions of composite operators in the planar limit of the 3D N=6 superconformal Chern–Simons gauge theory and derive functional relations for the jump discontinuities across the branch cuts in the complex rapidity plane. These relations encode the analytic structure of the Y functions and are extremely similar to the ones obtained for the previously-studied AdS5/CFT4 case. Together with the Y-system and more basic analyticity conditions, they are completely equivalent to the TBA equations. We expect these results to be useful to derive alternative nonlinear integral equations for the AdS4/CFT3 spectrum.

Quantisation of the effective string with TBA

In presence of a static pair of sources, the spectrum of low-lying states of whatever confining gauge theory in D space-time dimensions is described, at large source separations, by an effective string theory. In the far infrared the latter flows, in the static gauge, to a two-dimensional massless free-field theory. It is known that the Lorentz invariance of the gauge theory fixes uniquely the first few subleading corrections of this free-field limit. We point out that the first allowed correction - a quartic polynomial in the field derivatives - is exactly the composite field $ T\overline{T} $ , built with the chiral components,T and $ \overline{T} $ , of the energy-momentum tensor. This irrelevant perturbation is quantum integrable and yields, through the thermodynamic Bethe Ansatz (TBA), the energy levels of the string which exactly coincide with the Nambu-Goto spectrum. We obtain this way the results recently found by Dubovsky, Flauger and Gorbenko. This procedure easily generalizes to any two-dimensional CFT. It is known that the leading deviation of the Nambu-Goto spectrum comes from the boundary terms of the string action. We solve the TBA equations on an infinite strip, identify the relevant boundary parameter and verify that it modifies the string spectrum as expected.