⚛️ lattice

Solutions of Roy equations for variable pion mass

This paper solves the Roy equations for pion scattering at variable pion masses, utilizing chiral perturbation theory and lattice data to provide an accurate description of partial waves and resonance behavior that fully incorporates analyticity, unitarity, and crossing symmetry constraints.

Original authors: Gilberto Colangelo

Published 2026-10-14

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Gilberto Colangelo

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

To understand the story of how matter holds together, physicists often look at the simplest building blocks of the atomic nucleus: protons and neutrons. These particles are not solid spheres but are instead made of smaller, more fundamental pieces called quarks, which are glued together by a force carried by particles known as gluons. The behavior of these quarks and gluons is governed by a set of rules called quantum chromodynamics. While this theory is mathematically precise, it is notoriously difficult to use for making predictions about the everyday world because the equations become incredibly complex when the particles interact strongly. To bridge the gap between the fundamental rules and the physical world, scientists rely on a method called lattice quantum chromodynamics. Imagine a grid, like a three-dimensional chessboard, where the laws of physics are simulated on a computer. By placing the quarks on this grid, researchers can calculate how they move and interact. However, these simulations are incredibly expensive to run, and to make them manageable, scientists often have to use simplified versions of reality, such as making the particles heavier than they actually are in nature. This creates a puzzle: how can we trust results from a simulation where the particles are heavier than the ones we see in the real world?

A researcher has developed a new way to solve this puzzle by connecting the heavy, simulated world with the light, real world using a set of mathematical tools known as Roy equations. These equations act as a rigorous bridge, ensuring that the results of a simulation obey the fundamental laws of physics, such as the conservation of energy and the symmetry of time and space. The researcher used these equations to explore what happens to the interactions between pions—particles made of a quark and an antiquark—as their mass is changed. They started by solving the equations for pion masses lighter than those found in nature, all the way down to a theoretical limit where the mass is nearly zero. They found that even as the mass changed, the behavior of the particles remained smooth and predictable, following the patterns expected by theory. This confirmed that their mathematical framework was robust enough to handle the extreme conditions of the chiral limit, a state where the quarks become massless.

The researcher then turned their attention to the other side of the problem: what happens when the pion mass is heavier than in our real world? This is the exact scenario used in many current computer simulations. They applied their method to pion masses up to 240 MeV, a value significantly higher than the physical mass of about 140 MeV. When they compared their results to existing data from lattice simulations at a mass of 239 MeV, the two sets of numbers matched with remarkable precision. The researcher did not need to force the data to fit; the mathematical framework naturally produced the correct answer. This success suggests that the Roy equations can serve as a powerful tool for interpreting future lattice simulations, allowing scientists to extract accurate physical insights even when the simulations are run with heavier, unphysical particles.

The study also provided a deeper look at the properties of the lightest resonances, which are short-lived particles that appear when pions collide. Specifically, the researcher examined the sigma and rho particles. While the current work establishes the mathematical framework necessary to study these particles, a dedicated analysis of how the sigma and rho resonance parameters depend on the mass of the light quarks will be presented in a subsequent paper. This is a crucial step because the sigma particle, in particular, is a subject of long-standing debate in physics, with some theories suggesting it plays a special role in how the universe acquired its mass. The ability to study these particles across a wide range of masses, from the lightest possible values to the heavier ones used in simulations, opens the door to a more complete understanding of the strong force. The researcher confirmed that their approach works well across this entire range, providing a reliable map that connects the abstract world of computer simulations to the tangible reality of particle physics.

One of the most significant findings was the confirmation that the high-energy behavior of these particle interactions does not change drastically when the mass of the quarks is altered. The researcher had to make an assumption that the complex interactions occurring at higher energies would remain essentially the same, regardless of whether the particles were light or heavy. They tested this assumption by checking if the results still made sense when they varied the mass. The fact that the equations held together so well, even when the mass was changed by large amounts, validated their approach. It means that scientists can now use these equations to take the raw data from heavy-particle simulations and translate it into accurate predictions for the real world, without losing the rigorous constraints of physics along the way.

