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:1409.6660 · Journal of Physics A Mathematical and Theoretical

Non-equilibrium steady states in the Klein–Gordon theory

Benjamin Doyon, Andrew Lucas, Koenraad Schalm, M. J. Bhaseen

Published Feb 10, 2015• 111 citations
arXivHTMLPDF

Abstract

We construct non-equilibrium steady states in the Klein–Gordon theory in arbitrary space dimension d following a local quench. We consider the approach where two independently thermalized semi-infinite systems, with temperatures and , are connected along a -dimensional hypersurface. A current-carrying steady state, described by thermally distributed modes with temperatures and for left and right-moving modes, respectively, emerges at late times. The non-equilibrium density matrix is the exponential of a non-local conserved charge. We obtain exact results for the average energy current and the complete distribution of energy current fluctuations. The latter shows that the long-time energy transfer can be described by a continuum of independent Poisson processes, for which we provide the exact weights. We further describe the full time evolution of local observables following the quench. Averages of generic local observables, including the stress-energy tensor, approach the steady state with a power-law in time, where the exponent depends on the initial conditions at the connection hypersurface. We describe boundary conditions and special operators for which the steady state is reached instantaneously on the connection hypersurface. A semiclassical analysis of freely propagating modes yields the average energy current at large distances and late times. We conclude by comparing and contrasting our findings with results for interacting theories and provide an estimate for the timescale governing the crossover to hydrodynamics. As a modification of our Klein–Gordon analysis we also include exact results for free Dirac fermions.

DOI
10.1088/1751-8113/48/9/095002

Related papers

Non-equilibrium steady states for the interacting Klein-Gordon field in 1+3 dimensions

Thomas-Paul Hack, Rainer Verch

Non-equilibrium steady states (NESS) describe particularly simple and stationary non-equilibrium situations. A possibility to obtain such states is to consider the asymptotic evolution of two infinite heat baths brought into thermal contact. In this work we generalise corresponding results of Doyon~et.~al. (J.\ Phys.\ A 18 (2015) no.9) for free Klein-Gordon fields in several directions. Our analysis is carried out directly at the level of correlation functions and in the algebraic approach to QFT. We discuss non-trivial chemical potentials, condensates, inhomogeneous linear models and homogeneous interacting ones. We shall not consider a sharp contact at initial time, but a smooth transition region. As a consequence, the states we construct will be of Hadamard type, and thus sufficiently regular for the perturbative treatment of interacting models. Our analysis shows that perturbatively constructed interacting NESS display thermodynamic properties similar to the ones of the NESS in linear models. In particular, perturbation theory appears to be insufficient to describe full thermalisation in non-linear QFT. Notwithstanding, we find that the NESS for linear and interacting models is stable under small perturbations, which is one of the characteristic features of equilibrium states.

arXiv:1806.00504•Jun 1, 2018•2 citations

Non-equilibrium steady states in quantum critical systems with Lifshitz scaling

Daniel Fernández, Aruna Rajagopal, Lárus Thorlacius A

Abstract We study out-of-equilibrium energy transport in a quantum critical fluid with Lifshitz scaling symmetry following a local quench between two semi-infinite fluid reservoirs. The late time energy flow is universal and is accommodated via a steady state occupying an expanding central region between outgoing shock and rarefaction waves. We consider the admissibility and entropy conditions for the formation of such a non-equilibrium steady state for a general dynamical critical exponent z in arbitrary dimensions and solve the associated Riemann problem. The Lifshitz fluid with z = 2 can be obtained from a Galilean boost invariant field theory and the non-equilibrium steady state is identified as a boosted thermal state. A Lifshitz fluid with generic z is scale invariant but without boost symmetry and in this case the non-equilibrium steady state is genuinely non-thermal.

arXiv:1909.06377•Journal of High Energy Physics•Dec 1, 2019•12 citations

Quantum quench from a thermal tensor state: Boundary effects and generalized Gibbs ensemble

Mario Collura, Dragi Karevski

We consider a quantum quench in a noninteracting fermionic one-dimensional field theory. The system of size $L$ is initially prepared into two halves $\mathcal{L}$ $([\ensuremath{-}L/2,0])$ and $\mathcal{R}$ $([0,L/2])$, each of them thermalized at two different temperatures ${T}_{\mathcal{L}}$ and ${T}_{\mathcal{R}}$, respectively. At a given time, the two halves are joined together by a local coupling and the whole system is left to evolve unitarily. For an infinitely extended system $(L\ensuremath{\rightarrow}\ensuremath{\infty})$, we show that the time evolution of the particle and energy densities is well described via a hydrodynamic approach which allows us to evaluate the correspondent stationary currents. We show, in such a case, that the two-point correlation functions are deduced, at large times, from a simple nonequilibrium steady state. Otherwise, whenever the boundary conditions are retained (in a properly defined thermodynamic limit), any current is suppressed at large times, and the stationary state is described by a generalized Gibbs ensemble, which is diagonal and depends only on the post-quench mode occupation.

