⚛️ phenomenology

Vector mesons leading-twist longitudinal distribution amplitudes and related semi-leptonic decays within QCD sum rules

This paper employs an improved QCD sum rule scheme to calculate the leading-twist longitudinal distribution amplitudes and Gegenbauer moments for light vector mesons ( $\rho, K^\ast, \phi$ ) up to the tenth order, subsequently utilizing these results to determine their distribution behaviors and recalculate the transition form factors and branching ratios for related semi-leptonic $D$ and $D_s$ meson decays.

Original authors: Wan-Bing Luo, Ru-Meng Pan, Ya-Xiong Wang, Tao Zhong

Published 2026-07-21

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Wan-Bing Luo, Ru-Meng Pan, Ya-Xiong Wang, Tao Zhong

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is built from tiny, invisible Lego bricks called quarks. These bricks snap together to form larger structures called hadrons, like protons and neutrons, which make up the matter in our bodies and the stars above. But quarks are shy; they never hang out alone. They are always glued together by a mysterious force called the "strong interaction," which acts like an incredibly stretchy rubber band. To understand how these quark teams move and behave, physicists use a special map called a "distribution amplitude." Think of this map not as a picture of where the quarks are, but as a recipe that tells you how much "momentum" (or zooming power) each quark carries inside the particle. If the recipe is wrong, our predictions about how these particles crash and change into other things will be off, making it hard to test the fundamental rules of the universe.

For a long time, scientists have been trying to write down the perfect recipe for these quark teams, specifically for a group of particles called "vector mesons" (which include the $\rho$ , $K^*$ , and $\phi$ particles). The problem is that the math to figure out these recipes is incredibly messy, like trying to solve a puzzle where some pieces are missing and others are blurry. Different teams of scientists have been using different methods to guess the shape of these recipes, and they often get very different answers. This matters because these recipes are the starting point for calculating how heavy particles, like the "D meson," decay into lighter ones. If we get the recipe wrong, we might miss subtle clues about new physics hiding in the shadows of the standard model.

In this paper, the authors from Guizhou Minzu University decided to tackle this messy puzzle with a fresh approach. They focused on the "leading-twist longitudinal distribution amplitudes" for the $\rho$ , $K^*$ , and $\phi$ mesons. In simpler terms, they wanted to figure out the most important part of the momentum recipe for these specific particles. Instead of relying on the old, clunky ways of calculating these recipes, which often get stuck after just a few steps, they used a newer, sharper tool called "QCD sum rules" within a framework known as "background field theory."

Think of their method like trying to guess the shape of a hidden object by measuring its shadow from many different angles. Previous attempts often only looked at the shadow from one or two angles (low-order moments), which wasn't enough to see the whole picture. The authors, however, calculated the shadows from ten different angles (moments up to the tenth order). They found that the old way of stitching these shadows together was unreliable, so they used a clever trick: they took all ten measurements and used a computer to find the smoothest, most accurate curve that fit all of them at once. This is like having ten different sketches of a face and using a smart algorithm to blend them into one perfect portrait.

The result of this blending act is a set of highly detailed "recipes" for the $\rho$ , $K^*$ , and $\phi$ mesons. They found that the $\rho$ meson's recipe looks a bit like a smooth hill, while the $K^*$ meson's recipe is slightly lopsided, leaning a bit to one side because it contains a heavier "strange" quark. The $\phi$ meson's recipe is very symmetrical. These new recipes are more precise than many previous guesses, sitting comfortably between the predictions of other major theories like Lattice QCD and Dyson-Schwinger equations.

But the authors didn't stop at just drawing the map. They took their new, improved recipes and used them to predict what happens when heavy D-mesons decay into these lighter vector mesons. This is like using the new map to predict exactly how fast a car will go when it drives down a specific hill. They calculated the "branching ratios," which is just a fancy way of saying "how often this specific crash happens compared to all other possible crashes." Their predictions for how often a D-meson turns into a $\rho$ , $K^*$ , or $\phi$ particle match up very well with recent measurements from big experiments like BESIII and CLEO.

