Kaon leptonic and semileptonic decays with $N_f=2+1+1$ HISQ fermions

This paper presents progress in a correlated lattice-QCD analysis using $N_f=2+1+1$ HISQ fermions to address the CKM unitarity deficit by updating light-meson decay constants via staggered chiral perturbation theory and reanalyzing kaon semileptonic form factors to estimate their correlations.

The Gist

Imagine the Standard Model of physics as a giant, incredibly complex puzzle. For decades, scientists have been trying to fit all the pieces together to see if the picture is perfect. Recently, however, they noticed a tiny crack in the frame: the "first row" of the puzzle isn't quite adding up. This is known as the Cabibbo-Kobayashi-Maskawa (CKM) unitarity deficit, or more casually, the "Cabibbo Angle Anomaly."

This paper is a report from a team of scientists (the Fermilab Lattice and MILC Collaborations) who are trying to fix that crack by re-measuring two specific puzzle pieces with extreme precision.

Here is a breakdown of what they are doing, using simple analogies:

1. The Problem: A Wobbly Scale

In the world of subatomic particles, there are rules about how different types of "quarks" (the building blocks of matter) change into one another. One rule says that if you add up the probabilities of these changes in a specific way, the total must equal exactly 1.

Currently, when scientists measure this, the total comes up slightly short (like a scale that reads 0.99 instead of 1.00). This suggests either:

• There is a mistake in the measurements.
• There is a hidden piece of physics we don't understand yet (New Physics).

To know for sure, they need to measure the "ingredients" of the calculation with perfect precision. The two main ingredients they are re-measuring are:

1. How often a Kaon turns into a Pion (a specific type of particle decay).
2. How "heavy" or "strong" the Kaon and Pion are (their decay constants).

2. The Tool: A Digital Time Machine

You can't just put a Kaon on a scale in a lab; it exists for a fraction of a second and is too small to touch. Instead, these scientists use Lattice QCD (Quantum Chromodynamics).

Think of this as a giant 3D grid (like a digital chessboard) that represents space and time. They simulate the laws of physics on this grid, creating virtual universes where they can watch virtual Kaons and Pions interact.

• The "HISQ" Method: This is the specific "recipe" they use to simulate the quarks. It's like upgrading from a blurry, low-resolution camera to a 4K ultra-high-definition camera. It reduces the "pixelation" (errors) in their simulation.
• The "MILC" Configurations: These are the specific digital universes they built. They have created new, larger, and more detailed universes than before, including some that mimic the exact real-world conditions (physical quark masses) and others that are slightly different to test the rules of the universe.

3. The Strategy: The "Two-Step" Dance

The team is tackling this in two main phases, which they call a "correlated analysis."

Phase A: Re-measuring the "Shape" (The Form Factor)

They are looking at how a Kaon morphs into a Pion. Imagine a clay sculpture of a Kaon slowly turning into a Pion. The "form factor" is a measure of how the shape changes during that transformation.

• The Upgrade: In the past, they had to throw away some of their data because the computer simulations were too "noisy" (like trying to hear a whisper in a storm).
• The Fix: They developed a new mathematical trick (called "shrinkage") to clean up the noise. Now, they can keep all the data, even the messy parts. This is like using noise-canceling headphones to hear the whisper clearly.
• Result: They have a much clearer picture of the shape change, with smaller margins of error.

Phase B: Re-measuring the "Weight" (The Decay Constants)

They are also measuring the intrinsic "weight" or strength of the Kaon and Pion.

• The Old Way: Previously, they only looked at universes that perfectly matched our real world.
• The New Way: They are now using Staggered Chiral Perturbation Theory (SChPT). Think of this as a predictive map. Instead of just measuring the weight at one point, they measure it at many different "what-if" points (universes with slightly different quark masses) and use the map to predict the weight at the real-world point.
• Why this helps: By using the map, they can use more data points to draw a smoother line, making the final prediction much more accurate.

4. The Secret Sauce: Connecting the Dots

The most exciting part of this paper is the correlation.

