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:1711.00672 · Physical review. E

Relation between far-from-equilibrium dynamics and equilibrium correlation functions for binary operators

Jonas Richter, Robin Steinigeweg

Published Jan 11, 2019• 13 citations
arXivHTMLPDF

Abstract

Linear response theory (LRT) is one of the main approaches to the dynamics of quantum many-body systems. However, this approach has limitations and requires, e.g., that the initial state is (i) mixed and (ii) close to equilibrium. In this paper, we discuss these limitations and study the nonequilibrium dynamics for a certain class of properly prepared initial states. Specifically, we consider thermal states of the quantum system in the presence of an additional static force which, however, become nonequilibrium states when this static force is eventually removed. While for weak forces the relaxation dynamics is well captured by LRT, much less is known in the case of strong forces, i.e., initial states far away from equilibrium. Summarizing our main results, we unveil that, for high temperatures, the nonequilibrium dynamics of so-called binary operators is always generated by an equilibrium correlation function. In particular, this statement holds true for states in the far-from-equilibrium limit, i.e., outside the linear response regime. In addition, we confirm our analytical results by numerically studying the dynamics of local fermionic occupation numbers and local energy densities in the spin-1/2 Heisenberg chain. Remarkably, these simulations also provide evidence that our results qualitatively apply in a more general setting, e.g., in the anisotropic XXZ model where the local energy is a non-binary operator, as well as for a wider range of temperature. Furthermore, exploiting the concept of quantum typicality, all of our findings are not restricted to mixed states, but are valid for pure initial states as well.

DOI
10.1103/physreve.99.012114

Related papers

Dynamics of Pure States and Linear Response for Binary Operators

Jonas Richter, Robin Steinigeweg

We study static and dynamical expectation values for a class of pure states. This class mimics non-equilibrium density matrices via the concept of quantum typicality and allows for a controlled preparation of expectation values close to or far away from their equilibrium values at arbitrary temperatures. While these pure states do not require specific types of observables, we investigate so-called binary operators such as local fermionic occupation numbers or local energy densities in the isotropic Heisenberg spin-1/2 chain. For these two examples, we use numerical forward-propagation methods in real and imaginary time to demonstrate the accuracy of the pure-state approach within and outside the linear-response regime. At high temperatures, our numerical results show that the dynamics in all regimes is generated by a single correlation function. These numerical results are in excellent agreement with our analytical predictions for binary operators.

•Nov 2, 2017•1 citation

Relaxation of dynamically prepared out-of-equilibrium initial states within and beyond linear response theory

Jonas Richter, Mats H. Lamann, Christian Bartsch +2 more

We consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.

arXiv:1905.03292•Physical review. E•Sep 16, 2019•6 citations

Impact of eigenstate thermalization on the route to equilibrium

Jonas Richter, Jochen Gemmer B, Robin Steinigeweg

The eigenstate thermalization hypothesis (ETH) and the theory of linear response (LRT) are celebrated cornerstones of our understanding of the physics of many-body quantum systems out of equilibrium. While the ETH provides a generic mechanism of thermalization for states arbitrarily far from equilibrium, LRT extends the successful concepts of statistical mechanics to situations close to equilibrium. In our work, we connect these cornerstones to shed light on the route to equilibrium for a class of properly prepared states. We unveil that, if the off-diagonal part of the ETH applies, then the relaxation process can become independent of whether or not a state is close to equilibrium. Moreover, in this case, the dynamics is generated by a single correlation function, i.e., the relaxation function in the context of LRT. Our analytical arguments are illustrated by numerical results for idealized models of random-matrix type and more realistic models of interacting spins on a lattice. Remarkably, our arguments also apply to integrable quantum systems where the diagonal part of the ETH may break down.

arXiv:1805.11625•Physical review. E•May 24, 2019•31 citations

Spin- 1/2 XXZ chain coupled to two Lindblad baths: Constructing nonequilibrium steady states from equilibrium correlation functions

Tjark Heitmann, Jonas Richter, Fengping Jin +7 more

State-of-the-art approaches to extract transport coefficients of many-body quantum systems broadly fall into two categories: (i) they target the linear-response regime in terms of equilibrium correlation functions of the closed system; or (ii) they consider an open-system situation typically modeled by a Lindblad equation, where a nonequilibrium steady state emerges from driving the system at its boundaries. While quantitative agreement between (i) and (ii) has been found for selected model and parameter choices, also disagreement has been pointed out in the literature. Studying magnetization transport in the spin-$1/2$ XXZ chain, we here demonstrate that at weak driving, the nonequilibrium steady state in an open system, including its buildup in time, can remarkably be constructed just on the basis of correlation functions in the closed system. We numerically illustrate this direct correspondence of closed-system and open-system dynamics, and show that it allows the treatment of comparatively large open systems, usually only accessible to matrix product state simulations. We also point out potential pitfalls when extracting transport coefficients from nonequilibrium steady states in finite systems.

