An Update on the Isospin-Breaking Effects in the Pion Decay Constant with Staggered Quarks

This paper presents an update from the BMW Collaboration on their ongoing calculation of isospin-breaking effects in the pion decay constant using $N_f=2+1+1$ staggered quarks with near-physical pion masses and QED$_{\text{L}}$, including preliminary results for the axial-pseudoscalar correlator and future plans.

The Gist

Imagine you are trying to bake the perfect cake, but you need to know the exact weight of a single grain of sugar to make sure your recipe is perfect. In the world of particle physics, scientists are trying to do something similar: they are trying to measure the "weight" (or decay rate) of a tiny particle called a pion with extreme precision.

Why? Because this measurement helps solve a massive mystery in the universe: Why does the math of the universe sometimes seem to add up to 99% instead of 100%?

This paper is an update from a team of scientists (the BMW collaboration) who are building a super-accurate "kitchen scale" using giant supercomputers to weigh these pions. Here is the breakdown of their work in simple terms:

1. The Big Mystery: The "Missing" 1%

In the Standard Model of physics, there is a rulebook called the CKM matrix. It's like a ledger that tracks how particles change into one another. For the universe to be stable and logical, the numbers in the first row of this ledger must add up to exactly 1.0.

However, when scientists measure these numbers today, they add up to about 0.99. It's a tiny gap, but in physics, a missing 1% is like a hole in a dam—it suggests we are missing a piece of the puzzle. One of the main suspects for this missing piece is something called Isospin Breaking.

2. The "Twin" Problem: Isospin Breaking

Imagine you have a pair of identical twins. In a perfect world, they would be exactly the same. In the world of subatomic particles, the "up" quark and the "down" quark are like these twins. They are almost identical, but not quite.

• They have slightly different masses.
• They have different electric charges (one is positive, one is neutral-ish).

This tiny difference is called Isospin Breaking. Usually, scientists pretend the twins are identical to make the math easier. But to solve the "missing 1%" mystery, they can no longer ignore the differences. They have to calculate exactly how much the "charge" and "mass" differences mess up the pion's weight.

3. The Super-Computer Kitchen

To do this, the team uses Lattice QCD. Imagine space-time isn't a smooth sheet, but a giant 3D grid (like a giant egg carton). They place particles on the grid and simulate how they interact.

• The Staggered Quarks: This is the specific type of "flour" they use for their simulation. It's a sophisticated way of arranging the particles on the grid to keep the calculations efficient.
• The BMW Collaboration: This is the team of chefs (scientists) running the simulation. They are famous for making very precise measurements.

4. The Two Ingredients: Sea and Valence

When calculating the pion's weight, the scientists have to account for two types of "ingredients" inside the particle:

• The Valence Quarks (The Main Ingredients): These are the two quarks that make up the pion (like the flour and sugar in a cake). The team is calculating how the electric charge of these specific quarks changes the weight.
• The Sea Quarks (The Background Noise): Inside the pion, there is a chaotic "soup" of virtual particles popping in and out of existence. These are the "sea" quarks. The team also calculated how these background particles affect the weight.

The Analogy:
Think of the pion as a boat.

• The Valence quarks are the heavy cargo on the deck.
• The Sea quarks are the water splashing against the hull.
• Isospin Breaking is the fact that the cargo is slightly uneven, and the water is slightly salty on one side.
  The team is calculating exactly how much that unevenness and saltiness makes the boat sit lower in the water (changing the decay constant).

5. What They Found So Far

The paper presents a "preliminary" result. They have successfully built a new, highly accurate scale.

• The Result: They have calculated a value for the scale setting ($w_0$) which is crucial for converting their computer numbers into real-world units (like femtometers).
• The Precision: They are getting very close. The error margin is tiny (about 0.2%), but they need to get even smaller to solve the CKM matrix mystery.
• The Bottleneck: The biggest source of error right now is how they handle the "grid size" (continuum extrapolation). It's like trying to measure a curve with a ruler; if the ruler's marks are too far apart, your measurement is rough. They are currently running simulations on finer grids (smaller ruler marks) to get a sharper picture.

6. The Next Steps

The team is not done yet. They are:

1. Refining the Grid: Running simulations on even finer grids to reduce the "ruler" error.
2. Checking the Volume: Making sure their simulation box is big enough so the "boat" doesn't hit the walls of the virtual universe (finite volume effects).
3. Independent Checks: They are cross-checking their "Valence" calculations with another group's methods to ensure they aren't making a mistake.

The Bottom Line

This paper is a progress report from a team of physicists who are building the most precise "kitchen scale" in the universe. They are carefully accounting for the tiny differences between particle "twins" (isospin breaking) to ensure their measurements are perfect.

If they succeed, they might finally explain why the universe's math doesn't add up perfectly, potentially revealing new physics beyond what we currently know. It's a bit like finding a tiny crack in a perfect wall that suggests there's a whole new room behind it.

Technical Summary

1. Problem Statement

The primary motivation for this work is the unitarity violation of the first row of the CKM matrix ($|V_{ud}|^2 + |V_{us}|^2 + |V_{ub}|^2 = 1$), which remains an open problem in precision flavor physics.

• Context: At current levels of precision, Isospin-Breaking Effects (IBE) arising from the difference in up/down quark masses and electromagnetic interactions (QED) become significant.
• Specific Goal: To compute the pion decay constant ($F_\pi$) in full QCD+QED with high precision. This constant is crucial for:
  1. Determining the CKM element $|V_{ud}|$ from nuclear $\beta$-decay.
  2. Setting the physical scale (via the gradient-flow scale $w_0$) for the BMW collaboration's calculation of the muon anomalous magnetic moment ($a_\mu$).
