$B \to \pi$, $B_{(s)} \to D_{(s)}$ from 2+1+1 Flavor Lattice QCD

This paper presents a 2+1+1-flavor lattice QCD calculation of hadronic form factors for $B \to \pi$ and $B_{(s)} \to D_{(s)}$ semileptonic decays using highly improved staggered quarks on MILC ensembles, aiming for a percent-level determination of these form factors to enable high-precision measurements of $|V_{cb}|$.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks called quarks. These bricks snap together to form larger structures, like protons and neutrons, which make up the atoms in everything around us. But sometimes, these heavy Lego structures (called mesons) are unstable and want to break apart.

This paper is about a team of scientists (the Fermilab Lattice and MILC Collaborations) who are trying to predict exactly how these heavy Lego structures break apart, specifically when a "bottom" quark turns into a lighter "charm" or "up" quark.

Here is the breakdown of their work using simple analogies:

1. The Big Mystery: The "CKM" Puzzle

In the Standard Model of physics (our best rulebook for how the universe works), there are numbers called CKM matrix elements (specifically $|V_{cb}|$ and $|V_{ub}|$). Think of these as the "odds" or "probability" that a specific type of Lego break will happen.

• The Problem: Scientists have been measuring these odds in two different ways:
  1. Inclusive: Watching a huge pile of Legos break and counting the total debris.
  2. Exclusive: Watching one specific Lego break into one specific shape.
• The Conflict: The numbers from these two methods don't quite match. It's like if you counted the total number of red bricks in a pile, and it didn't match the number you got by counting every single red brick you saw in a specific box. This suggests we might be missing a piece of the puzzle, or perhaps there is "New Physics" (a hidden rule) we haven't discovered yet.

2. The Tool: The "Digital Time Machine"

To solve this, the scientists can't just watch real particles break apart easily; it's too messy and fast. Instead, they build a super-computer simulation of the universe.

• The Grid (Lattice): Imagine the empty space of the universe is a giant 3D grid, like a fish tank filled with water. The scientists place their Lego bricks on this grid.
• The Simulation: They run a movie of these bricks interacting. Because the math is incredibly hard, they need to use the world's fastest supercomputers (like Franklin, Perlmutter, and Frontier) to crunch the numbers.
• The "2+1+1" Flavor: This is just a fancy way of saying they included the right mix of Lego types in their simulation: two types of light bricks, one type of medium brick, and one type of heavy brick, to make the simulation look exactly like our real universe.

3. The Goal: Measuring the "Stretch" (Form Factors)

When a heavy meson breaks apart, it doesn't just snap; it stretches and changes shape as it transforms. The scientists need to calculate exactly how much it stretches. In physics, this is called a form factor.

• The Analogy: Imagine a rubber band being stretched. The "form factor" is the measurement of how much the rubber band stretches at different speeds.
• Why it matters: If you know exactly how the rubber band stretches (the theory), and you measure how fast it actually stretched in an experiment, you can calculate the "odds" (the CKM numbers) with extreme precision.

4. The New Tricks: Going "Physical" and "Blind"

This paper highlights two major improvements the team made:

• The "Physical Point" Ensemble: In the past, their simulations used Lego bricks that were slightly the wrong weight (too heavy or too light), and they had to guess how to fix it. In this new work, they built a simulation where the bricks are exactly the right weight (the "physical mass"). It's like finally baking a cake with the exact amount of flour and sugar the recipe calls for, instead of guessing and adjusting later. This makes the result much more accurate.
• The "Blind" Analysis: To make sure they didn't accidentally tweak their math to get the answer they wanted, they put a "blindfold" on their data. They multiplied the results by a secret random number. They did all their hard work and analysis without knowing the final number. Only at the very end did they remove the blindfold. This ensures the result is honest and unbiased.

5. The Result: A "Percent-Level" Precision

The team has successfully calculated these "stretching" numbers (form factors) with a precision of about 1%.

• Why 1% is a big deal: The experimental labs (like Belle II in Japan and LHCb in Europe) are now building machines that can measure these particle breaks with 1% precision.
• The Match: Previously, the scientists' predictions were too fuzzy (like a blurry photo) to compare with the sharp new photos from the experiments. Now, the scientists have sharpened their photo to match the experiments.

The Bottom Line

This paper is a progress report on building a perfectly calibrated ruler for the subatomic world. By simulating the universe on a supercomputer with perfect accuracy, the scientists are providing the "theoretical ruler" needed to measure the "odds" of particle decay.

If the new, super-precise measurements from the experiments still don't match this new, super-precise ruler, then we will know for sure that there is New Physics hiding in the shadows—something beyond our current understanding of the universe. If they do match, it confirms our current rulebook is correct. Either way, this work is a crucial step forward in understanding the fundamental laws of nature.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the long-standing discrepancy between inclusive and exclusive determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{ub}|$ and $|V_{cb}|$.

• Context: Exclusive semileptonic decays ($B \to \pi \ell \nu$ and $B_{(s)} \to D_{(s)} \ell \nu$) rely on non-perturbative hadronic form factors to extract these CKM elements from experimental data.
• Current Status: Current uncertainties are approximately 1.3% for $|V_{ub}|$ and 3.9% for $|V_{cb}|$.
• Goal: With experimental facilities like Belle II and LHCb aiming for percent-level precision in measuring decay rates, theoretical calculations must reduce uncertainties to the ~1% level to resolve tensions and test the Standard Model (including Lepton Flavor Universality ratios $R(D^{(*)})$).

2. Methodology