arXiv:2303.00430•Physical review. B./Physical review. B•Nov 27, 2023•6 citations

Quantum transport of strongly interacting fermions in one dimension far out of equilibrium

Jie Zou, Xiaopeng Li

In the study of quantum transport, much has been known about dynamics near thermal equilibrium. However, quantum transport far away from equilibrium is much less well understood---the linear response approximation does not hold for physics far out of equilibrium in general. In this work, motivated by recent cold atom experiments on probing quantum many-body dynamics of a one-dimensional XXZ spin chain, where a transition from ballistic to diffusive dynamics has been established by increasing the interaction strengths, we study the strong interaction limit of the one-dimensional spinless fermion model, which is dual to the XXZ spin chain. We develop a highly efficient computation algorithm for simulating the nonequilibrium dynamics of this system exactly, and examine the nonequilibrium dynamics starting from a density modulation quantum state. We find ballistic transport in this strongly correlated setting and show that a plane-wave description emerges at long-time evolution. We also observe a sharp distinction between transport velocities in short and long times as induced by interaction effects and provide a quantitative interpretation for the long-time transport velocity. We expect our results to shed light on the understanding of the dynamics of the XXZ spin chain in the strong interaction regime.

arXiv:2206.13023•Physical Review A•Nov 28, 2022•1 citation

Sudden removal of a static force in a disordered system: Induced dynamics, thermalization, and transport

Jonas Richter, J. Herbrych, Robin Steinigeweg

We study the real-time dynamics of local occupation numbers in a one-dimensional model of spinless fermions with a random on-site potential for a certain class of initial states. The latter are thermal (mixed or pure) states of the model in the presence of an additional static force but become nonequilibrium states after a sudden removal of this static force. For this class and high temperatures, we show that the induced dynamics is given by a single correlation function at equilibrium, independent of the initial expectation values being prepared close to equilibrium (by a weak static force) or far away from equilibrium (by a strong static force). Remarkably, this type of universality holds true in both the ergodic phase and the many-body localized regime. Moreover, it does not depend on the specific choice of a unit cell for the local density. We particularly discuss two important consequences. First, the long-time expectation value of the local density is uniquely determined by the fluctuations of its diagonal matrix elements in the energy eigenbasis. Thus, the validity of the eigenstate thermalization hypothesis is not only a sufficient but also a necessary condition for thermalization. Second, the real-time broadening of density profiles is always given by the current autocorrelation function at equilibrium via a generalized Einstein relation. In the context of transport, we discuss the influence of disorder for large particle-particle interactions, where normal diffusion is known to occur in the disorder-free case. Our results suggest that normal diffusion is stable against weak disorder, while they are consistent with anomalous diffusion for stronger disorder below the localization transition. Particularly, for weak disorder, Gaussian density profiles can be observed for single disorder realizations, which we demonstrate for finite lattices up to 31 sites.

