Continuum limit of a qubit-regularized SU(3) lattice gauge theory with glueballs

This paper demonstrates that a qubit-regularized SU(3) lattice gauge theory on a plaquette chain admits a continuum limit governed by a $\mathbb{Z}_3$ parafermion conformal field theory, yielding a massive one-dimensional model of glueballs with specific mass ratios and string tension values.

The Gist

Imagine the universe is built from tiny, invisible Lego bricks. Most of us know about the bricks that make up atoms (protons and neutrons), but there's a special "glue" that holds them together. In the world of physics, this glue is made of particles called gluons.

Here's the weird part: unlike the glue in your garage, which just sticks things together and stays put, these gluons are sticky to each other. They can clump together to form their own little balls of pure energy, even without any atoms attached. Physicists call these mysterious clumps glueballs.

The problem? Glueballs are incredibly hard to study. They are heavy, short-lived, and the math required to predict how they behave is so complex that even the world's fastest supercomputers struggle with it.

The New Approach: A "Toy" Universe

In this paper, a team of physicists (Rui Xian Siew, Shailesh Chandrasekharan, and Tanmoy Bhattacharya) decided to stop trying to simulate the whole, messy universe at once. Instead, they built a miniature, simplified model—a "toy universe"—to see if they could find glueballs there.

Think of it like this: If you want to understand how a real car engine works, you don't need to build a full-size Ferrari. You can build a tiny, working model out of wood and wire. If the little model revs up and makes the right sounds, you know you've captured the essential mechanics.

Their "toy universe" is a qubit-regularized lattice gauge theory.

• Qubits: These are the basic units of information for quantum computers (like 0s and 1s, but with a twist). The researchers designed their model so it could theoretically run on a future quantum computer.
• The "Plaquette Chain": Imagine a ladder. The sides of the ladder are the "links," and the rungs are the "plaquettes." They arranged their model as a long chain of these ladders. It's a one-dimensional strip, which is much easier to solve than the full 3D space we live in.

The Magic Trick: Tuning the Dial

The researchers didn't just build the model; they had to find the "sweet spot" where the model behaves like real physics.

1. The Critical Point: They adjusted the settings (couplings) of their model until it hit a "phase transition." Imagine a pot of water. If you turn the heat up slowly, it stays liquid. But at exactly 100°C, it hits a critical point and turns into steam. The researchers found a similar critical point in their math.
2. The "UV" and "IR":
   • UV (Ultraviolet): This represents the very small, high-energy scale. At this scale, their model behaves like a specific, beautiful mathematical pattern known as a Conformal Field Theory (think of it as a perfectly symmetrical, frictionless dance of particles).
   • IR (Infrared): This represents the large, low-energy scale (our everyday world). By adding a tiny "push" (a magnetic perturbation) to the system, they forced the model to break that perfect symmetry.

The Result: Glueballs Appear!

When they pushed the system into the "IR" regime, something magical happened. The particles in their toy model stopped behaving like massless waves and suddenly gained mass. They became heavy, stable particles.

Because the model was designed to mimic the rules of the strong nuclear force (SU(3) gauge theory), these new heavy particles are the toy versions of glueballs.

What Did They Measure?

Once they had these "toy glueballs," they measured two things to see if their model was realistic:

1. The Mass Ratio: They looked at two types of glueballs: one that is "even" (symmetric) and one that is "odd" (anti-symmetric) under a mirror test called "charge conjugation." They calculated the ratio of their masses.

   • The Result: The "odd" glueball is about 1.46 times heavier than the "even" one. This matches what we expect from more complex, traditional calculations.

2. The String Tension: In the real world, if you try to pull a quark (a particle inside a proton) away from another, the "glue" between them acts like a rubber band. The harder you pull, the tighter it gets. The energy required to stretch this band is called "string tension."

   • The Result: They calculated the ratio of this string tension to the glueball mass. The number they got (0.2648) is a precise fingerprint of their theory.

