⚛️ quantum physics

Steady-state phase transition in one-dimensional quantum contact process

This paper investigates the steady-state phase transition of the one-dimensional quantum contact process by combining mean-field approximations with linked-cluster expansions to reveal a discontinuous saddle-node bifurcation and a non-divergent correlation length, offering predictions testable in Rydberg atom quantum simulators.

Original authors: Lin Shang, Shuai Geng, Xingli Li, Jiasen Jin

Published 2026-02-06

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Lin Shang, Shuai Geng, Xingli Li, Jiasen Jin

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a long line of light switches, where each switch can be either OFF (empty) or ON (occupied). This is the setup for the "Quantum Contact Process" described in the paper.

In the real world, if you have a disease spreading, a rumor, or a fire, it usually needs a neighbor to spread to you. If you are healthy (OFF), you only get sick (ON) if someone next to you is already sick. If you are sick, you might recover on your own.

This paper studies what happens when these switches are quantum (they can be in a fuzzy state of being both ON and OFF at the same time) and when they are constantly being "reset" by a noisy environment (dissipation).

Here is the breakdown of their discovery using simple analogies:

1. The Two States: The "Dead Zone" vs. The "Party"

The researchers are looking for a tipping point where the system changes from a "Dead Zone" (where every switch is OFF and stays OFF) to a "Party" (where a steady number of switches stay ON).

The Absorbing Phase (Dead Zone): If the "spreading" force is too weak, the system eventually dies out completely. Everyone is OFF.
The Active Phase (Party): If the "spreading" force is strong enough, the system finds a way to keep a steady number of switches ON forever.

2. The "Ghost" Problem (Metastability)

When the researchers tried to simulate this on a computer, they ran into a tricky ghost.
Imagine you are trying to push a heavy boulder up a hill to get it to the top (the "Party" state).

The Trap: There is a small dip or a "ghost valley" just before the top. If you push the boulder, it might get stuck in this valley for a very long time. It looks like it's at the top, but it's actually stuck in a temporary holding pattern.
The Paper's Finding: They discovered that near the transition point, the system gets stuck in this "ghost valley" (metastable state) for a long time before finally falling back down to the "Dead Zone."
The Lesson: If you run a simulation for too short a time, you might think the system is stable when it's actually just pretending to be. You have to wait much longer to see the truth.

3. The "Magic Mirror" Solution

To fix the "ghost trap" problem, the authors invented a new way of looking at the system, which they call a self-consistent effective field.

The Analogy: Imagine trying to find the perfect balance point on a seesaw. Usually, you just sit on it and see where it settles. But if the seesaw is wobbly, it might get stuck in a weird spot.
The New Trick: Instead of just sitting on the seesaw, they built a "magic mirror" that tells the seesaw exactly what it should be doing based on what its neighbors are doing. This forces the system to ignore the "ghost valleys" and jump straight to the real, stable solutions.
The Result: This method allowed them to see the true shape of the transition clearly, without the system getting confused by the temporary traps.

4. The "Cliff" vs. The "Ramp" (Discontinuous Transition)

This is the most important discovery.

The Old Belief: Many scientists thought that as you turn up the "spreading" knob, the system would slowly, smoothly slide from "Dead" to "Party," like walking up a gentle ramp.
The Paper's Discovery: They found it's not a ramp; it's a cliff.
- As you turn the knob, nothing happens. The system stays dead.
- Suddenly, at a specific point, the system jumps instantly from "all OFF" to "some ON."
- This is called a saddle-node bifurcation. It's like a light switch that doesn't dim; it just snaps on.

5. The "Chain Reaction" Check (Correlations)

To make sure their "magic mirror" wasn't lying, they used a method called Linked-Cluster Expansion.

The Analogy: Imagine trying to predict the weather by looking at a single tree, then a small group of trees, then a whole forest.
The Finding: They checked if the "Party" state was caused by long-distance connections (like a signal traveling from one end of the line to the other). They found that the system mostly relies on neighbors (short-range connections).
The Proof: The "susceptibility" (how easily the system reacts to a nudge) didn't explode or go crazy at the transition point. If it were a smooth, continuous transition (like a ramp), the reaction would have gone to infinity. Since it didn't, it confirms the "Cliff" theory: the transition is sudden and discontinuous.

Summary

The paper argues that in this specific one-dimensional quantum system:

The system can get "stuck" in a fake state for a long time (metastability).
By using a clever new calculation method, they avoided this trap.
They proved that the system doesn't slowly wake up; it snaps from a dead state to an active state like a light switch flipping on.
This behavior is driven by local neighbors, not long-distance signals.

