$f_K/f_{\pi}$ in iso-symmetric QCD and the CKM matrix unitarity

This paper presents precise lattice QCD calculations of the ratio $f_K/f_{\pi}$ in the iso-symmetric limit using a combined Wilson unitary and mixed-action approach, which are subsequently utilized to determine $|V_{us}|/|V_{ud}|$ and test the unitarity of the first row of the CKM matrix after incorporating strong isospin-breaking and QED effects.

The Gist

Imagine the universe is built from a giant, invisible LEGO set. The smallest bricks in this set are called quarks. They snap together to form larger structures called protons and neutrons, which make up the atoms in everything around us.

However, these LEGO bricks don't just snap together randomly; they follow strict rules written in a "rulebook" called the Standard Model of physics. One of the most important pages in this rulebook is the CKM Matrix. Think of this matrix as a "translation table" or a "passport system" that tells us how likely it is for one type of quark to transform into another.

For this passport system to work correctly, the numbers in the first row of the table must add up to exactly 1. If they add up to 0.99 or 1.01, it means our rulebook is missing a page, or there's a secret new type of physics hiding in the shadows (what scientists call "New Physics").

The Problem: Measuring the "Translation"

To check if the numbers add up to 1, scientists need to measure how often a specific quark (the "strange" quark) turns into a "up" quark. This happens when a particle called a Kaon decays (falls apart) into a Pion.

The rate of this decay depends on two things:

1. The "passport" number (the CKM value we want to find).
2. The "stiffness" or "strength" of the particles involved. In physics, this is called the decay constant ($f_K$ for Kaons and $f_\pi$ for Pions).

To get the passport number, you have to divide the Kaon's strength by the Pion's strength. But here's the catch: Kaons and Pions are not simple objects. They are fuzzy clouds of quarks and gluons (the "glue" holding them together) that are constantly vibrating. You can't just weigh them on a scale; you have to simulate the entire universe inside a supercomputer to figure out their strength.

The Solution: A Digital Laboratory

This paper is about a team of scientists who built a digital laboratory (using Lattice QCD) to simulate these particles.

The Analogy of the Two Maps:
Imagine you are trying to draw a map of a city.

• Method A (Wilson Unitary): You draw the map using a grid of square tiles. It's accurate, but the corners of the tiles might make the roads look a little jagged.
• Method B (Mixed Action): You draw the map using square tiles for the background, but you use smooth, round stickers for the cars (the valence quarks) driving on them.

The authors of this paper used both methods simultaneously. They ran the simulation twice: once with Method A and once with Method B. By comparing the two, they could see where the "jagged corners" of the digital grid were distorting the picture. This allowed them to smooth out the errors and get a much more precise measurement of the "stiffness" ratio ($f_K/f_\pi$).

The Journey to the "Physical Point"

The simulations were run with different settings, like changing the size of the LEGO bricks or the weight of the particles. But the real world has specific weights for Pions and Kaons. The team had to use a mathematical "ramp" (called Chiral Perturbation Theory) to slide their simulation results from their "fake" settings to the Physical Point (the real-world values).

They also had to account for tiny "glitches" in their simulation:

• Finite Volume: Their digital universe was a box, not infinite. They had to calculate how the walls of the box squeezed the particles.
• Isospin Breaking: In the real world, the "up" and "down" quarks have slightly different weights, and electricity (QED) plays a role. The team calculated how to adjust their pure simulation to include these messy real-world effects.

The Big Result

After all this work, they calculated the ratio of the Kaon's strength to the Pion's strength with incredible precision:
Ratio = 1.1848 (with a tiny margin of error).

Using this number, they calculated the CKM "passport" value for the strange quark ($|V_{us}|$). When they added this to the known value for the up quark ($|V_{ud}|$), the sum was:
$|V_{ud}|^2 + |V_{us}|^2 = 0.9995$

What Does This Mean?

The result is 0.9995, which is incredibly close to 1.0.

• The Good News: The Standard Model is holding up! The "passport system" is working correctly. There is no massive evidence yet for "New Physics" hiding in this specific corner of the universe.
• The Catch: The uncertainty (the margin of error) in their measurement is still dominated by the difficulty of simulating the particles themselves. It's like trying to measure the height of a mountain with a ruler that has a slightly wobbly edge.

