Towards a formalism for $\pi\pi$ scattering from staggered lattice QCD

This paper addresses the challenges of extracting $\pi\pi$ scattering amplitudes from staggered lattice QCD by proposing two complementary approaches: calculating one-loop amplitudes using Rooted Staggered Chiral Perturbation Theory to verify quantization conditions, and generalizing the Lüscher formalism to explicitly incorporate taste-splitting and fourth-rooting effects.

The Gist

Imagine you are trying to understand how two billiard balls bounce off each other. In the real world, you can watch them collide in an infinite, open room. But in the world of Lattice QCD (a super-computer simulation of the strong nuclear force), we can't build an infinite room. Instead, we build a tiny, finite "box" (a lattice) to simulate these collisions.

The standard rulebook for translating what happens in this tiny box back to the real, infinite world is called the Lüscher formalism. It's like a magic translator that says, "If the balls have these specific energy levels in the box, then in the real world, they bounce off each other with this specific strength."

However, there's a problem. The specific type of computer code used by many physicists (called Staggered Fermions) is incredibly fast and efficient, like a sports car. But it has a weird quirk: it's built on a grid that distorts reality slightly.

Here is the breakdown of the paper's challenges and solutions, using simple analogies:

The Problem: The "Taste" Confusion and the "Root" Trick

1. The "Taste" Multiplication (The Clone Army)
In this specific computer code, every time you try to simulate a single particle (like a pion), the math accidentally creates four copies of it. These aren't real copies; they are mathematical artifacts called "tastes."

• Analogy: Imagine you order one pizza, but the delivery guy brings you four slightly different pizzas because he got confused by the address. One is pepperoni, one is cheese, one is veggie, and one is burnt. In the real world, you only want one pizza. In the simulation, you have a "multiplet" of pizzas with slightly different weights (masses).

2. The "Fourth Root" Trick (The Magic Eraser)
To fix the "four pizzas" problem, physicists use a mathematical trick called "taking the fourth root." They essentially say, "We will mathematically pretend that only 1/4th of these pizzas exist."

• The Catch: While this works great for calculating simple things, it breaks a fundamental law of physics called Unitarity. In physics, Unitarity means "probability is conserved" (everything that happens adds up to 100%). The "fourth root" trick is like a magic eraser that fixes the pizza count but accidentally deletes the rule that says "what goes in must come out." It creates a theory that is slightly "unreal" at the scale of the computer grid.

The Goal: The authors want to figure out how to use this fast, slightly "broken" code to accurately predict how particles scatter, without waiting for the computer grid to become infinitely small (which takes forever and costs too much money).

---

The Solution: Two New Approaches

The paper proposes two ways to fix the translator (the Lüscher formalism) so it works with this "broken" code.

Approach 1: The "Recipe Book" (Effective Theory)

Instead of trying to fix the computer code itself, the authors wrote a new "recipe book" (called Rooted Staggered Chiral Perturbation Theory) that describes exactly how these "taste" particles and the "fourth root" trick behave.

• The Analogy: Imagine you are a chef trying to bake a cake using a slightly broken oven that heats unevenly. Instead of fixing the oven, you write a new cookbook that says, "If the oven is 5 degrees hot on the left and 10 degrees on the right, here is exactly how the cake will rise."
• What they did: They calculated the "one-loop" (a specific level of complexity) recipes for how these particles bounce. This allows them to check if the standard translation rules (Lüscher formalism) still work when you add these "oven quirks."

Approach 2: Rewriting the Rulebook (Generalizing the Formalism)

This is the bigger, more ambitious idea. The authors suggest we need to rewrite the "Magic Translator" (Lüscher formalism) to explicitly account for the broken nature of the code. They propose three specific changes:

1. Multiple Pathways (Topologies):

   • Normal Physics: Particles usually bounce off each other in one main way.
   • Staggered Physics: Because of the "taste" clones, particles can bounce off each other in many different, weird ways simultaneously.
   • The Fix: The new rulebook must allow for multiple scattering paths to happen at once, not just one.