arXiv:1402.1944•Physical Review B•Jun 30, 2014•94 citations

Thermodynamic Bethe ansatz for non-equilibrium steady states: exact energy current and fluctuations in integrable QFT

Olalla A. Castro-Alvaredo, Yixiong Chen, Benjamin Doyon +1 more

We evaluate the exact energy current and scaled cumulant generating function (related to the large-deviation function) in non-equilibrium steady states with energy flow, in any integrable model of relativistic quantum field theory (IQFT) with diagonal scattering. Our derivations are based on various recent results of D. Bernard and B. Doyon. The steady states are built by connecting homogeneously two infinite halves of the system thermalized at different temperatures Tl, Tr, and waiting for a long time. We evaluate the current J(Tl, Tr) using the exact QFT density matrix describing these non-equilibrium steady states and using Al.B. Zamolodchikov’s method of the thermodynamic Bethe ansatz (TBA). The scaled cumulant generating function is obtained from the extended fluctuation relations which hold in integrable models. We verify our formula in particular by showing that the conformal field theory (CFT) result is obtained in the high-temperature limit. We analyze numerically our non-equilibrium steady-state TBA equations for three models: the sinh-Gordon model, the roaming trajectories model, and the sine-Gordon model at a particular reflectionless point. Based on the numerics, we conjecture that an infinite family of non-equilibrium c-functions, associated to the scaled cumulants, can be defined, which we interpret physically. We study the full scaled distribution function and find that it can be described by a set of independent Poisson processes. Finally, we show that the “additivity ” property of the current, which is known to hold in CFT and was proposed to hold more generally, does not hold in general IQFT, that is J(Tl, Tr) is not of the form f(Tl)−f(Tr).

•Aug 14, 2016•23 citations

Thermodynamic Bethe ansatz for non-equilibrium steady states: exact energy current and fluctuations in integrable QFT

Olalla Castro-Alvaredo, Yixiong Chen, Benjamin Doyon +1 more

We evaluate the exact energy current and scaled cumulant generating function (related to the large-deviation function) in non-equilibrium steady states with energy flow, in any integrable model of relativistic quantum field theory (IQFT) with diagonal scattering. Our derivations are based on various recent results of D. Bernard and B. Doyon. The steady states are built by connecting homogeneously two infinite halves of the system thermalized at different temperatures $T_l$, $T_r$, and waiting for a long time. We evaluate the current $J(T_l,T_r)$ using the exact QFT density matrix describing these non-equilibrium steady states and using Al.B. Zamolodchikov's method of the thermodynamic Bethe ansatz (TBA). The scaled cumulant generating function is obtained from the extended fluctuation relations which hold in integrable models. We verify our formula in particular by showing that the conformal field theory (CFT) result is obtained in the high-temperature limit. We analyze numerically our non-equilibrium steady-state TBA equations for three models: the sinh-Gordon model, the roaming trajectories model, and the sine-Gordon model at a particular reflectionless point. Based on the numerics, we conjecture that an infinite family of non-equilibrium $c$-functions, associated to the scaled cumulants, can be defined, which we interpret physically. We study the full scaled distribution function and find that it can be described by a set of independent Poisson processes. Finally, we show that the "additivity" property of the current, which is known to hold in CFT and was proposed to hold more generally, does not hold in general IQFT, that is $J(T_l,T_r)$ is not of the form $f(T_l)-f(T_r)$.

arXiv:1310.4779•Journal of Statistical Mechanics Theory and Experiment•Mar 20, 2014•37 citations

Non-equilibrium dynamics of interacting many-body quantum systems in one dimension

Bruno Bertini

