Untangling the heavy-flavor mess: status of the Fermilab-MILC calculation of the $B_{(s)}\to D^{(\ast)}_{(s)}\ell\nu$ form factors

This paper presents the status of Fermilab-MILC calculations for $B_{(s)}\to D^{(\ast)}_{(s)}\ell\nu$ form factors using seven physical-mass HISQ ensembles, aiming to resolve current unexplained tensions in heavy-to-light results and clarify the ambiguous situation in the $B\to D^{\ast}\ell\nu$ channel despite recent advances in heavy-to-heavy decay characterization.

The Gist

Imagine the universe is a giant, complex puzzle. For decades, scientists have been trying to solve it using a rulebook called the Standard Model. But recently, they've found some pieces that just don't seem to fit. These "glitches" are called B-anomalies, and they involve heavy particles (like the "B" meson) decaying into lighter ones.

This paper is a report from a team of scientists (the Fermilab Lattice and MILC collaborations) who are trying to fix the puzzle by recalculating the "instructions" for how these heavy particles break apart.

Here is the breakdown of their work in simple terms:

1. The Problem: Two Different Messes

The scientists are dealing with two different types of "messes" in their data:

• Mess #1: The "Heavy-to-Heavy" Tangle (The B to D problem)*

  • The Analogy: Imagine two heavy trucks (B-mesons) crashing and turning into slightly smaller trucks (D-mesons). Scientists have been arguing about exactly how much "fuel" (energy) is released in this crash.
  • The Issue: There are two ways to measure this. One way (inclusive) looks at the total fuel, and the other (exclusive) looks at the specific parts. For a long time, these two measurements disagreed. Recently, new computer simulations (lattice QCD) were made to help, but they are still slightly confused. They agree with each other reasonably well, but they don't quite match what the real-world experiments (like LHCb and Belle II) are seeing. It's like three different weather forecasters predicting rain, but the actual weather is sunny.
  • The Goal: The team wants to sharpen their computer models to see if the "rain" was just a calculation error or if it means there is a new, unknown force of nature at play.

• Mess #2: The "Heavy-to-Light" Tangle (The B to Pion problem)

  • The Analogy: Now imagine a heavy truck crashing and turning into a tiny bicycle (a pion). This is a much harder crash to simulate.
  • The Issue: This is the bigger mess. Different teams of scientists have been running these simulations, and they are getting completely different answers. It's like if Team A says the bicycle weighs 10 pounds, and Team B says it weighs 50 pounds, and they can't figure out why. This disagreement is so bad that it makes everyone lose confidence in the data.
  • The Goal: The team is trying to figure out why their calculations are so messy and fix the "glitch" in their code so everyone agrees on the weight of the bicycle.

2. The Solution: A New, Super-Powered Simulation

To fix these messes, the authors are running a massive, new computer simulation. Think of it as building a super-detailed virtual laboratory.

• The Tools: They are using seven different "ensembles" (sets of virtual universes) with different levels of detail. Some are blurry (coarse), and some are incredibly sharp (fine), ranging from a resolution of 0.15 down to 0.06.
• The Setup: They are simulating these crashes with the actual physical masses of the particles (not fake, approximated masses). They are using a specific method called the "Fermilab interpretation" to handle the heavy particles, which is like using a special lens to see the heavy trucks clearly without blurring the image.
• The Strategy: They are calculating everything at once. They aren't just looking at one crash; they are looking at all the heavy-to-heavy and heavy-to-light crashes simultaneously. This allows them to see how the errors in one calculation might affect the others (correlations).

3. The Current Status: "Blinded" but Promising

The team has finished the initial calculations, but they are currently "blinded."

• The Analogy: Imagine you are baking a cake to win a contest. You have mixed the ingredients and baked it, but you haven't tasted it yet because you don't want to accidentally change the recipe based on how it tastes. You are waiting to "unblind" the results to see the final numbers.
• The Results: So far, the preliminary (unseen) numbers look promising. They seem to be consistent with previous work, and the "heavy-to-light" calculations (the bicycle crash) are looking a bit more stable than before.

4. Why Does This Matter?

If the scientists can fix these calculations and the numbers still don't match the real-world experiments, it's a huge deal. It would mean the Standard Model is incomplete and that New Physics (something entirely new, like a new force or particle) exists.

However, if the calculations are just messy and wrong, then the "anomalies" disappear, and the Standard Model remains safe.

In summary: This paper is a progress report from a team of digital architects who are rebuilding the foundation of their virtual laboratory. They are trying to untangle a knot of confusing data to see if the universe is behaving exactly as we think it does, or if it's hiding a secret. They are almost ready to "unblind" their results, and the physics community is waiting with bated breath.

Technical Summary

1. Problem Statement

The paper addresses critical tensions in the flavor sector of the Standard Model (SM), specifically regarding "B-anomalies" and the determination of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Two primary issues are highlighted:

• Heavy-to-Heavy Decays ($|V_{cb}|$ and LFU): There is a persistent, mild tension (approx. $2\sigma$–$3\sigma$) between inclusive and exclusive determinations of the CKM matrix element $|V_{cb}|$. Furthermore, Lepton Flavor Universality (LFU) ratios, specifically $R(D)$ and $R(D^*)$, show a $\approx 3.5\sigma$ tension between experimental measurements and SM predictions. While recent lattice QCD (LQCD) calculations of $B \to D^*\ell\nu$ form factors have improved, they have not yet resolved these discrepancies.
• Heavy-to-Light Decays ($|V_{ub}|$ and Form Factor Inconsistencies): While the tension in $|V_{ub}|$ is less severe, a more fundamental problem exists: LQCD calculations of heavy-to-light form factors (specifically $B \to \pi\ell\nu$ and $B_s \to K\ell\nu$) from different collaborations (e.g., FNAL/MILC, RBC/UKQCD, HPQCD) are mutually incompatible. The FLAG 2024 review reports $\chi^2/\text{dof} > 3$ for combined fits, indicating a lack of consensus that undermines confidence in the data used for phenomenology.

2. Methodology

The Fermilab Lattice and MILC Collaborations are undertaking a comprehensive calculation to address these issues using a specific set of lattice ensembles and actions:

• Ensembles: The calculation utilizes seven $N_f = 2+1+1$ Highly Improved Staggered Quark (HISQ) gauge ensembles.
  • Lattice Spacings: Ranging from $0.15$ fm down to $0.06$ fm to facilitate a robust continuum limit extrapolation.
  • Quark Masses: More than half of the ensembles feature physical pion masses.
• Action and Interpretation:
  • Light/Strange Quarks: Simulated using the HISQ action.
  • Heavy Quarks ($b$ and $c$): Simulated using the Wilson-clover action with the Fermilab interpretation. This approach is designed to handle heavy quarks where the mass is comparable to the inverse lattice spacing ($am \sim 1$).
• Scope of Calculation: The project aims for a fully correlated analysis of:
  • Heavy-to-Heavy: $B_{(s)} \to D^{(*)}_{(s)}\ell\nu$ (four channels).
  • Heavy-to-Light: $B \to \pi\ell\nu$, $B_s \to K\ell\nu$, and $B \to K\ell\ell$ (three channels).
• Cross-Validation: Three ensembles are shared with a concurrent "HISQ-on-HISQ" calculation (where heavy quarks are also HISQ), allowing for a direct comparison of heavy-quark actions on the same problem to assess systematic errors.

3. Key Contributions and Status

The paper presents the current status of these calculations, noting that results are preliminary and blinded (masked to prevent bias during analysis).

• Heavy-to-Heavy Sector ($B \to D^{(*)}$):
  • Form factors for all heavy-to-heavy channels have been calculated.
  • Preliminary plots (Figs. 8–11) show the form factors ($h_+, h_-, h_{A1}, h_V$) across the recoil variable $w$.
  • The results appear consistent with existing calculations from JLQCD and HPQCD.
  • Notably, the $B_s$ decays are calculated with higher precision than $B$ decays due to better signal-to-noise ratios in the lattice data.
• Heavy-to-Light Sector ($B \to \pi, B_s \to K$):
  • The analysis is in a more advanced state than the heavy-to-heavy sector.
  • Preliminary results for the chiral-continuum extrapolation of $f_+$ and $f_0$ form factors are presented (Fig. 12).
  • The collaboration is working to finalize the analysis to resolve the discrepancies seen in previous FLAG reports.
• Systematic Error Mitigation:
  • The paper discusses the issue of taking the continuum limit in different bases (e.g., $f_\parallel, f_\perp$ vs. $f_+, f_0$).
  • Previous tests by the collaboration (Fig. 6) indicated that the choice of basis does not significantly affect results within $1\sigma$, suggesting that the large tensions seen in other collaborations may stem from other systematic sources or data handling, rather than the continuum limit procedure itself.

4. Results

• Current State: The results are currently blinded and preliminary. No final unblinded values for $|V_{cb}|$, $|V_{ub}|$, or $R(D^{(*)})$ are provided in this conference proceeding.
• Preliminary Observations:
  • The form factors for $B \to D^*\ell\nu$ show consistency with other LQCD groups, suggesting the "mess" in heavy-to-heavy calculations is one of precision rather than fundamental disagreement.
  • The heavy-to-light analysis is proceeding toward a resolution of the inter-collaboration tensions, with the chiral-continuum fits currently being finalized.

5. Significance

This work is significant for several reasons:

1. Resolving B-Anomalies: By providing high-precision, correlated form factors for both heavy-to-heavy and heavy-to-light decays, this calculation aims to definitively determine whether the observed tensions in $|V_{cb}|$ and LFU ratios are signs of New Physics or artifacts of theoretical uncertainties.
2. Addressing Lattice Inconsistencies: The heavy-to-light sector suffers from a lack of agreement between major LQCD groups. This calculation, using a distinct heavy-quark action (Fermilab-clover) on physical-mass ensembles, provides an independent check that could identify the source of the $\chi^2$ incompatibility reported by FLAG.
3. Methodological Benchmarking: By running two parallel calculations (Fermilab-clover vs. HISQ-on-HISQ) on shared ensembles, the collaboration creates a unique opportunity to benchmark heavy-quark discretization errors, a critical step for future precision flavor physics.
4. Future Impact: With experiments like LHCb and Belle II generating data at high rates, the community requires equally precise theoretical inputs. This work represents a concerted effort to match experimental precision, potentially clarifying the fate of the $B$-anomalies in the coming months.