Massive Schwinger Model in QED 1+1D

• The massive Schwinger model is QED in 1+1 dimensions with a finite fermion mass, exhibiting a rich vacuum structure, explicit chiral symmetry breaking, and CP-violating effects.
• It leverages bosonization, lattice simulations, and tensor-network methods to reveal nonperturbative dynamics and accurately predict spectral, chiral, and topological properties.
• Recent advancements using DMRG, TRG, and quantum simulation affirm its role as a benchmark for exploring real-time dynamics and effective field theories analogous to QCD.

The massive Schwinger model is quantum electrodynamics in $1+1$ dimensions with a nonzero fermion mass, formulated either for a single Dirac fermion or for multiple flavors depending on context. In continuum notation, its basic structure is that of QED$_{1+1}$ with a mass term $m\bar\psi\psi$, and, when present, a topological $\theta$-term proportional to $\epsilon^{\mu\nu}F_{\mu\nu}$ or equivalently to the electric field $E$ (Keegan, 2015, Zache et al., 2018, Grieninger et al., 2023). The model occupies a distinctive position in quantum field theory because it combines confinement, anomaly, nontrivial vacuum structure, and accessible nonperturbative control in a setting that is substantially simpler than QCD while still exhibiting closely analogous phenomena (Zache et al., 2018, Grieninger et al., 2023, Cruz et al., 2024). In current research it functions simultaneously as an exactly or semi-analytically tractable field theory, a benchmark for lattice and tensor-network methods, and a platform for studying real-time topology, entanglement, hadronization, and quantum criticality (Keegan, 2015, Florio et al., 2024, Cruz et al., 2024, Florio et al., 17 Jun 2025).

1. Definition and field-theoretic structure

The massive Schwinger model is obtained by adding a fermion mass to the Schwinger model, i.e. to QED in $1+1$ dimensions. In continuum form, representative Lagrangians used in the literature include

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$

for the one-flavor theory (Zache et al., 2018), and, for multi-flavor versions,

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$

for the two-flavor case (Albandea et al., 8 Jan 2025). Euclidean formulations with explicit $\theta$-terms are also standard,

$_{1+1}$0

in studies of finite temperature, density, and vacuum angle dependence (Grieninger et al., 2023).

The massless and massive cases differ qualitatively. The massless Schwinger model is exactly solvable and produces an anomalous gauge-boson mass, while the massive theory is not exactly solvable in general and develops a richer infrared structure, especially in the presence of a $_{1+1}$1-term or multiple fermion flavors (Keegan, 2015, Grieninger et al., 2023, Cruz et al., 2024). In $_{1+1}$2 dimensions the gauge coupling has dimension of mass, so ratios such as $_{1+1}$3 or $_{1+1}$4 control the dynamical regime (Zache et al., 2018, Cruz et al., 2024).

For the one-flavor model, the $_{1+1}$5-term can be written as

$_{1+1}$6

and, after a chiral rotation in temporal axial gauge, it appears as a complex phase in the mass term,

$_{1+1}$7

making explicit that $_{1+1}$8 cannot be rotated away once $_{1+1}$9 (Zache et al., 2018). This is the basis for the model’s use as a prototype of CP violation and topological vacuum structure.

For $m\bar\psi\psi$0, the model acquires an additional flavor structure. In the massless limit, bosonization yields one massive boson of mass

$m\bar\psi\psi$1

and $m\bar\psi\psi$2 massless bosons (Keegan, 2015). With nonzero mass, the theory is no longer exactly solvable, but strong-coupling bosonization, lattice gauge theory, tensor methods, and effective-field-theory constructions provide complementary access to its spectrum and vacuum structure (Keegan, 2015, Albandea et al., 8 Jan 2025).

2. Bosonization, anomaly, and low-energy descriptions

A central structural feature of the massive Schwinger model is its bosonized description. For the one-flavor theory at finite $m\bar\psi\psi$3, $m\bar\psi\psi$4, and $m\bar\psi\psi$5, the bosonized effective action takes the form

$m\bar\psi\psi$6

with

$m\bar\psi\psi$7

(Grieninger et al., 2023). Here $m\bar\psi\psi$8 is the anomalous Schwinger mass, while $m\bar\psi\psi$9 is a mass scale induced by explicit chiral symmetry breaking. In the strong-coupling regime $\theta$0, the low-energy theory is dominated by the anomaly-generated mass scale $\theta$1, and the physical pseudoscalar mass obeys

