Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1801.07646 · Physical review. B./Physical review. B

Time-dependent generalized Gibbs ensembles in open quantum systems

Florian Lange, Zala Lenarčič, Achim Rosch

Published Apr 26, 2018• 82 citations
arXivHTMLPDF

Abstract

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 break both 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 timescales 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 a small number of approximately conserved quantities, using the one-dimensional Heisenberg model with perturbations described by Lindblad operators as an example.

DOI
10.1103/physrevb.97.165138

Related papers

On generalized Gibbs ensembles with an infinite set of conserved charges

Balázs Pozsgay, Eric Vernier, Miklós Antal Werner

Abstract We revisit the question of whether and how the steady states arising after non-equilibrium time evolution in integrable models (and in particular in the XXZ spin chain) can be described by the so-called generalized Gibbs ensemble (GGE). Whereas it is known that the micro-canonical ensemble built on a complete set of charges correctly describes the long-time limit of local observables, it has been shown recently by Ilievski et al that the corresponding canonical ensemble is not well defined, and instead a different canonical ensemble was proposed in terms of particle occupation number operators. Here we provide an alternative construction by considering truncated GGEs (tGGEs) that include only a finite number of local and quasi-local conserved operators. It is shown that the tGGEs can approximate the steady states with arbitrary precision, i.e. all physical observables are exactly reproduced in the infinite truncation limit. We trace back the problems encountered in defining an untruncated GGE to the dependence of the associated Lagrange multipliers on the truncation index. Conversely, we show that this problem may be circumvented by considering a new set of (quasi)local charges which are linear combinations of the standard ones, and whose associated Lagrange multipliers are well-defined state functions. Our general arguments are applied to concrete quench situations in the XXZ chain, where the initial states are simple two-site or four-site product states. Depending on the quench we find that numerical results for the local correlators can be obtained with remarkable precision using truncated GGEs with only 10–100 charges.

arXiv:1703.09516•Journal of Statistical Mechanics Theory and Experiment•Sep 1, 2017•77 citations

Exact steady states for quantum quenches in integrable Heisenberg spin chains

Lorenzo Piroli, Eric Vernier, Pasquale Calabrese

What is the long-time behavior of a quantum system after it is brought out of equilibrium? In a generic case, it has been proposed that an effective thermalization occurs, and the local properties of the system can then be described by a thermal Gibbs ensemble. However, in integrable systems, with infinitely many conservation laws, the answer is far more complicated. It was shown in several special cases that the long-time steady state can be exactly described by the so called complete generalized Gibbs ensemble. The latter can be obtained by maximizing the entropy with constraints imposed by the conservation of energy, but also by additional carefully chosen subset of the integrals of motion. In this work, the authors analyze more general initial states, focusing on integrable XXZ Heisenberg spin chains. They find that the complete generalized Gibbs ensemble indeed correctly describes the system at long times after the quench.

arXiv:1606.00383•Physical review. B./Physical review. B•Aug 23, 2016•100 citations

Iterative Construction of Conserved Quantities in Dissipative Nearly Integrable Systems

Iris Ulčakar, Zala Lenarčič

Integrable systems offer rare examples of solvable many-body problems in the quantum world. Because of the fine-tuned structure, their realization in nature and experiment is never completely accurate, and therefore effects of integrability are observed only transiently. One way to overcome this limitation is to weakly couple nearly integrable systems to baths and driving: this will stabilize integrable effects up to arbitrary time and encode them in the stationary state approximated by a generalized Gibbs ensemble. However, the description of such driven dissipative nearly integrable models is challenging and no exact analytical methods have been proposed so far. Here, we develop an iterative scheme in which integrability breaking perturbations (baths) determine the conserved quantities that play the leading role in a highly efficient truncated generalized Gibbs ensemble description. Our scheme paves the way for easier calculations in thermodynamically large systems and can be used to construct unknown conserved quantities.

arXiv:2310.03809•Physical Review Letters•Jun 4, 2024•9 citations

Time evolution and thermodynamics for a nonequilibrium system in phase-space

Chen-Huan Wu

