Generalized Gibbs Ensemble Research Articles

Effective metal-insulator nonequilibrium quantum phase transition in the Su-Schrieffer-Heeger model

We consider the steady state behavior of observables in the Su-Schrieffer-Heeger model and in the one dimensional transverse field quantum Ising model after a sudden quantum quench of the parameter controlling the gap. In the thermodynamic limit, and for times $t\rightarrow\infty$, we find non-analyticities even in simple local observables as a function of the quench parameter, that is, a non-equilibrium quantum phase transition. We trace the appearance of this non-equilibrium quantum phase transition to an effective metal-insulator transition which occurs on the level of the generalized Gibbs Ensemble (describing the steady-state of the equilibrated system). Studying whether these transitions are robust, we find, in the paradigmatic case of the SSH model, that they persist for both quantum quench protocols of finite duration in time as well as thermal initial states, while they are washed out in the presence of fermion-fermion interactions and for finite system size.

Exact large-scale correlations in integrable systems out of equilibrium

Using the theory of generalized hydrodynamics (GHD), we derive exact Euler-scale dynamical two-point correlation functions of conserved densities and currents in inhomogeneous, non-stationary states of many-body integrable systems with weak space-time variations. This extends previous works to inhomogeneous and non-stationary situations. Using GHD projection operators, we further derive formulae for Euler-scale two-point functions of arbitrary local fields, purely from the data of their homogeneous one-point functions. These are new also in homogeneous generalized Gibbs ensembles. The technique is based on combining a fluctuation-dissipation principle along with the exact solution by characteristics of GHD, and gives a recursive procedure able to generate nn-point correlation functions. Owing to the universality of GHD, the results are expected to apply to quantum and classical integrable field theory such as the sinh-Gordon model and the Lieb-Liniger model, spin chains such as the XXZ and Hubbard models, and solvable classical gases such as the hard rod gas and soliton gases. In particular, we find Leclair-Mussardo-type infinite form-factor series in integrable quantum field theory, and exact Euler-scale two-point functions of exponential fields in the sinh-Gordon model and of powers of the density field in the Lieb-Liniger model. We also analyse correlations in the partitioning protocol, extract large-time asymptotics, and, in free models, derive all Euler-scale nn-point functions.

Quasilocal charges and the generalized Gibbs ensemble in the Lieb-Liniger model

We consider the construction of a generalized Gibbs ensemble composed of complete bases of conserved charges in the repulsive Lieb-Liniger model. We will show that it is possible to construct these bases with varying locality as well as demonstrating that such constructions are always possible provided one has in hand at least one complete basis set of charges. This procedure enables the construction of bases of charges that possess well defined, finite expectation values given an arbitrary initial state. We demonstrate the use of these charges in the context of two different quantum quenches: a quench where the strength of the interactions in a one-dimensional gas is switched suddenly from zero to some finite value and the release of a one dimensional cold atomic gas from a confining parabolic trap. While we focus on the Lieb-Liniger model in this paper, the principle of the construction of these charges applies to all integrable models, both in continuum and lattice form.

From the sinh-Gordon field theory to the one-dimensional Bose gas: exact local correlations and full counting statistics

We derive exact formulas for the expectation value of local observables in a one-dimensional gas of bosons with point-wise repulsive interactions (Lieb–Liniger model). Starting from a recently conjectured expression for the expectation value of vertex operators in the sinh-Gordon field theory, we derive explicit analytic expressions for the one-point K-body correlation functions in the Lieb–Liniger gas, for arbitrary integer K. These are valid for all excited states in the thermodynamic limit, including thermal states, generalized Gibbs ensembles and non-equilibrium steady states arising in transport settings. Our formulas display several physically interesting applications: most prominently, they allow us to compute the full counting statistics for the particle-number fluctuations in a short interval. Furthermore, combining our findings with the recently introduced generalized hydrodynamics, we are able to study multi-point correlation functions at the Eulerian scale in non-homogeneous settings. Our results complement previous studies in the literature and provide a full solution to the problem of computing one-point functions in the Lieb–Liniger model.

Squeezed thermal reservoir as a generalized equilibrium reservoir

We explore the perspective of considering the squeezed thermal reservoir as an equilibrium reservoir in a generalized Gibbs ensemble with two noncommuting conserved quantities. We outline the main properties of such a reservoir in terms of the exchange of energy, both heat and work, and entropy, giving some key examples to clarify its physical interpretation. This allows for a correct and insightful interpretation of all thermodynamical features of the squeezed thermal reservoir, as well as other similar nonthermal reservoirs, including the characterization of reversibility and the first and second laws of thermodynamics.