Why Does This Matter?

This paper is a breakthrough for three main reasons:

• Proof of Concept: It proves that you can build a simple, "qubit-friendly" model of the strong force that actually produces the right kind of physics (massive glueballs) when you zoom out.
• Quantum Computing Ready: Because the model is built on qubits, it's a perfect candidate for running on future quantum computers. This could eventually allow us to simulate the strong force in ways that are impossible for today's supercomputers.
• A New Map: It shows that even in a simplified, one-dimensional world, the complex dance of confinement (how quarks are trapped) and the emergence of mass can be understood. It suggests that the same "critical points" might exist in our full 3D universe, guiding us toward a deeper understanding of why the universe has mass at all.

In short: They built a tiny, simplified Lego version of the strong force, tuned it just right, and watched as "glueballs" spontaneously formed. It's a small step for a toy model, but a giant leap for understanding how the universe holds itself together.

Technical Summary

1. Problem Statement

The determination of the glueball spectrum (massive bound states of gluons) in SU(3) Yang-Mills theory is a fundamental challenge in non-perturbative quantum field theory. While standard lattice gauge theory (LGT) approaches discretize the theory on a lattice, they often suffer from infinite-dimensional local Hilbert spaces, making them difficult to simulate on quantum computers.

Qubit regularization aims to reformulate LGT using finite-dimensional local Hilbert spaces (qubits) to enable quantum simulation. However, a major open question is whether such simplified models can reproduce the correct continuum physics, specifically:

• Do they possess a relativistic quantum critical point?
• Can they flow to the asymptotically free fixed point of SU(3) Yang-Mills theory?
• Do they generate massive glueball excitations in the infrared (IR) without quarks?

This paper addresses these questions by investigating a specific qubit-regularized SU(3) lattice gauge theory on a plaquette chain (a quasi-one-dimensional geometry).

2. Methodology

Model Construction

The authors construct a model on a "plaquette chain," a ladder-like geometry where elementary plaquettes are formed by links along the chain ($\ell_c$) and rungs ($\ell_r$).

• Hilbert Space: The physical Hilbert space is constructed using the monomer-dimer tensor-network (MDTN) basis. Gauge degrees of freedom reside on links and are labeled by irreducible representations (irreps) of SU(3): $|1\rangle$ (singlet), $|3\rangle$ (triplet), and $|\bar{3}\rangle$ (anti-triplet).
• Hamiltonian: The system is governed by a local, gauge-invariant Hamiltonian ($H_{pc}$):
  $$H_{pc} = \kappa_c \sum_{\ell_c} \hat{E}_{\ell_c} + \kappa_r \sum_{\ell_r} \hat{E}_{\ell_r} - g \sum_P (\hat{U}_P + \hat{U}_P^\dagger)$$
  • $\hat{E}$ terms assign energy penalties to non-singlet states ($|3\rangle, |\bar{3}\rangle$).
  • $\hat{U}_P$ terms are plaquette operators that cyclically permute irreps ($1 \to 3 \to \bar{3} \to 1$) while preserving gauge invariance.
  • The model avoids the sign problem by restricting $g \geq 0$.

Mapping to the Quantum Clock Model

The authors demonstrate that within each topological sector of the Hilbert space (specifically the vacuum sector $H_I$), the SU(3) plaquette chain is equivalent to a one-dimensional three-state quantum clock model.

• The Hamiltonian maps to:
  $$H_{qcm} = -\sum_P \left( J \hat{Z}_P \hat{Z}_{P+1}^\dagger + g \hat{X}_P + h \hat{Z}_P + \text{h.c.} \right)$$
  where $J = \kappa_r/3$ and $h = 2\kappa_c$.