In short, the authors didn't just find a new number; they built a better, more reliable map for how quarks share their energy inside these particles. By using a method that gathers more data points than before, they reduced the guesswork. Their findings suggest that the current understanding of these particles is on the right track, and their new, sharper recipes provide a solid foundation for future experiments to hunt for even stranger physics. They didn't discover a new particle, but they gave the scientific community a much clearer lens through which to view the ones we already know.

Technical Summary: Vector Mesons Leading-Twist Longitudinal Distribution Amplitudes and Related Semi-Leptonic Decays within QCD Sum Rules

Problem Statement
The leading-twist longitudinal distribution amplitude (DA), $\phi^{\parallel}_{2;V}(x, \mu)$ , is a critical non-perturbative input for calculating exclusive processes involving light vector mesons ( $V = \rho, K^*, \phi$ ) within Quantum Chromodynamics (QCD). While various theoretical methods exist (e.g., Lattice QCD, Dyson-Schwinger equations, AdS/QCD), they often yield significantly different predictions for the moments and functional behaviors of these DAs. Furthermore, traditional QCD Sum Rule (QCD SR) approaches typically provide reliable results only for the lowest few Gegenbauer moments. Truncating the Gegenbauer series expansion with only a few low-order moments is often insufficient to accurately describe the full behavior of the DA, leading to discrepancies between theoretical predictions and increasingly precise experimental data from collaborations such as BESIII, CLEO, and BaBar regarding semi-leptonic decays $D_{(s)} \to V$ .

Methodology
This work adopts a specific research scheme for QCD SRs proposed by the authors in a previous study [34], utilizing the Background Field Theory (BFT) framework. The methodology proceeds in two main stages:

Calculation of $\xi$ -moments:
Instead of relying on the traditional extraction of Gegenbauer moments, the authors calculate the $\xi$ -moments, $\langle \xi^n \rangle^{\parallel}_{2;V}$ , up to the tenth order ( $n=1, \dots, 10$ ).
- A new sum rule formula is employed to address the issue where the zero-order moment sum rule cannot be normalized across the entire Borel parameter region due to the absence of high-dimensional vacuum condensate terms. The new formula (Eq. 12) normalizes the $n$ -th moment by the zeroth moment derived from the same sum rule.
- The calculation involves the Operator Product Expansion (OPE) in the deep Euclidean region, including contributions up to dimension-six vacuum condensates and quark mass corrections.
- The Borel windows are strictly determined to ensure continuum contributions remain below 45-50% and dimension-six condensate contributions below 5-8%.
Determination of DA Behavior:
Rather than converting high-order $\xi$ -moments into unreliable high-order Gegenbauer moments, the authors fit the calculated $\langle \xi^n \rangle^{\parallel}_{2;V}$ values directly using the least squares method.
- A normalized power-law parametrization (PLP) model (Eq. 14) is used to describe the DA behavior.
- The model parameters ( $\alpha^{\parallel}_{2;V}, \beta^{\parallel}_{2;V}$ ) are optimized to minimize the likelihood function $\chi^2$ , thereby determining the precise functional form of $\phi^{\parallel}_{2;V}(x, \mu)$ .
Application to Semi-Leptonic Decays:
The determined DAs are applied within the Light-Cone QCD Sum Rules (LCSRs) framework to calculate:
- Transition Form Factors (TFFs) $A_1(q^2), A_2(q^2), V(q^2)$ for $D \to \rho, K^*$ and $D_s \to \phi$ .
- The TFFs are extrapolated from the low/intermediate $q^2$ region (where LCSRs are valid) to the full physical region using the Simplified Series Expansion (SSE) method based on $z$ -expansion.
- Branching fractions and decay widths for $D^0 \to \rho^- \ell^+ \nu_\ell$ , $D^+ \to \bar{K}^{*0} \ell^+ \nu_\ell$ , and $D^+_s \to \phi \ell^+ \nu_\ell$ are computed.

