⚛️ lattice

Symplectic lattice gauge theories in the Grid framework: domain wall fermions and continuum extrapolations

This paper presents the first numerical lattice study of Sp(4) gauge theory with two fundamental Dirac fermions using the Möbius domain-wall fermion formulation within the Grid framework, demonstrating that this approach achieves order-of-magnitude gains in continuum extrapolation efficiency compared to Wilson fermions with moderate computational resources.

Original authors: Ed Bennett, Peter A. Boyle, Luigi Del Debbio, Niccolò Forzano, Ryan C. Hill, Deog Ki Hong, Jong-Wan Lee, C. -J. David Lin, Biagio Lucini, Maurizio Piai, Gianmarco Simonetti, Davide Vadacchino, Alexis

Published 2026-09-18

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ed Bennett, Peter A. Boyle, Luigi Del Debbio, Niccolò Forzano, Ryan C. Hill, Deog Ki Hong, Jong-Wan Lee, C. -J. David Lin, Biagio Lucini, Maurizio Piai, Gianmarco Simonetti, Davide Vadacchino, Alexis Verney-Provatas

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 universe as we understand it is built on a foundation of particles and forces, a framework known as the Standard Model. This model has been spectacularly successful, predicting the existence of particles like the Higgs boson before they were ever found in experiments. Yet, physicists know this picture is incomplete. It cannot explain why the universe is filled with invisible dark matter, nor does it account for the mysterious dark energy driving the cosmos apart. To fill these gaps, scientists propose new theories where the familiar particles are not fundamental, but are instead made of even smaller, tightly bound constituents held together by a new, powerful force. This idea, called composite dynamics, suggests that just as protons are made of quarks, the Higgs particle or dark matter candidates could be made of new, unknown particles interacting through a force far stronger than the ones we see today.

To test these ideas, researchers cannot simply build a machine to smash these hypothetical particles together, because the forces involved are too strong and the energies too high for current technology. Instead, they turn to a different kind of experiment: a virtual one. They create a digital grid, a four-dimensional lattice that mimics the fabric of space-time, and simulate the behavior of these new particles on a supercomputer. By watching how these virtual particles interact, move, and bind together on the grid, scientists can predict what the real universe might look like if these new theories are true. The challenge lies in the mathematics; the equations governing these strong interactions are incredibly difficult to solve, and the computer simulations must be precise enough to distinguish between a real physical signal and the noise of the digital approximation.

A team of researchers from the TELOS collaboration has taken a significant step forward in this virtual exploration by applying a new, more precise method to a specific theory that is a leading candidate for explaining both the Higgs boson and dark matter. This theory involves a group of particles called fermions, which are the building blocks of matter, interacting through a force described by a mathematical structure known as Sp(4). The team focused on a version of this theory with two types of these matter particles. Their goal was to see if they could simulate this system with enough accuracy to predict the masses and other properties of the composite particles that would form, specifically the lightest particles made of these fermions.

The researchers used a sophisticated technique called domain-wall fermions to perform their simulations. In simpler terms, this method adds an extra, invisible dimension to their digital grid. By doing this, they can separate the two "handedness" properties of the particles, which helps preserve a crucial symmetry of nature that is often lost in less precise simulations. This extra dimension allows the particles to behave more like they would in the real, continuous world, rather than being distorted by the blocky nature of the computer grid. The team implemented this method using a powerful software environment called Grid, running their calculations on supercomputers equipped with graphics processing units. They carefully tuned the settings of their simulation, adjusting parameters to ensure that the artificial effects of the extra dimension were minimized and that the results were as close to the true physical laws as possible.

Once the simulation was running smoothly, the team measured the properties of the lightest particles that emerged from the interactions. They calculated the masses of these particles and how strongly they decay, which are key numbers needed to compare the theory with real-world observations. They performed these measurements across several different versions of their digital grid, using different spacings to see how the results changed as the grid became finer. By comparing these results, they were able to project what the values would be if the grid were infinitely fine, effectively removing the digital artifacts to reveal the true physical prediction.

The results of this study show that the new method is highly effective. When the team compared their findings with previous studies of the same theory that used an older, less precise simulation technique, they found that their new approach reached the true physical answer much faster. The older method required a massive amount of computing power to get close to the correct answer, while the new method achieved similar precision with far fewer resources. This is a major breakthrough because it means that scientists can now explore these complex theories with much greater confidence and efficiency. The team found that the masses and decay rates of the particles in this theory are stable and well-defined, providing a solid benchmark for future research.

This work does not prove that this specific theory is the correct description of our universe, but it provides a reliable map of what the universe would look like if it were true. The researchers have demonstrated that with the right tools and careful tuning, it is possible to simulate these exotic, strongly interacting worlds with high precision. Their success opens the door to exploring even more complex versions of these theories, including those that might explain the nature of dark matter or the origin of the Higgs boson. By refining the way they simulate these forces, the team has laid the groundwork for a new generation of discoveries, bringing us closer to understanding the hidden layers of reality that lie beyond the Standard Model.

Technical Summary: Symplectic Lattice Gauge Theories in the Grid Framework: Domain Wall Fermions and Continuum Extrapolations