The authors suggest this could be tested in real life using Rydberg atoms (a type of atom used in quantum simulators), which act like these quantum switches.

Problem Statement
The paper investigates the steady-state phase transitions of the one-dimensional quantum contact process (QCP), a model describing non-equilibrium quantum many-body systems. While the classical contact process is a well-established model for phenomena like disease spreading and information dissemination, its quantum extension introduces complex interplay between coherent time evolution and dissipation. A central controversy in the field concerns the nature of the phase transition in the one-dimensional QCP: specifically, whether the transition from the absorbing phase (empty sites) to the active phase (occupied sites) is continuous or discontinuous. Previous studies have yielded conflicting results, with some suggesting continuous transitions and others hinting at discontinuities, leaving the order of the transition in one-dimensional systems unresolved.

Methodology
The authors employ a multi-faceted theoretical approach to resolve the nature of the transition:

Liouvillian Spectrum Analysis: The system's time evolution is described by a Lindblad master equation. The authors analyze the Liouvillian gap (the magnitude of the maximal negative real part of the Liouvillian eigenvalues) in the thermodynamic limit to identify phase transitions and metastability.
Mean-Field Approximations:
- Single-Site Mean-Field (MF): Used as a baseline to identify steady-state solutions and bifurcation behaviors.
- Cluster Mean-Field (CMF): To incorporate short-range correlations, the authors utilize the CMF method. A key innovation is the introduction of a novel self-consistent condition based on an effective field. This condition parameterizes the inter-cluster interaction and allows the authors to capture all possible steady-state solutions, including unstable ones, while preventing the system from being trapped in metastable states during numerical time evolution.
Numerical Linked-Cluster Expansion (LCE): To investigate long-range correlations and the behavior of the magnetic susceptibility ( $\chi$ ) in the thermodynamic limit, the authors implement a numerical LCE. This method sums contributions from clusters of increasing size to extrapolate results for an infinite system.

Key Results

Metastability and Liouvillian Gap: The analysis of the Liouvillian spectrum reveals a closing of the gap at $\Omega/\Gamma \approx 5.83$ in the thermodynamic limit. The spectrum structure indicates the existence of a metastable state near the transition point, characterized by a long-lived plateau in the time evolution of the population before decaying to the true absorbing state. This explains why previous simulations might have failed to observe the true steady state if the simulation time was insufficient.
Discontinuous Phase Transition: Both the single-site MF and the CMF analyses reveal a saddle-node bifurcation of the order parameter (steady-state population $\langle \bar{n} \rangle_{ss}$ ). As the control parameter $\Omega$ varies, the system exhibits a coexistence of absorbing and active states, with the active branch emerging via a saddle-node point ( $S_1$ ) rather than a supercritical pitchfork bifurcation. This behavior is characteristic of a discontinuous (first-order) phase transition.
Transition Point Extraction: Using the CMF method with the novel self-consistent condition, the authors extract the transition point for infinite-size systems. By extrapolating the cluster-size dependence of the transition point $\Omega_c(L)$ , they determine $\Omega_c/\Gamma \approx 6.04$ , which aligns well with the Liouvillian gap prediction of $\approx 5.83$ .
Correlation Length and Susceptibility: The LCE results show that in the absorbing phase, cluster weights decay exponentially, indicating short-range correlations. In the active phase, the decay follows a power law. Crucially, the extrapolated magnetic susceptibility $\chi$ exhibits a monotonic decrease across the transition point without diverging. This lack of divergence in susceptibility and correlation length further supports the conclusion that the transition is discontinuous, distinguishing it from continuous transitions where such quantities typically diverge.

Significance and Claims
The paper claims to resolve the controversy regarding the nature of the phase transition in the one-dimensional QCP by demonstrating that it is a discontinuous transition from the absorbing to the active phase. The authors emphasize that the observation of this discontinuity relies on correctly handling the system's metastability, which can trap numerical simulations in false steady states. The proposed self-consistent condition based on effective fields is presented as a robust tool to bypass these metastable interferences and capture the full landscape of steady-state solutions.

The results suggest that the steady-state properties are dominated by short-range correlations, validating the use of cluster mean-field approximations for this system. The authors note that these theoretical findings could be tested in quantum simulators using Rydberg atoms. They conclude by identifying the study of hybrid contact processes (combining quantum and classical branching/coagulation) and the extraction of critical exponents in such systems as promising future directions.