2. The "Weight" Adjustment (Rescaling):

   • The Issue: Remember the "fourth root" trick? It means some paths in the simulation are "fake" or "partial."
   • The Fix: The new rulebook includes a "weight factor" (like a volume knob). If a path involves a "taste loop" (a clone interacting with itself), we turn the volume down by a factor of 1/4. This mathematically cancels out the extra copies and restores the correct physics.

3. The Multi-Channel System (The Orchestra):

   • The Issue: In the real world, a pion is just a pion. In this simulation, a "pion" can be a "taste-I pion," a "taste-V pion," etc. They can mix and match.
   • The Fix: Instead of treating the collision as a solo act, the new rulebook treats it like a full orchestra. We have to track how a "Taste-I" pair turns into a "Taste-V" pair, and so on. It turns a simple equation into a giant matrix (a grid of numbers) that tracks all these different "flavors" of pions.

Why Does This Matter?

Currently, many of the most precise calculations in particle physics (like the magnetic moment of the muon) use this fast "staggered" code. But if we want to study scattering (how particles crash and bounce), the standard tools fail because of the "taste" and "root" issues.

This paper is the first step in building a universal translator that can take data from these fast, slightly "broken" simulations and translate them into accurate, real-world physics predictions. It's like teaching a translator who speaks "Staggered-Grid" to finally understand "Real-World Physics" without needing to wait for the computer to become perfect.

In short: They are building a new set of instructions to fix the "glitches" in the fastest supercomputer simulations, so we can finally understand how the smallest building blocks of the universe crash into each other.

Technical Summary

1. Problem Statement

The extraction of hadronic scattering amplitudes from Lattice QCD (LQCD) typically relies on the Lüscher formalism, which relates finite-volume energy spectra to infinite-volume scattering amplitudes via a quantization condition. This formalism assumes unitarity and analyticity of the $S$-matrix.

However, many modern LQCD simulations utilize rooted staggered fermions due to their computational efficiency. This approach introduces specific theoretical challenges that violate the standard assumptions of the Lüscher formalism:

• Taste Splitting: Staggered fermions possess four "tastes" (doublers) per flavor. In the meson sector, this results in a multiplet of pions with non-degenerate masses ($O(a^2)$ artifacts), rather than a single physical pion.
• Fourth-Rooting: To recover the correct number of flavors in the continuum limit, the fermion determinant is raised to the power of $1/4$. This procedure breaks unitarity at non-zero lattice spacing ($a \neq 0$) because the rooted determinant cannot be associated with a local fermionic action.
• Limitation of Existing Methods: Previous extensions of the Lüscher formalism to include $O(a^2)$ discretization effects (e.g., Ref. [4]) apply to Wilson, overlap, and domain-wall fermions but do not account for the specific non-local and unitarity-violating nature of rooted staggered artifacts.

Consequently, applying the standard Lüscher formalism directly to staggered data at non-zero lattice spacing yields incorrect scattering amplitudes.

2. Methodology

The authors propose two complementary approaches to resolve these issues, focusing initially on $\pi\pi$ scattering in the isospin-2 ($I=2$) channel (where unitarity violation in the $s$-channel is absent) before generalizing to other channels.

Approach A: Rooted Staggered Chiral Perturbation Theory (rS$\chi$PT)

The authors perform the first-ever calculation of one-loop $\pi\pi$ scattering amplitudes using rS$\chi$PT.

• Framework: They utilize the Lee-Sharpe Lagrangian generalized to 3 flavors, incorporating the taste-breaking potential ($U$) and the fourth-rooting prescription.
• Calculation: They compute amplitudes for the $I=2$ channel ($\pi^+\pi^+ \to \pi^+\pi^+$) at Next-to-Leading Order (NLO).
• Diagrammatic Analysis: They explicitly analyze the topology of Feynman diagrams:
  • $s$-channel: Contains a single topology but includes contributions from different pion tastes.
  • $t$-channel: Contains multiple topologies, including those with internal taste loops (requiring a $1/4$ rooting factor) and disconnected "hairpin" diagrams.