The authors employ Lattice Quantum Chromodynamics (Lattice QCD) to compute the required form factors with high precision.

A. Simulation Setup

• Ensembles: Calculations are performed on 2+1+1 flavor ensembles generated by the MILC Collaboration. These include up, down, strange, and charm sea quarks.
• Action: The Highly Improved Staggered Quark (HISQ) action is used for both valence and sea quarks. This action provides superior control over discretization errors, particularly for heavy quarks.
• Lattice Parameters:
  • Lattice Spacings ($a$): A range from $0.09$ fm down to a new finest spacing of $0.03$ fm.
  • Pion Masses: Includes ensembles with physical pion masses (specifically at $a \approx 0.04$ fm and $a \approx 0.03$ fm), eliminating the need for chiral extrapolation in those regions.
  • Heavy Quarks: Heavy quark masses ($m_h$) are varied from the charm mass up to the physical bottom mass ($m_b$). On the finest ensembles, the calculation reaches the physical $b$-quark mass directly without extrapolation.
• Blinding: To prevent bias, the analysis is blinded; three-point correlation functions are multiplied by a random factor in $[0.95, 1.05]$ until the analysis is finalized.

B. Correlation Functions and Form Factor Extraction

• Two-Point Functions: Used to extract energies and amplitudes of initial ($B$) and final ($\pi, D, D_s$) state hadrons. Bayesian fits with multiple states (including oscillating states inherent to staggered fermions) are used to control excited-state contamination.
• Three-Point Functions: Computed using scalar ($S$), temporal vector ($V_0$), and spatial vector ($V_i$) currents.
• Ratios: Specific ratios of three-point to two-point functions ($R_\parallel, R_\perp, R_0$) are constructed to plateau to the bare form factors ($f_\parallel, f_\perp, f_0$) as the source-sink separation increases.
• Renormalization: Non-perturbative renormalization is performed using the Partial Conservation of the Vector Current (PCVC) relation. A key advantage of the HISQ action is that the product of renormalization factors $Z_m Z_S = 1$, simplifying the extraction of $Z_{V_0}$ and $Z_{V_i}$.

C. Continuum and Chiral Extrapolation

• Fitting Model: For $B_{(s)} \to D_{(s)}$, the authors use a fit model based on Hard SU(2) Chiral Perturbation Theory ($\chi$PT).
  • Includes chiral logarithms up to Next-to-Leading Order (NLO).
  • Includes analytic terms up to Next-to-Next-to-Leading Order (NNLO) to account for mass mistunings, heavy-quark mass dependence, and discretization effects ($a^2, h^2, c^2$).
• $z$-Expansion: To compare with experimental data and ensure analyticity, the continuum results are reparameterized using the $z$-expansion (BCL parameterization). This maps the complex $q^2$ plane to the unit disk, allowing for a convergent series representation of the form factors.

3. Key Contributions

1. New Ensembles: The inclusion of ensembles with physical pion masses and the finest lattice spacing ($a \approx 0.03$ fm) significantly anchors the continuum extrapolation and reduces systematic uncertainties associated with chiral and discretization limits.
2. Direct $b$-Quark Calculation: On the finest lattices, the form factors are computed directly at the physical $b$-quark mass, removing the need for heavy-quark extrapolation in that specific regime.
3. Unified Analysis Strategy: A consistent analysis framework is applied to both $B \to \pi$ and $B_{(s)} \to D_{(s)}$ decays, utilizing joint correlated fits of two- and three-point functions to extract bare form factors.
4. Preliminary $B_{(s)} \to D_{(s)}$ Results: The paper presents the first continuum-limit results for the scalar and vector form factors of $B_s \to D_s$ using this specific 2+1+1 flavor HISQ setup.

4. Results

• Form Factor Behavior: The scalar form factors ($f_0$) are shown to be smooth, rising functions of $q^2$. The data spans a kinematic range of roughly $25 \text{ GeV}^2$.
• Continuum Fits: Preliminary chiral-continuum fits for $B_s \to D_s$ yield $\chi^2/\text{DOF} \approx 1$, indicating a good description of the data. The results are stable under variations of the fit function, prior choices, and the inclusion/exclusion of coarsest or finest ensembles.
• Uncertainty Budget:
  • The dominant uncertainty remains statistical, arising from the underlying lattice data on each ensemble.
  • For the $B_s \to D_s$ channel, the preliminary total uncertainty is less than 1% across the kinematic range.
  • The $z$-expansion parameterization further reduces overall uncertainties compared to the $\chi$PT parameterization due to its rapid convergence.
• Comparison: The continuum extrapolation results from $\chi$PT fits and the $z$-expansion show essentially identical central values, validating the robustness of the extraction.

5. Significance

• Precision Physics: This work represents a critical step toward achieving percent-level precision in theoretical predictions for exclusive $B$-meson decays. This precision is necessary to match upcoming experimental data from Belle II and LHCb.
• Resolving CKM Tensions: By reducing theoretical uncertainties, this calculation aims to clarify whether the discrepancy between inclusive and exclusive $|V_{ub}|$ and $|V_{cb}|$ measurements is due to unaccounted systematic errors or potential New Physics.
• Standard Model Tests: High-precision form factors are essential for testing Lepton Flavor Universality (via $R(D^{(*)})$) and searching for deviations from the Standard Model in rare decays.
• Future Outlook: The methodology established here is being extended to $B \to \pi$ and $B_s \to K$ decays, with the ultimate goal of a correlated analysis to extract the ratio $|V_{ub}|/|V_{cb}|$ directly from lattice data, minimizing experimental systematic errors.