Generalized Gibbs Ensemble (GGE)

• Generalized Gibbs Ensemble (GGE) is a framework in statistical mechanics that extends the conventional Gibbs ensemble by incorporating an infinite set of conserved quantities to describe post-quench stationary states.
• It leverages the fusion hierarchy of transfer matrices and the thermodynamic Bethe–Takahashi equations to accurately determine root densities in integrable lattice models.
• The formalism employs explicit string–charge relations and quasi-local charges to provide a complete, non-perturbative characterization of stationary states in SU(N)-invariant systems.

The Generalized Gibbs Ensemble (GGE) is a fundamental framework in the statistical mechanics of non-equilibrium integrable quantum systems. It extends the conventional Gibbs (thermal) ensemble by incorporating an infinite set of conserved quantities—typically both local and quasi-local integrals of motion—thereby allowing an accurate and complete description of the stationary state arising after quantum quenches or other nonequilibrium processes in integrable models. In lattice systems with higher internal symmetries, such as those exhibiting $SU(N)$ invariance and solvable by the nested Bethe Ansatz, the construction and operational characterization of the GGE requires a sophisticated algebraic machinery exploiting the fusion hierarchy of transfer matrices, their string-charge relations, and the thermodynamic Bethe–Takahashi equations (Fehér et al., 2019). This article presents a detailed account of the structure, formulation, and significance of the Generalized Gibbs Ensemble in nested Bethe Ansatz solvable models, emphasizing the role of quasi-local charges, string-charge relations, and the functional completeness of the GGE formalism.

1. Formal Definition and Variational Principle

In integrable quantum chains, the GGE is defined by the density matrix

$$\rho_{\mathrm{GGE}} = \frac{1}{Z} \exp\left( -\sum_{a=1}^\infty \lambda_a Q_a \right), \qquad Z = \mathrm{Tr} \exp\left(-\sum_{a=1}^\infty \lambda_a Q_a\right)$$

where $\{Q_a\}_{a=1}^\infty$ denotes a (typically infinite) family of mutually commuting conserved operators, with $Q_1 \equiv H$ being the Hamiltonian and the remainder generated by logarithmic derivatives of the fundamental transfer matrix. Each Lagrange multiplier $\lambda_a$ enforces the matching of the GGE expectation value of $Q_a$ to its initial value, $$\langle Q_a \rangle_{\textrm{initial}} = \mathrm{Tr} \left[ \rho_{\mathrm{GGE}} Q_a \right]$$ (Fehér et al., 2019). The GGE is thus the unique maximum entropy (Jaynes) state subject to the full set of constraints imposed by local and quasi-local conservation laws.

2. Complete Set of Charges via Fusion Hierarchy

A structurally complete set of conserved charges for $SU(N)$-invariant integrable lattice spin chains ($N \geq 3$) is obtained from the fusion hierarchy of transfer matrices. Specifically, the transfer matrices $t^{(a)}_m(u)$, defined for rows $a=1,\ldots,N-1$ and columns $m=1,2,\ldots$, satisfy the $SU(N)$ Hirota (T-system) functional relations: $$t^{(a)}_m(u+\tfrac{i}{2})\, t^{(a)}_m(u-\tfrac{i}{2}) = t^{(a)}_{m+1}(u)\, t^{(a)}_{m-1}(u) + t^{(a-1)}_m(u)\, t^{(a+1)}_m(u)$$ with appropriate boundary conditions. One defines the generating operators as

$X^{(a)}_m(u) = -i\, \partial_u \ln\, t^{(a)}_m(u)$

which, for $|{\rm Im}\,u| < 1/2$ (within the physical strip), are quasi-local: they admit decompositions into densities with Hilbert–Schmidt norm growing linearly with system length $L$, and local overlaps remaining finite as $L \to \infty$. The quasi-locality of the $X$-operators stems from local inversion relations satisfied by the fused Lax operators. This construction ensures that all relevant root distributions characterizing thermodynamic macrostates are fully fixed by the infinite set of (quasi-)local charges so constructed (Fehér et al., 2019).

3. Thermodynamic Bethe–Takahashi Equations and Root Densities

Bethe eigenstates in $SU(N)$ spin chains are parametrized by $N-1$ families ("levels") of Bethe rapidities $\{\lambda_j^{(a)}\}$, each corresponding to a nesting level. In the thermodynamic limit, these become continuous string densities $\rho_n^{(a)}(\lambda)$, indexed by level $a$ and string length $n$, with hole densities $\rho_{h,n}^{(a)}(\lambda)$ and total densities $\rho_{t,n}^{(a)} = \rho_n^{(a)} + \rho_{h,n}^{(a)}$.