• Goal: To derive explicit expressions for the amplitudes that include $O(a^2)$ corrections and to identify how taste splitting and rooting factors manifest in the perturbative expansion.

Approach B: Generalizing the Lüscher Formalism

The authors propose a modification of the Lüscher quantization condition to incorporate staggered artifacts directly. They adapt the "skeleton expansion" derivation (originally used for discretization effects) with three specific modifications:

1. Multiple $s$-channel Topologies: Unlike the continuum case where a single topology dominates, the staggered theory requires summing over multiple diagram topologies in the $s$-channel to account for different taste configurations.
2. Rescaling Factors ($\alpha_n$): To implement the fourth-rooting trick perturbatively, they introduce a rescaling factor $\alpha_n$ for each $s$-channel diagram.
   • If a diagram contains an internal taste loop, $\alpha_n = 1/4$.
   • Otherwise, $\alpha_n = 1$.
   • This effectively rescales the propagators in the skeleton expansion, breaking unitarity in a controlled manner consistent with the rooted theory.
3. Multichannel System: Since staggered pions exist in a multiplet of tastes ($I, A, V, P, T$), the scattering of two pions with taste $\xi$ can couple to states with different tastes $\xi'$. The formalism is extended to a multichannel scattering matrix ($M_{\xi \to \xi'}$) where total taste quantum numbers are conserved.

3. Key Contributions and Results

• First rS$\chi$PT Amplitudes: The paper provides the first explicit one-loop expressions for $\pi\pi$ scattering amplitudes in the $I=2$ channel within the rooted staggered framework. These expressions include:
  • Dependence on taste-breaking Low-Energy Constants (LECs).
  • Explicit $O(a^2)$ corrections arising from single-trace and double-trace operators.
  • The distinction between connected and disconnected (hairpin) propagators.
• Modified Quantization Condition: The authors derive a generalized quantization condition:
  $$\det \left[ \mathcal{M}^{-1}(E) - \mathcal{F}_{\text{tot}}(E, L) \right] = 0$$
  Where $\mathcal{F}_{\text{tot}}$ is a sum over all $N$ $s$-channel topologies, weighted by the rooting factors $\alpha_n$. This condition allows for the extraction of amplitudes directly from staggered lattice data without first taking the continuum limit.
• Handling Unitarity Violation: The formalism explicitly incorporates the non-unitary nature of the rooted theory by allowing the rescaling factors $\alpha_n$ to differ from unity, thereby modeling the specific way unitarity is broken at $O(a^2)$.
• Multichannel Matrix Structure: The authors define the full scattering matrix for the 5 taste channels, highlighting the increased complexity (number of unknowns) compared to single-channel scattering, which necessitates the use of model-dependent Ansätze or large datasets to constrain.

4. Significance and Future Outlook

• Bridging the Gap: This work provides the necessary theoretical framework to utilize the vast amount of existing staggered lattice data (e.g., from the FNAL/MILC HISQ ensembles) for scattering studies. Currently, such data is often discarded for scattering analyses due to the lack of a formalism to handle taste splitting and rooting.
• Validation Strategy: The authors plan to use their rS$\chi$PT results to explicitly test the validity of their proposed quantization condition. By comparing the energy levels predicted by the modified formalism against the amplitudes calculated in rS$\chi$PT, they can verify the consistency of the approach.
• Continuum Limit Strategy: An alternative future path involves calculating energy levels at fixed volume, taking the continuum limit ($a \to 0$) where artifacts vanish and unitarity is restored, and then applying the standard Lüscher formalism. However, the proposed direct approach allows for analysis at non-zero lattice spacing, potentially reducing computational costs.
• Broader Impact: Successfully implementing this formalism will enable precise determinations of scattering lengths, resonance properties, and other hadronic observables using the most computationally efficient fermion action available, significantly advancing the field of Lattice QCD phenomenology.