In this thesis we study three examples of interacting many-body systems undergoing a non equilibrium time evolution. Firstly we consider the time evolution in an integrable system: the sine-Gordon field theory in the repulsive regime. We will focus on the one point function of the semi-local vertex operator eiβφ(x)/2 on a specific class of initial states. By analytical means we show that the expectation value considered decays exponentially to zero at late times and we determine the decay time. The method employed is based on a form-factor expansion and uses the "Representative Eigenstate Approach" of Ref. [73] (a.k.a. "Quench Action"). In a second example we study the time evolution in models close to "special" integrable points characterised by hidden symmetries generating infinitely many local conservation laws that do not commute with one another, in addition to the infinite commuting family implied by integrability. We observe that both in the case where the perturbation breaks the integrability and when it breaks only the additional symmetries maintaining integrability, the local observables show a crossover behaviour from an initial to a final quasi stationary plateau. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin-1/2 chain with additional perturbations that break integrability. Finally, we study the effects of integrability breaking perturbations on the non-equilibrium evolution of more general many-particle quantum systems, where the unperturbed integrable model is generic. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalisations of quantum Boltzmann equations. We benchmark our method against time dependent density matrix renormalisation group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalisation plateaux for local observables on all accessible time scales. Increasing the strength of the integrability breaking term induces a "drift" away from the prethermalisation plateaux towards thermal behaviour. We identify a time scale characterising this crossover.

