Skip to content
Quantum.world

Research

Classical Cellular Automaton for Measurement-Only Entanglement Transitions

Overview Research area: Quantum physics — monitored/measurement-only quantum systems, entanglement phase transitions, and classical stochastic cellular automata. Technical level: Intermediate. The pap

Classical Cellular Automaton for Measurement-Only Entanglement Transitions
arXiv
2609.28743
Published
2026-09-23
Authors
Will Holdhusen, Bailey M. McAmis, Armin Rahmani

AI summary

Overview

Research area: Quantum physics — monitored/measurement-only quantum systems, entanglement phase transitions, and classical stochastic cellular automata. Technical level: Intermediate. The paper is written for readers comfortable with qubit chains, Bell states, entanglement entropy, and stabilizer circuits, though the core idea (entanglement counted as pairings of sites) is easy to state. Scope: The paper introduces a one-dimensional measurement-only qubit model with competing single-site and power-law long-range Bell measurements, maps it exactly onto a classical stochastic matching automaton, and derives the volume-law, fractal, and area-law entanglement phases analytically and numerically.

What This Paper Is About

In ordinary monitored quantum circuits, unitary gates build entanglement while measurements destroy it. The authors ask what happens when measurements are the only ingredient — no unitary evolution at all — and when the Bell measurements are long-ranged, with pair separations drawn from a power-law distribution. Their goal is to turn this quantum dynamics into a purely classical pairing process that can be solved analytically in some limits and simulated far more cheaply than the quantum circuit itself.

Key Contributions

  1. An exact classical mapping. Every trajectory of the quantum model remains a tensor product of Bell pairs and unpaired qubits, so the full entanglement structure is encoded in a classical matching of the sites. This matching evolves as a long-range generalized stochastic cellular automaton with pairing moves.
  2. Two analytically solvable limits. The authors derive steady-state entanglement exactly in the all-to-all limit (α = 0) and in the dilute-pair limit (p → 1⁻, with p = 1 − ε and ε ≪ 1).
  3. A complete phase diagram. They identify volume-law, fractal (sub-volume-law), and area-law entanglement regimes separated by transitions at α* = 1 and α_c = 2, and confirm numerically that all three regimes persist for finite p and intermediate α.
  4. A postselection-free probe. Because Bell measurement outcomes only relabel Bell states and never change the matching, all subsystem entropies are outcome independent — removing the exponential postselection cost that complicates experimental studies of measurement-induced entanglement.

Main Findings

  • Elementary automaton moves: A Bell measurement on an unpaired qubit plus one member of a pair replaces the bond (b, c) with (a, b), leaving c unpaired. A Bell measurement on one qubit from each of two pairs performs entanglement swapping, (a, c), (b, d) → (a, b), (c, d). The pairing update is identical for all four Bell outcomes.
  • Entropy equals crossing bonds: With entropy defined as S_A(N_A) = −Tr(ρ_A log₂ ρ_A), each Bell pair connecting subsystem A to its complement contributes exactly one unit (using natural logarithms, ln 2).
  • All-to-all limit (α = 0): The unpaired-qubit density obeys dn₁/dt = (2/N)[p(1 − n₁) − (1 − p)n₁²], with steady state n₊ = (−p + √(4p − 3p²))/(2(1 − p)) and relaxation time τ = N/(2√(4p − 3p²)). At p = 0 the density decays algebraically as n₁(t) = 1/(1 + 2t/N).
  • Volume law at α = 0: The entropy integrates to S_A(N_A, t) = N n_A(1 − n_A)(1 − n₁(t)), giving volume-law scaling for every 0 ≤ p < 1 except the trivial p = 1 case.
  • Dilute-pair limit: The steady-state pair-length distribution P(r) approaches the input distribution P_in(r). Simulations at p = 0.95 and α = 0.5 (400 qubits, 10⁵ trials) confirm this, while p = 0.05 shows a clear enhancement of the long tail from entanglement swapping.
  • Dilute-limit scaling: Half-chain entropy follows S_A(N/2) ∼ N^β, with β = 1 for α < 1 (volume law), β = 2 − α for 1 < α < 2 (fractal), and β = 0 for α > 2 (area law).
  • Boundary scaling: The transitions sit at α* = 1 and α_c = 2, with S ∼ N/log N at α* and S ∼ log N at α_c.
  • Numerical persistence: Heatmap exponents were extracted from finite-size fits over system sizes from 400 to 12800 qubits, averaged over 10⁴ trials, and show the dilute-limit regimes surviving at finite p. Phase boundaries depart appreciably from dilute-limit values only for p ≲ 0.1.
  • Nearest-neighbor connection: In the limit α → ∞, the p = 0 dynamics are related to Majorana loop models with logarithmic entanglement scaling.
  • Efficient simulation: Each update of the matching costs 𝒪(1) time, and adequate convergence was reached with 10³ to 10⁶ realizations.