The Takeaway

This paper is a triumph of precision engineering in theoretical physics. The authors didn't just measure a number; they built a dual-lens camera (using two different simulation methods) to get the sharpest possible picture of how quarks interact.

While they didn't find a "new universe" today, they tightened the screws on our current understanding. If there is a new particle or force hiding in the shadows, it's going to be very hard to find because the "old physics" is now measured with such extreme accuracy. The next step is to build even bigger digital universes and run the simulations longer to shrink that tiny margin of error even further.

Technical Summary

1. Problem Statement

The primary objective of this work is to provide a high-precision determination of the ratio of kaon to pion decay constants, $f_K/f_\pi$, within the iso-symmetric limit of pure QCD (isoQCD) with $N_f = 2+1$ flavors. This quantity is critical for testing the Standard Model (SM), specifically the unitarity of the first row of the Cabibbo-Kobayashi-Maskawa (CKM) matrix:
$$|V_{ud}|^2 + |V_{us}|^2 + |V_{ub}|^2 = 1$$
While $|V_{ud}|$ is known with high precision from super-allowed nuclear $\beta$-decays and $|V_{ub}|$ is negligible, $|V_{us}|$ is extracted from leptonic meson decays ($K \to \ell \nu$ and $\pi \to \ell \nu$). The extraction relies on the experimental ratio of decay rates combined with the theoretical ratio $f_K/f_\pi$. A discrepancy in this unitarity test could signal New Physics. The challenge lies in reducing the lattice QCD uncertainty of $f_K/f_\pi$, which currently dominates the error budget in CKM unitarity tests.

2. Methodology

Lattice Setup and Ensembles

• Ensembles: The authors utilize $N_f = 2+1$ CLS (Coordinated Lattice Simulations) ensembles generated with the Lüscher-Weisz gauge action and non-perturbatively $O(a)$-improved Wilson fermions.
• Regularizations: A unique feature of this study is the combination of two distinct lattice regularizations to control continuum extrapolation:
  1. Wilson Unitary: Wilson fermions in both the sea and valence sectors.
  2. Mixed-Action: Wilson fermions in the sea and Wilson Twisted Mass (Wtm) fermions in the valence sector. The valence parameters are tuned to match sea and valence pion/kaon masses, ensuring unitarity in the continuum limit.
• Scale Setting: Unlike many studies that use gradient flow scales ($t_0, w_0$), this work sets the scale using the physical pion decay constant ($f_\pi$). This avoids systematic uncertainties associated with the continuum extrapolation of the scale itself, as the physical value of $f_\pi$ is fixed by the "Edinburgh Consensus" and FLAG 24.
• Physical Point: The chiral interpolation is performed to the physical point defined by:
  • $m_\pi^{ph} = 135$ MeV
  • $m_K^{ph} = 494.6$ MeV
  • $f_\pi^{ph} = 130.5$ MeV

Analysis Strategy

• Finite Volume Corrections: Corrections are applied using Next-to-Leading Order (NLO) Chiral Perturbation Theory ($\chi$PT). The authors find these corrections are subleading (few per mille) compared to statistical errors.
• Chiral-Continuum Extrapolation:
  • Data is interpolated using $\chi$PT-inspired formulae (NLO) and Taylor expansions in terms of dimensionless variables $\rho_2 = m_\pi^2/f_\pi^2$ and $\rho_4 = (m_K^2 + m_\pi^2/2)/f_\pi^2$.
  • Cutoff effects ($O(a^2)$) are modeled by expanding Low Energy Constants (LECs) with terms proportional to $(a f_\pi^{sym})^2$.
  • Model Averaging: To assess systematic uncertainty, the authors perform a model average over various fit functions (different $\chi$PT forms, Taylor expansions), renormalization scales ($\mu = f_\pi$ vs. $f_{\pi K}$), and data cuts (removing coarsest lattices or heavy pions). The final error includes the spread between these models.
