SU(3) Gauge Theory on Small Lattices

• The paper formulates SU(3) gauge theory with precise Schrödinger Functional boundary conditions and O(a) corrections to control finite-volume renormalization effects.
• It introduces minimal qutrit truncations and optimized quantum simulation circuits, efficiently capturing essential Yang–Mills dynamics with resource-aware strategies.
• Gradient flow techniques and analytic methods are applied to extract thermodynamic observables and hybrid potentials, advancing nonperturbative QCD insights.

SU(3) gauge theory on small lattices investigates the nonperturbative dynamics, renormalization, and quantum simulation properties of the pure Yang–Mills sector in QCD discretized on finite volumes. Small lattices—where “small” means low spatial and temporal extent in units of the lattice spacing—offer a controlled environment for studying boundary effects, strong-coupling truncations, quantum circuit construction, the equation of state at phase transitions, and precision extraction of potential-related observables. Fundamental developments include the formulation of renormalized couplings under finite-volume Schrödinger functional (SF) boundary conditions, the use of irreducible electric-field truncations (especially in “qutrit” bases), the application of gradient-flow methods for thermodynamic observables, and efficient resource allocation for quantum simulation. Modern approaches integrate explicit handling of O(a) and O(a²) cutoff effects, mapping between continuum and lattice effective field theories, and the transfer of gauge invariance to qudit or qubit architectures for NISQ devices.

1. Formulation and Boundary Conditions of SU(3) on Small Lattices

The standard approach for discretizing SU(3) gauge theory on a lattice employs the Wilson plaquette action, typically in four Euclidean dimensions. Small lattices require careful attention to boundary conditions and artefacts:

• Schrödinger Functional (SF) Setup: The SF framework defines a finite cylinder with spatial extent $L$ ($L/a\in 2\mathbb{N}$ even) and temporal extent $T$ ($T/a\in 2\mathbb{N}+1$ odd), subject to periodic spatial and Dirichlet temporal boundaries. The prescription $T = L + s\,a$ (with $s = \pm1$) allows for the implementation of spatially constant Abelian background fields interpolated exactly on the lattice, minimizing interpolation artefacts (Pérez-Rubio et al., 2011).
• O(a) Boundary Counterterm: Artefacts linear in the lattice spacing, $O(a)$, can be removed by augmenting the action with a weight factor $c_t$ on time-like plaquettes at the Dirichlet boundaries. The tree-level coefficient is $c_t^{(0)} = 2/(2+s)$, with $s=(T-L)/a$ (Pérez-Rubio et al., 2011). Precise one-loop corrections to $c_t$ are extracted from a perturbative expansion of the lattice SF coupling, ensuring cancellation of all $O(a)$ effects.
• Step Scaling Function: The response of the renormalized coupling to scale changes is captured by $$\Sigma(u,a/L)=\bar{g}^2_{\rm SF}(2L)|_{\bar{g}^2_{\rm SF}(L)=u}$$, with the continuum limit $\sigma(u)=\lim_{a\to 0}\Sigma(u,a/L)$. For sufficiently large $L/a\geq 8$, cutoff-induced deviations $\delta(a/L)$ are within $0.2$–$0.5\%$, demonstrating high precision in finite-volume couplings (Pérez-Rubio et al., 2011).

2. Hamiltonian Truncations and Qutrit Basis on Small Lattices

Efficient analysis and simulation of SU(3) gauge theory on small lattices frequently relies on optimal truncations of the local Hilbert space, designed to capture the essential gauge dynamics while minimizing resource demands:

• Minimal Irrep (“Qutrit”) Truncation: The minimal electric-basis truncation retains only the irreps $\mathbf{1}$, $\mathbf{3}$, and $\mathbf{\overline{3}}$ for each link (Chen et al., 15 Jan 2026). The link basis is then a $\mathbb{Z}_3$ clock, where local gauge-invariant fusion at vertices restricts states to singlets of SU(3).
• Kogut–Susskind Hamiltonian: In this basis, the Hamiltonian is written as