Quantum quench and thermalization of one-dimensional Fermi gas via phase-space hydrodynamics

By exploring a phase space hydrodynamics description of one-dimensional free Fermi gas, we discuss how systems settle down to steady states described by the generalized Gibbs ensembles through quantum quenches. We investigate time evolutions of the Fermions which are trapped in external potentials or a circle for a variety of initial conditions and quench protocols. We analytically compute local observables such as particle density and show that they always exhibit power law relaxation at late times. We find a simple rule which determines the power law exponent. Our findings are, in principle, observable in experiments in an one dimensional free Fermi gas or Tonk's gas (Bose gas with infinite repulsion).

Machine learning algorithms based on generalized Gibbs ensembles

Machine learning algorithms often take inspiration from the established results and knowledge from statistical physics. A prototypical example is the Boltzmann machine algorithm for supervised learning, which utilizes knowledge of classical thermal partition functions and the Boltzmann distribution. Recently, a quantum version of the Boltzmann machine was introduced by Amin et al, however, non-commutativity of quantum operators renders the training process by minimizing a cost function inefficient. Recent advances in the study of non-equilibrium quantum integrable systems, which never thermalize, have lead to the exploration of a wider class of statistical ensembles. These systems may be described by the so-called generalized Gibbs ensemble (GGE), which incorporates a number of ‘effective temperatures’. We propose that these GGEs can be successfully applied as the basis of a Boltzmann-machine–like learning algorithm, which operates by learning the optimal values of effective temperatures. We show that the GGE algorithm is an optimal quantum Boltzmann machine: it is the only quantum machine that circumvents the quantum training-process problem. We apply a simplified version of the GGE algorithm, where quantum effects are suppressed, to the classification of handwritten digits in the MNIST database. While lower error rates can be found with other state-of-the-art algorithms, we find that our algorithm reaches relatively low error rates while learning a much smaller number of parameters than would be needed in a traditional Boltzmann machine, thereby reducing computational cost.

Quantum dynamics with stochastic reset

We study non-equilibrium dynamics of integrable and non-integrable closed quantum systems whose unitary evolution is interrupted with stochastic resets, characterized by a reset rate $r$, that project the system to its initial state. We show that the steady state density matrix of a non-integrable system, averaged over the reset distribution, retains its off-diagonal elements for any finite $r$. Consequently a generic observable $\hat O$, whose expectation value receives contribution from these off-diagonal elements, never thermalizes under such dynamics for any finite $r$. We demonstrate this phenomenon by exact numerical studies of experimentally realizable models of ultracold bosonic atoms in a tilted optical lattice. For integrable Dirac-like fermionic models driven periodically between such resets, the reset-averaged steady state is found to be described by a family of generalized Gibbs ensembles (GGE s) characterized by $r$. We also study the spread of particle density of a non-interacting one-dimensional fermionic chain, starting from an initial state where all fermions occupy the left half of the sample, while the right half is empty. When driven by resetting dynamics, the density profile approaches at long times to a nonequilibrium stationary profile that we compute exactly. We suggest concrete experiments that can possibly test our theory.

Fluctuation relations for heat exchange in the generalized Gibbs ensemble

In this work, we investigate the heat exchange between two quantum systems whose initial equilibrium states are described by the generalized Gibbs ensemble. First, we generalize the fluctuation relations for heat exchange discovered by Jarzynski and W\'ojcik to quantum systems prepared in the equilibrium states described by the generalized Gibbs ensemble at different generalized temperatures. Second, we extend the connections between heat exchange and R\'enyi divergences to quantum systems with very general initial conditions.These relations are applicable for quantum systems with conserved quantities and are universally valid for quantum systems in the integrable and chaotic regimes.

Universal spectral correlations in the chaotic wave function and the development of quantum chaos