$\theta$2

for small $\theta$3 (Grieninger et al., 2023).

For the two-flavor theory in the strong-coupling, small-mass regime, integrating out the heavy singlet yields a sine-Gordon description for the light field. The bosonized Lagrangian quoted for the two-flavor model is

$\theta$4

with $\theta$5 and $\theta$6 (Albandea et al., 8 Jan 2025). In the limit $\theta$7, the heavy $\theta$8 field can be integrated out, leaving an effective sine-Gordon theory for $\theta$9 (Albandea et al., 8 Jan 2025). This yields exact or semi-exact statements about the light spectrum, including

$\epsilon^{\mu\nu}F_{\mu\nu}$0

in the strong-coupling regime (Albandea et al., 8 Jan 2025).

Bosonization also clarifies the role of the anomaly and of $\epsilon^{\mu\nu}F_{\mu\nu}$1. In the one-flavor model, the $\epsilon^{\mu\nu}F_{\mu\nu}$2-dependence is encoded in the phase shift of the cosine term and in the vacuum structure of the scalar potential (Grieninger et al., 2023). In the two-flavor theory, the anomaly makes the singlet $\epsilon^{\mu\nu}F_{\mu\nu}$3-like mode heavy, while the light sector exhibits nontrivial mass scaling but no spontaneous continuous chiral symmetry breaking in the strict massless limit, consistent with the Mermin–Wagner–Coleman theorem as discussed in the two-flavor analysis (Albandea et al., 8 Jan 2025). A plausible implication is that the massive model’s utility derives precisely from this coexistence of anomalous singlet dynamics and a nontrivial light sector that remains tractable.

3. Vacuum angle, topology, and CP structure

The $\epsilon^{\mu\nu}F_{\mu\nu}$4-dependence of the massive Schwinger model is one of its defining characteristics. In $\epsilon^{\mu\nu}F_{\mu\nu}$5 dimensions the topological term is local and linear in the electric field, and the total electric flux serves as the topological charge (Zache et al., 2018). In the one-flavor theory, $\epsilon^{\mu\nu}F_{\mu\nu}$6 labels distinct vacuum sectors and becomes a physical parameter whenever $\epsilon^{\mu\nu}F_{\mu\nu}$7, producing genuine CP-violating effects away from $\epsilon^{\mu\nu}F_{\mu\nu}$8 (Zache et al., 2018).

For the one-flavor massive theory, Euclidean lattice studies with quantized topological charge reconstruct the full $\epsilon^{\mu\nu}F_{\mu\nu}$9-dependence of the order parameter and find that the results at $E$0 are compatible with Coleman’s conjectured phase diagram (Azcoiti et al., 2017). In that framework, the topological charge density behaves as an order parameter for CP at $E$1, and the extracted behavior distinguishes a CP-broken regime at large fermion mass from a CP-symmetric regime at small mass (Azcoiti et al., 2017). The same study emphasizes exact $E$2-periodicity through a lattice definition of integer-valued topological charge (Azcoiti et al., 2017).

For $E$3, recent Grassmann tensor renormalization group studies with staggered fermions and a $E$4-periodic logarithmic $E$5-term compute the free energy in the thermodynamic limit and analyze its $E$6-dependence across a broad mass range (Kanno et al., 2024, Kanno et al., 23 Jan 2025). In the large-mass limit, the free energy approaches the Maxwell-theory form

$E$7

which has a cusp at $E$8, indicating a first-order transition and two degenerate vacua (Kanno et al., 23 Jan 2025). The tensor-network calculation reproduces this limit and shows twofold ground-state degeneracy at $E$9 in the large-mass regime (Kanno et al., 23 Jan 2025).

The same lattice tensor studies also report that, at finite lattice spacing, the $1+1$0 lattice theory can display a phase structure differing from continuum expectations in the small-mass regime (Kanno et al., 2024, Kanno et al., 23 Jan 2025). In particular, the free energy at $1+1$1 still shows visible $1+1$2-dependence at moderate $1+1$3, even though the continuum massless $1+1$4 theory should be $1+1$5-independent (Kanno et al., 23 Jan 2025). This suggests that the finite-$1+1$6 lattice phase structure is strongly affected by lattice artifacts. The authors explicitly interpret the apparent small-mass deviations from continuum behavior as finite-$1+1$7 effects rather than as evidence against the continuum picture (Kanno et al., 23 Jan 2025).

