Quantum predator-prey cycles in dissipative Rydberg array

Technical Summary: Quantum Predator-Prey Cycles in Dissipative Rydberg Arrays

Problem Statement
The predator-prey cycle is a paradigmatic example of self-organized population oscillations in far-from-equilibrium systems, classically described by deterministic models like the Lotka–Volterra equations. While spatial stochastic lattice models have revealed rich phenomena such as pattern formation and noise-induced transitions, empirical validation in natural ecosystems remains challenging due to the inability to isolate and control microscopic parameters (e.g., reproduction rates, interaction ranges, and intrinsic noise). Furthermore, while classical Lotka–Volterra dynamics have been observed in thermal Rydberg vapors, the role of genuine quantum effects—specifically fluctuations and nonperturbative coherence—in shaping such dynamics remains an open frontier. There is a need for a clean, tunable platform to explore these dynamics in the quantum regime, particularly to understand how interaction ranges and quantum fluctuations stabilize oscillations.

Methodology
The authors propose a quantum analogue of predator-prey dynamics using a tunable two-dimensional (2D) array of Rydberg atoms. The system is modeled as a spin-1 dissipative Rydberg gas where a ground state $|g\rangle$ is coupled to two Rydberg states, $|r\rangle$ (prey) and $|s\rangle$ (predator), via laser fields with Rabi frequencies $\Omega_{r,s}$ and detunings $\Delta_{r,s}$. The dynamics are governed by a quantum master equation incorporating coherent Hamiltonian evolution (including strong, long-range Rydberg-Rydberg interactions) and dissipative decay to the ground state.

To analyze the system, the authors employ a multi-faceted approach:

1. Mean-Field (MF) Analysis: They assume a factorized density matrix to derive equations of motion for expectation values. Stability is assessed via the eigenvalues of the Jacobian matrix to distinguish between stationary states, limit cycles (LC), and quasicycles.
2. Open-System Discrete Truncated Wigner Approximation (OSDTWA): To capture quantum noise and go beyond MF theory, they utilize OSDTWA, which combines the semiclassical truncated Wigner framework with the method of quantum jumps. This allows for large-scale numerical simulations of the stochastic spin-1 system, which is computationally intractable via exact diagonalization due to the exponential growth of the Liouville-space dimension ($O(9^N)$).
3. Correlation Analysis: The study employs two-time auto-correlation and cross-correlation functions, as well as Fourier spectra, to characterize oscillation frequencies, phase shifts, and coherence lifetimes.

Key Results

• Emergence of Quantum Limit Cycles: The system exhibits stable predator-prey cycles where the populations of $|r\rangle$ and $|s\rangle$ oscillate out-of-phase. This dynamics is driven by asymmetric effective interactions induced by a relative detuning ($\Delta_r - \Delta_s$). The cycle evolves through three phases: prey accumulation, predator growth, and mutual suppression.
• Role of Quantum Coherence: The authors demonstrate that these limit cycles and the associated time-translation symmetry breaking are intrinsically quantum phenomena. A rate-equation analysis (which neglects quantum coherence) yields only stable fixed points, indicating that the oscillatory behavior resides in the full eight-dimensional phase space driven by nonperturbative quantum coherence.
• Quasicycles and Noise Amplification: In regimes where the MF fixed point is linearly stable, the system exhibits "quasicycles"—slow oscillatory modes driven by the amplification of intrinsic quantum noise from quantum jumps. These quasicycles possess an intrinsic frequency distinct from classical Gaussian-noise-induced oscillations.
• System Size and Synchronization:
  • Small Systems: Under finite-range van der Waals (vdW) interactions, the system maintains robust global synchronization.
  • Large Systems: As system size increases, local quantum jumps induce significant desynchronization, suppressing global oscillations in the average populations. However, the predator-prey relationship persists locally. The system fragments into spatio-temporal clusters, each sustaining local cycles.
  • Scaling: The amplitude of the quasicycles scales inversely with the square root of the system size ($\delta n \propto 1/\sqrt{N}$), confirming their origin in finite-size noise.
• Interaction Range Effects: The study contrasts short-range vdW interactions with all-to-all coupling. All-to-all coupling stabilizes global oscillations against quantum jumps, whereas finite-range interactions lead to the loss of global coherence while preserving local synchronization.

Significance and Claims
The paper claims to extend the study of predator-prey models into the quantum realm, offering a fully quantum theoretical investigation of these dynamics in the coherent regime of Rydberg arrays. Key contributions include:

• Demonstration of Quantum-Driven Symmetry Breaking: Showing that nonperturbative quantum coherence is indispensable for breaking time-translation symmetry in this context, a mechanism distinct from stochastic contact processes in classical systems.
• Identification of Local vs. Global Dynamics: Revealing that in large systems with finite-range interactions, global observables may appear as a stochastic steady state, while spatially resolved observables and correlation functions uncover hidden, locally synchronized predator-prey cycles.
• Experimental Feasibility: Positioning two-component Rydberg atomic lattices as a highly controllable platform for exploring self-organizing phenomena in nonequilibrium systems. The authors note that the predicted microsecond timescales and specific parameter dependencies (laser detuning, Rabi frequencies) are resolvable in current Rydberg atom array experiments.
• Distinction from Dissipative Time Crystals: The work distinguishes these cycles from dissipative time crystals by noting that while both involve spontaneous time-translation symmetry breaking, the latter typically exhibits global phase coherence, whereas the finite-range predator-prey cycles here exhibit only local synchronization due to quantum-jump-induced desynchronization.

The authors conclude that their work advances quantum simulation strategies by leveraging engineered many-body nonequilibrium effects to probe mechanisms that are difficult to isolate in natural ecosystems.