Nonequilibrium Steady State Research Articles

Emerging new phases in correlated Mott insulatorCa2RuO4.

The Mott insulator Ca2RuO4is a paradigmatic example among transition
metal oxides, where the interplay of charge, spin, orbital, and lattice degrees of freedom leads to competing quantum phases. In this paper, we focus on and review some key aspects, from the underlying physical framework and its basic properties, to recent theoretical efforts that aim to trigger unconventional quantum ground states, using several external parameters and stimuli. Using first-principle calculations, we demonstrate that Ca2RuO4shows a spin splitting in the reciprocal space, and identify it as an altermagnetic candidate material. The non relativistic spin-splitting has an orbital selective nature, dictated by the local crystallographic symmetry. Next, we consider two routes that may trigger exotic quantum states. The first one corresponds to transition metal substitution of the 4d4Ru with isovalent 3d3ions. This substitutional doping may alter the spin-orbital correlations favoring the emergence of negative thermal expansion. The second route explores fledgling states arising in a nonequilibrium steady state under the influence of an applied electric field. We show that the electric field can directly affect the orbital density, eventually leading to strong orbital fluctuations and the suppression of orbital imbalance, which may, in turn, reduce antiferromagnetism. These aspects suggest possible practical applications, as its unique properties may open up possibilities for augmenting existing technologies, surpassing the limitations of conventional materials.

Minimum time connection between non-equilibrium steady states: the Brownian gyrator

Minimum time connection between non-equilibrium steady states: the Brownian gyrator

Relaxed hydrodynamic theory of electrically driven nonequilibrium steady states

The capability of hydrodynamics to accurately describe slow and long-wavelength fluctuations around nonequilibrium steady states (NESS), characterized by a stationary flow of energy or matter in the presence of a driving force, remains an open question. In this study, we explicitly construct a hydrodynamic description of electrically driven nonequilibrium charged steady states in the limit in which the relaxation of the first nonhydrodynamic excitation is parametrically slow. Our approach involves introducing gapped modes and extending the effective description into a relaxed hydrodynamic theory (RHT). Leveraging the gauge-gravity duality as a tool for controlled computations within nonequilibrium systems, we establish an ultraviolet complete model for these NESS that confirms the validity of our RHT. In summary, our findings provide a concrete realization of the validity of hydrodynamics beyond thermal equilibrium, offering valuable insights into the dynamics of nonequilibrium systems. Published by the American Physical Society 2024

On the nonequilibrium dynamics of gravitational algebras

Abstract We explore nonequilibrium features of certain operator algebras which appear in quantum gravity. The algebra of observables in a black hole background is a Type II ∞ von Neumann algebra. We discuss how this algebra can be coupled to the algebra of observable of an infinite reservoir within the canonical ensemble, aiming to induce nonequilibrium dynamics. The resulting dynamics can lead the system towards a nonequilibrium steady state which can be characterized through modular theory. Within this framework we address the definition of entropy production and its relationship to relative entropy, alongside exploring other applications.

Stable-to-unstable transition in quantum friction

We investigate the frictional force arising from quantum fluctuations when two dissipative metallic plates are set in a shear motion. While early studies showed that the electromagnetic fields in the quantum friction setup reach nonequilibrium steady states, yielding a time-independent force, other works have demonstrated the failure to attain steady states, leading to instability and time-varying friction under sufficiently low-loss conditions. Here, we develop a fully quantum-mechanical theory without perturbative approximations, and we unveil the transition from stable to unstable regimes of the quantum friction setup. Due to the relative motion of the plates, their electromagnetic response may be active in some conditions, resulting in optical gain. We prove that the standard fluctuation-dissipation leads to inconsistent results when applied to our system, and, in particular, it predicts a vanishing frictional force. Using a modified fluctuation-dissipation relation tailored for gain media, we calculate the frictional force in terms of the system Green's function, thereby recovering early works on quantum friction. Remarkably, we also find that the frictional force diverges to infinity as the relative velocity of the plates approaches a threshold. This threshold is determined by the damping strength and the distance between the metal surfaces. Beyond this critical velocity, the system exhibits instability, akin to the behavior of a laser cavity, where no steady state exists. In such a scenario, the frictional force escalates exponentially. Our findings pave the way for experimental exploration of the frictional force in proximity to this critical regime. Published by the American Physical Society 2024

Nonequilibrium steady state full counting statistics in the noncrossing approximation.

