First-Order Dissipative Phase Transitions

• First-order dissipative phase transitions are non-equilibrium phenomena characterized by abrupt discontinuities in order parameters due to engineered dissipation in open many-body systems.
• These transitions exhibit hallmark features such as bistability, hysteresis, and phase coexistence, modeled effectively with Lindblad master equations and observed in systems like Kerr resonators and spin chains.
• Critical slowing down occurs as the Liouvillian spectral gap closes exponentially with system size, offering robust platforms for quantum simulations and precision sensing applications.

A first-order dissipative phase transition is a non-equilibrium phenomenon in open quantum or classical many-body systems where the steady-state of the driven-dissipative dynamics displays a discontinuous change of some order parameter as a system or control parameter is varied. Unlike their equilibrium counterparts, dissipative phase transitions depend not only on the underlying Hamiltonian but also crucially on the nature of the coupling to the environment and the details of dissipation, as captured by Lindblad or more general master equations. Such transitions manifest as abrupt jumps in observables, bistability, phase coexistence, hysteresis, and the exponential closing of the Liouvillian spectral gap. These features have been extensively characterized in spin chains, bosonic cavities, solid-state platforms, and models of collective spin systems, and their understanding is supported by rigorous analytic solutions, semiclassical theory, and numerically exact methods.

1. Definition and Characteristic Phenomena

A first-order dissipative phase transition (DPT) occurs when the nonequilibrium steady-state of an open, driven quantum system undergoes a discontinuous change in a macroscopic observable (order parameter) as control parameters (e.g., drive strength, dissipation rate, or external field) are tuned. Hallmarks include:

• Bistability and Hysteresis: The steady state exhibits two or more macroscopically distinct stable solutions in a finite parameter region. Sweeping the control parameter up or down leads to hysteretic loops due to the metastability of the phases (Casteels et al., 2016, Röhrle et al., 2023).
• Phase Coexistence: Within the coexistence regime, the system can randomly switch between phases, and the steady-state can be decomposed as a mixture of the two, e.g., $\rho_\mathrm{ss} = \frac12(\rho^+ + \rho^-)$ in the thermodynamic limit (Minganti et al., 2018).
• Discontinuity in Observables: The order parameter (e.g., photon number, magnetization, condensate density) jumps discontinuously at the critical point (Dagvadorj et al., 2021, Beaulieu et al., 2023, Raghunandan et al., 2017).
• Exponential Closing of Liouvillian Gap: The slowest decay rate (Liouvillian spectral gap) vanishes exponentially with system size, leading to extremely long-lived metastable states (critical slowing down) (Casteels et al., 2016, Minganti et al., 2018, Beaulieu et al., 2023).
• Metastability and Critical Fluctuations: For finite-size systems, long-lived metastable states and large critical fluctuations are pronounced in the vicinity of the transition (Meglio et al., 2020, Röhrle et al., 2023).

2. Theoretical Framework and Model Examples

The prototypical setting for analyzing first-order dissipative phase transitions is an open many-body system governed by a Lindblad master equation,

$$\frac{d\rho}{dt} = -i[H, \rho] + \sum_j u_j \mathcal{D}[L_j]\rho \,,$$

with $$\mathcal{D}[L_j]\rho = L_j \rho L_j^\dagger - \tfrac12\{L_j^\dagger L_j, \rho\}$$. The specific structure of $H$ and the dissipators $L_j$ determines the phase behavior.

Key illustrative models include:

• Driven-dissipative Kerr Resonators: Single or two-photon driven cavities (e.g., with Hamiltonian $$H = -\Delta a^\dagger a + (U/2)a^{\dagger2} a^2 + F(a^\dagger + a)$$), displaying S-shaped multistability curves, optical bistability, and first-order DPTs (Casteels et al., 2016, Beaulieu et al., 2023).
• Spin Chains and Quantum Ising Models: 1D or 2D Ising chains with transverse field and local quantum baths can display continuous-to-first-order transition tuning as a function of dissipation strength or bath quantum nature (Meglio et al., 2020, 2207.13782, Jin et al., 2018).
• Fully Connected p-Spin Models: The canonical mean-field limit, where first-order transitions arise due to multistability in semiclassical equations; discontinuities persist in the thermodynamic limit (Wang et al., 2020, Nava et al., 2019).
• Central Spin Systems: Collective dissipative transitions, e.g., in quantum dots or NV centers coupled to nuclear spins, with discontinuous jumps in nuclear polarization and optical pumping observables (Kessler et al., 2012, Raghunandan et al., 2017).
• Rabi and Dicke-Ising Models: The interplay of drive, dissipation, and light-matter coupling in cavity QED systems can stabilize bistable nonequilibrium steady states, giving rise to first-order DPTs with multicritical and tetracritical points (Lyu et al., 2023, Li et al., 2023, Wang et al., 10 Feb 2026).

3. Dynamics and Critical Behavior

At the first-order dissipative critical point, the competition between unitary dynamics and dissipation produces unique scaling phenomena:

• Scaling of Dissipation with Energy Gap: For example, in the quantum Ising chain at a first-order transition, the energy gap $\Delta_L$ between nearly degenerate ground states closes exponentially with system size. Only when the Lindblad dissipation strength $u$ is tuned as $u \sim \Delta_L \sim e^{-cL}$ is a universal dynamical scaling regime accessed (Meglio et al., 2020).
• Universal Scaling Functions: The time evolution of observables under matched scaling variables (rescaled time $\theta = \Delta_L t$, dissipation $\gamma = u/\Delta_L$) collapses onto universal scaling functions, with nontrivial dependence on the global vs. local nature of the dissipation (Meglio et al., 2020).
• Metastable Timescales: The switching time between metastable branches scales as $\tau_\mathrm{switch} \sim e^{+cL}$, with critical slowing down observed as the Liouvillian gap closes (Casteels et al., 2016, Beaulieu et al., 2023, Minganti et al., 2023).
• Experimental Observables: Hysteresis area versus sweep time, critical exponents for phase boundaries and multicritical points, spectral signatures (e.g., bimodal Wigner/Husimi functions), discontinuity in correlation length and density (Röhrle et al., 2023, Lyu et al., 2023, Dagvadorj et al., 2021).