The integrable system is constrained strictly by the conservation law during the time evolution, and the prethermal state from the nearly integrable system is also constrained by the conserved parameters (the constants of motion) with the corresponding generalized Gibbs ensemble (GGE), which is indubitably a powerful tool in the prediction of the relaxation dynamics. For stochastic evolution dynamics with considerable noise, the two-point correlation of local operators (like the density of kinks or transverse magnetization correlators), which do not exhibit the thermal features, display the behaviors of nonthermalization and an asymptotic GGE. In fact it is an asymptotic quasi-steady state with an infinite temperature, therefore the required distance to the nonthermal steady state is in an infinite time average. In this paper, we unambiguously investigate the relaxation of a nonequilibrium system in a canonical ensemble for integrable and nonintegrable systems. Temporal behavior of the many-body quantum system and the corresponding linear-coupling between the harmonic oscillators are discussed. The matrix-method in entropy ensemble is utilized to discuss the boundary and the diagonalization algebraically. The approximation results for nonintegrable system under the considerable perturbations are also presented.

Canadian Journal of Physics•Sep 28, 2018•6 citations

Eigenstate Gibbs ensemble in integrable quantum systems

Sourav Nandy, Arnab Sen, Arnab Das +1 more

The eigenstate thermalization hypothesis conjectures that for a thermodynamically large system in one of its energy eigenstates, the reduced density matrix describing any finite subsystem is determined solely by a set of relevant conserved quantities. In a chaotic quantum system, only the energy is expected to play that role and hence eigenstates appear locally thermal. Integrable systems, on the other hand, possess an extensive number of such conserved quantities and therefore the reduced density matrix requires specification of all the corresponding parameters (generalized Gibbs ensemble). However, here we show by unbiased statistical sampling of the individual eigenstates with a given finite energy density that the local description of an overwhelming majority of these states of even such an integrable system is actually Gibbs-like, i.e., requires only the energy density of the eigenstate. Rare eigenstates that cannot be represented by the Gibbs ensemble can also be sampled efficiently by our method and their local properties are then shown to be described by appropriately truncated generalized Gibbs ensembles. We further show that the presence of these rare eigenstates differentiates the model from the chaotic case and leads to the system being described by a generalized Gibbs ensemble at long time under a unitary dynamics following a sudden quench, even when the initial state is a typical (Gibbs-like) eigenstate of the prequench Hamiltonian.

arXiv:1605.09225•Physical review. B./Physical review. B•Dec 21, 2016•45 citations

Perturbative approach to weakly driven many-particle systems in the presence of approximate conservation laws

Zala Lenarčič, Florian Lange, Achim Rosch

We develop a Liouville perturbation theory for weakly driven and weakly open quantum systems in situations when the unperturbed system has a number of conservations laws. If the perturbation violates the conservation laws, it drives the system to a new steady state which can be approximately but efficiently described by a (generalized) Gibbs ensemble characterized by one Lagrange parameter for each conservation law. The value of those has to be determined from rate equations for conserved quantities. Remarkably, even weak perturbations can lead to large responses of conserved quantities. We present a perturbative expansion of the steady state density matrix; first we give the condition that fixes the zeroth-order expression (Lagrange parameters) and then determine the higher-order corrections via projections of the Liouvillian. The formalism can be applied to a wide range of problems including two-temperature models for electron-phonon systems, Bose condensates of excitons or photons, or weakly perturbed integrable models. We test our formalism by studying interacting fermions coupled to nonthermal reservoirs, approximately described by a Boltzmann equation.