$$H = \frac{g^2\Sigma}{n_\ell}\sum_{\rm links}\hat{E}^2 + \frac{1}{g^2\Sigma}\sum_{\rm plaq}(n_\ell-(U_P+U_P^\dagger))$$

where explicit matrix elements for clock-shift operators are provided, and all nonzero transitions can be efficiently computed (Chen et al., 15 Jan 2026).

• Extensions to SU($N_c$): Truncation strategies generalize naturally by retaining totally antisymmetric irreps, allowing for efficient “qudit” encodings with basis states $|k\rangle$ ($k\in \mathbb{Z}_{N_c}$), and recoupling governed by leading-order $\mathbb{Z}_{N_c}$ clock models.

3. Quantum Simulation Circuits and Resource-Efficient Configurations

Quantum simulation of SU(3) gauge theory on small lattices leverages basis truncations and circuit-level optimizations:

• Electric-Basis Encodings: Electric basis truncations (e.g., $T_1$ for $\mathbf{1}$, $\mathbf{3}$, $\mathbf{\overline{3}}$; $T_2$ for up to two-index irreps) admit encoding in small numbers of qubits per link, with additional registers for site singlet multiplicities (Balaji et al., 11 Mar 2025, Ciavarella et al., 14 Mar 2025).
• Controlled-Plaquette Circuit Designs: Time evolution under $H_B$ is implemented via multi-qudit Givens rotations, with optimization via control pruning (minimizing the number of conditioned qubits) and fusion (aggregating identical gate actions) (Balaji et al., 11 Mar 2025).
• Perturbation Theory and Ansatz Circuits: Strong-coupling perturbative expansions yield compact ground-state ansatz circuits, e.g., single-step “EM” and multi-step “EMEM” Trotter sequences, achieving fidelities $\gtrsim$99% in truncated models on $2\times 2$ and $2\times 2\times 2$ lattices (Balaji et al., 30 Sep 2025).
• Gate and Qubit Resource Scaling: For $T_1$ truncation, one finds per-plaquette counts of $\sim$20–134 CNOTs and $\sim$17–157 $R_z$ gates, with circuit depths scaling as $O(V\log V)$ in volume $V$. Full $2\times 2\times 2$ cubes require $\approx 20$–$48$ qubits and $10^5$–$10^6$ gates per Trotter step in NISQ architectures (Balaji et al., 11 Mar 2025, Ciavarella et al., 14 Mar 2025).

4. Thermodynamics and Equation of State at Phase Transitions

On small lattices, precise determination of nonzero-temperature behavior and first-order deconfinement dynamics is feasible:

• Gradient Flow Methods: The Yang–Mills gradient flow is applied to gauge configurations, evolving link variables $U_\mu(x)$ to flowed fields $V_t(x,\mu)$ over radius $\sqrt{8t}$ and constructing ultraviolet-finite operators for the energy-momentum tensor (Shirogane et al., 2018).
• Determination of Latent Heat and Pressure Gap: Near the critical temperature $T_c$, phase-separated averages of the flowed energy–momentum tensor are computed to yield the latent heat $$\Delta\epsilon/T^4=1.10(5)_{\rm stat}(8)_{\rm syst}$$ and $\Delta p/T^4=0$, with systematic uncertainties well controlled via two complementary extrapolation schemes (Shirogane et al., 2018).
• Finite-Volume and Continuum Corrections: Finite-size effects are suppressed, with $10\%$ amplitude decrease from $R=6$ to $R=8$, and discretization errors of $O(a^2)$ removed via straightforward linear fits in $a^2$. The protocols are robust to changes in fit windows and operator discretizations.

5. Precision Extraction of Hybrid Potentials at Small Separation

Small lattice volumes and fine lattice spacings enable extraction of hybrid static potentials with sub-percent precision:

• Computation of $V_{\Pi_u}(r)$ and $V_{\Sigma_u^-}(r)$: Correlation functions with explicit quantum-number projectors are measured at $a=0.04$–$0.093$ fm, for quark–antiquark separations as small as $r=0.08$ fm (Schlosser et al., 2021, Schlosser et al., 2021).
• Removal of Discretization Effects: Potentials are corrected for $O(a)$ self-energy and $O(a^2)$ shifts via tree-level improved Coulomb subtraction and ensemble-specific additive constants.
• Parameterizations and Born–Oppenheimer Applications: Fitted hybrid potentials yield accurate input for BO calculations of heavy hybrid meson masses, e.g., $M_{b\bar{b}}^{\Pi_u}=10.78(10)$ GeV and $M_{b\bar{b}}^{\Sigma_u^-}=10.82(12)$ GeV, matching earlier lattice predictions (Schlosser et al., 2021).
• Systematic Error Control: Topological freezing, finite-volume effects, and glueball decay contamination are shown to be negligible below $r\lesssim 0.1$ fm, validating the continuum extrapolations.

6. Prepotential and Loop-String-Hadron Formulations on Small Lattices

Alternative approaches to gauge-invariant subspaces have been developed for small lattices:

• Prepotential Construction: SU(3) harmonic oscillator “prepotentials” $a^\dagger_\alpha$ and $b^{\dagger\,\alpha}$, subject to Sp(2,$\mathbb{R}$) constraints, generate local gauge-invariant Hilbert spaces and simplify Mandelstam constraints at the site level, efficiently constructing loop states and organizing basis truncations in a manifestly local manner (0909.2394).
• Local Loop-String-Hadron (LSH) Variables: In $d=1$, the on-site LSH basis employs two bosonic flux modes and three fermionic string-end (hadronic) modes per site, eliminating all nonabelian redundancy and leaving only two abelian Gauss laws per link. The LSH basis yields exponential savings in Hilbert-space dimension compared to naive gauge-fixed fermion bases, and is fully extensible to higher dimensions with local interactions (Kadam et al., 2022).

7. Confinement, Duality, and Analytic Results on Small $R^3\times S^1$ Lattices

Compactification on small circles reveals analytically tractable nonperturbative features:

• Double-Trace Deformations and Adjoint Fermions: Center symmetry in SU(3) is enforced at weak coupling by adding double-trace Polyakov loop terms or periodic adjoint fermions, leading to a confining phase at small circle circumference $L$ (Ogilvie, 2014).
• Monopole Gas and Abelian Duality: Nonperturbative dynamics are governed by three fundamental monopole constituents of the finite-temperature instanton (“caloron”), enabling analytic calculation of the string tension in the semiclassical regime:

$$\sigma = \frac{4g}{\pi L} \sqrt{\xi\,e^{-8\pi^2/(3g^2)}\,\frac{\sqrt{3}}{2}}$$

and its lattice counterpart via dual sine-Gordon theory, both maintaining $Z(3)$ symmetry and a continuously connected confining phase (Ogilvie, 2014).

• Numerical Scaling: Lattice simulations confirm the smooth scaling of the string tension and the preservation of center symmetry for $N_t\leq 6$, and show the analytic regime is connected to the large-$L$ bulk.

---

In summary, the study of SU(3) gauge theory on small lattices synthesizes refinements in boundary-condition implementation, irrep truncations, quantum circuit optimization, thermodynamic measurement, and analytic nonperturbative field theory, establishing the theoretical and practical foundation for precision lattice computations and scalable quantum simulation of Yang–Mills dynamics (Pérez-Rubio et al., 2011, Chen et al., 15 Jan 2026, Balaji et al., 30 Sep 2025, Shirogane et al., 2018, Kadam et al., 2022, Schlosser et al., 2021, Balaji et al., 11 Mar 2025, Ciavarella et al., 14 Mar 2025, Ogilvie, 2014, 0909.2394, Ciavarella et al., 2021).