The work also involved a careful check of the results using specific mathematical relationships known as sum rules. These rules act as a consistency check, ensuring that the calculated properties of the particles at the very lowest energies match the predictions of a different theoretical framework called chiral perturbation theory. The researcher found that their results agreed perfectly with these predictions across the entire range of masses they studied. This agreement is a strong indicator that the method is not just a numerical trick, but a genuine reflection of how nature works. It confirms that the underlying assumptions about how the particles behave are correct and that the mathematical tools are being applied properly.

Looking ahead, this research paves the way for a more systematic analysis of lattice data. Instead of treating each simulation as an isolated event, scientists can now use these equations to weave together results from different mass values into a coherent picture. This will be particularly useful for studying the sigma and rho particles, whose properties are sensitive to the mass of the quarks. By having a reliable method to connect the heavy simulation world with the light physical world, researchers can finally extract precise information about these elusive particles. The study demonstrates that it is possible to vary the mass of the quarks within a rigorous mathematical framework without losing the essential features of the theory. This is a significant step forward, offering a new path to understanding the fundamental forces that shape our universe.

The researcher emphasized that while their method works well for the range of masses they tested, there are limits to how far it can be pushed. As the mass of the pion increases beyond a certain point, the assumptions about the stability of the high-energy input might begin to break down. However, for the masses currently used in lattice simulations, the method is robust and reliable. The researcher plans to apply this approach in future work to provide a detailed study of how the sigma and rho resonances depend on the mass of the light quarks. This will allow for a more precise determination of their properties, which has been a goal of the field for many years.

In the end, this paper is about building a bridge. It connects the heavy, computationally manageable world of lattice simulations with the light, physical world of real particles. By using the Roy equations as a sturdy structure, the researcher has shown that it is possible to cross this gap without losing the fundamental laws of physics. The results are not just a collection of numbers, but a validated method that can be used to unlock the secrets of the strong force. The work stands as a proof of principle that dispersive approaches, which rely on fundamental principles like analyticity and unitarity, can be successfully extended to unphysical masses. This opens up new possibilities for exploring the behavior of matter under conditions that are difficult to reach directly, bringing us closer to a complete understanding of the subatomic world.

Technical Summary: Solutions of Roy equations for variable pion mass

Problem Statement
The Roy equations provide a rigorous framework for determining the $\pi\pi$ scattering amplitude by combining analyticity, unitarity, and crossing symmetry. Historically, these equations have been solved for physical quark masses, relying on phenomenological input for high-energy regions and chiral perturbation theory ( $\chi$ PT) for subtraction constants (scattering lengths). A significant limitation of this approach is its applicability to lattice QCD calculations performed at unphysical quark masses. While $\chi$ PT can rigorously control the dependence on light-quark masses, it is limited to low-energy regimes and requires low-energy constants (LECs) that are difficult to determine phenomenologically. Conversely, dispersive analyses like the Roy equations extend the region of validity but were previously restricted to the physical point because the high-energy phenomenological input was assumed to be mass-dependent. This paper addresses the limitation of connecting dispersive analyses with lattice simulations at unphysical pion masses, specifically aiming to solve the Roy equations for pion masses ( $M_\pi$ ) both below and above the physical value.

Methodology
The author extends the numerical solution of the Roy equations for the S and P waves to variable pion masses in the range $0 \le M_\pi \le 240$ MeV. The methodology relies on a central assumption: the phenomenological input for the imaginary parts of the S and P waves above a matching point $s_A = 0.84 \text{ GeV}^2$ (and for all other waves from threshold to infinity) remains essentially unchanged as light-quark masses vary.

Key methodological components include:

