Determination of heavy meson light-cone distribution amplitudes: theoretical framework and lattice simulations

This paper presents a first-principles, continuum-limit lattice QCD determination of heavy meson light-cone distribution amplitudes (LCDAs) using the HQLaMET framework across multiple ensembles, yielding precise QCD and HQET results with total uncertainties below 30% that are consistent with experimental constraints.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks called quarks. When these bricks stick together, they form larger structures called "mesons," which are like small, unstable Lego towers. Some of these towers are "heavy" because they contain a massive brick (a heavy quark), while others are light.

For decades, physicists have been trying to understand exactly how the tiny bricks inside these heavy towers are arranged and how they move. This arrangement is described by something called a Light-Cone Distribution Amplitude (LCDA). Think of the LCDA as a "blueprint" or a "map" that tells you the probability of finding a specific piece of the tower at a specific speed or position while the whole thing is zooming past you.

Knowing this blueprint is crucial. It helps scientists predict how these heavy towers will break apart (decay) and interact with other particles. However, for a long time, this blueprint was missing. Physicists had to guess what it looked like using models, and different guesses led to very different predictions, creating a lot of uncertainty in their calculations.

The Problem: A Broken Compass

The main reason this blueprint was so hard to find is that the heavy towers behave in a tricky way. When you try to look at them using the standard tools of physics (called Lattice QCD), the math gets stuck. It's like trying to take a photo of a speeding car with a camera that only works for stationary objects. The standard method involves looking at a "cusp" (a sharp corner) in the math, which causes the calculation to blow up and become meaningless. This is known as the "cusp divergence."

The Solution: A New Way to Look

The authors of this paper, a large collaboration of scientists, developed a new strategy to fix this. They used a clever two-step approach they call HQLaMET (Heavy-Quark Large-Momentum Effective Theory).

Here is the analogy for their method:

1. The "Quasi" Photo: Instead of trying to take a picture of the tower while it's moving at the speed of light (which is impossible in their computer simulations), they take a picture of the tower while it's moving very fast, but not quite at light speed. This gives them a "blurred" but usable picture called a "quasi-distribution."
2. The "Sharpening" Filter: Once they have this fast-moving picture, they use a mathematical "filter" (called matching) to sharpen it. This filter removes the blur caused by the speed and translates the "quasi" picture into the real, light-speed blueprint they were looking for.

What They Did

To make this work, the team didn't just run one simulation. They ran six different simulations on supercomputers.

• They used different sizes of "pixels" (lattice spacings) to make sure their picture wasn't just a result of low resolution.
• They used different weights for the "light" bricks (pion masses) to ensure the results worked even when the bricks were at their natural, physical weight.
• They used special tricks to make the signal clearer, like "smearing" the connections between the bricks to reduce static noise.

They focused on a specific heavy tower called the D meson (made of a charm quark and a light quark). By analyzing this, they could map out the entire blueprint of how the light quark moves inside the heavy tower.

The Results

The team successfully produced the first "first-principles" (meaning calculated from the basic laws of physics without guessing) maps for these heavy mesons.

• The Shape: They found that the light quark inside the D meson isn't evenly spread out. Instead, it tends to cluster in a specific region, peaking at about 20-30% of the total speed, and then trails off.
• The Precision: Their map has an uncertainty of less than 30% in the most important areas. This is a huge improvement over previous guesses.
• The Check: To make sure they didn't make a mistake, they used a completely different method (calculating specific "moments" or averages) to double-check their work. The two methods agreed perfectly, confirming their results are solid.

Why It Matters (According to the Paper)

The paper states that these new blueprints are essential for the next generation of physics experiments. Specifically, they help scientists calculate the "inverse moment" (a specific number that summarizes the shape of the map) with high precision.