The partially decoupled Thermodynamic Bethe–Takahashi equations are: $$\rho_{t,n}^{(a)}(\lambda) = \delta_{n,1}\, \delta_{a,1}\, s(\lambda) + s \ast \left[ \rho_{h,n-1}^{(a)} + \rho_{h,n+1}^{(a)} \right](\lambda) + s \ast \left[ \rho_n^{(a-1)} + \rho_n^{(a+1)} \right](\lambda)$$ with $s(\lambda) = 1/2\cosh(\pi \lambda)$, where the convolution $s\ast f (\lambda)$ is defined as usual. The hierarchy is completed by the boundary conventions $\rho_{h, 0}^{(a)} = 0$. In Fourier space, $s(\lambda)$ has transform $\hat{s}(k) = 1/(2\cosh(k/2))$. These equations and their associated root densities provide the complete thermodynamic characterization of stationary states.

4. String–Charge Relations and Information-Completeness

A central result is the explicit relation between the $X$-operators and the string distributions. On a generic Bethe eigenstate, one finds

$$X^{(a)}_m(u) \approx \sum_{j}\left[ \frac{1}{u - \lambda_j^{(a)} - mi/2} - \frac{1}{u - \lambda_j^{(a)} + mi/2} \right]$$

which can be further represented as

$$X^{(a)}_m(u) \approx 2\pi L \sum_{n=1}^{\infty} \int d\lambda\, \rho_n^{(a)}(\lambda) \sum_{\ell=1}^{\min(m,n)} a_{|n-m|-1 + 2\ell}(u-\lambda)$$

with $a_p(\lambda) = p/[2\pi(\lambda^2 + (p/2)^2)]$.

Inversion of these expressions, typically by convolution and Fourier techniques, yields the string–charge relations: $$\rho_m^{(a)}(\lambda) = \frac{ X_m^{(a)}(\lambda + i/2) + X_m^{(a)}(\lambda - i/2) - X_{m-1}^{(a)}(\lambda) - X_{m+1}^{(a)}(\lambda) }{2\pi L}$$ Alternatively, the hole densities are expressed as

$$\rho_{h,m}^{(a)}(\lambda) = a_m(\lambda) - \frac{ X_m^{(a)}(\lambda + i/2)+ X_m^{(a)}(\lambda - i/2) - X_m^{(a+1)}(\lambda) - X_m^{(a-1)}(\lambda) }{2\pi L}$$

These relations, together with the TBA, provide an invertible, complete mapping between the infinite set of quasi-local charges and the full set of thermodynamic macrostates.

5. Characterization of the GGE Steady State and Quench Dynamics

The information-complete description offered by the GGE in the fusion hierarchy framework ensures that for any global quench in an integrable $SU(N)$ spin chain, the late-time stationary state is given by the unique root distributions consistent with the initial values of all $X^{(a)}_m(u)$. The GGE is operationalized by enforcing

$$\langle X^{(a)}_m(u) \rangle_{\mathrm{GGE}} = \langle X^{(a)}_m(u) \rangle_{\mathrm{initial}}$$

over the physical strip $|{\rm Im}\,u| < 1/2$, yielding the corresponding set of root densities $\{\rho_m^{(a)}(\lambda)\}$ via the string–charge relations. All local observables and correlation functions, as well as the entire dynamics of relaxation following the quench, are thereby determined (Fehér et al., 2019).

6. General Properties and Broader Significance

The GGE constructed from the fusion hierarchy:

• Involves both local and quasi-local charges, the latter being crucial in Bethe Ansatz models where strictly local charges alone are insufficient for completeness.
• Supplies a non-perturbative and fully explicit procedure for characterizing stationary states in $SU(N)$ spin chains and, by extension, other nested Bethe Ansatz systems.
• Justifies, at the algebraic and spectral level, the accuracy of GGE predictions for late-time observables following unitary dynamics from arbitrary initial conditions.
• Resolves the ambiguities in selecting the relevant set of conservation laws for the GGE, by relating them directly to the transfer-matrix fusion hierarchy and its functional relations (Fehér et al., 2019).

7. Schematic Overview of GGE Construction in $SU(N)$-Invariant Chains

Step	Object/Concept	Formal Structure/Role
1. Symmetry	$SU(N)$, Nested Bethe Ansatz	Bethe rapidities per nesting level
2. Charges	Transfer-matrix fusion hierarchy $t_m^{(a)}$	Quasi-local charges: $X^{(a)}_m(u)$
3. Thermodynamics	TBA equations for $\rho_n^{(a)}(\lambda)$	Complete thermodynamic macrostate
4. String–Charge	Invertible relations $X^{(a)}_m \rightarrow \rho_n^{(a)}$	Information completeness of the GGE
5. GGE State	$\rho_{\mathrm{GGE}}$ with all $\lambda_a$	Stationary state with prescribed initial $\langle Q_a\rangle$

This hierarchical methodology encapsulates a fully explicit, non-perturbative description of stationary (post-quench) states in $SU(N)$ integrable lattice models, providing a robust and general framework for the generalized thermalization and equilibration of such systems (Fehér et al., 2019).