Imagine you are trying to guess the weight of a watermelon and the weight of a cantaloupe. If you measure them separately, you might make a small mistake on each. But if you realize that both fruits were grown in the same garden with the same soil and rain, you know that if your soil measurement was slightly off, it affected both fruits in the same way.

The scientists are now calculating the Kaon-to-Pion shape change and the Kaon/Pion weights together. They are mapping out how the errors in one measurement affect the other. By understanding this "dance" between the two measurements, they can cancel out some of the uncertainty.

5. The Goal: Solving the Anomaly

By combining these two super-precise measurements and understanding how they relate to each other, the team hopes to:

1. Pinpoint the error: Determine if the "missing 0.01" in the puzzle is just a measurement mistake that can be fixed.
2. Find New Physics: If the measurements are perfect and the puzzle still doesn't add up, it proves that there is a new, unknown force or particle in the universe that the Standard Model doesn't know about.

Summary

In short, this paper is about a team of scientists upgrading their digital microscope, cleaning up their data, and using a clever mathematical map to measure two tiny particles with unprecedented precision. They aren't just measuring the particles; they are measuring how the errors in those measurements are linked. This "correlated" approach is their best shot at solving one of the biggest mysteries in modern physics: why the universe's fundamental rules seem to have a tiny, unexplained glitch.

Technical Summary

1. Problem Statement

The Standard Model (SM) predicts that the Cabibbo-Kobayashi-Maskawa (CKM) matrix is unitary. However, current experimental data reveals a deficit in the unitarity relation for the first row, known as the Cabibbo Angle Anomaly:
$$|V_{ud}|^2 + |V_{us}|^2 + |V_{ub}|^2 - 1 \neq 0$$
To test this relation with high precision, the values of $|V_{ud}|$ and $|V_{us}|$ must be determined accurately.

• $|V_{ud}|$ is typically derived from superallowed beta decays but suffers from nuclear systematic uncertainties.
• $|V_{us}|$ is derived from kaon semileptonic decays ($K \to \pi \ell \nu$, or $K_{\ell3}$) and the ratio of kaon-to-pion leptonic decay widths ($K_{\ell2}/\pi_{\ell2}$).

These extractions rely heavily on non-perturbative lattice QCD inputs:

1. The vector form factor at zero momentum transfer, $f_+^{K\pi}(0)$.
2. The ratio of decay constants, $f_K/f_\pi$.

The Challenge: Previous lattice calculations treated these inputs independently. To fully reduce the uncertainty in the CKM unitarity test, it is crucial to estimate the correlations between $f_+^{K\pi}(0)$ and $f_K/f_\pi$. Furthermore, previous analyses often relied on data only at physical quark masses or required "data thinning" (discarding data points) to handle statistical noise in covariance matrices.

2. Methodology

The authors present a correlated analysis using Highly Improved Staggered Quarks (HISQ) on $N_f = 2+1+1$ MILC gauge configurations (including up, down, strange, and charm sea quarks).

A. Lattice Setup and Data Generation

• Ensembles: The analysis utilizes a wide range of lattice spacings ($a \approx 0.03$ fm to $0.15$ fm) and physical/lighter-than-physical quark masses.
• New Data: The study includes two new physical quark mass ensembles ($a \approx 0.12$ fm and $0.09$ fm) and increased statistics for the $a \approx 0.06$ fm ensemble (increasing configurations from 692 to 844).
• Kinematics: Partially twisted boundary conditions are used to access the specific kinematic point $q^2=0$ directly in three-point correlation functions.

B. Correlator Analysis and Covariance Handling

• Shrinkage Method: A key methodological improvement is the use of the shrinkage method to construct the covariance matrix. This addresses the issue of small eigenvalues in the covariance matrix that previously forced the authors to shorten fit ranges or discard data ("thinning").
• Result: This allows the inclusion of data from additional source-sink separations ($T$) without instability, effectively increasing statistics and reducing uncertainties without data loss.
• Fit Stability: The analysis employs a conservative fit function with $N_{exp} = 3+3$ exponentials (including opposite-parity states) and varies the minimum time slice ($t_{min}$) to ensure stability.