Quantum transport is often characterized not just by mean observables like the particle or energy current but by their fluctuations and higher moments, which can act as detailed probes of the physical mechanisms at play. However, relatively few theoretical methods are able to access the full counting statistics (FCS) of transport processes through electronic junctions in strongly correlated regimes. While most experiments are concerned with steady state properties, most accurate theoretical methods rely on computationally expensive propagation from a tractable initial state. Here, we propose a simple approach for computing the FCS through a junction directly at the steady state, utilizing the propagator noncrossing approximation. Compared to time propagation, our method offers reduced computational cost at the same level of approximation, but the idea can also be used within other approximations or as a basis for numerically exact techniques. We demonstrate the method's capabilities by investigating the impact of lead dimensionality on electronic transport in the nonequilibrium Anderson impurity model at the onset of Kondo physics. Our results reveal a distinct signature of one dimensional leads in the noise and Fano factor not present for other dimensionalities, showing the potential of FCS measurements as a probe of the environment surrounding a quantum dot.

Partial stochastic resetting with refractory periods

Abstract The effect of refractory periods in partial resetting processes is studied. Under Poissonian partial resets, a state variable jumps to a value closer to the origin by a fixed fraction at constant rate, $x\to a x$. Following each reset, a stationary refractory period of arbitrary duration takes place. We derive an exact closed-form expression for the propagator in Fourier-Laplace space, which shows rich dynamical features such as connections not only to other resetting schemes but also to intermittent motion. For diffusive processes, we use the propagator to derive exact expressions for time dependent moments of $x$ at all orders. At late times the system reaches a non-equilibrium steady state which takes the form of a mixture distribution that splits the system into two subpopulations; trajectories that at any given time in the stationary regime find themselves in the freely evolving phase, and those that are in the refractory phase. In contrast to conventional resetting, partial resets give rise to non-trivial steady states even for the refractory subpopulation. Moments and cumulants associated with the steady state density are studied, and we show that a universal optimum for the kurtosis can be found as a function of mean refractory time, determined solely by the strength of the resetting and the mean inter-reset time. The presented results could be of relevance to growth-collapse processes with periods of inactivity following a collapse.

Generalized hydrodynamics and approach to generalized Gibbs equilibrium for a classical harmonic chain

Abstract We study the evolution of a classical harmonic chain with nearest-neighbor interactions starting from domain wall initial conditions. The initial state is taken to be either a product of two Gibbs ensembles (GEs) with unequal temperatures on the two halves of the chain or a product of two generalized Gibbs ensembles (GGEs) with different parameters in the two halves. For this system, we construct the Wigner function and demonstrate that its evolution defines the generalized hydrodynamics (GHD) describing the evolution of the conserved quantities. We solve the GHD for both finite and infinite chains and compute the evolution of conserved densities and currents. For a finite chain with fixed boundaries, we show that these quantities relax as ∼ 1 / t to their respective steady-state values given by the final expected GE or GGE state, depending on the initial conditions. Exact expressions for the Lagrange multipliers of the final expected GGE state are obtained in terms of the steady-state densities. In the case of an infinite chain, we find that the conserved densities and currents at any finite time exhibit ballistic scaling while, at infinite time, any finite segment of the system can be described by a current-carrying non-equilibrium steady state (NESS). We compute the scaling functions analytically and show that the relaxation to the NESS occurs as ∼ 1 / t for the densities and as ∼ 1 / t 2 for the currents. We compare the analytical results from hydrodynamics with those from exact microscopic numerics and find excellent agreement.

Nonequilibrium Steady States in the Weakly-Coupled XXZ Model

Nonequilibrium Steady States in the Weakly-Coupled XXZ Model

A local–global principle for nonequilibrium steady states

The global steady state of a system in thermal equilibrium exponentially favors configurations with lesser energy. This principle is a powerful explanation of self-organization because energy is a local property of configurations. For nonequilibrium systems, there is no such property for which an analogous principle holds, hence no common explanation of the diverse forms of self-organization they exhibit. However, a flurry of recent empirical results has shown that a local property of configurations called "rattling" predicts the steady states of some nonequilibrium systems, leading to claims of a far-reaching principle of nonequilibrium self-organization. But for which nonequilibrium systems is rattling accurate, and why? We develop a theory of rattling in terms of Markov processes that gives simple and precise answers to these key questions. Our results show that rattling predicts a broader class of nonequilibrium steady states than has been claimed and for different reasons than have been suggested. Its predictions hold to an extent determined by the relative variance of, and correlation between, the local and global "parts" of a steady state. We show how these quantities characterize the local-global relationships of various random walks on random graphs, spin-glass dynamics, and models of animal collective behavior. Surprisingly, we find that the core idea of rattling is so general as to apply to equilibrium and nonequilibrium systems alike.

Energy Aware Technology Mapping of Genetic Logic Circuits.

