Determination of $B$-meson distribution amplitudes from $B\to π,K,D$ transition form factors

This paper presents a global fit of $B\to \pi, K, D$ transition form factors using lattice QCD data, light-cone sum rules, and experimental measurements to constrain the inverse moment $\lambda_B$ of the $B$-meson light-cone distribution amplitude and determine the CKM matrix element $|V_{\text{ub}}|$.

The Gist

The Big Picture: Measuring the "Shape" of a Heavy Particle

Imagine the B-meson as a heavy, complex delivery truck driving through a busy city. Inside this truck, there is a heavy driver (the bottom quark) and a light passenger (a "spectator" quark) who is bouncing around in the back.

Physicists want to know exactly how that light passenger is moving. Is it sitting still in the corner? Is it bouncing wildly all over the place? This "movement pattern" is called the Light-Cone Distribution Amplitude (LCDA). It's like a map showing the probability of finding the passenger at any specific spot inside the truck.

The most important number on this map is called $\lambda_B$ (lambda-B). Think of $\lambda_B$ as the "average bounce factor."

• A low $\lambda_B$ means the passenger is mostly huddled near the driver (low momentum).
• A high $\lambda_B$ means the passenger is bouncing around wildly (high momentum).

Knowing this number is crucial because it helps physicists calculate how fast the truck can turn into other vehicles (decays) and helps them measure the fundamental rules of the universe (specifically a number called $|V_{ub}|$).

The Problem: We Didn't Have a Good Map

For a long time, scientists had two ways to guess this "bounce factor," but both were flawed:

1. Theoretical Guesses (QCD Sum Rules): Like trying to guess the passenger's speed by listening to the engine noise. It's useful, but the engine is noisy, and the guesses vary wildly (some say 300, some say 400).
2. Computer Simulations (Lattice QCD): Like trying to film the passenger with a super-fast camera. This is very accurate, but the camera can only film when the truck is moving slowly (low recoil). It can't film the truck when it's speeding up or turning sharply (high recoil).

Because of this gap, scientists couldn't get a precise, single number for the "bounce factor."

The Solution: A Global "Fit"

The authors of this paper decided to play a game of Jigsaw Puzzle. They didn't just look at one piece; they gathered every available piece from different sources to force the picture to make sense.

They combined three types of data:

1. The "Slow Motion" Photos: High-precision data from computer simulations (Lattice QCD) showing how the B-meson turns into pions, kaons, and D-mesons when moving slowly.
2. The "Speeding" Photos: Experimental data from real-world particle colliders (BaBar, Belle, Belle II) showing how often these decays happen when the B-meson is moving fast.
3. The "Theoretical Bridge": A mathematical formula (Light-Cone Sum Rules) that connects the slow photos to the fast photos, using the "bounce factor" ($\lambda_B$) as the key variable.

The Method: Tuning the Radio

Imagine you are trying to tune a radio to a specific station, but the signal is fuzzy.

• The Radio Station is the true value of the "bounce factor" ($\lambda_B$).
• The Static is the uncertainty in our models.
• The Knob is the parameter $\lambda_B$.

The authors took all their data points (the slow photos and the fast photos) and turned the knob ($\lambda_B$) until the theoretical curve perfectly matched all the data points at once. This is called a Global Fit.

They also had to account for the "shape" of the passenger's movement, which they modeled with a three-parameter recipe. They tested thousands of different recipes to see which one made the radio signal the clearest.

The Results: A Clearer Signal

After running this massive global fit, they found:

1. The Bounce Factor ($\lambda_B$): They determined the value to be approximately 217 MeV.

   • Note: This is lower than many previous guesses (which were often around 300–400). Why? Because their new math included a subtle correction (called "Next-to-Leading Power") that previous studies missed. It's like realizing the passenger was actually sitting slightly closer to the driver than we thought.
   • They also found a range: 208 to 324 MeV, acknowledging that our model of the passenger's shape isn't perfect yet.

2. The Universal Constant ($|V_{ub}|$): By pinning down the bounce factor, they could also measure a fundamental constant of nature called $|V_{ub}|$ with high precision: 3.68. This number tells us how likely a bottom quark is to turn into an up quark. Their result matches other major studies, giving physicists more confidence in the Standard Model.

The Takeaway

This paper didn't just guess the value; it forced the value to be consistent with everything we know about how B-mesons behave.

• Before: Scientists had a blurry picture where the "bounce factor" could be almost anything between 300 and 400.
• Now: By combining computer simulations, real-world experiments, and better math, they narrowed it down to a much tighter range around 217.

While there is still some uncertainty (because we don't perfectly know the "shape" of the passenger's movement yet), this is the most precise and comprehensive determination of this number to date. It's like finally getting a high-definition map of the inside of the B-meson truck, which helps us understand the fundamental rules of the universe a little better.

Technical Summary

Technical Summary: Determination of B-meson Distribution Amplitudes from $B \to \pi, K, D$ Transition Form Factors

Problem Statement
The precision of theoretical predictions for exclusive hadronic processes, particularly semi-leptonic $B$-meson decays, is currently limited by the incomplete knowledge of the $B$-meson light-cone distribution amplitudes (LCDAs). These LCDAs encode the momentum distribution of the light spectator parton and are indispensable inputs for Light-Cone Sum Rules (LCSR) calculations at large recoil. A critical parameter in modeling the leading-twist $B$-meson LCDA is its inverse moment, $\lambda_B(\mu)$. While various approaches—including QCD sum rules, lattice QCD, and indirect extractions from $B \to \gamma \ell \nu$ or other form factors—have provided estimates for $\lambda_B$, results vary significantly (ranging roughly from 200 to 500 MeV) and often suffer from large systematic uncertainties or model dependence. Furthermore, a global analysis utilizing all available $B$-to-pseudoscalar transition form factors from lattice QCD combined with experimental semi-leptonic decay data has been missing.