The dynamical role of $1+1$8 extends beyond equilibrium. Quenches of $1+1$9 in the one-flavor massive model generate real-time dynamical quantum phase transitions, with the Loschmidt amplitude

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$0

and a gauge-invariant two-time fermion correlator serving as central observables (Zache et al., 2018). For strong quenches $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$1, the theory exhibits zeros of the mode-resolved correlator, nonanalyticities in the rate function, and quantized jumps of a dynamical topological order parameter built from winding numbers of the correlator phase (Zache et al., 2018). This establishes the massive Schwinger model as a concrete setting where topological vacuum structure survives in controlled real-time dynamics.

4. Spectrum, chiral observables, and quantum criticality

The massive Schwinger model supports several distinct spectral regimes, depending on flavor number, mass ratio, and vacuum angle. In the one-flavor case at $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$2, varying $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$3 drives a second-order phase transition between a symmetric phase and a symmetry-broken phase. Bosonization at $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$4 leads to a scalar potential with two degenerate minima above a critical mass, and the transition belongs to the 2D Ising universality class (Cruz et al., 2024, Jentsch et al., 2021).

A recent high-precision DMRG study of the one-flavor model at $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$5 determines the critical mass as

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$6

excluding the conjecture that $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$7 exactly (Cruz et al., 2024). That study uses four finite-size criticality criteria based on gap scaling and entanglement entropy, and finds perfect agreement among them (Cruz et al., 2024). The same work confirms the Ising CFT interpretation through gap ratios and entanglement scaling with central charge $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$8 (Cruz et al., 2024). Earlier NISQ-based quantum computation on very small lattices gave a hardware estimate $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \bar\psi(i\gamma^\mu D_\mu - m)\psi$$9, consistent with the classical value $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$0 within the quoted uncertainties (Thompson et al., 2021).

The nonperturbative functional renormalization group reaches the same qualitative picture from the bosonized side. It finds that the phase transition driven by $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$1 belongs to the 2D Ising universality class and estimates the critical ratio as

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$2

in good agreement with DMRG (Jentsch et al., 2021). The FRG analysis further computes temperature and $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$3-dependence of the chiral density, electric field, and entropy density, and concludes that screening of fractional charges and deconfinement occur only at infinite temperature (Jentsch et al., 2021).

In the two-flavor theory at $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$4, the strong-coupling low-energy spectrum is described by the sine-Gordon sector. Lattice simulations confirm the exact sine-Gordon prediction

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$5

and discriminate against semiclassical and WKB coefficients, which are slightly larger (Albandea et al., 8 Jan 2025). They also confirm the exact ratio

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$6

through the low-energy effective theory based on a nonlinear $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$7-model plus a dilaton (Albandea et al., 8 Jan 2025). For the decay constant, the same study finds lattice values consistent with

$$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$8

and inconsistent with identifying $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \sum_{i=1}^2 \bar{\psi}_i \left( i\!\not\!\partial - g \!\not\!A - m_i \right) \psi_i$$9 with the Witten–Veneziano value $\theta$0 (Albandea et al., 8 Jan 2025).

The chiral condensate shows nontrivial mass scaling in both one- and two-flavor settings. For $\theta$1, the condensate scales as

$\theta$2

so it vanishes as $\theta$3 in the chiral limit (Albandea et al., 8 Jan 2025). In the one-flavor thermal and density-dependent bosonized theory, the finite-temperature chiral condensate enters through a temperature-dependent boson mass $\theta$4, with the cosine contribution melting at high temperature (Grieninger et al., 2023). These results collectively indicate that the massive Schwinger model encodes explicit symmetry breaking in a way that preserves delicate anomalous and infrared features rather than reducing to a trivial mass deformation.

5. Lattice, tensor-network, and spectral methods

The massive Schwinger model has become a standard testbed for nonperturbative numerical methods. On the lattice, staggered fermions, Wilson fermions, and momentum-space constructions have all been used, depending on the observable of interest (Keegan, 2015, Thompson et al., 2021, Cruz et al., 2024, Albandea et al., 8 Jan 2025). Open-boundary Kogut–Susskind Hamiltonians are particularly useful because Gauss’s law can be used to integrate out the gauge field, leaving long-range spin Hamiltonians suitable for DMRG or MPS simulations (Cruz et al., 2024, Florio et al., 2024, Florio et al., 17 Jun 2025).

