Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. In this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.
Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. In this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.
Traditional inelastic neutron scattering (INS) characterizes excitations-such as phonons and magnons-in condensed matter systems at thermodynamic equilibrium. However, the most intriguing and puzzling many-body effects in open quantum systems often emerge from dissipative dynamics that are inherently out of equilibrium. Here, we use a combination of laser pumping and INS to experimentally observe long-lived nonequilibrium magnons in a two-dimensional (2D) square-lattice Heisenberg antiferromagnet. These nonequilibrium magnons manifest themselves as a violation of detailed balance in the dynamic structure factor and reach steady states under periodic driving, analogous to nonequilibrium steady states in driven dissipative systems. Furthermore, we show that the violation of detailed balance reflects the quantum-mechanical nature of the underlying dynamical system, where out-of-time-ordered correlations of creation and annihilation operators do not satisfy commutation relations. The in operando INS technique developed here provides a new approach to studying nonequilibrium magnons in prototypical 2D quantum magnets and can be extended to other systems, including one-dimensional spin chains and topological many-body spin systems, where nonequilibrium effects are widespread and rich in discovery potential.
A deeper theoretical understanding of driven-dissipative interacting systems and their nonequilibrium phase transitions is essential both to advance our fundamental physics understanding and to harness technological opportunities arising from optically controlled quantum many-body states. This paper provides a numerical study of dynamical phases and the transitions between them in the nonequilibrium steady state of the prototypical two-dimensional Heisenberg antiferromagnet with drive and dissipation. We demonstrate a nonthermal transition that is characterized by a qualitative change in the magnon distribution from subthermal at low drive to a generalized Bose-Einstein form including a nonvanishing condensate fraction at high drive. A finite-size analysis reveals static and dynamical critical scaling at the transition, with a discontinuous slope of the magnon number versus driving field strength and critical slowing down at the transition point. Implications for experiments on quantum materials and polariton condensates are discussed.
Metastability, i.e., partial relaxation to long-lived, quasi-stationary states before true asymptotic equilibrium sets in, emerges ubiquitously in classical and quantum dynamical systems as a result of timescales separation. In open quantum systems, an intrinsically nonequilibrium analogue, dynamical metastability, can originate from the spectral geometry of a non-Hermitian operator. In noninteracting models, this mechanism produces boundary-sensitive anomalous relaxation, transient amplification, and topologically mandated long-lived edge modes, all of which are enhanced as system size grows. Here we extend dynamical metastability into the nonlinear, interacting regime and identify magnetic heterostructures as a natural platform for its exploration. We introduce an interacting spin Lindbladian whose linearized magnon dynamics map onto a dynamically metastable Hatano-Nelson chain, and show that dynamical metastability in the noninteracting limit seeds genuinely nonlinear phenomena, including size-dependent spin dipping and anomalous attraction to unstable equilibria. Long-lived edge states associated to topologically mandated Dirac bosons persist under nonlinearities and disorder. We further analyze the magnetization dynamics in magnetic multilayers within the classical Landau-Lifshitz-Gilbert-Slonczewski framework, identifying Dzyaloshinskii-Moriya interaction, nonlocal damping, and spin-transfer torque as control parameters governing bulk-boundary stability mismatch and band topology. While all the distinctive dynamical phenomena previously identified reappear in this experimentally relevant setting, the LLGS framework also supports multistability and limit cycles that are absent in the quantum model. Our results constitute the first systematic study of dynamical metastability in nonlinear dynamics, directly relevant to spin-torque oscillator arrays, magnonic devices, and beyond.
We investigate nonequilibrium steady states of driven-dissipative ideal quantum gases of both bosons and fermions. We focus on systems of sharp particle number that are driven out of equilibrium either by the coupling to several heat baths of different temperature or by time-periodic driving in combination with the coupling to a heat bath. Within the framework of (Floquet-)Born-Markov theory, several analytical and numerical methods are described in detail. This includes a mean-field theory in terms of occupation numbers, an augmented mean-field theory taking into account also nontrivial two-particle correlations, and quantum-jump-type Monte Carlo simulations. For the case of the ideal Fermi gas, these methods are applied to simple lattice models and the possibility of achieving exotic states via bath engineering is pointed out. The largest part of this work is devoted to bosonic quantum gases and the phenomenon of Bose selection, a nonequilibrium generalization of Bose condensation, where multiple single-particle states are selected to acquire a large occupation [Phys. Rev. Lett. 111, 240405 (2013)]. In this context, among others, we provide a theory for transitions where the set of selected states changes, describe an efficient algorithm for finding the set of selected states, investigate beyond-mean-field effects, and identify the dominant mechanisms for heat transport in the Bose-selected state.
We study three aspects of quantum nonequilibrium dynamics; (1) transport of conserved quantities, (2) entanglement spreading, and (3) construction of local operators which slowly relax to thermal equilibrium. Motivated by recent progresses in ultracold atom experiments, we first analyze transport phenomena of a population imbalanced two-component fermi gas with arbitrary strength of inter-species interaction in three dimension. Using the Boltzmann kinetic equation at dilute regime, we obtain the transport coefficients of linear responses to gradients of temperature and chemical potential imbalance. We identify the magneto-caloric effects, and determine how these effects depend on interaction strength and population imbalance. Then, we propose an experimental protocol to observe these effects in an experiment with ultracold atoms. Next, we study entanglement spreading in a one-dimensional quantum Ising chain with longitudinal and transverse fields, which is diffusive and nonintegrable. Fully diagonalizing the Hamiltonian matrix, we explicitly show that the entanglement spreading is ballistic, thus faster than diffusive transport of conserved quantities. We provide a local spreading picture of entanglement entropy in terms of logarithmic negativity. Then, we discuss a role of energy conservation in entanglement spreading through analyzing a Floquet system. Lastly, we construct local operators that relax slowly to equilibrium with the same one-dimensional Hamiltonian. By showing that the Hamiltonian satisfies the eigenstate thermalization hypothesis, we first conclude that this model relaxes local operators to thermal equilibrium. Then, we systematically construct local operators that relax slower than conventional diffusive modes.
We investigate the relaxation dynamics of two distinct non-equilibrium processes: relaxation of three-dimensional antiferromagnetic lattice spin models with Heisenberg interaction following a critical quench, and a one-dimensional exclusion process inspired by the gear-like motion of molecular motors. In a system of three-dimensional Heisenberg antiferromagnets the non-conserved staggered magnetization components couple non-trivially to the conserved magnetization densities inducing fully reversible terms that enter the Langevin dynamic equation. We simulate the exact microscopic dynamics of such a system of antiferromagnets by employing a hybrid simulation algorithm that combines the reversible spin precession implemented by the fourth-order Runge-Kutta integration method with the standard relaxational dynamics at finite temperatures using Monte Carlo updates. We characterize the dynamic universality class of this system by probing the early temporal window where the system exhibits aging scaling properties. We also verify an earlier renormalization group prediction that the temporal decay exponent in the two-time spin autocorrelation function exhibits non-universality, specifically it depends on the width of the initial spin orientation distribution. We employ a similar numerical technique to study the critical dynamics of an anisotropic Heisenberg antiferromagnet in the presence of an external field. The phase diagram of this system exhibits two critical lines that meet at a bicritical point. We study the aging scaling dynamics for the model C critical line, probe the model F critical line by investigating the system size dependence of the characteristic spin-wave frequencies near criticality, and measure the dynamic critical exponents for the order parameter including its aging scaling at the bicritical point. We introduce a one-dimensional non-equilibrium lattice gas model representing the processive motion of dynein molecular motors over the microtubule. We study both dynamical and stationary state properties for the model consisting of hardcore particles hopping on the lattice with variable step sizes. We find that the stationary state gap-distribution exhibits striking peaks around gap sizes that are multiples of the maximum step size, for both open and periodic boundary conditions, and verify this using a mean-field calculation. For open boundary conditions, we observe intriguing damped oscillator-like distribution of particles over the lattice with a periodicity equal to the maximum step size. To characterize transient dynamics, we measure the mean square displacement that shows weak superdiffusive growth with exponent γ≈ 1.34 for periodic boundary and ballistic growth ( γ≈ 2) for open boundary conditions at early times. We also study the effect of Langmuir dynamics on the density profile.