We investigate the appearance of quantum chaos in a single many-body wave function by analyzing the statistical properties of the eigenvalues of its reduced density matrix $\rho_A$ of a spatial subsystem A. We find that the spectrum of $\rho_A$ is described by a so-called Wishart random matrix, which exhibits universal spectral correlations between eigenvalues separated by distances ranging from one up to many mean level spacings. We use these universal spectral characteristics of $\rho_A$ as a definition of chaos in the wave function. A simple and precise characterization of such correlations is a segment of linear growth at sufficiently long times, recently called the "ramp", of the spectral form factor. Specifically, numerical results for the spectral form factor of $\rho_A$ of generic non-integrable systems, such as one-dimensional quantum Ising and Floquet spin models, are found to exhibit an universal ramp identical to that appearing for a random pure state, in which $\rho_A$ is the Wishart random matrix. In addition, we study the development of chaos in the wave function by letting an initial direct product state evolve under the unitary time evolution. We find that the universal spectral correlations as manifested by the ramp set in as soon as the entanglement entropy begins to grow, and first develop for the eigenvalues at the top of the spectrum of $\rho_A$, subsequently spreading over the entire spectrum at later times. Finally, we study a prethermalized regime described by a generalized Gibbs ensemble, which develops in a rapidly driven Floquet model at intermediate times. We find that the prethermalized regime exhibits no chaos, as evidenced by the absence of a ramp in the spectral form factor of $\rho_A$, while the universal spectral correlations start to develop when the prethermalized regime finally relaxes at late times to the fully thermalized chaotic regime.

Generalized hydrodynamics of classical integrable field theory: the sinh-Gordon model

Using generalized hydrodynamics (GHD), we develop the Euler hydrodynamics of classical integrable field theory. Classical field GHD is based on a known formalism for Gibbs ensembles of classical fields, that resembles the thermodynamic Bethe ansatz of quantum models, which we extend to generalized Gibbs ensembles (GGEs). In general, GHD must take into account both solitonic and radiative modes of classical fields. We observe that the quasi-particle formulation of GHD remains valid for radiative modes, even though these do not display particle-like properties in their precise dynamics. We point out that because of a UV catastrophe similar to that of black body radiation, radiative modes suffer from divergences that restrict the set of finite-average observables; this set is larger for GGEs with higher conserved charges. We concentrate on the sinh-Gordon model, which only has radiative modes, and study transport in the domain-wall initial problem as well as Euler-scale correlations in GGEs. We confirm a variety of exact GHD predictions, including those coming from hydrodynamic projection theory, by comparing with Metropolis numerical evaluations.

Quenched dynamics of classical isolated systems: the spherical spin model with two-body random interactions or the Neumann integrable model

We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N − 1 integrals of motion, notably, their scaling with N, and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs–Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.

Heating in Integrable Time-Periodic Systems.

We investigate a heating phenomenon in periodically driven integrable systems that can be mapped to free-fermion models. We find that heating to the high-temperature state, which is a typical scenario in nonintegrable systems, can also appear in integrable time-periodic systems; the amount of energy absorption rises drastically near a frequency threshold where the Floquet-Magnus expansion diverges. As the driving period increases, we also observe that the effective temperatures of the generalized Gibbs ensemble for conserved quantities go to infinity. By the use of the scaling analysis, we reveal that, in the limit of infinite system size and driving period, the steady state after a long time is equivalent to the infinite-temperature state. We obtain the asymptotic behavior L^{-1} and T^{-2} as to how the steady state approaches the infinite-temperature state as the system size L and the driving period T increase.

Relations between dissipated work and Rényi divergences in the generalized Gibbs ensemble

In this work, we show that the dissipation in a many-body system under an arbitrary non-equilibrium process is related to the R\'{e}nyi divergences between two states along the forward and reversed dynamics under very general family of initial conditions. This relation generalizes the links between dissipated work and Renyi divergences to quantum systems with conserved quantities whose equilibrium state is described by the generalized Gibbs ensemble. The relation is applicable for quantum systems with conserved quantities and can be applied to protocols driving the system between integrable and chaotic regimes. We demonstrate our ideas by considering the one-dimensional transverse quantum Ising model which is driven out of equilibrium by the instantaneous switching of the transverse magnetic field.

Time-dependent generalized Gibbs ensembles in open quantum systems

Generalized Gibbs ensembles have been used as powerful tools to describe the steady state of integrable many-particle quantum systems after a sudden change of the Hamiltonian. Here we demonstrate numerically, that they can be used for a much broader class of problems. We consider integrable systems in the presence of weak perturbations which both break integrability and drive the system to a state far from equilibrium. Under these conditions, we show that the steady state and the time-evolution on long time-scales can be accurately described by a (truncated) generalized Gibbs ensemble with time-dependent Lagrange parameters, determined from simple rate equations. We compare the numerically exact time evolutions of density matrices for small systems with a theory based on block-diagonal density matrices (diagonal ensemble) and a time-dependent generalized Gibbs ensemble containing only small number of approximately conserved quantities, using the one-dimensional Heisenberg model with perturbations described by Lindblad operators as an example.