• Challenge: While isosymmetric QCD calculations are mature, the full QCD+QED calculation requires separating and computing sea-quark and valence-quark electromagnetic contributions, including non-factorizable diagrams and finite-volume effects.

2. Methodology

The calculation employs $N_f = 2+1+1$ Staggered Quarks with near-physical pion masses. The methodology combines results from the BMW collaboration (isosymmetric QCD + sea QED) and the RM123-Southampton collaboration (valence QED).

A. Mixed Determination Strategy

The authors utilize a hybrid approach to separate the contributions:

1. Parametrization: The pion decay constant ($F_\pi$) scaled by the gradient flow scale ($w_0$) is parametrized as a function of quark masses and electric charges ($e_v, e_s$):
   $$w_0 F_\pi = A + B F_\pi^{-2} M_\pi^2 + C F_\pi^{-2} (\dots) + E e_v^2 + F e_v e_s + G e_s^2$$
   • Coefficients $B, C$ are mass derivatives computed in QCD.
   • Coefficients $E, F, G$ represent valence-valence, sea-valence, and sea-sea electromagnetic derivatives.
2. Scheme Consistency: The calculation combines the BMW isosymmetric QCD results with the RM123 valence QED results. To ensure consistency, the authors map the RM123 results (defined in the GRS scheme) to the BMW framework.
   • The relation used is: $[w_0 F_\pi]_{\text{QCD+QED}} = [w_0 F_\pi]_{\text{QCD+seaQED}} \sqrt{1 + [\delta R_\pi]_{\text{GRS}}}$.
   • Here, $\delta R_\pi$ represents the isospin-breaking correction to the decay rate.

B. Computational Details

• Ensembles: Simulations use $N_f=2+1+1$ Staggered fermions with 4 levels of stout smearing (and 4 levels of HEX smearing for finer lattices).
• QED Treatment:
  • Sea Quarks: Computed using QED$_L$ (finite volume QED with zero-mode removal) in Coulomb gauge.
  • Valence Quarks: Computed using sequential propagators with 16 photon sources and 3 U(1) random sources.
  • Finite Volume: The study explicitly addresses finite-volume effects, noting that QED$_L$ introduces power-law corrections of order $O(1/L^2)$.
• Scale Setting: The scale is set using the gradient flow scale $w_0$, computed via both Wilson-flow and Zeuthen-flow to control continuum extrapolation systematics.

3. Key Contributions

1. Isosymmetric QCD Value: A high-precision determination of the isosymmetric quantity $w_0 f_\pi$ in the FLAG scheme, achieving 0.23% precision. The error budget is dominated by the continuum extrapolation.
2. Sea Quark IBE: A direct computation of the sea-sea electromagnetic derivative using BMW ensembles. The result shows the sea quark IBE contributes approximately 0.1% to the total $w_0 F_\pi$ value.
3. Valence Quark IBE (Preliminary):
   • Presentation of preliminary results for the factorizable valence-valence diagrams.
   • Demonstration of clear plateaus in the effective derivative at large time separations, indicating low excited-state contamination.
   • Analysis of volume dependence across four different volumes ($L \approx 3.2$ to $8.4$ fm) and two lattice spacings ($a \approx 0.13$ fm and $0.095$ fm).
4. Negligible Terms: The sea-valence IBE is estimated to be $\sim 20\%$ of the sea-sea term and is currently neglected in the final error budget.

4. Results

The paper presents a preliminary value for the pion decay constant in full QCD+QED and the derived scale $w_0$:

• Pion Decay Constant ($w_0 F_\pi$):
  $$[w_0 F_\pi]_{\text{QCD+QED}} = 0.11527(15)_{\text{stat}}(26)_{\text{syst}}[31]_{\text{total}}$$

  • The systematic error is dominated by the QCD continuum extrapolation (23) and the valence QED contribution (11).

• Gradient Flow Scale ($w_0$):
  Combining the above with the experimental value of $F_\pi = 131.711(45)$ MeV:
  $$[w_0]_{\text{QCD+QED}} = 0.17270(22)_{\text{stat}}(40)_{\text{syst}}[46]_{\text{total}} \text{ fm}$$

  • Error Breakdown: The total uncertainty is 0.27%. The dominant systematic error (0.20%) comes from the QCD continuum extrapolation. The QED contribution to the error is smaller but non-negligible, driven by the valence QED determination.

5. Significance and Future Outlook

• Impact on Muon $g-2$: This result provides a more rigorous scale setting ($w_0$) for the BMW collaboration's calculation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment, reducing theoretical uncertainties in the Standard Model prediction.
• CKM Unitarity: By refining the determination of $F_\pi$ with full QED effects, the calculation aids in resolving the tension in the first row of the CKM matrix.
• Next Steps:
  • Finer Lattices: Generating data at finer lattice spacings to improve the continuum extrapolation of the isosymmetric part.
  • Independent Valence Calculation: Performing a fully independent determination of valence effects, including the non-factorizable diagrams (involving leptons) and the renormalization of the electroweak Hamiltonian.
  • Finite Volume Study: Conducting a detailed study of finite-volume effects using volumes up to 10.8 fm to address the poor convergence of the QED$_L$ series observed in previous works.

In summary, this work represents a critical step toward a sub-percent precision determination of the pion decay constant in full QCD+QED, bridging the gap between pure QCD simulations and the physical reality of electromagnetic interactions.