Energy and its dissipation are fundamental to all living systems, including cells. Insufficient abundance of energy carriers─as caused by the additional burden of artificial genetic circuits─shifts a cell's priority to survival, also impairing the functionality of the genetic circuit. Moreover, recent works have shown the importance of energy expenditure in information transmission. Despite living organisms being non-equilibrium systems, non-equilibrium models capable of accounting for energy dissipation and non-equilibrium response curves are not yet employed in genetic design automation (GDA) software. To this end, we introduce Energy Aware Technology Mapping, the automated design of genetic logic circuits with respect to energy efficiency and functionality. The basis for this is an energy aware non-equilibrium steady state model of gene expression, capturing characteristics like energy dissipation─which we link to the entropy production rate─and transcriptional bursting, relevant to eukaryotes as well as prokaryotes. Our evaluation shows that a genetic logic circuit's functional performance and energy efficiency are disjoint optimization goals. For our benchmark, energy efficiency improves by 37.2% on average when comparing to functionally optimized variants. We discover a linear increase in energy expenditure and overall protein expression with the circuit size, where Energy Aware Technology Mapping allows for designing genetic logic circuits with the energetic costs of circuits that are one to two gates smaller. Structural variants improve this further, while results show the Pareto dominance among structures of a single Boolean function. By incorporating energy demand into the design, Energy Aware Technology Mapping enables energy efficiency by design. This extends current GDA tools and complements approaches coping with burden in vivo.

Soluble Model of a Nonequilibrium Steady State: The Van Kampen Objection and Other Lessons.

A simple model of charge transport is provided by a classical particle in a random potential and a dissipative coupling to the environment. The corresponding nonequilibrium steady state (NESS) can be determined analytically when both the disorder and dissipation are weak. We use it to illuminate some foundational issues in nonequilibrium statistical mechanics. We show that at nonlinear level dissipative response sensitively depends on the system-environment coupling, and the range of validity of the linear response theory is set by this coupling. This validates the Van Kampen objection. We also show that the principle of minimum entropy production does not determine the NESS beyond linear order in the electric field, while entropy maximization fails to produce the correct NESS already at linear order.

An optimization-based equilibrium measure describes fixed points of non-equilibrium dynamics: application to edge of chaos

Abstract Understanding neural dynamics is a central topic in machine learning, non-linear physics and neuroscience. However, the dynamics is non-linear, stochastic and particularly non-gradient, i.e., the driving force can not be written as gradient of a potential. These features make analytic studies very challenging. The common tool is the path integral approach or dynamical mean-field theory, but the drawback is that one has to solve the integro-differential or dynamical mean-field equations, which is computationally expensive and has no closed form solutions in general. From the aspect of associated Fokker-Planck equation, the steady state solution is generally unknown. Here, we treat searching for the steady states as an optimization problem, and construct an approximate potential related to the speed of the dynamics, and find that searching for the ground state of this potential is equivalent to running an approximate stochastic gradient dynamics or Langevin dynamics. Only in the zero temperature limit, the distribution of the original steady states can be achieved. The resultant stationary state of the dynamics follows exactly the canonical Boltzmann measure. Within this framework, the quenched disorder intrinsic in the neural networks can be averaged out by applying the replica method, which leads naturally to order parameters for the non-equilibrium steady states. Our theory reproduces the well-known result of edge-of-chaos, and further the order parameters characterizing the continuous transition are derived, and the order parameters are explained as fluctuations and responses of the steady states. Our method thus opens the door to analytically study the steady state landscape of the deterministic or stochastic high dimensional dynamics.

Localization in Boundary-Driven Lattice Models

Several systems display an equilibrium condensation transition, where a finite fraction of a conserved quantity is spatially localized. The presence of two conservation laws may induce the emergence of such transition in an out-of-equilibrium setup, where boundaries are attached to different and subcritical heat baths. We study this phenomenon in a class of stochastic lattice models, where the local energy is a general convex function of the local mass, mass and energy being both globally conserved in the isolated system. We obtain exact results for the nonequilibrium steady state (spatial profiles, mass and energy currents, Onsager coefficients) and we highlight important differences between equilibrium and out-of-equilibrium condensation.

Energy Aware Technology Mapping of Genetic Logic Circuits.