This number is a key ingredient in predicting how B mesons (another type of heavy tower) decay. Since B meson decays are used to test the Standard Model of physics and look for "new physics" (things we haven't discovered yet), having a precise blueprint for the D meson helps remove the "guesswork" from these tests.

In short, the paper claims to have solved a decades-old puzzle by building a new, more reliable camera and a better way to develop the photos, giving physicists their first clear, model-free look at the internal structure of heavy mesons.

Technical Summary

1. Problem Statement

Heavy meson weak decays (e.g., $B \to \pi\pi$, $B \to K^*\ell\ell$) are crucial for testing the Standard Model and searching for new physics. Theoretical predictions for these processes rely on factorization, which separates short-distance perturbative physics from long-distance non-perturbative hadronic dynamics. The latter is encoded in Light-Cone Distribution Amplitudes (LCDAs).

• The Challenge: While LCDAs for light mesons are well-studied, heavy meson LCDAs (both in full QCD and Heavy Quark Effective Theory, HQET) have historically been inaccessible via first-principles lattice QCD.
• Specific Obstacles:
  • HQET Limitation: The defining operator for HQET LCDAs involves a cusp between a time-like heavy quark field and a light-like Wilson line. This induces a cusp divergence, preventing a well-defined local limit and rendering standard Operator Product Expansion (OPE) methods (which rely on local moments) inapplicable.
  • Model Dependence: Previous determinations relied on phenomenological models, leading to large systematic uncertainties (often the dominant error source) in form factor predictions and branching ratios.
  • Lattice Limitations: Conventional lattice QCD works in Euclidean space and cannot directly access light-like (Minkowski) correlations.

2. Methodology: Heavy-Quark Large-Momentum Effective Theory (HQLaMET)

The authors employ a refined framework called HQLaMET, which reorders the scale separation to bypass the cusp divergence issue.

• Scale Hierarchy: Instead of integrating out the heavy mass $m_Q$ first (which leads to HQET and the cusp problem), the framework starts with full QCD fields where the hierarchy is $\Lambda_{\text{QCD}} \ll m_Q \ll P^z$ (large hadron momentum).
  1. Step 1 (LaMET): Integrate out the large momentum $P^z$ to match Euclidean equal-time correlators (quasi-DAs) to QCD LCDAs.
  2. Step 2 (HQET): Integrate out the heavy mass $m_Q$ to extract the HQET LCDA from the QCD LCDA.
• Simulation Setup:
  • Ensembles: Six $N_f=2+1$ gauge ensembles generated by the CLQCD collaboration with lattice spacings $a \in [0.0519, 0.1053]$ fm and pion masses $m_\pi \in [135.5, 317.2]$ MeV.
  • Signal Enhancement: Used momentum-smeared sources, Hypercubic (HYP) smearing for Wilson lines, and optimized interpolating operators to significantly improve the signal-to-noise ratio for nonlocal correlators at large spatial separations and high boosts ($P^z$ up to 3.5 GeV).
• Renormalization & Matching:
  • Hybrid Renormalization: Applied a hybrid scheme to handle the linear divergence of the Wilson line and logarithmic UV divergences, separating short-distance (perturbative) and long-distance (non-perturbative) regions.
  • LaMET Matching: Matched the renormalized quasi-DA to the QCD LCDA using Next-to-Leading Order (NLO) kernels with renormalon resummation to stabilize endpoint behavior.
  • Peak-and-Tail Factorization: For HQET LCDAs, the distribution is split into a peak region (non-perturbative, derived from lattice QCD) and a tail region (perturbative, calculated via HQET). These are merged using a model-independent Laguerre polynomial parametrization.
• Cross-Validation: The QCD LCDA moments derived from LaMET were cross-checked against direct lattice calculations of local twist-two operators using the OPE method in the RI/SMOM scheme.

3. Key Contributions

1. First-Principles Determination: Achieved the first continuum-limit, first-principles determination of heavy meson LCDAs in both full QCD and HQET, removing the reliance on model ansätze.
2. Controlled Extrapolations: Overcame the "single-lattice-spacing" limitation of previous studies by performing controlled continuum, chiral, and infinite-momentum extrapolations to the physical point.
3. Technical Innovations:
   • Demonstrated the viability of HQLaMET for heavy mesons, effectively bypassing the cusp divergence obstruction.
   • Implemented a rigorous cross-check between LaMET (distribution) and OPE (moments), confirming that power corrections are under control.
   • Developed a robust "peak-and-tail" construction for HQET LCDAs that combines non-perturbative lattice data with perturbative QCD.

4. Key Results

• QCD LCDA of the D Meson:
  • The distribution $\phi(y, \mu)$ peaks at a light-quark momentum fraction $y \approx 0.2 - 0.3$.
  • Total uncertainties are below 30% in the physically relevant region $0.1 < y < 0.9$ at $\mu = m_D$.
  • The first two Gegenbauer moments were extracted: $a_1 \approx -0.45$ and $a_2 \approx 0.15$.
• HQET LCDA and Moments:
  • Constructed a continuous HQET LCDA $\phi_+(\omega, \mu)$ over the full $\omega$ range.
  • Inverse Moment ($\lambda_B$): Determined at $\mu = 1$ GeV to be $\lambda_B = 0.340(20)$ GeV. This is consistent with experimental constraints ($\lambda_B > 0.24$ GeV) and recent phenomenological determinations.
  • Inverse-Logarithmic Moment ($\sigma_B^{(1)}$): Determined as $\sigma_B^{(1)} = 1.685(63)$.
• Validation: The Gegenbauer moments derived from the LaMET-reconstructed distribution agree with those calculated directly from OPE local operators within uncertainties, validating the control of finite-$P^z$ power corrections.

5. Significance

• Reduction of Theoretical Uncertainties: The inverse moment $\lambda_B$ is a dominant source of uncertainty in predictions for $B$-meson decay form factors (e.g., $B \to K^*$) and branching ratios. This work reduces the uncertainty on $\lambda_B$ significantly compared to previous model-dependent estimates.
• Phenomenological Impact: The results provide robust, first-principles inputs for next-generation heavy flavor physics, enabling more precise tests of the Standard Model and tighter constraints on New Physics in rare decays like $B \to K^*\ell\ell$.
• Framework Advancement: The successful application of HQLaMET establishes a generalizable pathway for calculating other non-perturbative quantities involving heavy quarks and light-like operators that were previously inaccessible to lattice QCD.

In summary, this paper represents a major breakthrough in lattice QCD, transforming the determination of heavy meson LCDAs from a model-dependent exercise into a rigorous, first-principles calculation with quantified uncertainties, directly addressing a critical bottleneck in heavy flavor phenomenology.