A Precision Test of First Row CKM Unitarity from Lattice QCD

This paper reviews the Fermilab Lattice and MILC collaborations' efforts to perform a correlated lattice QCD analysis using Highly Improved Staggered Quarks and Staggered Chiral Perturbation Theory to precisely determine the CKM matrix elements $|V_{ud}|$ and $|V_{us}|$ without nuclear inputs, thereby testing the first-row unitarity of the CKM matrix and probing physics beyond the Standard Model.

The Gist

Imagine the Standard Model of physics as a giant, incredibly complex jigsaw puzzle that scientists have been assembling for decades. This puzzle explains how the universe works at its most fundamental level. One of the most important pieces in this puzzle is a set of numbers called the CKM matrix. Think of these numbers as a "rulebook" that dictates how different types of particles (quarks) can change into one another.

According to the rules of the puzzle, if you add up the squares of the first three numbers in this rulebook, they should equal exactly 1. It's like a mathematical law of conservation: the total probability of all possible outcomes must add up to 100%.

The Problem: A Missing Piece

Recently, scientists have measured these numbers with extreme precision, and they found a problem. When they add up the first two numbers (which represent the most common particle changes), the total is slightly less than 1. It's like trying to close a suitcase, but there's a tiny gap left over.

This gap is called the Cabibbo Angle Anomaly. It's a big deal because:

1. It might mean the puzzle is broken: The "rulebook" might be wrong, suggesting there are new, undiscovered particles or forces (Physics Beyond the Standard Model) that we haven't found yet.
2. It might mean our measurements are slightly off: We might just be bad at measuring the pieces.

The Challenge: The "Nuclear Noise"

To fix this, scientists need to measure these numbers again, but with a twist. The most common way to measure one of these numbers involves looking at atomic nuclei (the core of atoms). However, nuclei are messy, chaotic places. It's like trying to hear a whisper in a crowded, noisy stadium. The "noise" from the nuclear structure makes it hard to know if the gap is real or just a measurement error.

The goal of this paper is to measure these numbers without looking at the noisy stadium. They want to listen to the "whisper" in a soundproof room.

The Solution: A Supercomputer Kitchen

The authors (from Fermilab and MILC collaborations) are using a technique called Lattice QCD. Imagine the universe not as a smooth fabric, but as a giant 3D grid (like a chessboard made of space and time). They use supercomputers to simulate how particles move on this grid.

Instead of measuring real atoms, they create "virtual atoms" on the computer. They calculate two specific things:

1. How fast pions and kaons decay: Think of these particles as unstable soap bubbles. The authors calculate exactly how long it takes for them to pop (decay) and what they turn into.
2. The "Decay Constants": These are like the "size" or "strength" of the bubbles.

By combining the experimental data (how fast they actually pop in the real world) with their supercomputer calculations (how big the bubbles are), they can extract the rulebook numbers ($|V_{ud}|$ and $|V_{us}|$) without ever touching a messy atomic nucleus.

The New Strategy: Cooking a Two-Part Meal

In the past, scientists calculated these two numbers separately, like cooking a soup and a salad in different pots without talking to each other. This meant they couldn't see if a mistake in the soup affected the salad.

This paper describes a new, smarter strategy:

• The Correlated Kitchen: They are now cooking both the "soup" (decay constants) and the "salad" (form factors) in the same pot, using the same ingredients and the same recipe.
• The Recipe (SChPT): They use a mathematical tool called "Staggered Chiral Perturbation Theory." Think of this as a master recipe book that helps them adjust for the fact that their computer grid isn't perfect (it's a bit "pixelated") and that they sometimes have to simulate particles with the wrong weight (like using heavy flour instead of light flour).
• The "Priors": They use the results from the first part of the calculation to help guide the second part. It's like using the taste of the soup to adjust the seasoning of the salad, ensuring both dishes taste perfect and consistent.

Why This Matters

By doing this "correlated analysis," they are reducing the uncertainty. They are essentially saying, "We know exactly how these two numbers relate to each other, so we can be much more confident in the final result."

If their new, ultra-precise numbers still show a gap (a deficit) in the rulebook, it will be a very strong signal that new physics exists. It would be like finding a crack in the foundation of a building that we thought was perfectly solid. This would force physicists to rewrite the rulebook and discover new particles or forces that explain the missing piece.

In short: These scientists are using supercomputers to simulate the universe's most basic rules, cleaning up the "noise" from real-world experiments, and checking if the universe's math actually adds up. If it doesn't, we might be on the verge of a massive discovery.

Technical Summary

1. Problem Statement

The central issue addressed is the Cabibbo Angle Anomaly (CAA), a persistent deficit in the unitarity relation of the first row of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In the Standard Model (SM), the CKM matrix must be unitary, implying:
$$\Delta \equiv |V_{ud}|^2 + |V_{us}|^2 + |V_{ub}|^2 - 1 = 0$$
Current experimental and theoretical determinations show a 2–3$\sigma$ deficit ($\Delta \neq 0$). Since $|V_{ub}|$ is negligible at current precision levels, the anomaly is driven by the values of $|V_{ud}|$ and $|V_{us}|$.