For the two-flavor massive Schwinger model without a $\theta$5-term, lattice studies of the Dirac spectral density use Wilson fermions and the mode number of the hermitian Dirac operator to extract a scale-dependent mass anomalous dimension $\theta$6 (Keegan, 2015). In the $\theta$7 model, analytic arguments imply an infrared value

$\theta$8

and lattice fits to the mode number and spectral density show $\theta$9 running from $_{1+1}$00 in the ultraviolet toward $_{1+1}$01 in the infrared, with the smallest accessible eigenvalues within roughly $_{1+1}$02 of the expected value (Keegan, 2015). The quenched theory, by contrast, yields an effective infrared $_{1+1}$03, corresponding to nonzero chiral condensate in the fit-form language (Keegan, 2015). This study established the model as a benchmark for systematic errors in the mode-number method.

Tensor-network methods have been especially influential for $_{1+1}$04-dependent and real-time problems. Grassmann TRG formulations with a $_{1+1}$05-periodic logarithmic $_{1+1}$06-term avoid sign problems and enable direct access to the thermodynamic limit for the free energy and topological observables (Kanno et al., 2024, Kanno et al., 23 Jan 2025). DMRG with up to $_{1+1}$07 staggered-fermion sites has been used to determine critical properties with five-digit precision near $_{1+1}$08 (Cruz et al., 2024). MPS time evolution via TDVP has made possible direct simulations of jet-like real-time dynamics, entanglement growth, and emergent hydrodynamics (Florio et al., 2024, Florio et al., 17 Jun 2025).

Momentum-space truncations offer an alternative route, particularly for NISQ-era quantum computation. In a Hamiltonian formulation on a circle, solving Gauss’s law leaves only a gauge zero mode plus fermionic modes, and large gauge transformations define explicit $_{1+1}$09-sector basis states (Thompson et al., 2021). Truncating to the most relevant basis states allowed a three-qubit VQE implementation of the $_{1+1}$10 transition with results consistent with the known critical mass (Thompson et al., 2021). This suggests that the model’s combination of gauge structure and low-dimensionality makes it unusually well-suited for symmetry-adapted quantum simulation.

6. Real-time dynamics, entanglement, hadronization, and thermalization

The massive Schwinger model has become a major venue for studying real-time nonequilibrium gauge dynamics. In jet-like setups, external sources create back-to-back charges that mimic quark jets, producing a flux tube that fragments through pair production (Florio et al., 2024, Florio et al., 17 Jun 2025). Using a staggered-fermion Hamiltonian with open boundaries and source-dependent external electric fields, simulations track the full quantum evolution from the interacting vacuum (Florio et al., 2024).

These studies show that the entanglement entropy between left and right halves grows rapidly in time, driven by an increasing number of reduced-density-matrix eigenvalues with comparable weight (Florio et al., 2024). At early times the dominant entanglement eigenvectors are close to simple fermionic Fock states, while at stronger coupling and later times they evolve into extended superpositions interpreted as meson-like bound states (Florio et al., 2024). This provides a dynamical picture of hadronization: the relevant entangled degrees of freedom interpolate from parton-like excitations to bound states as the flux tube breaks (Florio et al., 2024).

Local observables in the central rapidity region approach approximately constant values at late times, suggesting local equilibration (Florio et al., 2024). A later study makes this statement quantitative by comparing central reduced density matrices and local observables to true thermal states of the same Hamiltonian (Florio et al., 17 Jun 2025). It finds that the temperatures extracted from the chiral condensate, nearest-neighbor condensate correlator, kinetic energy density, entanglement entropy density, and Hilbert–Schmidt overlap with Gibbs states all converge at late times, with especially good agreement near $_{1+1}$11 (Florio et al., 17 Jun 2025). The same work shows that the entanglement entropy crosses over from an area law to a volume law, and that the entanglement entropy density matches the Gibbs entropy density of the equilibrium theory when expressed at the corresponding effective temperature (Florio et al., 17 Jun 2025). This supports the interpretation that quantum entanglement drives the emergence of local thermal behavior in the closed field-theoretic system.

