Quantum to classical relaxation dynamics of the dissipative Rydberg gas

Using the truncated Wigner approximation to study large two-dimensional systems, this paper demonstrates that quantum kinetically constrained dynamics and pronounced relaxation slowdowns persist in dissipative Rydberg gases even in the weakly dissipative regime where coherent and dissipative processes compete.

The Gist

The Big Picture: A Dance of Atoms

Imagine a crowded dance floor filled with atoms. These aren't just ordinary atoms; they are Rydberg atoms, which are like atoms that have been stretched out and are wearing giant, fluffy hats. Because of these hats, they are very sensitive to each other. If one atom puts on its hat (gets excited), it creates a "personal space bubble" (called a blockade radius) where no other atom nearby is allowed to put on a hat at the same time.

The scientists in this paper wanted to watch what happens when these atoms try to relax from a chaotic state back to a calm, resting state. But there's a twist: the atoms are dancing in two different worlds at once:

1. The Quantum World: Where atoms can be in two places at once, dance in perfect sync, and do weird, wave-like things.
2. The Classical (Noisy) World: Where atoms are just regular particles bumping into things, losing energy, and getting confused by the environment (like wind or noise).

The paper asks: What happens when the "Quantum Dance" and the "Noisy Wind" are fighting for control?

The Problem: The "Glass" Trap

In the past, scientists knew that if the "Noisy Wind" (dissipation) was very strong, the atoms would get stuck in a glassy state.

• The Analogy: Imagine trying to move through a crowded room where everyone is holding hands. If you try to move, you can't because your neighbors are blocking you. You get stuck in a jam.
• In physics terms, this is called kinetic constraints. The atoms want to relax, but the rules of the game (the Rydberg blockade) say, "You can't move unless your neighbor moves first," but your neighbor is also stuck. This leads to a very slow, sluggish relaxation, like honey dripping.

The big mystery was: What happens if the "Quantum Dance" is strong? Does the atoms' ability to be "quantum" (coherent) help them break out of this jam, or does it make the jam even weirder?

The Tool: The "Truncated Wigner Approximation" (TWA)

To study this, the scientists needed to simulate millions of atoms.

• The Problem: Simulating quantum atoms is like trying to calculate the path of every single grain of sand in a beach storm. The math gets so huge that even the world's fastest supercomputers crash.
• The Solution: They used a method called TWA.
• The Analogy: Instead of tracking every single grain of sand perfectly, imagine taking a photo of the beach and then creating thousands of slightly different "what-if" versions of that photo. You run a simulation for each version, and then you take the average.
  • This method treats the atoms like little classical balls (easy to calculate) but adds a little bit of "quantum fuzziness" (noise) to the start to make it feel real. It's a clever shortcut that lets them study huge systems (like a 2D grid of 150x150 atoms) that were previously impossible to simulate.

The Experiments: Two Starting Points

The scientists started the atoms in two different "dance formations" to see how they relaxed:

1. The "All Asleep" Start (Fully Polarized State)

• The Setup: Every atom starts in the ground state (no hats).
• What Happened:
  • At first, the atoms start putting on hats (getting excited) because of a laser.
  • But as soon as a few get hats, the "personal space bubbles" kick in. The atoms get stuck in a traffic jam.
  • The Result: The system hits a plateau. It stops relaxing for a while. It's like a traffic light turning red; cars (atoms) are stuck waiting for a gap to open up.
  • The Surprise: In 2D (a flat grid), this jam was even deeper and lasted longer than in 1D (a line). The atoms were so constrained by their neighbors that they couldn't move at all for a long time.

2. The "Chess Board" Start (Néel State / Quantum Scars)

• The Setup: The atoms start in a perfect checkerboard pattern (Up, Down, Up, Down). This is a special "Quantum Scar" state, meaning it's a very organized pattern that usually oscillates (wiggles back and forth) for a long time without settling down.
• What Happened:
  • Instead of hitting a flat plateau, the magnetization (the average "hat-ness") dipped down and then started wiggling (oscillating).
  • The Result: The quantum nature kept the atoms "alive" and dancing for a while, fighting against the noise trying to stop them. Eventually, the noise won, and they settled down, but the journey was full of swings and oscillations rather than a simple stop.

The Key Findings