Problem Statement
The paper addresses the numerical study of the $Sp(4)$ gauge theory coupled to $N_f=2$ Dirac fermions in the fundamental representation. This theory serves as a minimal benchmark for Beyond the Standard Model (BSM) scenarios, specifically Composite Higgs Models (CHMs) and Strongly Interacting Massive Particle (SIMP) dark matter models. While previous lattice studies of this theory utilized Wilson fermions, these formulations suffer from $O(a)$ discretization errors and the explicit breaking of chiral symmetry, necessitating extensive extrapolations to the continuum limit and the massless limit. The authors aim to demonstrate that the Möbius Domain Wall Fermion (MDWF) formulation, implemented within the Grid software framework, offers a superior approach by providing automatic $O(a^2)$ improvement and better control over chiral symmetry, thereby reducing computational costs required to reach the continuum limit.

Methodology
The study employs the MDWF formulation, which extends the four-dimensional lattice into a fifth dimension to separate chiral modes and suppress their overlap, effectively realizing an approximate chiral symmetry. The implementation utilizes the Grid software suite on GPU-based supercomputing architectures.

Key methodological steps include:

Algorithm Optimization: The authors perform extensive scans of the MDWF algorithmic parameters—specifically the Möbius scaling parameter $\alpha$ , the fifth-dimensional scale $a_5/a$ , the domain-wall height $am_5$ , and the Pauli-Villars mass $am_{PV}$ —to minimize the residual mass ( $am_{res}$ ). The residual mass quantifies the explicit breaking of chiral symmetry due to the finite extent of the fifth dimension ( $N_5$ ).
Ensemble Generation: Gauge configurations are generated for three lattice couplings ( $\beta = 7.4, 7.5, 7.6$ ) and several bare fermion masses ( $am_0$ ). The fifth dimension is fixed at $N_5=8$ for the main production runs, with tuning scans performed to optimize $\alpha$ for each $\beta$ .
Scale Setting and Topology: The Wilson flow is used to set the lattice scale via the $w_0$ parameter. The topological charge is monitored to ensure the absence of topological freezing and to verify that the ensembles are well-thermalized.
Spectroscopy: Two-point correlation functions are computed for pseudoscalar (PS), vector (V), and axial-vector (AV) mesons. Masses and decay constants are extracted using multi-exponential fits and the Hadrons analysis framework.
Continuum Extrapolation: Results are extrapolated to the continuum limit ( $a \to 0$ ) using an ansatz inspired by Wilson Chiral Perturbation Theory (W $\chi$ PT), adapted for the $O(a^2)$ scaling expected of MDWF.

Key Contributions

First MDWF Study of $Sp(4)$ : This work presents the first numerical lattice study of the $Sp(4)$ gauge theory with $N_f=2$ using domain-wall fermions.
Parameter Optimization: The paper provides a detailed optimization of MDWF parameters for this specific gauge group, identifying optimal values for $\alpha$ that minimize the residual mass while maintaining computational efficiency.
Implementation Validation: The authors validate their implementation through gauge invariance tests and consistency checks, such as verifying the relation between the PCAC mass and the sum of the bare mass and residual mass ( $am_{PCAC} \approx am_0 + am_{res}$ ).
Comparative Analysis: The study directly compares MDWF results with previous Wilson fermion results for the same theory, highlighting the differences in discretization artifacts.

Results

Residual Mass Suppression: Through the tuning of $\alpha$ , the authors successfully suppressed the residual mass to values significantly smaller than the input bare fermion masses (e.g., $am_{res} \approx 0.0015$ for $am_0 = 0.02$ ), ensuring good chiral symmetry properties.
Spectroscopy: Precise measurements of the lightest flavored meson masses ( $m_{PS}, m_V$ ) and the pseudoscalar decay constant ( $f_{PS}$ ) were obtained across the generated ensembles.
Continuum Extrapolation:
- MDWF: The MDWF data exhibits a leading discretization dependence of $O(a^2)$ . The linear $O(a)$ term is consistent with zero, allowing for a robust continuum extrapolation using fewer lattice spacings.
- Wilson Comparison: In contrast, the Wilson fermion data from previous studies showed a clear $O(a)$ dependence, requiring more complex fitting ansätze and more data points to achieve a reliable continuum limit.
- Agreement: The continuum extrapolated values for masses and decay constants obtained via MDWF are compatible with those derived from Wilson fermions within statistical uncertainties.
Computational Efficiency: The authors demonstrate that with moderate computational resources, the MDWF formulation yields order-of-magnitude gains in the approach to the continuum limit compared to the unimproved Wilson formulation in the region of physical parameter space relevant to phenomenology.

Significance and Claims
The paper claims that the adoption of the MDWF formulation, combined with the Grid software environment, enables high-precision lattice studies of new strong dynamics theories with significantly reduced systematic errors from discretization. The authors assert that their results supersede the systematic precision of existing measurements in the literature for the $Sp(4)$ , $N_f=2$ theory.

The significance of this work is twofold:

Phenomenological: It provides the most precise available measurements of the properties (masses and decay constants) of the lightest spin-0 and spin-1 flavored mesons in this theory, which are critical inputs for building Composite Higgs and Dark Matter models.
Methodological: It demonstrates that the MDWF formulation can be effectively deployed for symplectic gauge theories, offering a pathway to high-precision continuum limits without the prohibitive computational costs often associated with chiral fermions. The authors note that while their current study uses moderately large fermion masses, the methodology sets the stage for future investigations into lighter mass regimes and more complex theories (e.g., those including fermions in the antisymmetric representation) where chiral symmetry is essential.

The paper concludes modestly, stating that while the results are promising, the full potential of this technology for BSM lattice physics remains to be fully explored in more challenging regions of parameter space and for more complex gauge groups.

More like this