Hydrodynamic behavior has also been observed. Tensor-network simulations of a localized excitation atop the vacuum in the one-flavor model show that, for light fermions $_{1+1}$12, the energy density, fluid velocity, and bulk pressure evolve consistently with Bjorken flow (Shao et al., 13 Sep 2025). The fluid velocity approaches $_{1+1}$13, the energy density develops a rapidity plateau, and the proper-time dependence follows the Bjorken scaling law

$_{1+1}$14

with the speed of sound $_{1+1}$15 obtained from the equation of state of the same lattice Hamiltonian (Shao et al., 13 Sep 2025). The hydrodynamic onset time aligns with the thermalization time of the quantum distribution function (Shao et al., 13 Sep 2025). By contrast, in the heavier-mass regime $_{1+1}$16, Bjorken-like hydrodynamics does not emerge (Shao et al., 13 Sep 2025).

A complementary bosonized treatment of a receding charge pair in Bjorken coordinates studies pair production and string breaking through a coherent electric field plus Gaussian fluctuations (Batini et al., 2024). The electric field shows damped oscillations, and the asymptotic total particle density per rapidity interval at large mass fits a Boltzmann factor, allowing the extraction of an effective temperature proportional to $_{1+1}$17, with $_{1+1}$18 the string tension (Batini et al., 2024). This has been interpreted as an analog of a hadronization temperature in QCD (Batini et al., 2024). A plausible implication is that, in low-dimensional confining systems, thermal-like spectra may emerge from string breaking and entanglement even before conventional kinetic equilibration becomes the appropriate description.

7. Entanglement, finite density, condensed-matter analogs, and broader significance

Spatial entanglement in the massive Schwinger model has recently been analyzed in the strong-coupling regime at finite temperature, density, and $_{1+1}$19-angle using bosonization (Grieninger et al., 2023). The UV-finite entropic function

$_{1+1}$20

interpolates between a short-distance CFT-like regime and a large-interval regime governed by the boson mass $_{1+1}$21 (Grieninger et al., 2023). For $_{1+1}$22, the entropic function behaves as

$_{1+1}$23

so the entanglement is controlled by the effective correlation length $_{1+1}$24 (Grieninger et al., 2023). At finite density, the bosonized vacuum forms a chiral density wave,

$_{1+1}$25

but the leading-order entropic function remains independent of $_{1+1}$26, because the light bosonic excitation mass does not depend on density at that order (Grieninger et al., 2023).

The massive Schwinger model also supports condensed-matter analogies. An extended multi-species version describes carbon nanotubes when electric flux is confined to one dimension, yielding an effective multi-flavor massive Schwinger model with linearly confining potential $_{1+1}$27 (Oka, 2010). In that context, excitons play the role of mesons, and the light-front bound-state equation is a multi-species extension of the ’t Hooft–Bergknoff equation (Oka, 2010). The work further argues that a dark exciton obeys a Gell-Mann–Oakes–Renner-like relation and that nonlinear transport can probe Coleman’s half-asymptotic state (Oka, 2010). This suggests that some aspects of confinement and anomaly physics in the massive Schwinger model are not confined to high-energy analogies alone.

More broadly, the model’s contemporary role is methodological as much as conceptual. It has become a proving ground for DMRG, MPS real-time evolution, exact diagonalization, FRG, Grassmann TRG, and NISQ-oriented variational algorithms (Zache et al., 2018, Thompson et al., 2021, Cruz et al., 2024, Kanno et al., 2024, Florio et al., 17 Jun 2025). It is used both to validate numerical methods against known analytic limits and to explore regions where no closed-form solution exists, such as real-time topological transitions, finite-density entanglement, and isospin breaking in the two-flavor theory (Zache et al., 2018, Grieninger et al., 2023, Albandea et al., 8 Jan 2025).

A recurring theme across these directions is that the adjective “massive” does not merely designate a trivial deformation of an exactly solvable model. Rather, it marks the point at which the Schwinger model acquires the full complexity of explicit chiral symmetry breaking, $_{1+1}$28-dependent vacuum structure, nontrivial low-energy spectroscopy, and rich real-time dynamics. The available evidence indicates that this complexity remains sufficiently constrained to allow precision tests of nonperturbative ideas, yet sufficiently rich to serve as a credible analog of more complicated confining gauge theories (Keegan, 2015, Zache et al., 2018, Cruz et al., 2024, Florio et al., 17 Jun 2025).