• Current Limitations: The most precise determination of $|V_{ud}|$ relies on superallowed nuclear beta decays, which are dominated by theoretical uncertainties related to nuclear structure. Similarly, $|V_{us}|$ is often extracted from semileptonic kaon decays ($K_{\ell3}$), requiring non-perturbative inputs.
• The Goal: To perform a rigorous test of first-row unitarity without relying on nuclear inputs. This requires a precise, correlated determination of:
  1. The ratio of decay constants $f_K/f_\pi$ (from leptonic decays $K_{\ell2}/\pi_{\ell2}$).
  2. The semileptonic form factor $f_+(0)$ (from $K_{\ell3}$ decays).
  3. The correlation between these two lattice QCD inputs to minimize systematic errors in global fits.

2. Methodology

The Fermilab Lattice and MILC collaborations employ Lattice QCD to compute the necessary non-perturbative inputs from first principles.

• Ensembles: The analysis utilizes $N_f = 2 + 1 + 1$ gauge-field configurations generated by the MILC collaboration, which include up, down, strange, and charm quarks in the sea. These ensembles cover multiple lattice spacings, quark masses, and physical volumes.
• Chiral-Continuum Extrapolation: To extrapolate lattice data to the physical point (continuum limit, infinite volume, and physical quark masses), the authors use Staggered Chiral Perturbation Theory (SChPT).
  • The fit function is constructed from Next-to-Leading Order (NLO) SChPT, Next-to-Next-to-Leading Order (NNLO) ChPT, and additional terms parameterizing higher-order discretization effects.
  • This approach allows for the systematic inclusion of discretization errors, chiral effects, finite volume corrections, and strong-isospin breaking.
• Two-Step Analysis Strategy for $f_K/f_\pi$:
  1. Step 1: Restrict the analysis to ensembles with lighter-than-physical strange quark masses where SChPT converges well. This yields a correlated distribution of SChPT Low-Energy Constants (LECs).
  2. Step 2: Use the LEC distributions from Step 1 as priors in a global fit including all available data (including heavier masses). Additional higher-order terms are added to the fit function to account for convergence issues at the strange quark mass.
• Correlated Analysis: The same SChPT LECs appear in the chiral expressions for both the decay constants ($f_K/f_\pi$) and the semileptonic form factor ($f_+(0)$). The authors use the posterior distributions of LECs from the leptonic analysis as priors for the semileptonic analysis to propagate correlations.
• Statistical Treatment: A resampling analysis is performed using synthetic samples drawn from the data covariance matrix to estimate parameter correlations.

3. Key Contributions

• First Correlated Analysis: This work presents the first effort to simultaneously determine $f_K/f_\pi$ and $f_+(0)$ with a rigorous treatment of their statistical correlations using a unified SChPT framework.
• Improved Precision via SChPT: By moving away from a reliance solely on physical mass ensembles (as in previous work) to a full SChPT-based fit function, the collaboration can include unphysical quark mass data. This significantly improves statistical precision and allows for better control of systematic uncertainties.
• Robust LEC Determination: The two-step strategy allows for a precise determination of SChPT Low-Energy Constants (LECs) even in regimes where ChPT convergence is poor (near the physical strange mass), by constraining them with data from the lighter mass region.
• Propagation of Correlations: The methodology successfully propagates correlations between the leptonic and semileptonic sectors, demonstrating non-zero correlations between $f_+(0)$ and the LECs.

4. Preliminary Results

• LECs: The analysis yields precise determinations of specific LEC combinations with strong correlations between them (visualized in the paper's correlation matrix). The central values and correlations remain consistent between the two-step analysis phases.
• $f_K/f_\pi$: A precise value for the ratio of decay constants is obtained, with systematic uncertainties better controlled than in previous iterations.
• $f_+(0)$: The semileptonic form factor at zero momentum transfer is determined with improved precision. Crucially, the analysis confirms non-zero correlations between $f_+(0)$ and the LECs, validating the resampling strategy.
• Status: While the central values are robust, the authors note that the increased number of parameters in the fit function makes the estimation of correlations in the resampling analysis more challenging and is an ongoing area of optimization.

5. Significance and Outlook

• Testing New Physics: The correlated inputs ($f_K/f_\pi$ and $f_+(0)$) are essential for testing the SM without nuclear model dependencies. A precise, correlated determination reduces the uncertainty in the unitarity test, sharpening the search for Beyond the Standard Model (BSM) physics.
• SMEFT Implications: Recent global fits within the Standard Model Effective Field Theory (SMEFT) suggest that tensions in the first row of the CKM matrix could be resolved by operators modifying right-handed charged currents (e.g., in vector-like quark theories). The sensitivity of these global fits to new physics is directly dependent on the precision and correlation structure of the SM inputs.
• Future Impact: By providing a rigorous, correlated lattice input, this work enhances the ability to distinguish between statistical fluctuations and genuine signs of new physics in the flavor sector. It sets the stage for a definitive resolution of the Cabibbo Angle Anomaly.