Framework Definition: The solution is sought between the threshold and $s_A$ . The subtraction constants (scattering lengths $a_0^0, a_0^2$ ) are fixed using $\chi$ PT predictions up to $O(p^6)$ .
Input Flexibility: To accommodate the constraints of the Roy equations as $M_\pi$ decreases, the author introduces "hyperparameters" to allow for slight adjustments in the input at $s_A$ . Specifically, parameters are introduced to control the value and the derivative of the imaginary parts of the S-waves at $s_A$ , and a parameter $x_2^0$ controls the interpolation of the S2 wave between $s_A$ and $s_2 = (1.6 \text{ GeV})^2$ .
Threshold and Kinematic Shifts: The position of the $K\bar{K}$ threshold is shifted according to $\chi$ PT predictions for the charged kaon mass dependence on $M_\pi$ . A kinematic mapping is applied to the imaginary part of the S0 wave around this threshold to account for the $f_0(980)$ resonance behavior.
Parametrization: A Chebyshev polynomial parametrization is adopted for the phase shifts, replacing the Schenk parametrization used in previous works, to improve numerical efficiency and solution quality.
Validation: The solutions are tested against Wanders' sum rules to evaluate threshold parameters (scattering lengths and effective ranges) and compare them with $\chi$ PT predictions. Additionally, solutions for $M_\pi > M_\pi^{phys}$ are compared with existing lattice QCD data at $M_\pi = 239$ MeV.

Key Contributions and Results

Solutions for Light Quark Masses: The Roy equations were successfully solved for pion masses ranging from the physical value down to 5 MeV (near the chiral limit). The solutions maintain excellent quality, with the residual error ( $Q^2$ ) increasing only by two orders of magnitude from the physical point to 5 MeV, remaining well within acceptable limits.
Chiral Behavior Verification: The threshold parameters (S, P, and D waves) calculated via Wanders' sum rules were found to agree excellently with $\chi$ PT predictions across the entire range $0 \le M_\pi \le M_\pi^{phys}$ . This confirms that the mass dependence of the dispersive integrals is consistent with chiral symmetry constraints. A minor discrepancy was noted in the chiral limit of the parameter $b_1^1$ , attributed to a small non-vanishing chiral logarithm in the dispersive solution compared to the strictly vanishing one in leading-order $\chi$ PT; this was shown to be a result of numerical precision rather than a fundamental flaw.
Extension to Heavy Pion Masses: The framework was extended to pion masses up to 240 MeV. Solutions were found to be of high quality.
Comparison with Lattice Data: At $M_\pi = 239$ MeV, the central Roy solution (using physical input above $s_A$ ) showed good agreement with lattice data from the Hadron Spectrum Collaboration, particularly for the S-waves. Local discrepancies in the P-wave were reduced by allowing the input parameters at $s_A$ to vary slightly, demonstrating that the approach can be adapted to fit lattice data without free parameters in the low-energy region.
Mass Dependence of Imaginary Parts: The imaginary parts of the partial waves were found to depend smoothly on $M_\pi$ , well-described by a quadratic polynomial in $M_\pi^2$ (linear in light-quark masses) across the energy region between threshold and $s_A$ .

Significance
The paper establishes a proof of principle that Roy equations can be solved for variable light-quark masses, effectively bridging the gap between dispersive analyses and lattice QCD simulations at unphysical masses. By demonstrating that the high-energy phenomenological input can be treated as largely independent of the light-quark masses (within uncertainties), the work enables the application of rigorous analyticity and unitarity constraints to lattice data. This paves the way for a more accurate study of the dependence of resonances (specifically the $\sigma$ and $\rho$ ) on light-quark masses and allows for the determination of pole positions in the complex plane for unphysical pion masses without relying solely on model-dependent unitarization methods like the Inverse Amplitude Method (IAM). The author notes that a dedicated study of the $\sigma$ and $\rho$ resonances using this framework will be presented in future work.

More like this