1. The "Glass" is Real: Even with quantum effects, the atoms still get stuck in kinetic jams. The "blockade" creates a rule where you can't move unless your neighbor moves, leading to slow relaxation.
2. Quantum vs. Noise:
   • If the noise is loud (strong dissipation), the atoms behave like classical particles and get stuck in a slow, glassy jam.
   • If the noise is quiet (weak dissipation), the quantum "dance" takes over. The atoms oscillate and wiggle before finally getting stuck.
3. The 2D Difference: In two dimensions, the jam is much more complex. The atoms aren't just blocked by their immediate neighbors; the "blockade" ripples out further, creating a more complicated traffic jam than in a simple line.

The Conclusion: Why This Matters

This paper is like a map for understanding how complex systems (from quantum computers to materials) relax when they are noisy.

• For Quantum Computers: If you want to build a quantum computer, you need to know how long your atoms can stay "quantum" before the noise ruins the party. This paper shows that even with noise, there are specific patterns (like the checkerboard) that can survive for a while, but eventually, the "traffic jams" (kinetic constraints) will slow everything down.
• The Takeaway: Nature has a way of getting stuck. Whether it's atoms in a lab or people in a crowded room, if everyone is waiting for someone else to move first, the whole system slows down to a crawl. The scientists found that even when the atoms are "quantum," they still get stuck in these jams, but the path to getting stuck looks very different depending on how noisy the room is.

In short: The scientists used a clever math trick to watch a huge crowd of quantum atoms. They found that when the atoms try to relax, they often get stuck in a "traffic jam" caused by their own rules. Sometimes they wiggle and dance before getting stuck, and sometimes they just freeze immediately, depending on how much "noise" is in the room.

Technical Summary

1. Problem Statement

The paper addresses the relaxation dynamics of open quantum many-body systems, specifically Rydberg gases, where coherent driving and dissipation compete.

• The Context: Rydberg atoms exhibit strong, tunable long-range interactions leading to the Rydberg blockade, which prevents simultaneous excitations within a characteristic radius. This imposes kinetic constraints on the system's dynamics.
• The Gap:
  • In the strongly dissipative limit (dephasing rate $\gamma \gg$ Rabi frequency $\Omega$), quantum coherences are suppressed, and the system is well-described by effective classical rate equations exhibiting glassy, kinetically constrained relaxation.
  • In the weakly dissipative limit ($\gamma \sim \Omega$), quantum coherences become relevant. While small 1D systems have been studied, the fate of kinetic constraints in large systems and higher dimensions (2D) remains unexplored due to the exponential growth of the Hilbert space and the computational cost of tensor-network methods for open, long-range interacting systems.
• The Core Question: How do kinetic constraints induced by the Rydberg blockade manifest in the relaxation dynamics when coherent and dissipative processes are comparable in large 1D and 2D systems?

2. Methodology

To overcome the limitations of exact diagonalization and tensor networks for large 2D systems, the authors employ the Truncated Wigner Approximation (TWA).

• Model: A system of $N$ atoms on a $d=1$ or $d=2$ lattice, modeled as two-level systems (ground $|\downarrow\rangle$, Rydberg $|\uparrow\rangle$).
  • Hamiltonian: Includes a transverse field drive ($\Omega$) and long-range van der Waals interactions ($V/r^6$).
  • Dissipation: Modeled via local dephasing (jump operators $\hat{L}_i = \sqrt{\gamma}\hat{n}_i$).
  • Master Equation: The evolution is governed by a Lindblad master equation.
• Truncated Wigner Approximation (TWA):
  • Maps the quantum master equation onto a stochastic evolution of classical phase-space variables (spin vectors $\mathbf{s}_i$).
  • Initial State Sampling: Uses a discrete sampling scheme (DTWA) to reproduce the first and second moments of quantum operators for specific initial states:
    1. Fully Polarized State: All spins down ($|\Psi_\downarrow\rangle$).
    2. Néel State: Alternating up/down spins ($|\Psi_{N\acute{e}el}\rangle$), which belongs to a manifold of quantum many-body scars.
  • Dynamics: The classical spins evolve according to stochastic differential equations (Stratonovich interpretation) incorporating deterministic precession (from $H$) and stochastic noise (from dissipation).
  • Observables: Expectation values (e.g., magnetization $m_z$) and correlation functions are obtained by averaging over an ensemble of classical trajectories.
• Benchmarking: The TWA results are validated against:
  • Exact quantum-jump Monte Carlo simulations (for small systems, $L=10, A=32$).
  • Kinetic Monte Carlo (kMC) based on classical rate equations (valid only in the strong dissipation limit).