arXiv:1808.00497•Physical review. B./Physical review. B•Oct 8, 2018•11 citations
Linear response theory (LRT) is one of the main approaches to the dynamics of quantum many-body systems. However, this approach has limitations and requires, e.g., that the initial state is (i) mixed and (ii) close to equilibrium. In this paper, we discuss these limitations and study the nonequilibrium dynamics for a certain class of properly prepared initial states. Specifically, we consider thermal states of the quantum system in the presence of an additional static force which, however, become nonequilibrium states when this static force is eventually removed. While for weak forces the relaxation dynamics is well captured by LRT, much less is known in the case of strong forces, i.e., initial states far away from equilibrium. Summarizing our main results, we unveil that, for high temperatures, the nonequilibrium dynamics of so-called binary operators is always generated by an equilibrium correlation function. In particular, this statement holds true for states in the far-from-equilibrium limit, i.e., outside the linear response regime. In addition, we confirm our analytical results by numerically studying the dynamics of local fermionic occupation numbers and local energy densities in the spin-1/2 Heisenberg chain. Remarkably, these simulations also provide evidence that our results qualitatively apply in a more general setting, e.g., in the anisotropic XXZ model where the local energy is a non-binary operator, as well as for a wider range of temperature. Furthermore, exploiting the concept of quantum typicality, all of our findings are not restricted to mixed states, but are valid for pure initial states as well.
We study static and dynamical expectation values for a class of pure states. This class mimics non-equilibrium density matrices via the concept of quantum typicality and allows for a controlled preparation of expectation values close to or far away from their equilibrium values at arbitrary temperatures. While these pure states do not require specific types of observables, we investigate so-called binary operators such as local fermionic occupation numbers or local energy densities in the isotropic Heisenberg spin-1/2 chain. For these two examples, we use numerical forward-propagation methods in real and imaginary time to demonstrate the accuracy of the pure-state approach within and outside the linear-response regime. At high temperatures, our numerical results show that the dynamics in all regimes is generated by a single correlation function. These numerical results are in excellent agreement with our analytical predictions for binary operators.
We consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.
The eigenstate thermalization hypothesis (ETH) and the theory of linear response (LRT) are celebrated cornerstones of our understanding of the physics of many-body quantum systems out of equilibrium. While the ETH provides a generic mechanism of thermalization for states arbitrarily far from equilibrium, LRT extends the successful concepts of statistical mechanics to situations close to equilibrium. In our work, we connect these cornerstones to shed light on the route to equilibrium for a class of properly prepared states. We unveil that, if the off-diagonal part of the ETH applies, then the relaxation process can become independent of whether or not a state is close to equilibrium. Moreover, in this case, the dynamics is generated by a single correlation function, i.e., the relaxation function in the context of LRT. Our analytical arguments are illustrated by numerical results for idealized models of random-matrix type and more realistic models of interacting spins on a lattice. Remarkably, our arguments also apply to integrable quantum systems where the diagonal part of the ETH may break down.
State-of-the-art approaches to extract transport coefficients of many-body quantum systems broadly fall into two categories: (i) they target the linear-response regime in terms of equilibrium correlation functions of the closed system; or (ii) they consider an open-system situation typically modeled by a Lindblad equation, where a nonequilibrium steady state emerges from driving the system at its boundaries. While quantitative agreement between (i) and (ii) has been found for selected model and parameter choices, also disagreement has been pointed out in the literature. Studying magnetization transport in the spin-1/21/21/2 XXZ chain, we here demonstrate that at weak driving, the nonequilibrium steady state in an open system, including its buildup in time, can remarkably be constructed just on the basis of correlation functions in the closed system. We numerically illustrate this direct correspondence of closed-system and open-system dynamics, and show that it allows the treatment of comparatively large open systems, usually only accessible to matrix product state simulations. We also point out potential pitfalls when extracting transport coefficients from nonequilibrium steady states in finite systems.
In the study of quantum transport, much has been known about dynamics near thermal equilibrium. However, quantum transport far away from equilibrium is much less well understood---the linear response approximation does not hold for physics far out of equilibrium in general. In this work, motivated by recent cold atom experiments on probing quantum many-body dynamics of a one-dimensional XXZ spin chain, where a transition from ballistic to diffusive dynamics has been established by increasing the interaction strengths, we study the strong interaction limit of the one-dimensional spinless fermion model, which is dual to the XXZ spin chain. We develop a highly efficient computation algorithm for simulating the nonequilibrium dynamics of this system exactly, and examine the nonequilibrium dynamics starting from a density modulation quantum state. We find ballistic transport in this strongly correlated setting and show that a plane-wave description emerges at long-time evolution. We also observe a sharp distinction between transport velocities in short and long times as induced by interaction effects and provide a quantitative interpretation for the long-time transport velocity. We expect our results to shed light on the understanding of the dynamics of the XXZ spin chain in the strong interaction regime.
We study the real-time dynamics of local occupation numbers in a one-dimensional model of spinless fermions with a random on-site potential for a certain class of initial states. The latter are thermal (mixed or pure) states of the model in the presence of an additional static force but become nonequilibrium states after a sudden removal of this static force. For this class and high temperatures, we show that the induced dynamics is given by a single correlation function at equilibrium, independent of the initial expectation values being prepared close to equilibrium (by a weak static force) or far away from equilibrium (by a strong static force). Remarkably, this type of universality holds true in both the ergodic phase and the many-body localized regime. Moreover, it does not depend on the specific choice of a unit cell for the local density. We particularly discuss two important consequences. First, the long-time expectation value of the local density is uniquely determined by the fluctuations of its diagonal matrix elements in the energy eigenbasis. Thus, the validity of the eigenstate thermalization hypothesis is not only a sufficient but also a necessary condition for thermalization. Second, the real-time broadening of density profiles is always given by the current autocorrelation function at equilibrium via a generalized Einstein relation. In the context of transport, we discuss the influence of disorder for large particle-particle interactions, where normal diffusion is known to occur in the disorder-free case. Our results suggest that normal diffusion is stable against weak disorder, while they are consistent with anomalous diffusion for stronger disorder below the localization transition. Particularly, for weak disorder, Gaussian density profiles can be observed for single disorder realizations, which we demonstrate for finite lattices up to 31 sites.