Key Results

Moments and Parameters: The authors provide precise values for the first ten $\xi$ $ξ$ -moments and the corresponding Gegenbauer moments ( $a_1, a_2$ $a_{1}, a_{2}$ ) at scales $\mu = 1, 1.4, 2$ $μ = 1, 1.4, 2$ GeV.
- For $\rho$ : $\langle \xi^2 \rangle^{\parallel}_{2;\rho} = 0.225^{+0.013}_{-0.012}$ at $\mu=1$ GeV.
- For $K^*$ : $\langle \xi^1 \rangle^{\parallel}_{2;K^*} = -0.0228^{+0.0042}_{-0.0040}$ and $\langle \xi^2 \rangle^{\parallel}_{2;K^*} = 0.217^{+0.007}_{-0.007}$ at $\mu=1$ GeV.
- For $\phi$ : $\langle \xi^2 \rangle^{\parallel}_{2;\phi} = 0.209^{+0.020}_{-0.020}$ at $\mu=1$ GeV.
DA Behaviors: The fitted PLP model parameters yield DA curves that:
- For $\rho$ : Show a peak region higher than Dyson-Schwinger Equation (DSE) predictions but lower than the asymptotic form.
- For $K^*$ : Exhibit a slight asymmetry with a peak shifted to $x < 0.5$ due to $s$ -quark mass effects (SU(3) flavor symmetry breaking), consistent with Lattice QCD (LQCD) and asymptotic expectations.
- For $\phi$ : Show behaviors very close to LQCD and the asymptotic form, differing from earlier QCD SR predictions.
- All three DAs tend toward the asymptotic form as the scale $\mu$ increases.
Transition Form Factors (TFFs): At the large recoil point ( $q^2=0$ ), the calculated TFFs (e.g., $A_1^{D\to\rho}(0) = 0.528^{+0.039}_{-0.036}$ ) are consistent with experimental data from CLEO and theoretical predictions from LQCD and other models within error margins. The ratios $r_V = V(0)/A_1(0)$ and $r_2 = A_2(0)/A_1(0)$ agree with recent BESIII measurements.
Branching Fractions: The calculated branching fractions for $D^0 \to \rho^- e^+ \nu_e$ ( $1.877^{+0.271}_{-0.241} \times 10^{-3}$ ), $D^+ \to \bar{K}^{*0} e^+ \nu_e$ ( $5.558^{+0.712}_{-0.635} \times 10^{-2}$ ), and $D^+_s \to \phi e^+ \nu_e$ ( $2.501^{+0.329}_{-0.240} \times 10^{-2}$ ) are in agreement with PDG averages and recent BESIII/CLEO measurements within uncertainties.

Significance and Claims
The paper claims that the adopted scheme—calculating high-order $\xi$ -moments and fitting them directly to a phenomenological model—offers a more reliable determination of DA behavior than traditional methods relying on truncated Gegenbauer expansions. By avoiding the unreliability of high-order Gegenbauer moments and utilizing a larger set of constraints (up to the tenth $\xi$ -moment), the authors provide DAs that are consistent with diverse theoretical approaches (LQCD, DSE) and experimental data.

The study demonstrates that these improved DAs can successfully resolve discrepancies in semi-leptonic decay calculations. The resulting branching fractions and TFF ratios align with high-precision experimental data from BESIII and CLEO, suggesting that the improved treatment of non-perturbative inputs is crucial for testing the Standard Model and exploring potential new physics in heavy meson decays. The work emphasizes that the specific treatment of non-perturbative QCD dynamics, particularly through the BFTSR framework with the new sum rule formula, is effective in reducing theoretical uncertainties in these exclusive processes.

More like this