Simulating the Generalized Gibbs Ensemble (GGE): a Hilbert space Monte Carlo approach

Abstract

By combining classical Monte Carlo and Bethe ansatz techniques we devise a numerical method to construct the Truncated Generalized Gibbs Ensemble (TGGE) for the spin-1/2 isotropic Heisenberg ($XXX$) chain. The key idea is to sample the Hilbert space of the model with the appropriate GGE probability measure. The method can be extended to other integrable systems, such as the Lieb-Liniger model. We benchmark the approach focusing on GGE expectation values of several local observables. As finite-size effects decay exponentially with system size, moderately large chains are sufficient to extract thermodynamic quantities. The Monte Carlo results are in agreement with both the Thermodynamic Bethe Ansatz (TBA) and the Quantum Transfer Matrix approach (QTM). Remarkably, it is possible to extract in a simple way the steady-state Bethe-Gaudin-Takahashi (BGT) roots distributions, which encode complete information about the GGE expectation values in the thermodynamic limit. Finally, it is straightforward to simulate extensions of the GGE, in which, besides the local integral of motion (local charges), one includes arbitrary functions of the BGT roots. As an example, we include in the GGE the first non-trivial quasi-local integral of motion.

6 pages, 3 figures

Cited by in corpus (8)

• Quench dynamics and relaxation in isolated integrable quantum spin chains
• Entanglement and thermodynamics after a quantum quench in integrable systems
• Quantum quenches in 1+1 dimensional conformal field theories
• Quench action and Renyi entropies in integrable systems
• Entanglement spreading in non-equilibrium integrable systems
• Rényi entropies after releasing the Néel state in the XXZ spin-chain
• Quasi-local conserved charges and spin transport in spin-$1$ integrable chains
• The quench action approach in finite integrable spin chains