Bootstrapping Open Quantum Many-body Systems with Absorbing Phase Transitions

This paper introduces a systematic bootstrap method that leverages the positivity of density matrices and steady-state conditions to analyze open quantum many-body systems on infinite lattices, successfully deriving bounds on critical parameters and spectral gaps for the quantum contact process exhibiting absorbing phase transitions.

The Gist

Imagine you are trying to understand a massive, chaotic city where people are constantly arriving, leaving, and changing their minds. This city represents an open quantum system—a world of tiny particles that are constantly interacting with their environment (like the air or heat), causing them to lose energy and information.

In the past, physicists had a hard time predicting how this city would behave in the long run because the rules were messy and the math was incredibly difficult. They couldn't just "solve" the equations to find the answer.

This paper introduces a new, clever detective method called Bootstrapping to figure out the city's secrets without needing to solve the whole puzzle at once.

The Problem: The "Absorbing" City

The specific city the authors studied is called the Quantum Contact Process. Imagine a grid of light switches (the particles).

• The "Off" State (Absorbing State): If all switches are off, they stay off forever. Nothing can turn them on because the "off" state is a trap. Once everyone is asleep, the party is over.
• The "On" State (Active Phase): If you have enough energy (a parameter called $\Omega$), the switches can flip on and off, keeping the party going.

The big question is: How much energy is needed to keep the party going? There is a "critical point" (a tipping point). Below it, the city goes silent. Above it, the city buzzes with activity.

The Solution: The "Bootstrap" Method

Usually, to predict the future of a system, you need to know the exact starting point and run the simulation forward. But for infinite quantum systems, this is impossible.

Instead, the authors use a method called Bootstrapping. Think of it like this:

Imagine you are trying to guess the height of a giant, invisible mountain. You can't see the top, but you know a few rules:

1. The mountain must be solid (it can't be made of negative rock).
2. The mountain must be stable (it doesn't change shape over time).
3. The mountain must look the same if you walk a few steps to the left or right (symmetry).

You don't need to see the whole mountain. You just take a small piece of it, apply those rules, and ask: "What is the tallest the top of this piece could possibly be?" and "What is the shortest it could be?"

By tightening these rules and looking at slightly larger pieces, you get a "box" that the mountain's height must fit inside. Even if you don't know the exact height, you can say, "It's definitely between 1,000 and 2,000 feet."

In this paper, the "rules" are:

• Positivity: The probability of finding a particle in a certain state can't be negative (you can't have -50% chance of being alive).
• Steady State: The system isn't changing anymore; it's settled.

What Did They Find?

Using this "mathematical squeeze," the authors managed to trap the answers to three big questions:

1. The Tipping Point: They calculated a strict lower limit for the energy needed to keep the party going. They proved that the critical energy is at least 2.85. (Previous guesses were around 6, so this is a huge improvement in precision, even if it's not the final number yet).
2. The "Ghost" State: In the active phase, there's a special "ghost" state that represents the difference between the noisy party and the silent trap. They figured out exactly how much the "noise" in the system can vary compared to the "silence."
3. The Speed of Silence: In the silent phase, if you start with a noisy system, how fast does it die down? They calculated the "spectral gap," which is basically the speed limit of how fast the system forgets its past and settles into silence.

Why Is This a Big Deal?

Think of it like trying to navigate a foggy forest.

• Old methods: You try to map the whole forest at once. If the fog is too thick, you get lost.
• This new method: You don't need the whole map. You just need to know that the trees are solid and the ground is flat. By checking small patches of the forest, you can draw a fence around the impossible areas. You know you can't walk there, so you must be walking here.

The authors showed that you can build a rigorous "fence" around the behavior of these complex quantum systems without ever needing to solve the impossible, infinite equations.

The Catch (and the Future)

The method is powerful, but it's computationally expensive. It's like trying to solve a Sudoku puzzle where the grid gets bigger every time you add a new rule. The authors are currently limited to looking at small "patches" of the city (subsystems of size 8).

To get the exact answer, they need to look at larger patches. But they've proven the method works. It's a new flashlight in a very dark room, and it's showing us that we can find rigorous answers to problems that were previously thought to be unsolvable.

In short: They built a mathematical "safety net" that catches the possible behaviors of quantum systems, proving exactly where the line is between a dead, silent system and a lively, active one.

Technical Summary

1. Problem Statement

The paper addresses the challenge of rigorously analyzing open quantum many-body systems governed by Lindblad master equations on infinite lattices. Unlike closed systems described by Hermitian Hamiltonians, open systems undergo non-unitary evolution due to environmental coupling, leading to quantum dissipation.

Specifically, the authors focus on systems exhibiting absorbing phase transitions, a universality class characterized by:

• Absorbing State: A pure state (e.g., the "all-off" state) that is a steady state for all parameter values.
• Phase Transition: A transition between an absorbing phase (where the absorbing state is the unique steady state) and an active phase (where non-trivial, genuinely mixed steady states exist).
• The Challenge: Analytical solutions for these transitions are rare, even in 1D. Existing numerical methods (like Tensor Networks or TEBD) often lack rigorous error bounds or struggle with the infinite lattice limit. The absence of Hermitian generators and standard symmetries (Landau paradigm) makes deriving rigorous theoretical results difficult.

2. Methodology: The Bootstrap Framework

The authors adapt the quantum bootstrap method, traditionally used for closed systems, to open quantum systems. The core idea is to constrain the space of possible physical observables using the fundamental mathematical properties of the system, rather than solving the dynamics directly.

Core Constraints

The method relies on two defining properties of steady states $\rho$ in Lindblad dynamics:

1. Positivity: The density matrix must be positive semidefinite ($\rho \succeq 0$).
2. Stationarity: The state must be invariant under the Lindblad super-operator ($\mathcal{L}(\rho) = 0$).

Implementation Steps

1. Truncation Scheme: The infinite lattice is approximated by considering local observables supported on a finite subsystem of size $N$.
2. Moment Matrices: The authors construct a matrix $M$ of expectation values $\langle O_i^\dagger O_j \rangle$. The positivity of the density matrix implies $M \succeq 0$.
3. Linear Constraints:
   • Translation Invariance: $\langle T(O) \rangle = \langle O \rangle$.
   • Equations of Motion: $\langle \hat{\mathcal{L}}(O) \rangle = 0$ for all local operators $O$.
4. Optimization: The problem is formulated as a Semidefinite Programming (SDP) task. The goal is to minimize or maximize a target observable (e.g., magnetization) subject to the positivity and stationarity constraints.

Key Technical Innovations

• Reduced Density Matrix (RDM) Basis: To improve computational efficiency, the authors utilize a basis corresponding to reduced density matrices. This reduces the dimension of the positivity matrix from $4^N$ to $2^N$, allowing for larger subsystem sizes ($N$) to be tested.
• Handling the Active Phase: In the supercritical phase, the system has a non-trivial steady state $\rho_1$. The authors define a deviation matrix $D = \rho_1 - \rho_0$ (where $\rho_0$ is the absorbing state). Since $\rho_0$ is a steady state, $D$ satisfies the same linear constraints but with $\langle 1 \rangle_D = 0$. They bootstrap the ratio of deviations $r(O) = \delta\langle O \rangle / \delta\langle Z_1 \rangle$ to handle the scale invariance of $D$.
• Spectral Gap Bootstrapping: For the absorbing phase, the authors analyze the late-time relaxation behavior $\langle O \rangle(t) \sim e^{-\Delta t}$. They formulate a feasibility problem to find the Liouvillian spectral gap $\Delta$ by checking if a moment matrix constructed from the decay modes is positive semidefinite.

3. Key Contributions

• Systematic Framework: Established the first systematic bootstrap framework specifically for open quantum systems with absorbing phase transitions on infinite lattices.
• Rigorous Bounds: Provided rigorous, non-perturbative lower and upper bounds on physical quantities without relying on finite-size scaling extrapolations or uncontrolled approximations.
• Extension to Non-Equilibrium: Successfully extended the bootstrap philosophy (combining positivity and dynamical constraints) from equilibrium/stochastic systems to quantum Lindblad dynamics.
• New Constraints for Active Phases: Developed specific constraints for the "deviation" matrix $D$ to probe the non-trivial steady state in the active phase.

4. Key Results

The authors applied their method to the Quantum Contact Process (QCP), a 1D chain of qubits with Hamiltonian coupling $\Omega$ and dissipation rate $\gamma=1$.

• Critical Coupling ($\Omega^*$):
  • By bootstrapping the steady-state magnetization $\langle Z_1 \rangle$, they derived a rigorous lower bound for the critical coupling.
  • At subsystem size $N=8$, they found $\Omega^*_{LB} \approx 2.85$.
  • While this is lower than previous estimates ($\Omega^* \approx 6$) from TEBD, the bounds have not yet converged, indicating that higher $N$ will tighten the estimate.
• Active Phase Observables:
  • They obtained bounds on the ratio of deviations $r(Z_1 Z_2) = \delta\langle Z_1 Z_2 \rangle / \delta\langle Z_1 \rangle$.
  • The results confirmed the consistency of the method, as the critical point derived from the active phase constraints matched the steady-state bounds.
• Liouvillian Spectral Gap ($\Delta$):
  • In the absorbing phase, they computed bounds on the spectral gap $\Delta$ (the rate of relaxation to the absorbing state).
  • At $\Omega=2.0$, for $N=8$, the gap was bounded between $0.187$ and $0.629$.
  • These bounds are consistent with, though currently looser than, estimates from Matrix Product States (MPS) and TEBD ($\Delta \approx 0.316$).
• Convergence Behavior: The results show that bounds tighten as $N$ increases, but convergence is slow. The authors note that the number of variables grows exponentially, posing a computational bottleneck.

5. Significance and Outlook

• Rigorous Validation: This work provides a mathematically rigorous alternative to numerical simulations for open quantum systems, offering guaranteed error bounds rather than heuristic estimates.
• Universality Class Insights: It offers a new tool to study the Directed Percolation (DP) universality class in the quantum regime, a problem that has resisted analytical solution.
• Future Directions:
  • Coarse-Graining: The authors suggest that standard coarse-graining (effective for 1D closed systems) may not work for open systems due to potential volume-law entanglement in steady states. Identifying the "most relevant" subset of constraints is a key open question.
  • Quantum Circuits: The framework is naturally extendable to discrete-time open systems (quantum circuits).
  • Information-Theoretic Quantities: Future work aims to bootstrap quantities like entropy and free energy directly, moving beyond local operator expectation values.

In summary, the paper demonstrates that combining density matrix positivity with steady-state conditions creates a powerful, systematic tool for deriving rigorous bounds in open quantum many-body physics, successfully applied to the critical behavior of the Quantum Contact Process.