C. Chiral-Continuum Analysis Strategy

The authors employ Staggered Chiral Perturbation Theory (SChPT) to guide the fits. The strategy is divided into two main components:

1. Semileptonic Form Factor ($f_+^{K\pi}(0)$):

   • The fit function combines Next-to-Leading Order (NLO) Partially Quenched SChPT (PQSChPT) with Next-to-Next-to-Leading Order (NNLO) continuum ChPT.
   • It includes terms for discretization effects ($a^2$) and higher-order chiral corrections.
   • Systematics: The authors are implementing Bayesian Model Averaging (BMA) to account for systematic uncertainties arising from fit function variations and data subset selection.

2. Decay Constants ($f_K/f_\pi$) and Correlations:

   • Two-Step Approach:
     • Step 1: Fit light-meson data (using ensembles with lighter-than-physical strange quark masses) to determine Low Energy Constants (LECs) of the SChPT Lagrangian (e.g., $B_0$, $L_{1-8}$, hairpin parameters).
     • Step 2: Use the LECs from Step 1 as priors to fit the full dataset (including physical masses) to extract $f_K/f_\pi$.
   • Correlation Extraction: By treating the LECs as fit parameters in the semileptonic analysis (rather than fixed values), the authors can directly compute the covariance matrix between $f_+^{K\pi}(0)$ and $f_K/f_\pi$.

3. Key Contributions

• Correlated Analysis: This is the first step toward a fully correlated determination of the two key lattice inputs ($f_+^{K\pi}(0)$ and $f_K/f_\pi$) using the same SChPT framework and LECs.
• Improved Statistical Precision: The use of the shrinkage method eliminates the need for data thinning, allowing the inclusion of more source-sink separations and resulting in reduced statistical uncertainties.
• New Ensembles: Inclusion of new physical mass ensembles at $a \approx 0.12$ fm and $0.09$ fm improves the control over the continuum extrapolation.
• Systematic Error Framework: The introduction of Bayesian Model Averaging (BMA) provides a more robust estimation of systematic errors compared to previous single-model fits.
• Scale Setting: A transition from the $r_1$ scale to the gradient-flow scale $w_0$ (calibrated via the Omega baryon mass) for lattice scale setting.

4. Preliminary Results

• Form Factor ($f_+^{K\pi}(0)$): The reanalysis shows good agreement with previous results but with a slight reduction in statistical uncertainty.
  • Note: The $a \approx 0.12$ fm physical ensemble showed unstable results and was temporarily dropped from the preliminary chiral-continuum fit, though it will be included in systematic error estimates.
• Decay Constants ($f_K/f_\pi$): The two-step SChPT analysis successfully extracts LECs and provides a precise value for the ratio.
• Correlations: Preliminary correlation matrices (Fig. 4 and Fig. 5) demonstrate strong correlations between the form factor and specific LECs (including hairpin parameters). This confirms that the resampling technique effectively propagates correlations through the analysis.
• Consistency: The central values obtained when treating LECs as fixed parameters versus fit parameters are consistent, validating the methodology.

5. Significance

This work represents a critical advancement in the precision testing of the Standard Model:

1. Reducing CKM Uncertainty: By estimating the correlations between $f_+^{K\pi}(0)$ and $f_K/f_\pi$, the total uncertainty in the first-row CKM unitarity test can be significantly reduced.
2. Robustness: The methodological improvements (shrinkage, BMA, SChPT-guided fits) set a new standard for handling lattice QCD data, particularly regarding the treatment of covariance matrices and systematic errors.
3. Phenomenological Impact: The extraction of SChPT Low Energy Constants (LECs) provides essential inputs for other phenomenological studies beyond just the CKM matrix.

In summary, the paper outlines a rigorous, correlated lattice QCD framework that leverages advanced statistical techniques and improved datasets to resolve the Cabibbo Angle Anomaly with higher precision than previously possible.