Energy and its dissipation are fundamental to all living systems, including cells. Insufficient abundance of energy carriers -as caused by the additional burden of artificial genetic circuits- shifts a cell's priority to survival, also impairing the functionality of the genetic circuit. Moreover, recent works have shown the importance of energy expenditure in information transmission. Despite living organisms being non-equilibrium systems, non-equilibrium models capable of accounting for energy dissipation and non-equilibrium response curves are not yet employed in genetic design automation (GDA) software. To this end, we introduce Energy Aware Technology Mapping, the automated design of genetic logic circuits with respect to energy efficiency and functionality. The basis for this is an energy aware non-equilibrium steady state (NESS) model of gene expression, capturing characteristics like energy dissipation -which we link to the entropy production rate- and transcriptional bursting, relevant to eukaryotes as well as prokaryotes. Our evaluation shows that a genetic logic circuit's functional performance and energy efficiency are disjoint optimization goals. For our benchmark, energy efficiency improves by 37.2% on average when comparing to functionally optimized variants. We discover a linear increase in energy expenditure and overall protein expression with the circuit size, where Energy Aware Technology Mapping allows for designing genetic logic circuits with the energy efficiency of circuits that are one to two gates smaller. Structural variants improve this further, while results show the Pareto dominance among structures of a single Boolean function. By incorporating energy demand into the design, Energy Aware Technology Mapping enables energy efficiency by design. This extends current GDA tools and complements approaches coping with burden in vivo.

Tensor-network-based variational Monte Carlo approach to the non-equilibrium steady state of open quantum systems

We introduce a novel method of efficiently simulating the non-equilibrium steady state of large many-body open quantum systems with highly non-local interactions, based on a variational Monte Carlo optimization of a matrix product operator ansatz. Our approach outperforms and offers several advantages over comparable algorithms, such as an improved scaling of the computational cost with respect to the bond dimension for periodic systems. We showcase the versatility of our approach by studying the phase diagrams and correlation functions of the dissipative quantum Ising model with collective dephasing and long-ranged power law interactions for spin chains of up to N=100 spins.

Semi-infinite Simple Exclusion Process: From Current Fluctuations to Target Survival.

The symmetric simple exclusion process (SEP), where diffusive particles cannot overtake each other, is a paradigmatic model of transport in the single-file geometry. In this model, the study of currents has attracted a lot of attention, but so far most results are restricted to two geometries: (i)a finite system between two reservoirs, which does not conserve the number of particles but reaches a nonequilibrium steady state, and (ii)an infinite system which conserves the number of particles but never reaches a steady state. Here, we obtain an expression for the full cumulant generating function of the integrated current in the important intermediate situation of a semi-infinite system connected to a reservoir, which does not conserve the number of particles and never reaches a steady state. This results from the determination of the full spatial structure of the correlations, which we infer to obey the very same closed equation recently obtained in the infinite geometry and argue to be exact. Besides their intrinsic interest, these results allow us to solve two open problems: the survival probability of a fixed target in the SEP, and the statistics of the number of particles injected by a localized source.

Steady-state properties of multi-orbital systems using quantum Monte Carlo.

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by signproblems. A dynamical signproblem makes it exponentially difficult to simulate long times. A multi-orbital signproblem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both signproblems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Long-range velocity correlations from active dopants

Active matter systems display collective behaviors that are impossible in thermodynamic equilibrium. One such feature, observed in in dense active matter systems is the appearance of long-range velocity correlations without explicit aligning interaction. However, the conditions for the appearance of these correlations remain largely unexplored. Here we show that such long-range velocity correlations can also be generated in a dense athermal passive system by the inclusion of a very small fraction of active Brownian particles. We develop a continuum theory to explain the emergence of velocity correlations generated via such active dopants. We validate the predictions for the effects of magnitude and persistence time of the active force and the area fractions of active and passive particles using extensive Brownian dynamics simulation of a canonical active-passive mixture. Our work decouples the roles that density and activity play in generating long-range velocity correlations in such exotic non-equilibrium steady states.

Computationally feasible bounds for the free energy of nonequilibrium steady states, applied to simple models of heat conduction

In this paper, we study computationally feasible bounds for relative free energies between two many-particle systems. Specifically, we consider systems out of equilibrium that admit a nonequilibrium steady state that is reached asymptotically in the long-time limit. The bounds that we suggest are based on the well-known Bogoliubov inequality and variants of Gibbs' and Donsker–Varadhan variational principles. As a general paradigm, we consider systems of oscillators coupled to heat baths at different temperatures. For such systems, we define the free energy of the system relative to any given reference system in terms of the Kullback–Leibler divergence between steady states. By employing a two-sided Bogoliubov inequality and a mean-variance approximation of the free energy (or cumulant generating function), we can efficiently estimate the free energy cost needed in passing from the reference system to the system out of equilibrium (characterised by a temperature gradient). A specific test case to validate our bounds are harmonic oscillator chains with ends that are coupled to Langevin thermostats at different temperatures; such a system is simple enough to allow for analytic calculations and general enough to be used as a prototype to estimate, e.g. heat fluxes or interface effects in a larger class of nonequilibrium particle systems.