3. Key Contributions

1. Scalability: Successfully applied TWA to study relaxation dynamics in large 1D ($L=100$) and 2D ($A=152$) Rydberg arrays, a regime inaccessible to exact quantum simulations.
2. Regime Mapping: Systematically explored the crossover from coherent-dominated to dissipation-dominated dynamics, identifying how kinetic constraints evolve as the ratio $\gamma/\Omega$ changes.
3. Initial State Dependence: Contrasted the relaxation behavior of a trivial initial state (polarized) against a non-trivial quantum scar state (Néel), revealing distinct dynamical signatures.
4. Dimensionality Effects: Provided the first systematic analysis of kinetic constraints in 2D Rydberg gases under weak dissipation, comparing them to 1D counterparts.

4. Key Results

A. Benchmarking and Validity

• In the strongly dissipative limit ($\gamma/\Omega = 10$), TWA, exact dynamics, and classical kMC agree perfectly, confirming the validity of the classical rate equation description.
• In the weakly dissipative limit ($\gamma/\Omega = 0.1, 0.5$), TWA captures coherent oscillations and qualitative features missed by kMC. While TWA slightly underestimates the amplitude of intermediate oscillations, it correctly predicts the long-time relaxation behavior and the emergence of kinetic constraints.

B. Relaxation Dynamics: Fully Polarized Initial State

• Emergence of Plateaus: As interaction strength $V/\Omega$ increases, the magnetization relaxation exhibits a transient plateau at intermediate times.
  • 1D: The plateau occurs at $m_z \approx -0.45$, corresponding to a "hard-dimer" constraint (nearest-neighbor exclusion).
  • 2D: A plateau emerges at $m_z \approx -0.62$. This value is lower than the theoretical prediction for a simple "hard-square" model ($m_z \approx -0.54$), suggesting that coherent dynamics in 2D impose effective constraints extending beyond nearest neighbors.
• Timescales:
  • Plateau formation time ($\tau_{plateau}$): Scales linearly with $\gamma$ ($\tau \sim \gamma$). Lower dissipation accelerates the approach to the plateau via coherent dynamics.
  • Final relaxation time ($\tau_{final}$): Scales inversely with $\gamma$ ($\tau \sim 1/\gamma$). The approach to the stationary state is governed by dissipative spin flips.
• Correlations: The Mandel $Q$-parameter shows the buildup of strong anticorrelations (negative $Q$) during the plateau, which decay slowly as the system approaches the uncorrelated stationary state.

C. Relaxation Dynamics: Néel Initial State (Quantum Scars)

• Oscillations vs. Plateaus: Unlike the polarized state, the Néel state does not form a plateau. Instead, it exhibits pronounced coherent oscillations at early times due to the quantum scar manifold.
• Intermediate Minimum: In the weakly dissipative regime, the magnetization develops a deep intermediate minimum before relaxing.
• Persistence of Constraints: Even after the initial Néel order is lost, the system remains in a transient regime with strong anticorrelations ($Q \approx -0.5$) for a long duration, indicating that blockade-induced constraints continue to restrict the accessible configuration space.

D. Spatial and Temporal Correlations

• Spatial: Anticorrelations are a hallmark of the kinetic constraints. In 2D, these correlations persist longer and are more complex than in 1D.
• Temporal: Two-time autocorrelation functions reveal that as dissipation decreases, the system retains memory of its initial configuration for longer times, with relaxation rates evolving dynamically (waiting-time dependence).

5. Significance and Outlook

• Theoretical Impact: The work demonstrates that kinetic constraints are robust even when quantum coherence is significant. It reveals that in 2D, coherent dynamics can generate effective constraints that are more restrictive than simple nearest-neighbor exclusion models predict.
• Methodological Advance: It validates the TWA as a powerful tool for studying open quantum systems in the "intermediate" regime where neither pure quantum nor pure classical descriptions are sufficient, particularly for 2D geometries.
• Experimental Relevance: The predicted slow, constraint-dominated relaxation and transient plateaus are observable in current Rydberg array experiments. However, the authors note that realizing the long-time dynamics requires platforms with extended coherence times, such as proposed Rydberg-ion systems, to avoid heating and loss processes that obscure these effects.
• Future Directions: The paper suggests that while TWA captures the onset of constraints, more advanced methods (e.g., cluster TWA, BBGKY hierarchies, or tensor networks) may be needed to fully describe the constraint-dominated regime in 2D at very strong interactions.