4. Physical Mechanisms and Universality Classes

First-order dissipative transitions originate from the interplay between multistability of the underlying coherent dynamics and structure of the dissipation:

• Bistability Mechanism: The system admits multiple locally stable steady states for the same parameters, with the ultimate steady state determined by stochastic switching or initial conditions when bistability is present (Casteels et al., 2016, Raghunandan et al., 2017, Röhrle et al., 2023).
• Microscopic Origin: In spin models, dissipation can convert continuous transitions to first order by inducing effective symmetry breaking fields through quantum bath fluctuations (2207.13782). In the p-spin models, the structure of the semiclassical dynamical equations determines whether first-order transitions are accessible depending on the order $p$ and symmetry.
• Universality: The character (first vs. second order) depends on system symmetries, spatial dimensionality, connectivity, and the competition of dissipation and interactions. For instance, the first-order transition appears in 2D driven polariton systems but not in 1D, due to the divergence (or lack thereof) of the nucleation barrier with system size (Li et al., 2021).

5. Mathematical and Spectral Signatures

The Liouvillian superoperator formalism provides unambiguous spectral criteria for first-order dissipative transitions:

• Liouvillian Gap: The spectral gap $\Delta = -\mathrm{Re}(\lambda_1)$, where $\lambda_1$ is the subleading eigenvalue of the Liouvillian, closes only at the critical point and does so exponentially in system size for first-order transitions (Casteels et al., 2016, Minganti et al., 2018).
• Metastability and Mixture Steady State: At the critical point, for $N\to\infty$, the steady state becomes an equal mixture of two distinct pure steady states: $\rho_{\rm ss}(\zeta_c) = \frac12 (\rho^+ + \rho^-)$ (Minganti et al., 2018).
• Switching Rates and Arrhenius Barrier: The rate of noise-induced switching between coexisting steady states is exponentially suppressed as the free energy barrier height increases, $\tau \sim e^{\beta \Delta f}$ (Nava et al., 2019).
• Observables: Discontinuity and bimodality in order parameters, ensemble-averaged expectation values, and dynamical correlators provide direct experimental access to the transition (Röhrle et al., 2023, Dagvadorj et al., 2021).

6. Experimental Realizations and Applications

First-order dissipative phase transitions have been realized or are experimentally feasible in a diverse range of platforms:

• Superconducting Circuits: Two-photon driven Kerr resonators exhibit clear signatures including hysteresis, phase coexistence, and critical slowing down, enabling robust bosonic qubits (“Kerr cats”) (Beaulieu et al., 2023).
• Quantum Optical Systems: Semiconductor microcavities, optical bistability in photonic lattices, polariton condensates, and arrays of nitrogen-vacancy centers in diamond employ interaction-induced first-order DPTs as a resource for enhanced quantum sensing, logic elements, and robust state preparation (Dagvadorj et al., 2021, Raghunandan et al., 2017, Li et al., 2021).
• Spin Qubits and Cold Atoms: Engineered dissipative environments in Rydberg atom arrays, trapped ions, and open central-spin systems provide platforms for studying both fundamental critical behavior and applications in quantum metrology (Kessler et al., 2012, Yan et al., 2024, Lyu et al., 2023).

Effective quantum simulation of first-order dissipative phase transitions is enabled by the ability to tune system size, coupling strengths, and dissipation rates to access the critical scaling regime (e.g., $u \sim \Delta_L$), which is often within reach in systems of modest size ($L \sim 10-20$) (Meglio et al., 2020). Criticality-enhanced sensitivity and the robustness of first-order transitions against disorder and decoherence are key features exploited in sensing platforms (Raghunandan et al., 2017).

7. Broader Context and Outlook

First-order dissipative phase transitions fundamentally extend the concept of phase transitions beyond equilibrium thermodynamics, bringing in new universality classes, critical phenomena, and technological opportunities:

• Nonequilibrium Effects: The universality and spectral theory extend to models involving strong symmetries ($Z_n$), higher-order nonlinearities, and even non-Hermitian dynamical generators, leading to multicritical and mixed-order transitions, as in open Rabi and Dicke-Ising models (Lyu et al., 2023, Wang et al., 10 Feb 2026, Yan et al., 2024).
• Connection to Classical Kinetics: Analogies to classical nucleation and first-order kinetics emerge in the scaling of nucleation barriers and the dimensional dependence of the transition (Li et al., 2021).
• Nontrivial Dynamics and Out-of-Equilibrium Phenomena: Rich dynamical signatures such as dynamical hysteresis, metastability, Mpemba-like memory effects, and global-to-local phase separation have been characterized in the context of Lindblad-driven fully-connected and locally-coupled systems (Nava et al., 2019).
• Extensions to Cosmology and Field Theory: Dissipative effects during cosmological first-order transitions imprint on gravitational wave spectra, demonstrating the universality of dissipative first-order mechanisms even in early-universe contexts (Guo, 2023).

The ongoing development of first-order dissipative phase transition theory integrates quantum optics, many-body physics, nonequilibrium statistical mechanics, and quantum information, with important implications for experimental quantum simulation and critical sensing technologies.