• Isospin and QED Corrections:
  • Strong isospin-breaking ($\delta_{SU(2)}$) and QED effects ($\delta_{EM}$) are added after the isoQCD calculation.
  • The authors adopt a conservative approach for these corrections, assigning large uncertainties to higher-order terms, noting that the final result is still dominated by the isoQCD lattice uncertainty.

3. Key Contributions

1. Dual Regularization Approach: The simultaneous use of Wilson unitary and Mixed-Action (Wtm) setups allows for a robust control of the continuum limit extrapolation. The agreement between these two distinct discretizations validates the removal of lattice artifacts.
2. Scale Setting via $f_\pi$: By using $f_\pi$ to set the scale, the authors eliminate the systematic uncertainty inherent in extrapolating theoretical scales (like $t_0$) to the continuum, a common source of error in other lattice determinations.
3. Comprehensive Systematic Error Analysis: The paper employs a rigorous model-averaging technique to quantify the uncertainty arising from the choice of chiral fit functions and data cuts, rather than relying on a single "best" fit.
4. CKM Unitarity Test: The work provides a complete pipeline from raw lattice data to a test of the CKM unitarity constraint, including all necessary isospin-breaking and QED corrections.

4. Key Results

Decay Constant Ratio

The final result for the ratio in isoQCD (pure QCD, isospin symmetric) is:
$$\frac{f_K}{f_\pi} = 1.1872(59)_{\text{stat}}(84)_{\chi-\text{cont}}$$
Total uncertainty: 1.03% (dominated by the chiral-continuum extrapolation systematic).

Including Isospin and QED

After applying strong isospin-breaking ($\delta_{SU(2)} = -0.004(4)$) and QED corrections ($\delta_{EM} = -0.0072(14)$), the physical ratio becomes:
$$\frac{f_{K^\pm}}{f_{\pi^\pm}} = 1.1848(59)_{\text{stat}}(84)_{\chi-\text{cont}}(24)_{SU(2)}$$

CKM Matrix Elements

Using the experimental decay rate ratio and the calculated $f_K/f_\pi$:

• Ratio: $|V_{us}|/|V_{ud}| = 0.2330(11)_{\text{stat}}(17)_{\chi-\text{cont}}(5)_{SU(2)}(2)_{\text{EM}}(3)_{\text{exp}}$
• Absolute Value: $|V_{us}| = 0.2269(13)_{\text{stat}}(15)_{\chi-\text{cont}}(4)_{SU(2)}(2)_{\text{EM}}(2)_{\beta-\text{dec}}(2)_{\text{exp}}$

Unitarity Test

The sum of the squares of the first row elements is:
$$|V_{ud}|^2 + |V_{us}|^2 = 0.9995(6)_{\text{stat}}(7)_{\chi-\text{cont}}(2)_{SU(2)}(1)_{\text{EM}}(7)_{\beta-\text{dec}}(1)_{\text{exp}}$$
Result: The result is $0.9995(11)$, which is consistent with the Standard Model unitarity constraint of 1.0 within approximately 0.5$\sigma$.

5. Significance and Conclusion

• Consistency with SM: The result confirms the unitarity of the CKM matrix to high precision, finding no evidence for New Physics in the first row.
• Error Budget: The study highlights that the uncertainty in $|V_{us}|$ is currently completely dominated by the lattice determination of $f_K/f_\pi$ (specifically the chiral-continuum extrapolation), rather than by the isospin-breaking or QED corrections.
• Future Outlook: The authors identify that further improvements require:
  1. Increased statistics on ensembles close to the physical point (specifically ensemble E250).
  2. Simulations at finer lattice spacings ($a < 0.04$ fm) at the physical point.
  3. Inclusion of Next-to-Next-to-Leading Order (NNLO) $\chi$PT terms to better constrain the chiral dependence.
• Benchmark: The value $f_K/f_\pi = 1.1872(103)$ serves as a precise, unambiguous benchmark for other lattice groups working in isoQCD, provided they adopt the same physical point definition.

In summary, this paper represents a state-of-the-art lattice QCD calculation that significantly tightens the constraints on CKM unitarity, demonstrating that the current tension (if any) is not driven by strong interaction uncertainties but is limited by statistical and extrapolation systematics in the lattice data itself.