arXiv:1706.05700•Physical review. B./Physical review. B•Jan 3, 2018•54 citations
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 break both 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 timescales 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 a small number of approximately conserved quantities, using the one-dimensional Heisenberg model with perturbations described by Lindblad operators as an example.
Abstract We revisit the question of whether and how the steady states arising after non-equilibrium time evolution in integrable models (and in particular in the XXZ spin chain) can be described by the so-called generalized Gibbs ensemble (GGE). Whereas it is known that the micro-canonical ensemble built on a complete set of charges correctly describes the long-time limit of local observables, it has been shown recently by Ilievski et al that the corresponding canonical ensemble is not well defined, and instead a different canonical ensemble was proposed in terms of particle occupation number operators. Here we provide an alternative construction by considering truncated GGEs (tGGEs) that include only a finite number of local and quasi-local conserved operators. It is shown that the tGGEs can approximate the steady states with arbitrary precision, i.e. all physical observables are exactly reproduced in the infinite truncation limit. We trace back the problems encountered in defining an untruncated GGE to the dependence of the associated Lagrange multipliers on the truncation index. Conversely, we show that this problem may be circumvented by considering a new set of (quasi)local charges which are linear combinations of the standard ones, and whose associated Lagrange multipliers are well-defined state functions. Our general arguments are applied to concrete quench situations in the XXZ chain, where the initial states are simple two-site or four-site product states. Depending on the quench we find that numerical results for the local correlators can be obtained with remarkable precision using truncated GGEs with only 10–100 charges.
What is the long-time behavior of a quantum system after it is brought out of equilibrium? In a generic case, it has been proposed that an effective thermalization occurs, and the local properties of the system can then be described by a thermal Gibbs ensemble. However, in integrable systems, with infinitely many conservation laws, the answer is far more complicated. It was shown in several special cases that the long-time steady state can be exactly described by the so called complete generalized Gibbs ensemble. The latter can be obtained by maximizing the entropy with constraints imposed by the conservation of energy, but also by additional carefully chosen subset of the integrals of motion. In this work, the authors analyze more general initial states, focusing on integrable XXZ Heisenberg spin chains. They find that the complete generalized Gibbs ensemble indeed correctly describes the system at long times after the quench.
Integrable systems offer rare examples of solvable many-body problems in the quantum world. Because of the fine-tuned structure, their realization in nature and experiment is never completely accurate, and therefore effects of integrability are observed only transiently. One way to overcome this limitation is to weakly couple nearly integrable systems to baths and driving: this will stabilize integrable effects up to arbitrary time and encode them in the stationary state approximated by a generalized Gibbs ensemble. However, the description of such driven dissipative nearly integrable models is challenging and no exact analytical methods have been proposed so far. Here, we develop an iterative scheme in which integrability breaking perturbations (baths) determine the conserved quantities that play the leading role in a highly efficient truncated generalized Gibbs ensemble description. Our scheme paves the way for easier calculations in thermodynamically large systems and can be used to construct unknown conserved quantities.
The integrable system is constrained strictly by the conservation law during the time evolution, and the prethermal state from the nearly integrable system is also constrained by the conserved parameters (the constants of motion) with the corresponding generalized Gibbs ensemble (GGE), which is indubitably a powerful tool in the prediction of the relaxation dynamics. For stochastic evolution dynamics with considerable noise, the two-point correlation of local operators (like the density of kinks or transverse magnetization correlators), which do not exhibit the thermal features, display the behaviors of nonthermalization and an asymptotic GGE. In fact it is an asymptotic quasi-steady state with an infinite temperature, therefore the required distance to the nonthermal steady state is in an infinite time average. In this paper, we unambiguously investigate the relaxation of a nonequilibrium system in a canonical ensemble for integrable and nonintegrable systems. Temporal behavior of the many-body quantum system and the corresponding linear-coupling between the harmonic oscillators are discussed. The matrix-method in entropy ensemble is utilized to discuss the boundary and the diagonalization algebraically. The approximation results for nonintegrable system under the considerable perturbations are also presented.
The eigenstate thermalization hypothesis conjectures that for a thermodynamically large system in one of its energy eigenstates, the reduced density matrix describing any finite subsystem is determined solely by a set of relevant conserved quantities. In a chaotic quantum system, only the energy is expected to play that role and hence eigenstates appear locally thermal. Integrable systems, on the other hand, possess an extensive number of such conserved quantities and therefore the reduced density matrix requires specification of all the corresponding parameters (generalized Gibbs ensemble). However, here we show by unbiased statistical sampling of the individual eigenstates with a given finite energy density that the local description of an overwhelming majority of these states of even such an integrable system is actually Gibbs-like, i.e., requires only the energy density of the eigenstate. Rare eigenstates that cannot be represented by the Gibbs ensemble can also be sampled efficiently by our method and their local properties are then shown to be described by appropriately truncated generalized Gibbs ensembles. We further show that the presence of these rare eigenstates differentiates the model from the chaotic case and leads to the system being described by a generalized Gibbs ensemble at long time under a unitary dynamics following a sudden quench, even when the initial state is a typical (Gibbs-like) eigenstate of the prequench Hamiltonian.
We develop a Liouville perturbation theory for weakly driven and weakly open quantum systems in situations when the unperturbed system has a number of conservations laws. If the perturbation violates the conservation laws, it drives the system to a new steady state which can be approximately but efficiently described by a (generalized) Gibbs ensemble characterized by one Lagrange parameter for each conservation law. The value of those has to be determined from rate equations for conserved quantities. Remarkably, even weak perturbations can lead to large responses of conserved quantities. We present a perturbative expansion of the steady state density matrix; first we give the condition that fixes the zeroth-order expression (Lagrange parameters) and then determine the higher-order corrections via projections of the Liouvillian. The formalism can be applied to a wide range of problems including two-temperature models for electron-phonon systems, Bose condensates of excitons or photons, or weakly perturbed integrable models. We test our formalism by studying interacting fermions coupled to nonthermal reservoirs, approximately described by a Boltzmann equation.