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.
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 reexamine the emergence of a universal nonequilibrium steady state following a local quench between quantum critical heat baths in spatial dimensions greater than one. We show that energy transport proceeds by the formation of an instantaneous shock wave and a broadening rarefaction wave on either side of the interface, and not by two shock waves as previously proposed. For small temperature differences the universal steady state energy currents of the two-shock and rarefaction-shock solutions coincide. Over a broad range of parameters, the difference in the energy flow across the interface between these two solutions is at the level of 2%. The properties of the energy flow remain fully universal and independent of the microscopic theory. We briefly discuss the width of the shock wave in a viscous fluid, the effects of momentum relaxation, and the generalization to charged fluids.
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.
In these notes I explain how to describe one-dimensional quantum systems that are simultaneously near to, but not exactly at, a critical point, and in a far-from-equilibrium steady state. This description uses a density matrix on scattering states (of the type of Hershfield's density matrix), or equivalently a Gibbs-like ensemble of scattering states. The steady state I am considering is one where there is a steady flow of energy along the chain, coming from the steady draining / filling of two far-away reservoirs put at different temperatures. The context I am using is that of massive relativistic quantum field theory, which is the framework for describing the region near quantum critical points in any universality class with translation invariance and with dynamical exponent z equal to 1. I show, in this completely general setup, that a particular steady-state density matrix occurs naturally from the physically motivated Keldysh formulation of the steady state. In this formulation the steady state occurs as a result of a large-time evolution from an initial state where two halves of the system are separately thermalized. I also show how this suggests a particular dependence (a "factorization") of the average current on the left and right temperatures. The idea of this density matrix was proposed already in a recent publication with my collaborator Denis Bernard, where we studied it in the context of conformal field theory.
We study the dynamics of the entanglement in one-dimensional critical quantum systems after a local quench in which two independently thermalized semi-infinite halves are joined to form a homogeneous infinite system and left to evolve unitarily. We show that under certain conditions a nonequilibrium steady state (NESS) is reached instantaneously as soon as the entanglement interval is within the light cone emanating from the contact point. In this steady state, the exact expressions for the entanglement entropy and the logarithmic negativity are in agreement with the steady state density matrix being a boosted thermal state, as expected. We derive various general identities: relating the negativity after the quench with unequal left and right initial temperatures with that where the left and right temperatures are equal; and relating these with the negativity in equilibrium thermal states. In certain regimes the resulting expressions can be analytically evaluated. Immediately after the interval intersects the light cone, we find logarithmic growth. For a very long interval, we find that the negativity approaches a plateau after sufficiently long times, different from its NESS value. The NESS value is reached instantly as soon as the entire interval is contained in the light cone. This provides a theoretical framework explaining recently obtained numerical results.
Universality and scaling laws are hallmarks of equilibrium phase transitions and critical phenomena. However, extending these concepts to non-equilibrium systems is an outstanding challenge. Despite recent progress in the study of dynamical phases, the universality classes and scaling laws for non-equilibrium phenomena are far less understood than those in equilibrium. In this work, using a trapped-ion quantum simulator with single-spin resolution, we investigate the non-equilibrium nature of critical fluctuations following a quantum quench to the critical point. We probe the scaling of spin fluctuations after a series of quenches to the critical Hamiltonian of a long-range Ising model. With systems of up to 50 spins, we show that the amplitude and timescale of the post-quench fluctuations scale with system size with distinct universal critical exponents, depending on the quench protocol. While a generic quench can lead to thermal critical behavior, we find that a second quench from one critical state to another (i.e. a double quench) results in a new universal non-equilibrium behavior, identified by a set of critical exponents distinct from their equilibrium counterparts. Our results demonstrate the ability of quantum simulators to explore universal scaling beyond equilibrium.
Using holographic duality, we investigate thermalization process when two finite-size quantum critical systems are brought into thermal contact along a perfectly transmitting interface. Through real-time simulations of gravitational dynamics, which are spatially inhomogeneous and anisotropic and are confined within two dynamical bulk branes, we identify three distinct thermalization patterns governed by the energy imbalance (temperature difference) and system size. For systems with large size and small energy imbalance, we observe recurrent cycles of formation and collapse of non-equilibrium steady states (NESS). Under large energy imbalance, shock waves persist for a prolonged period with sustained boundary reflections, while rarefaction waves rapidly homogenize. When the system size is sufficiently small, dissipation dominates and leads to oscillatory decay without sustained NESS or shock structure. In sharp contrast to diffusive systems, we uncover that wave-propagated energy transfer together with boundary reflections enables nearly complete energy swapping between subsystems during thermalization. Our results reveal rich thermalization dynamics in finite-size quantum critical systems across spatial scales and energy gradient regimes.