•Jan 1, 2015•1 citation
We construct non-equilibrium steady states in the Klein–Gordon theory in arbitrary space dimension d following a local quench. We consider the approach where two independently thermalized semi-infinite systems, with temperatures and , are connected along a -dimensional hypersurface. A current-carrying steady state, described by thermally distributed modes with temperatures and for left and right-moving modes, respectively, emerges at late times. The non-equilibrium density matrix is the exponential of a non-local conserved charge. We obtain exact results for the average energy current and the complete distribution of energy current fluctuations. The latter shows that the long-time energy transfer can be described by a continuum of independent Poisson processes, for which we provide the exact weights. We further describe the full time evolution of local observables following the quench. Averages of generic local observables, including the stress-energy tensor, approach the steady state with a power-law in time, where the exponent depends on the initial conditions at the connection hypersurface. We describe boundary conditions and special operators for which the steady state is reached instantaneously on the connection hypersurface. A semiclassical analysis of freely propagating modes yields the average energy current at large distances and late times. We conclude by comparing and contrasting our findings with results for interacting theories and provide an estimate for the timescale governing the crossover to hydrodynamics. As a modification of our Klein–Gordon analysis we also include exact results for free Dirac fermions.
Non-equilibrium steady states (NESS) describe particularly simple and stationary non-equilibrium situations. A possibility to obtain such states is to consider the asymptotic evolution of two infinite heat baths brought into thermal contact. In this work we generalise corresponding results of Doyon~et.~al. (J.\ Phys.\ A 18 (2015) no.9) for free Klein-Gordon fields in several directions. Our analysis is carried out directly at the level of correlation functions and in the algebraic approach to QFT. We discuss non-trivial chemical potentials, condensates, inhomogeneous linear models and homogeneous interacting ones. We shall not consider a sharp contact at initial time, but a smooth transition region. As a consequence, the states we construct will be of Hadamard type, and thus sufficiently regular for the perturbative treatment of interacting models. Our analysis shows that perturbatively constructed interacting NESS display thermodynamic properties similar to the ones of the NESS in linear models. In particular, perturbation theory appears to be insufficient to describe full thermalisation in non-linear QFT. Notwithstanding, we find that the NESS for linear and interacting models is stable under small perturbations, which is one of the characteristic features of equilibrium states.
Abstract We study out-of-equilibrium energy transport in a quantum critical fluid with Lifshitz scaling symmetry following a local quench between two semi-infinite fluid reservoirs. The late time energy flow is universal and is accommodated via a steady state occupying an expanding central region between outgoing shock and rarefaction waves. We consider the admissibility and entropy conditions for the formation of such a non-equilibrium steady state for a general dynamical critical exponent z in arbitrary dimensions and solve the associated Riemann problem. The Lifshitz fluid with z = 2 can be obtained from a Galilean boost invariant field theory and the non-equilibrium steady state is identified as a boosted thermal state. A Lifshitz fluid with generic z is scale invariant but without boost symmetry and in this case the non-equilibrium steady state is genuinely non-thermal.
We consider a quantum quench in a noninteracting fermionic one-dimensional field theory. The system of size LLL is initially prepared into two halves L\mathcal{L}L $([\ensuremath{-}L/2,0])$ and R\mathcal{R}R ([0,L/2])([0,L/2])([0,L/2]), each of them thermalized at two different temperatures TL{T}_{\mathcal{L}}TL​ and TR{T}_{\mathcal{R}}TR​, respectively. At a given time, the two halves are joined together by a local coupling and the whole system is left to evolve unitarily. For an infinitely extended system $(L\ensuremath{\rightarrow}\ensuremath{\infty})$, we show that the time evolution of the particle and energy densities is well described via a hydrodynamic approach which allows us to evaluate the correspondent stationary currents. We show, in such a case, that the two-point correlation functions are deduced, at large times, from a simple nonequilibrium steady state. Otherwise, whenever the boundary conditions are retained (in a properly defined thermodynamic limit), any current is suppressed at large times, and the stationary state is described by a generalized Gibbs ensemble, which is diagonal and depends only on the post-quench mode occupation.
We evaluate the exact energy current and scaled cumulant generating function (related to the large-deviation function) in non-equilibrium steady states with energy flow, in any integrable model of relativistic quantum field theory (IQFT) with diagonal scattering. Our derivations are based on various recent results of D. Bernard and B. Doyon. The steady states are built by connecting homogeneously two infinite halves of the system thermalized at different temperatures Tl, Tr, and waiting for a long time. We evaluate the current J(Tl, Tr) using the exact QFT density matrix describing these non-equilibrium steady states and using Al.B. Zamolodchikov’s method of the thermodynamic Bethe ansatz (TBA). The scaled cumulant generating function is obtained from the extended fluctuation relations which hold in integrable models. We verify our formula in particular by showing that the conformal field theory (CFT) result is obtained in the high-temperature limit. We analyze numerically our non-equilibrium steady-state TBA equations for three models: the sinh-Gordon model, the roaming trajectories model, and the sine-Gordon model at a particular reflectionless point. Based on the numerics, we conjecture that an infinite family of non-equilibrium c-functions, associated to the scaled cumulants, can be defined, which we interpret physically. We study the full scaled distribution function and find that it can be described by a set of independent Poisson processes. Finally, we show that the “additivity ” property of the current, which is known to hold in CFT and was proposed to hold more generally, does not hold in general IQFT, that is J(Tl, Tr) is not of the form f(Tl)−f(Tr).
We evaluate the exact energy current and scaled cumulant generating function (related to the large-deviation function) in non-equilibrium steady states with energy flow, in any integrable model of relativistic quantum field theory (IQFT) with diagonal scattering. Our derivations are based on various recent results of D. Bernard and B. Doyon. The steady states are built by connecting homogeneously two infinite halves of the system thermalized at different temperatures TlT_lTl​, TrT_rTr​, and waiting for a long time. We evaluate the current J(Tl,Tr)J(T_l,T_r)J(Tl​,Tr​) using the exact QFT density matrix describing these non-equilibrium steady states and using Al.B. Zamolodchikov's method of the thermodynamic Bethe ansatz (TBA). The scaled cumulant generating function is obtained from the extended fluctuation relations which hold in integrable models. We verify our formula in particular by showing that the conformal field theory (CFT) result is obtained in the high-temperature limit. We analyze numerically our non-equilibrium steady-state TBA equations for three models: the sinh-Gordon model, the roaming trajectories model, and the sine-Gordon model at a particular reflectionless point. Based on the numerics, we conjecture that an infinite family of non-equilibrium ccc-functions, associated to the scaled cumulants, can be defined, which we interpret physically. We study the full scaled distribution function and find that it can be described by a set of independent Poisson processes. Finally, we show that the "additivity" property of the current, which is known to hold in CFT and was proposed to hold more generally, does not hold in general IQFT, that is J(Tl,Tr)J(T_l,T_r)J(Tl​,Tr​) is not of the form f(Tl)−f(Tr)f(T_l)-f(T_r)f(Tl​)−f(Tr​).
In this thesis we study three examples of interacting many-body systems undergoing a non equilibrium time evolution. Firstly we consider the time evolution in an integrable system: the sine-Gordon field theory in the repulsive regime. We will focus on the one point function of the semi-local vertex operator eiβφ(x)/2 on a specific class of initial states. By analytical means we show that the expectation value considered decays exponentially to zero at late times and we determine the decay time. The method employed is based on a form-factor expansion and uses the "Representative Eigenstate Approach" of Ref. [73] (a.k.a. "Quench Action"). In a second example we study the time evolution in models close to "special" integrable points characterised by hidden symmetries generating infinitely many local conservation laws that do not commute with one another, in addition to the infinite commuting family implied by integrability. We observe that both in the case where the perturbation breaks the integrability and when it breaks only the additional symmetries maintaining integrability, the local observables show a crossover behaviour from an initial to a final quasi stationary plateau. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin-1/2 chain with additional perturbations that break integrability. Finally, we study the effects of integrability breaking perturbations on the non-equilibrium evolution of more general many-particle quantum systems, where the unperturbed integrable model is generic. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalisations of quantum Boltzmann equations. We benchmark our method against time dependent density matrix renormalisation group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalisation plateaux for local observables on all accessible time scales. Increasing the strength of the integrability breaking term induces a "drift" away from the prethermalisation plateaux towards thermal behaviour. We identify a time scale characterising this crossover.