Switching Dynamics of Metastable Open Quantum Systems

1. Problem Statement

The paper addresses a fundamental paradox in driven-dissipative quantum many-body systems:

• The Paradox: Markovian open quantum systems governed by the Lindblad equation are theoretically expected to relax to a unique stationary state as time approaches infinity. However, experimental and theoretical studies of systems like Rydberg atom ensembles exhibit collective quantum jumps (stochastic switching) between distinct metastable states (e.g., "bright" and "dark" states) that persist indefinitely.
• The Conflict: There is a conceptual gap between two facets of quantum metastability:
  1. Spectrum-level metastability: Characterized by a small Liouvillian spectral gap (slow relaxation of the density matrix).
  2. Trajectory-level metastability: Characterized by stochastic switching between disjoint basins of attraction in quantum trajectories.
• The Question: How do these two levels relate in finite systems? Specifically, does the small spectral gap always imply slow relaxation, and how do quantum fluctuations at zero temperature drive switching between metastable states in a manner analogous to thermal activation in classical systems?

2. Methodology

The authors employ a multi-faceted approach combining analytical theory, spectral analysis, and numerical simulations on a model of $N$ interacting Rydberg atoms:

• Physical Model: A system of $N$ two-level atoms driven by a laser (Rabi frequency $\Omega$, detuning $\Delta$) with long-range interactions (approximated as all-to-all coupling $V$) and dissipative decay (rate $\gamma$).
• Mean-Field (MF) Analysis: Derivation of coupled nonlinear equations for collective spin variables to identify fixed points and phase diagrams (monostable vs. bistable regimes).
• Spectral Decomposition:
  • Analysis of the Liouvillian superoperator eigenvalues ($\lambda_l$) to identify the spectral gap ($-\text{Re}[\lambda_1]$).
  • Decomposition of the stationary state into disjoint metastable subspaces ($\hat{\rho}_+$ and $\hat{\rho}_-$) using the trace-preserving symmetry of the Lindbladian.
• Large Deviation (LD) Theory:
  • Calculation of the scaled cumulant generating function (SCGF) for photon emission rates to detect dynamical phase transitions (kinks in the LD function).
• Semiclassical Instanton Approach:
  • Mapping the quantum dynamics to a classical stochastic process via quantum Langevin equations.
  • Using the Martin-Siggia-Rose (MSR) path integral formalism to derive the effective action.
  • Calculating optimal transition paths (instantons) and the quasipotential (effective energy barrier $\Phi$) to estimate switching rates.
• Quantum-Jump Monte Carlo Simulations (QJMC):
  • Direct simulation of individual quantum trajectories to measure waiting times between collective jumps and verify theoretical predictions for various system sizes ($N$) and detunings ($\Delta$).

3. Key Contributions

1. Distinction of Metastability Levels: The paper rigorously distinguishes between spectrum-level metastability (a small Liouvillian gap leading to slow relaxation only if specific initial conditions are met) and trajectory-level metastability (stochastic switching that erases memory of initial conditions and limits relaxation).
2. Quantum Arrhenius Law: The authors establish that switching rates between quantum metastable states follow an Arrhenius-like law: $T^{-1} \propto e^{-\Phi N}$, where $\Phi$ is the effective energy barrier and the inverse system size ($1/N$) acts as the nonequilibrium analog of temperature.
3. Connection to Quasipotential: They extend the concept of the quasipotential functional (typically used in classical non-equilibrium statistical mechanics) to the quantum realm, linking the probability of rare fluctuations to the steady-state occupation probabilities of the metastable states.
4. Resolution of the "Unique Stationary State" Paradox: The paper demonstrates that the unique stationary state in finite bistable systems is a statistical mixture of two disjoint metastable states, sustained by exponentially rare collective quantum jumps.

4. Key Results

• Spectral Gap Scaling: In the bistable regime, the Liouvillian spectral gap decreases exponentially with system size ($N$), consistent with the exponential scaling of switching times. In monostable regimes, the gap remains finite or scales differently.
• Disjoint Subspaces: The two low-lying eigenmodes of the Liouvillian correspond to two disjoint density matrices ($\hat{\rho}_+$ and $\hat{\rho}_-$), which map directly to the "bright" (high excitation) and "dark" (low excitation) states predicted by Mean-Field theory.
• Switching Dynamics:
  • Waiting Times: Quantum-jump simulations confirm that the mean waiting time for collective jumps scales exponentially with $N$ ($T \sim e^{\Phi N}$).
  • Barrier Dependence: The effective energy barriers ($\Phi_b, \Phi_d$) for escaping the bright and dark states depend on the detuning $\Delta$. As $\Delta$ increases, the barrier for the dark state increases, making it more stable, while the barrier for the bright state decreases.
  • Initial Condition Independence: Unlike systems with only a small spectral gap (where relaxation speed depends on initial state overlap), the presence of stochastic switching ensures that the system eventually loses memory of its initial state, and the relaxation time is strictly limited by the switching rate.
• Validation: The effective energy barriers calculated via the instanton approach (semiclassical path integral) match those extracted from quantum-jump simulations and spectral decomposition, validating the semiclassical approximation for these quantum systems.

5. Significance

• Theoretical Insight: The work provides a unified framework connecting classical large deviation theory with quantum open system dynamics. It clarifies that quantum bistability is not merely a finite-size artifact but a robust phenomenon driven by quantum fluctuations acting as "noise" at zero temperature.
• Experimental Relevance: The findings are directly applicable to current platforms like Rydberg atom arrays and circuit QED, where collective quantum jumps are observed. It offers a method to characterize phase transitions and metastability without requiring full diagonalization of the Liouvillian (which is computationally prohibitive for large $N$).
• Practical Applications:
  • Quantum Metrology & Information: Understanding the switching dynamics is crucial for designing robust quantum memories or sensors that utilize bistability.
  • Stochastic Resonance: The exponential size scaling suggests that system size can be tuned to control switching frequencies, potentially enhancing signal detection via stochastic resonance.
  • Distinguishing Bistability vs. Metastability: The paper provides a clear criterion to distinguish true bistability (where both states are occupied in the steady state) from mere metastability (where the system is trapped in one state), which is vital for interpreting experimental data in finite-size quantum simulators.

In summary, the paper elucidates the mechanism of quantum switching in dissipative systems, demonstrating that quantum fluctuations induce an Arrhenius-type activation process where the system size plays the role of inverse temperature, thereby bridging the gap between spectral properties and trajectory-level dynamics.