Methodology in Plain English

The authors place N qubits on an open one-dimensional chain, all starting in the state |00⋯0⟩. At each step they flip a weighted coin: with probability p they measure a single site in the computational basis, and with probability 1 − p they perform a Bell measurement on a pair of sites whose separation r is chosen with weight proportional to r^(−α). For numerical work they sample separations from P_in(r) ∝ (N − r)r^(−α) using a Walker–Vose alias table, which accounts for the N − r available pairs at each length.

The key simplification is that this dynamics never creates anything more complicated than Bell pairs. So instead of tracking the exponentially large quantum state, they track a vector of N integers recording each site's partner (or marking it unpaired), updating it in constant time. Entanglement entropy for a subsystem is then just a count of pairs crossing the cut. To reach steady states, simulations for p > 0 were run for T = 200N timesteps, well beyond the relaxation time in Eq. (5). To locate the transitions, they fit the entropy to the scaling form expected at α = 1 (A·N/log N + B·N/(log N)²) or α = 2 (A log N + B) on a grid of α values with spacing Δα = 0.01, choosing the value that maximizes R².

Why This Matters

The work supplies a rare measurement-only model whose entanglement transition can be handled both analytically and numerically, and it does so without the postselection bottleneck that has blocked experimental tests of measurement-induced entanglement. It also gives a concrete geometric picture — measurements as pairing moves on a classical matching — for how entanglement is generated and spread.

  • Experimental probing of monitored entanglement: Since entropies depend only on which sites were measured, not on the Bell outcomes, no exponentially unlikely measurement record needs to be repeated.
  • Efficient classical benchmarking: The automaton replaces a stabilizer circuit with an 𝒪(1)-per-update matching simulation, useful for validating hardware and algorithms.
  • Quantum teleportation and Bell-measurement protocols: Bell measurements are central to teleportation, and the paper discusses realizing them via ancilla-assisted stabilizer measurements in superconducting circuits or entangling gates plus local readout in trapped-ion teleportation experiments.
  • Theory of long-range entanglement structure: The model provides a controlled setting connecting power-law measurement range to entanglement scaling, including fractal sub-volume-law phases.

Industry relevance: The paper does not report industrial deployments or commercial benchmarks; its relevance is to quantum-hardware groups needing postselection-free diagnostics and to teams building efficient classical simulators for measurement-based protocols.

Future Directions

  • Higher dimensions and other geometries: The construction is stated to extend naturally to higher dimensions, other lattice geometries, and irregular spatial arrangements, all unexplored here.
  • Other pair-selection distributions: The authors propose exploring distributions beyond the r^(−α) power law, as well as spatially disordered measurement rates.
  • The unresolved p = 0 line: Except for the α = 0 and α → ∞ limits, the asymptotic entanglement scaling along p = 0 remains open; apparent area-law behavior in finite-time simulations may be an artifact of slow relaxation, and the Supplemental Material suggests logarithmic or small-exponent fractal scaling at longer times. Longer times and larger sizes are needed.
  • Experimental realization: Turning the protocol into an actual hardware demonstration using the ancilla-assisted stabilizer or trapped-ion routes the authors cite.

Target Audience

Condensed-matter and quantum-information theorists working on measurement-induced phases and monitored circuits; researchers interested in classical stochastic and cellular-automaton mappings of quantum dynamics; and experimentalists seeking postselection-free signatures of measurement-generated entanglement. Graduate students with a background in entanglement entropy and Bell states will find the model accessible, and readers focused on numerical methods will benefit from the 𝒪(1)-update matching algorithm.

Authors’ abstract

We introduce a generalized classical long-range stochastic cellular automaton that exactly captures the entanglement dynamics and transitions of a measurement-only monitored quantum system. The corresponding quantum model consists of a one-dimensional qubit chain subject to competing single-site and long-range Bell measurements, with separations drawn from a power-law distribution. At every step of each trajectory, the state remains a tensor product of Bell pairs and unpaired qubits, with bipartite entanglement entropy fully encoded in a classical matching of the sites. This exact representation allows us to solve the dynamics analytically in several limits and to determine the steady-state scaling of the entanglement entropy. We find volume-law, fractal, and area-law regimes. Numerical results show that these regimes persist beyond the analytically solvable limit. Because Bell measurement outcomes change only the Bell-state labels and leave the matching unchanged, all subsystem entropies are outcome independent, avoiding the exponential cost of trajectory postselection often encountered in experimental studies of measurement-induced entanglement dynamics.

Read the original paper