⚛️ quantum physics

Quantum simulation of gauge theories on dynamical spacetimes via Floquet-induced matrix models

This paper introduces a Floquet-based quantum simulation framework that utilizes large- $N$ matrix models and randomized benchmarking to efficiently simulate gauge theories on dynamical spacetimes with exponentially fewer qubits than conventional lattice methods, thereby preserving continuous symmetries and enabling the study of phenomena like deconfinement on expanding cosmological backgrounds.

Original authors: Samuel Buckley-Bonanno, Noah Eckstein, Susanne F. Yelin

Published 2026-07-07

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Samuel Buckley-Bonanno, Noah Eckstein, Susanne F. Yelin

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

The Big Problem: Trying to Map a Shapeshifting World with a Grid

Imagine you are trying to draw a map of a city. If the city is static (buildings never move), you can use a standard grid paper. You draw squares for blocks and lines for streets. This works great for a normal city.

However, imagine a city where the ground itself is stretching, shrinking, and warping like a rubber sheet (this is what happens in the universe during cosmic expansion or near black holes). If you try to draw this on a fixed grid, the grid breaks. You'd have to constantly rip up squares and add new ones, which ruins the rules of the map. In physics, this is the problem with simulating gauge theories (the rules governing forces like electromagnetism and the strong nuclear force) on dynamical spacetimes (moving, changing universes). Traditional computer simulations use a fixed "lattice" (grid), and they struggle when the universe itself changes shape.

The New Solution: A Jigsaw Puzzle of Moving Parts

The authors propose a completely different way to build the map. Instead of a grid of points, they use a set of matrices (think of them as complex, multi-layered spreadsheets of numbers).

In this new system:

Space isn't a grid of dots; it's defined by how these spreadsheets interact with each other.
The "shape" of the universe is determined by the mathematical "commutation" (the order in which you multiply these spreadsheets). If you multiply Spreadsheet A then B, you get a different result than B then A. This difference is the curvature of space.
No fixed grid: Because the geometry comes from the math of the numbers themselves, the "universe" can stretch or shrink without needing to add or remove grid squares. It's like a hologram that changes shape without breaking the projector.

The Magic Trick: The "Floquet" Dance

The paper's biggest breakthrough is figuring out how to make a quantum computer actually calculate with these matrices.

Normally, to simulate a system, you have to track every single particle. But here, the authors use a technique called Floquet engineering.

The Analogy: Imagine you want to know how heavy a box is, but you can't lift it. Instead, you shake the box back and forth in a specific rhythm. The way the box wobbles tells you its weight.
The Method: The researchers take a quantum computer and apply a rapid, repeating sequence of "shakes" (operations) to a random set of quantum states. They call this a Floquet sequence.
The Result: By measuring how much the quantum state "wobbles" or changes after this dance (a measurement called fidelity or a Loschmidt echo), they can directly read out the "weight" of that specific universe configuration.

Essentially, they found a way to turn the complex math of the universe's shape into a simple "yes/no" or "strong/weak" signal that a quantum computer can measure easily.

Why This is a Game-Changer

Fewer Qubits: Traditional methods require a massive number of quantum bits (qubits) to represent space. This new method uses exponentially fewer. It's like needing a single high-resolution camera instead of a million low-resolution ones to take a picture.
Handling the "Moving" Universe: Because they aren't using a fixed grid, they can simulate universes that are expanding (like our own Big Bang cosmology) without breaking the laws of physics (unitarity).
Real Results: The team didn't just write theory; they ran simulations. They successfully modeled a specific type of force field (SU(2) gauge field) on two types of backgrounds:
- A flat, boring space.
- An expanding, curved space (like a cosmological model).
- They observed a "phase transition" (a sudden change in the state of the matter), proving the system works even when the background is stretching.

The Bottom Line

The paper introduces a new "language" for quantum computers to speak about the universe. Instead of forcing the universe into a rigid grid (which breaks when the universe moves), they use a flexible system of interacting numbers (matrices). By making these numbers "dance" in a specific rhythm, they can extract the physics of forces and curved space directly, opening the door to simulating the dynamic, expanding universe in ways that were previously impossible.

Technical Summary: Quantum Simulation of Gauge Theories on Dynamical Spacetimes via Floquet-Induced Matrix Models

Problem Statement
The quantum simulation of gauge theories faces a fundamental conflict between continuous spacetime symmetries and the discrete resources of conventional computation. Standard lattice gauge theory resolves this by fixing a spatial grid, sacrificing exact translational and rotational symmetry for computational tractability. While successful for static backgrounds, this approach fails for dynamical spacetimes (e.g., cosmological expansion or gravitational horizons) because adding or removing lattice sites violates unitarity. Furthermore, lattice discretization is inherently a classical treatment of position, even when on-site states are quantum. Existing proposals often retain the lattice framework, limiting their ability to simulate curved or time-dependent geometries covariantly.

Methodology
The authors propose a framework based on large- $N$ Yang-Mills matrix models, where spacetime geometry and gauge fields are encoded in the commutation structure of a set of $N \times N$ Hermitian matrices $\{X_a\}$ , rather than on a lattice of points. In this formalism:

Geometry: The classical continuum emerges smoothly as $N \to \infty$ . The geometry is determined by the dynamical Lie algebra generated by $\{X_a\}$ under repeated commutation. Derivatives are replaced by commutators (e.g., $[X_a, \cdot]$ ), and the Laplacian is defined by nested commutators.
Fields: Scalar fields are introduced as additional matrices, while Yang-Mills fields arise as fluctuations around block-diagonal matrix backgrounds with internal unitary symmetry.
Floquet Protocol: The core innovation is a Floquet-induced simulation of the Euclidean path integral. The authors show that the Euclidean path-integral weights ( $e^{-S[X]}$ $e^{- S [X]}$ ) of a matrix model can be reproduced by the ensemble-averaged fidelities of Haar-random states evolved under periodic sequences of matrix operators.
- A four-step symmetric sequence $F_2(X_a, X_b; \Delta t)$ is constructed, which, via the Baker-Campbell-Hausdorff formula, isolates the commutator $i[X_a, X_b]$ at leading order in the time step $\Delta t$ .
- When applied to a Haar-random ensemble, the average fidelity (Loschmidt echo) scales as $\exp(-\text{tr}(H_{\text{eff}}^2))$ , effectively reproducing the Boltzmann factor $e^{-S[X]}$ where the action $S[X]$ is built from squared commutator terms.
Quantum Encoding: The $N \times N$ matrices are encoded as Hamiltonian operators acting on $2 \log_2 N$ qubits. This requires exponentially fewer qubits than canonical quantization approaches, which typically require $\Lambda N^2 d$ qubits for truncated Fock spaces.
Sampling: To evaluate operator product expansions, the protocol uses controlled Floquet sequences on entangled states. A register of $n$ control qubits activates $2^n$ distinct matrix states in parallel. Post-selecting on the target qubit remaining stationary yields a probability distribution proportional to the Euclidean weights $e^{-S[X]}$ , allowing for the sampling of the path integral measure.
Observables: Local field observables are extracted via quasi-coherent states $|\vec{x}\rangle$ , which are ground states of quadratic Hamiltonians localized around classical coordinates. These states allow the reconstruction of local field expectation values and Wilson loops, trading lattice discreteness for quantum uncertainty in spatial information.

Key Contributions

Floquet-Path Integral Correspondence: The paper establishes that the Euclidean path-integral weights of large- $N$ matrix models are directly accessible through the average fidelities of Haar-random states under periodic Floquet driving. This provides a "quantum-native" method to sample the partition function without explicit state preparation of the Boltzmann distribution.
Exponential Qubit Efficiency: By treating matrices as operators rather than states, the protocol reduces the qubit requirement from polynomial in $N$ (canonical quantization) to logarithmic ( $2 \log_2 N$ ).
Dynamical Spacetime Simulation: The framework naturally accommodates curved and expanding spacetimes (e.g., FLRW metrics) by adjusting the commutation structure of $\{X_a\}$ without adding or removing degrees of freedom, preserving unitarity and continuous symmetries.
Experimental Compatibility: The protocol relies on standard randomized benchmarking techniques (Loschmidt echo measurements) and does not require complex state preparation, making it directly implementable on programmable quantum platforms.

Results

Fidelity-Weight Validation: Numerical simulations confirm that the measured average fidelities of Haar-random states under Floquet evolution converge to the exact Euclidean weights $e^{-S[X]}$ as the Floquet frequency increases (high-frequency/weak-coupling limit).
Parallelized Sampling: The authors demonstrate a globally controlled quantum circuit that samples an ensemble of matrix states with probabilities proportional to their Euclidean weights, successfully reproducing the exponential weighting of the path integral.
Deconfinement Transition: The framework was used to simulate an $SU(2)$ gauge field on both flat Euclidean space and a negatively curved, expanding FLRW background. In both cases, the authors identified a deconfinement phase transition characterized by a drop in Wilson-loop expectation values as the coupling constant increases. The transition showed minimal hysteresis, consistent with a second-order phase transition.
Cosmological Application: The simulation on an expanding background demonstrates that field configurations at different points in cosmological expansion can be extracted by varying the time parameter used to generate quasi-coherent states, maintaining unitary evolution throughout.

Significance and Claims
The paper claims that this framework opens quantum simulation to field and spacetime dynamics that are beyond the reach of conventional lattice methods. By avoiding fixed spacetime lattices, the approach preserves continuous symmetries and unitarity on dynamical geometries.

Novelty in Cosmology: The authors state this is the first identification of a gauge-field deconfinement transition on an expanding cosmological background within a quantum-simulable encoding, a task previously inaccessible to lattice gauge theory due to its inability to represent time-dependent metrics covariantly.
Theoretical Implications: The work suggests that geometric structures may emerge spontaneously in periodically driven quantum systems, as the coherent error of any periodic drive takes the form of a Yang-Mills action.
Future Directions: The authors identify several extensions, including the introduction of fermionic fields (via Jordan-Wigner or physical platform statistics) to access supersymmetric models like IKKT, experimental encoding of the Lorentzian sector via non-Hermitian operations, and the derivation of matrix states for arbitrary metrics (e.g., Schwarzschild). They posit that this offers a new route to probing the relationship between operator-algebraic encodings of bulk geometry and entanglement-based tensor network encodings.

The paper concludes that while the formalism is complex, its language is deeply compatible with quantum information tools, offering a new pathway to understand high-energy physics and cosmology through experimental quantum platforms.