• Critical Point: The model is critical at $J=g=1$ and $h=0$. This point belongs to the universality class of the 2D three-state Potts model, described by a $Z_3$ parafermion Conformal Field Theory (CFT) with central charge $c=4/5$.
• Continuum Limit: By tuning the magnetic field $h$ (equivalent to $\kappa_c$) to be small but non-zero, the system is driven away from the critical point into a massive phase. This $h$ acts as a relevant perturbation, generating a mass gap and flowing to a massive continuum Quantum Field Theory (QFT) in the IR.

Numerical Simulation

The authors employ Density Matrix Renormalization Group (DMRG) to compute the low-energy spectrum of the mapped quantum clock model.

• Scaling Variable: They analyze the system using the dimensionless scaling variable $\mu = h^{15/28}L$, where $L$ is the number of plaquettes.
• Extrapolation: Calculations are performed for lattice sizes $L = 32, 48, 64, 80, 96$ and various values of $\mu$. Results are extrapolated to the thermodynamic limit ($L \to \infty$) and the continuum limit ($h \to 0$).

3. Key Contributions

1. First Demonstration of Massive Continuum Limit: This is the first proof that a simple qubit-regularized SU(3) LGT admits a massive continuum limit with relativistic excitations.
2. Identification of UV Fixed Point: The work explicitly identifies the $Z_3$ parafermion CFT as the ultraviolet (UV) fixed point governing the short-distance physics of the qubit-regularized theory.
3. Relativistic Dispersion: The authors verify that the massive excitations in the IR obey a relativistic dispersion relation, confirming the emergence of a Lorentz-invariant QFT from the discrete lattice model.
4. Precise Non-Perturbative Ratios: The study computes precise ratios of physical observables that are independent of the lattice scale, providing benchmark data for future theoretical and experimental (quantum simulation) comparisons.

4. Key Results

Glueball Mass Ratio

The authors compute the ratio of the masses of the lightest stable particles with opposite charge conjugation ($C$-parity):

• $m_-$: Mass of the lightest $C$-odd state.
• $m_+$: Mass of the lightest $C$-even state.
• Result: In the IR limit, the ratio converges to:
  $$m_- / m_+ = 1.459(2)$$
  This result is consistent with previous numerical estimates using the Truncated Conformal Space Approach (TCSA) but provides higher precision.

String Tension and Mass Scale

The authors calculate the string tension $\sigma$ between a static quark and antiquark and compare it to the lightest glueball mass $m_+$:

• Result:
  $$\sqrt{\sigma} / m_+ = 0.2648(2)$$
  This ratio characterizes the confinement scale relative to the particle mass.

Relativistic Dispersion

By analyzing the energy spectrum of single-particle excitations with non-zero lattice momentum, the authors confirmed the dispersion relation:
$$(\epsilon_{j,n})^2 - (m_j)^2 = 2(1 - \cos p_n)\zeta^2$$
The extracted constant $\zeta = 2.5907(2)$ confirms the relativistic nature of the emergent particles.

5. Significance

• Toy Model for QCD: The plaquette chain serves as a minimal, solvable toy model for strong interactions without quarks. It demonstrates that complex continuum physics (confinement, massive glueballs) can emerge from simple qubit-based Hamiltonians.
• Quantum Computing Roadmap: The results validate the "qubit regularization" approach. Since the model is free of sign problems and uses finite-dimensional Hilbert spaces, it is a prime candidate for implementation on near-term quantum computers to simulate non-Abelian gauge theories.
• Universality and Confinement: The work links the confinement-deconfinement transition in the gauge theory to the ordering transition in the quantum clock model. It suggests that similar critical points and massive continuum limits likely exist in higher dimensions (3+1D), offering a pathway to study full SU(3) Yang-Mills theory via qubit regularization.
• Benchmarking: The precise numerical values for $m_-/m_+$ and $\sqrt{\sigma}/m_+$ provide essential benchmarks for testing future analytical methods and quantum simulations of gauge theories.

In conclusion, the paper successfully bridges the gap between abstract qubit-regularized models and physical continuum gauge theories, proving that massive glueball-like excitations can be generated and studied in a controlled, non-perturbative setting.