Methodology
The authors propose a global fit strategy to simultaneously constrain the inverse moment $\lambda_B$, the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{ub}|$, and the parameters of the $B$-meson LCDA model. The methodology integrates three distinct data sources:

1. Lattice QCD Data: High-precision form factors ($f_+^P, f_0^P, f_T^P$) in the small recoil (large $q^2$) region for $B \to \pi$, $B \to K$, and $B \to D$ transitions, compiled from the RBC/UKQCD, FNAL/MILC, and HPQCD collaborations.
2. Experimental Data: $q^2$-binned partial branching fractions for $B \to \pi \ell \nu$ from BaBar, Belle, and Belle II collaborations.
3. Theoretical Framework (LCSR): In the large recoil (small $q^2$) region, form factors are calculated using LCSR with $B$-meson LCDAs as input. These calculations include Next-to-Leading Logarithmic (NLL) accuracy at leading power and Next-to-Leading Power (NLP) corrections.

The analysis employs a three-parameter ansatz for the $B$-meson LCDA, $\phi_+^B(\omega)$, expressed via a hypergeometric function. This model is characterized by the inverse moment $\lambda_B$ and two inverse logarithmic moments, $\hat{\sigma}_1$ and $\hat{\sigma}_2$. The $q^2$ dependence of the form factors is parameterized using the Bourrely-Caprini-Lellouch (BCL) expansion with truncation order $N=3$.

A global $\chi^2$ function is constructed by summing contributions from:

• $\chi^2_{LQCD}$: Comparing BCL-parameterized form factors against lattice data points.
• $\chi^2_{LCSR}$: Comparing BCL parameters at $q^2=0$ against LCSR predictions, where the LCSR values are treated as functions of the fit parameters ($\lambda_B, \hat{\sigma}_1, \hat{\sigma}_2$) rather than fixed inputs.
• $\chi^2_{exp}$: Comparing theoretical branching fractions against experimental $B \to \pi \ell \nu$ data.

The fit minimizes the total $\chi^2$ with respect to 20 free parameters: $\lambda_B$, $|V_{ub}|$, 18 BCL coefficients, and the logarithmic moments (with $\hat{\sigma}_2$ treated as a nuisance parameter in the primary fit).

Key Contributions

• Global Fit Strategy: The paper presents the first global analysis that simultaneously utilizes lattice QCD results for $B \to \pi, K, D$ and experimental $B \to \pi \ell \nu$ data to constrain $\lambda_B$. This approach treats the LCSR predictions as parameters to be determined, rather than fixed inputs, allowing for a self-consistent extraction of $\lambda_B$.
• Inclusion of NLP Corrections: The LCSR inputs incorporate NLP corrections, which the authors note reduce the magnitude of the form factors by approximately 30% compared to leading-power-only calculations, significantly impacting the extracted value of $\lambda_B$.
• Model-Dependent Uncertainty Quantification: The study explicitly quantifies the uncertainty arising from the model dependence of the LCDA by varying the inverse logarithmic moments $\hat{\sigma}_1$ and $\hat{\sigma}_2$ within physically motivated ranges.
• Joint Confidence Regions: The authors perform a profile $\chi^2$ analysis to map the joint confidence regions for $\lambda_B$ and $\hat{\sigma}_1$, revealing a strong correlation between these parameters.

Results
The global fit yields the following central values and uncertainties:

• Inverse Moment: $\lambda_B = 217(19)^{+82}_{-17}$ MeV.
  • The first uncertainty is statistical (from input parameters).
  • The second uncertainty is systematic, derived by varying $\hat{\sigma}_1 \in [-0.7, 0.7]$ and $\hat{\sigma}_2 \in [-6, 6]$ while constraining $\lambda_B > 200$ MeV.
• CKM Element: $|V_{ub}| = 3.68(13)^{+0}_{-1} \times 10^{-3}$.
  • The result is consistent with values from exclusive decays in the Particle Data Group and previous global fits.
• Joint Constraints: When $\lambda_B$ and $\hat{\sigma}_1$ are fitted simultaneously, the 68.3% confidence intervals are:
  • $\lambda_B = [208, 324]$ MeV
  • $\hat{\sigma}_1 = [-0.7, 0.27]$

The analysis indicates that while the determination of $\lambda_B$ remains subject to sizable model-dependent uncertainties, the extracted $|V_{ub}|$ is robust and largely insensitive to the specific modeling of the $B$-meson LCDA. The inclusion of $B \to K$ and $B \to D$ lattice data significantly improves the precision of the $B \to \pi$ form factor at $q^2=0$, which is crucial for the extraction of $|V_{ub}|$.

Significance
The paper claims that this work provides a more reliable and comprehensive determination of the $B$-meson LCDA inverse moment $\lambda_B$ by leveraging the complementarity between lattice QCD (small recoil) and LCSR (large recoil). By demonstrating that $|V_{ub}|$ can be precisely constrained even with the current uncertainties in $\lambda_B$, the study reinforces the reliability of exclusive decay determinations of CKM matrix elements. The authors note that the lower central value of $\lambda_B$ compared to many previous theoretical estimates is a direct consequence of including NLP corrections in the LCSR framework. They conclude that while the current theoretical precision is limited by the model dependence of the LCDAs, the methodology established here offers a robust framework for future improvements as lattice data and theoretical calculations (e.g., higher-twist contributions) become more precise.