Thermodynamic Bethe Ansatz (TBA)

• Thermodynamic Bethe Ansatz is a framework based on non-linear integral equations that describe thermodynamic properties and finite-size spectra in quantum integrable systems.
• It leverages variational principles and model-dependent kernel functions from two-body scattering to relate pseudoenergy functions with continuous Bethe root densities.
• The TBA extends to generalized Gibbs ensembles, linking infinite conserved charges with quantum integrable models and establishing connections to quantum field theory and gauge/string dualities.

The Thermodynamic Bethe Ansatz (TBA) is a non-linear integral equation framework central to the computation of thermodynamic properties, exact finite-size spectra, and the study of generalized ensembles in quantum integrable models. Originally developed for describing equilibrium states in integrable quantum systems, the TBA framework has evolved into a unifying approach with deep connections to quantum field theory, statistical mechanics, quantum spectral theory, and gauge/string dualities.

1. Conceptual Foundation and General Structure

At its core, the TBA formalism replaces discrete Bethe quantization conditions for finite systems with integral equations for continuous root densities in the thermodynamic limit. In an integrable system, the spectrum is encoded in a set of Bethe rapidities, whose distribution determines physical properties. The TBA arises from a variational principle: one minimizes the generalized free energy functional

$$G[\rho, \rho^h] = \sum_n \beta_n Q_n[\rho] - S[\rho, \rho^h]$$

where $\rho(\lambda)$ and $\rho^h(\lambda)$ denote the densities of particles and holes at rapidity $\lambda$, $\{Q_n\}$ are commuting conserved charges each with chemical potential $\beta_n$, and the entropy $S[\rho, \rho^h]$ emerges from the combinatorial counting of Bethe states (Mossel et al., 2012, Tongeren, 2016).

Minimization under Bethe-Yang constraints yields a non-linear integral equation for the pseudoenergy (or "Y-function"):

$$\varepsilon(\lambda) = \varepsilon_0(\lambda) - \int K(\lambda-\mu) \ln[1 + e^{-\varepsilon(\mu)}]\, d\mu$$

where the driving term $$\varepsilon_0(\lambda) = \sum_n \beta_n q_n(\lambda)$$ encodes information about the conserved charges and external fields, while the kernel $K$ is model-dependent, determined by the two-body S-matrix.

This pseudoenergy is directly related to root densities and occupation numbers via

$$\vartheta(\lambda) = \frac{1}{1 + e^{\varepsilon(\lambda)}} = \frac{\rho(\lambda)}{\rho(\lambda) + \rho^h(\lambda)}$$

with observables and free energies given as explicit functionals thereof.

2. Extensions to Generalized Gibbs Ensembles

Integrable models generically possess an infinite set of local and quasi-local conserved charges. The natural statistical ensemble is therefore the Generalized Gibbs Ensemble (GGE) with density matrix

$$\hat\rho_{\rm GGE} = \mathcal Z_{\rm GGE}^{-1} \exp\left( -\sum_{n=0}^\infty \beta_n \hat Q_n \right)$$

where all $\beta_n$ are fixed by $\langle \hat Q_m\rangle$ constraints (Mossel et al., 2012). The TBA formalism generalizes by encoding all these constraints in the driving term, and the existence and uniqueness of solutions are ensured under mild growth conditions on $\varepsilon_0(\lambda)$.

In the context of the Lieb-Liniger model (delta-interacting 1D bosons), for example,

$$\varepsilon_0(\lambda) = \sum_{n=0}^\infty \beta_n \lambda^n$$

and the Bethe integral equation connects all physical root densities smoothly to the set of chemical potentials, establishing a bijective correspondence between nonequilibrium stationary states and the set $\{\beta_n\}$.

3. Functional and Algebraic Structures: The Y-System and Kernel Properties

TBA integral equations possess deep functional relations known as Y-systems, which take the form of discrete difference equations governing the analytic structure of the solutions. For many lattice and continuum models, especially those based on SU(N) or related symmetries, the Y-system is local:

$$Y_a(\theta + i\eta)\, Y_a(\theta - i\eta) = \prod_b [1 + Y_b(\theta)]^{I_{ab}}$$

where $I_{ab}$ is the incidence matrix of the associated Dynkin diagram and $\eta$ encodes the fusion spectral parameter (Tongeren, 2016, 0902.3930, 0903.0141). The analytic data (zero/pole positions, asymptotics) supplement these relations to reconstruct the full pseudoenergy functions.

The explicit kernels $K_{ab}(\lambda)$, encoding two-body scattering or fusion interaction, are determined by the logarithmic derivatives of the S-matrix. For strings/bound-state sectors in spin chains (e.g., XXX/XXZ, Inozemtsev elliptic), kernel hierarchies emerge and the Y-systems acquire an infinite-tower structure (Klabbers, 2016, Pavlis et al., 2019).

4. Representation-theoretic and Gauge-theoretic Connections

A broad class of quantum integrable systems admits a "ODE/IM correspondence," where TBA equations solve the Riemann-Hilbert problem for WKB periods of corresponding quantum differential equations (Imaizumi, 2020, Ito et al., 2019, Suzuki, 2015, Masoero, 2010). For instance, the TBA for the Mathieu equation precisely computes quantum periods, effective central charges