Collective spin systems -- driven-dissipative spin ensembles whose unitary and dissipative processes act uniformly on all spins, allowing their dynamics to be described as a single macrospin -- play an important role in atomic and solid-state physics. Their intrinsic nonlinearity gives rise to multiple long-lived metastable states that ultimately relax to a unique most probable state. This dominant state can change with a control parameter, leading to first-order phase transitions. We develop a real-time instanton approach based on quantum quasiprobability dynamics that captures the stationary state in the large-spin limit and the asymptotic scaling of relaxation rates. We further show that these features are not accurately described by the previously applied semiclassical Wigner approach due to its neglect of non-Gaussian fluctuations.
Comments
9 pages, 4 figures + supplemental material. Companion paper added at arXiv:2610.21684
Collective spin systems -- driven-dissipative spin ensembles whose unitary and dissipative processes act uniformly on all spins, allowing their dynamics to be described as a single macrospin -- play an important role in atomic and solid-state physics. Their intrinsic nonlinearity gives rise to multiple long-lived metastable states that ultimately relax to a unique most probable state. This dominant state can change with a control parameter, leading to first-order phase transitions. We develop a real-time instanton approach based on quantum quasiprobability dynamics that captures the stationary state in the large-spin limit and the asymptotic scaling of relaxation rates. We further show that these features are not accurately described by the previously applied semiclassical Wigner approach due to its neglect of non-Gaussian fluctuations.
Driven-dissipative collective spin systems can host several long-lived metastable states; rare switching between them governs slow relaxation and first-order dissipative phase transitions associated with changes in the dominant stationary state. In the companion Letter [arXiv:2604.15091], we introduced a quantum instanton framework based on exact quasiprobability dynamics to describe this switching and illustrated it with a minimal model possessing a simplifying symmetry. Here we develop its numerical implementation, combining collocation and shooting methods, and apply it to two microscopically motivated systems: the Lipkin-Meshkov-Glick model with purely collective dissipation and the cooperative resonance fluorescence with collective and local dissipation. By determining the leading exponential scaling of switching rates between metastable attractors, we predict the stationary magnetization, the first-order dissipative phase-transition point, and the exponential scaling of slow relaxation rates. These predictions agree with quantum master-equation calculations and an exact analytical result for the second model without dephasing. We also examine how the choice of quasiprobability representation affects numerical convergence and find that the alternative semiclassical Wigner approach, based on Gaussian truncation, has varying accuracy across models and parameter regimes. These results substantiate the framework's applicability to experimentally relevant collective spin systems with both collective and local dissipation.
Metastability in many-body quantum systems and its associated exponentially-long timescales have been the subject of considerable recent interest. Here, we focus on a class of driven-dissipative many-body open quantum systems described by a Lindbladian having hidden time-reversal symmetry (a form of quantum detailed balance). Examples include boundary-driven interacting spin chains, bosonic lattice models and driven-dissipative collective spin models. We suggest that for such systems, slow timescales in the vicinity of a dissipative first-order phase transition can be analytically predicted using a special purification of the non-equilibrium steady state. We show the accuracy of our conjecture through detailed studies of a dissipative transverse-field Ising model with collective and local decay, and a driven-dissipative nonlinear cavity model. Our results allow quantitative insights into metastability and slow dynamics for a range of systems, including cases where semiclassical or path-integral instanton approaches are intractable.
We apply a recently developed theory for metastability in open quantum systems to a one-dimensional dissipative quantum Ising model. Earlier results suggest this model features either a nonequilibrium phase transition or a smooth but sharp crossover, where the stationary state changes from paramagnetic to ferromagnetic, accompanied by strongly intermittent emission dynamics characteristic of first-order coexistence between dynamical phases. We show that for a range of parameters close to this transition or crossover point the dynamics of the finite system displays pronounced metastability, i.e., the system relaxes first to long-lived metastable states before eventual relaxation to the true stationary state. From the spectral properties of the quantum master operator we characterize the low-dimensional manifold of metastable states, which are shown to be probability mixtures of two, paramagnetic and ferromagnetic, metastable phases. We also show that for long times the dynamics can be approximated by a classical stochastic dynamics between the metastable phases that is directly related to the intermittent dynamics observed in quantum trajectories and thus the dynamical phases.
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 present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.
We propose a semiclassical framework for solving open quantum dynamics in driven-dissipative spin systems. Our method consists of generalized spin-wave approximations tailored to describing quantum trajectories unravelled from the master equation, and generically applies to regimes beyond the reach of conventional spin-wave theories, including short-range interactions and local quantum jumps, enabling the efficient simulation of large-scale interacting spins. We illustrate the versatility of our framework by studying a variable-range driven-dissipative Ising model on a 2D lattice. When the dissipation acts along the drive axis, we find a continuous phase transition breaking the Z2 symmetry, and demonstrate that the interaction range, when tuned from fully-connected to nearest-neighbour, profoundly alters the universality class of the criticality. With the dissipation along the interaction axis, we show the emergence of a first-order transition. Demonstrated with both state-diffusion and quantum-jump types of trajectory dynamics, our framework provides a powerful toolbox for the efficient semiclassical description of non-equilibrium dynamics and many-body phases in spin systems.