Thermalization dynamics of two correlated bosonic quantum wires after a split

Coherently splitting a one-dimensional Bose gas provides an attractive, experimentally estab- lished platform to investigate many-body quantum dynamics. At short enough times, the dynamics is dominated by the dephasing of single quasi-particles, and well described by the relaxation to- wards a generalized Gibbs ensemble corresponding to the free Luttinger theory. At later times on the other hand, the approach to a thermal Gibbs ensemble is expected for a generic, interacting quantum system. Here, we go one step beyond the quadratic Luttinger theory and include the lead- ing phonon-phonon interactions. By applying kinetic theory and non-equilibrium Dyson-Schwinger equations, we analyze the full relaxation dynamics beyond dephasing and determine the asymptotic thermalization process in the two-wire system for a symmetric splitting protocol. The major ob- servables are the different phonon occupation functions and the experimentally accessible coherence factor, as well as the phase correlations between the two wires. We demonstrate that, depending on the splitting protocol, the presence of phonon collisions can have significant influence on the asymptotic evolution of these observables, which makes the corresponding thermalization dynamics experimentally accessible.

Quantum quench in a harmonically trapped one-dimensional Bose gas

We study the non-equilibrium dynamics of a one-dimensional Bose gas trapped by a harmonic potential for a quench from zero to infinite interaction. The different thermodynamic limits required for the equilibrium pre- and post-quench Hamiltonians are the origin of a few unexpected phenomena that have no counterparts in the translational invariant setting. We find that the dynamics is perfectly periodic with breathing time related to the strength of the trapping potential. For very short times, we observe a sudden expansion leading to an extreme dilution of the gas and to the emergence of slowly decaying tails in the density profile. The haste of the expansion induces a undertow effect with a pronounced local minimum of the density at the center of the trap. At half period there is a refocusing phenomenon characterized by a sharp central peak of the density, juxtaposed to algebraically decaying tails. We finally show that the time-averaged density is correctly captured by a generalized Gibbs ensemble built with the conserved mode occupations.

Double light-cone dynamics establish thermal states in integrable 1D Bose gases

We theoretically investigate the non-equilibrium dynamics in a quenched pair of one-dimensional Bose gases with density imbalance. We describe the system using its low-energy effective theory, the Luttinger liquid model. In this framework the system shows strictly integrable relaxation dynamics via dephasing of its approximate many-body eigenstates. In the balanced case, this leads to the well-known light-cone-like establishment of a prethermalized state, which can be described by a generalized Gibbs ensemble. In the imbalanced case the integrable dephasing leads to a state that, counter-intuitively, closely resembles a thermal equilibrium state. The approach to this state is characterized by two separate light-cone dynamics with distinct characteristic velocities. This behavior is a result of the fact that in the imbalanced case observables are not aligned with the conserved quantities of the integrable system. We discuss a concrete experimental realization to study this effect using matterwave interferometry and many-body revivals on an atom chip.

Conservation laws shape dissipation

Starting from the most general formulation of stochastic thermodynamics—i.e. a thermodynamically consistent nonautonomous stochastic dynamics describing systems in contact with several reservoirs—we define a procedure to identify the conservative and the minimal set of nonconservative contributions in the entropy production. The former is expressed as the difference between changes caused by time-dependent drivings and a generalized potential difference. The latter is a sum over the minimal set of flux-force contributions controlling the dissipative flows across the system. When the system is initially prepared at equilibrium (e.g. by turning off drivings and forces), a finite-time detailed fluctuation theorem holds for the different contributions. Our approach relies on identifying the complete set of conserved quantities and can be viewed as the extension of the theory of generalized Gibbs ensembles to nonequilibrium situations.

Generalized thermalization for integrable system under quantum quench.

We investigate equilibration and generalized thermalization of the quantum Harmonic chain under local quantum quench. The quench action we consider is connecting two disjoint harmonic chains of different sizes and the system jumps between two integrable settings. We verify the validity of the generalized Gibbs ensemble description for this infinite-dimensional Hilbert space system and also identify equilibration between the subsystems as in classical systems. Using Bogoliubov transformations, we show that the eigenstates of the system prior to the quench evolve toward the Gibbs Generalized Ensemble description. Eigenstates that are more delocalized (in the sense of inverse participation ratio) prior to the quench, tend to equilibrate more rapidly. Further, through the phase space properties of a generalized Gibbs ensemble and the strength of stimulated emission, we identify the necessary criterion on the initial states for such relaxation at late times and also find out the states that would potentially not be described by the generalized Gibbs ensemble description.