⚛️ phenomenology

Pseudoscalar meson decay constants and distribution amplitudes up to twist-4 in the light-front quark model

This paper utilizes the light-front quark model with the Bakamjian-Thomas construction to demonstrate the frame, component, and process independence of pseudoscalar meson decay constants and to analyze their twist-4 distribution amplitudes and moments for both $1S$ and $2S$ heavy meson states.

Original authors: Ahmad Jafar Arifi, Ho-Meoyng Choi, Chueng-Ryong Ji

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ahmad Jafar Arifi, Ho-Meoyng Choi, Chueng-Ryong Ji

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

Inside the heart of every atom, protons and neutrons are built from smaller, more fundamental particles called quarks. These quarks are never found alone; they are permanently bound together by the strong force, a powerful interaction that acts like a cosmic glue. When a quark and an antiquark pair up, they form a short-lived particle known as a meson. To understand how the universe works at its most basic level, physicists must map out exactly how these quarks move and share energy within their tiny prison. They do this by studying "distribution amplitudes," which are essentially probability maps showing where a quark is likely to be found and how fast it is moving relative to its partner. These maps are crucial for predicting how particles behave in high-energy collisions, such as those created in massive particle accelerators around the world.

A team of researchers has recently refined the way these maps are calculated for heavy mesons, which contain at least one very heavy quark. Using a theoretical framework called the light-front quark model, they examined the internal structure of these particles in their ground state and in a higher-energy excited state. Their goal was to resolve a long-standing puzzle: when physicists calculate the same physical property using different mathematical tools, they often get slightly different answers. This inconsistency suggests that something is missing from the calculation, much like a scale that gives different weights depending on which side you stand on. The researchers set out to prove that their method could produce a single, consistent answer regardless of the mathematical path taken.

The team focused on a specific property called the decay constant, which measures how quickly a meson breaks apart into other particles. In their study, they calculated this value using three distinct mathematical approaches, each relying on a different type of interaction operator. In previous attempts, these three methods often yielded conflicting results, creating confusion about the true nature of the particle. However, by applying a specific consistency rule to their model—essentially ensuring that the mass used in the calculation reflects the internal energy of the quark pair rather than the total mass of the moving particle—they found that all three methods converged on the exact same number. This result confirms that their model is robust and that the decay constant is a true, unchanging property of the meson, independent of how it is measured.

Beyond confirming the decay constant, the researchers also mapped out the distribution amplitudes up to a high level of detail, known as twist-four. This level of detail accounts for complex motions, such as the quarks moving sideways or interacting with extra particles, which are usually ignored in simpler models. They discovered that for heavier mesons, the probability maps shift significantly, with the lighter quark carrying most of the forward momentum. They also found that the excited state of the meson, which is like a vibrating version of the ground state, has a distinct internal structure with a node, or a point where the probability of finding a quark drops to zero. This nodal structure creates a double-humped shape in the probability map, a clear signature of the particle's excited nature.

The study also revealed how the internal spin of the quarks, known as helicity, contributes to the overall behavior of the meson. When the researchers looked at the "good" components of the particle's motion, only the standard spin configurations appeared. However, when they examined the more complex "bad" components, contributions from higher, more unusual spin states emerged. These unusual states are usually suppressed in heavy particles but become more prominent in the excited states and at lower transverse momenta. By carefully accounting for these contributions, the team ensured that their final results remained consistent and did not depend on the specific angle or frame of reference used for the calculation.

The researchers provided detailed numerical predictions for the decay constants of various heavy mesons, including those containing charm and bottom quarks. For the ground state of the D meson, they calculated a value of 208 MeV, which aligns closely with experimental measurements of 206.7 MeV. For the excited D meson, they predicted a value of 110, a figure for which no experimental data currently exists. Similar predictions were made for other mesons like the B and Bs families, offering a comprehensive set of benchmarks for future experiments. These numbers are not just abstract figures; they serve as essential inputs for calculating how these particles interact in the real world, helping to refine our understanding of the strong force.

Ultimately, this work demonstrates that a consistent theoretical framework can successfully describe the complex internal lives of heavy mesons without contradictions. By proving that different mathematical routes lead to the same physical destination, the researchers have strengthened the reliability of the light-front quark model. Their findings provide a clearer picture of how quarks organize themselves inside these particles, offering a solid foundation for interpreting data from current and future particle physics experiments. The ability to accurately predict the behavior of these particles, especially in their excited states, opens new doors for exploring the fundamental forces that shape our universe.

Technical Summary: Pseudoscalar Meson Decay Constants and Distribution Amplitudes up to Twist-4 in the Light-Front Quark Model

Problem Statement
The distribution amplitudes (DAs) of mesons are essential non-perturbative inputs for understanding light-cone dominated processes via collinear factorization. While leading-twist DAs are well-studied, higher-twist DAs (twist-3 and twist-4) are crucial for precision analyses in modern high-energy experiments (e.g., KEKII, LHC, EIC) but remain challenging to compute. In the Light-Front Quark Model (LFQM), a significant theoretical hurdle exists: the lack of manifest covariance often leads to frame-dependent results and ambiguities regarding Light-Front (LF) zero-mode contributions and instantaneous terms. Specifically, calculating higher-twist DAs requires utilizing "bad" current components (e.g., $J^-$ ), which are sensitive to these zero modes. Previous studies have shown that without a rigorous treatment, decay constants derived from different current operators ( $\Gamma$ ) or different current components can yield inconsistent results, violating process independence and rotational invariance.

Methodology
The authors employ a Light-Front Quark Model (LFQM) based on the Bakamjian-Thomas (BT) construction. In this framework, the meson state is constructed from non-interacting constituent quark and antiquark representations, while interactions are incorporated into the mass operator ( $M = M_0 + V_{Q\bar{Q}}$ ).

Wave Functions: The model utilizes radial wave functions (Gaussian type) mixed for $1S$ and $2S$ states, combined with interaction-independent Melosh-transformed spin-orbit wave functions.
Self-Consistency Condition: A critical methodological feature is the enforcement of the "Type II" link between the manifestly covariant Bethe-Salpeter (BS) model and the LFQM. This involves replacing the physical meson mass $M$ with the invariant mass $M_0$ in the integrands of matrix element calculations. This replacement effectively includes LF zero-mode contributions and restores Lorentz and rotational invariance.
Operators: The study calculates matrix elements $\langle 0 | \bar{q} \Gamma q | P \rangle$ using three distinct current operators: axial-vector ( $\Gamma_A = \gamma_\mu \gamma_5$ ), pseudoscalar ( $\Gamma_P = i\gamma_5$ ), and pseudotensor ( $\Gamma_T = \sigma_{\mu\nu} \gamma_5$ ).
Frame Independence: Calculations are performed in a general reference frame where the transverse momentum $P_\perp \neq 0$ , allowing for a rigorous test of boost and rotational invariance.
Twist-4 Analysis: The study explicitly derives the twist-4 distribution amplitude from the minus component of the axial-vector current, a component often difficult to handle due to zero-mode issues.

Key Contributions and Results

Process-Independent Decay Constants: The authors demonstrate that the pseudoscalar meson decay constants ( $f_P$ ) derived from the three distinct current operators ( $\Gamma_A, \Gamma_P, \Gamma_T$ ) are numerically identical within their LFQM framework. This equality holds regardless of the specific current component used (plus, transverse, or minus).
- This consistency is achieved only when the self-consistency condition ( $M \to M_0$ ) is imposed. Without this, results differ between current components.
- For the pseudotensor current, the inclusion of the factor $(1 - \rho_+)$ in the definition of the nonlocal matrix element is shown to be essential for obtaining the same decay constant as the local operators.
Lorentz and Rotational Invariance: The study explicitly proves that the calculated decay constants are independent of the frame's transverse momentum ( $P_\perp$ ) and invariant under rotations.
- While individual helicity contributions to the decay constant (particularly from the "bad" minus current) depend on $P_\perp$ , their sum is strictly $P_\perp$ -independent.
- The analysis of the wave function in 3D momentum space confirms that the decay constant is rotationally invariant, with odd functions of $k_z$ integrating to zero.
Helicity Contributions: The paper provides a quantitative breakdown of helicity contributions ( $\uparrow\downarrow, \downarrow\uparrow$ vs. $\uparrow\uparrow, \downarrow\downarrow$ ).
- "Good" currents (plus component) receive contributions only from ordinary helicities.
- "Bad" currents (minus component) involve higher helicity contributions, which are suppressed as quark mass increases but are non-negligible, particularly for radially excited ( $2S$ ) states.
- The higher helicity contributions are more pronounced in the low $P_\perp$ region.
Distribution Amplitudes (DAs) up to Twist-4:
- The authors derive twist-2, twist-3, and twist-4 DAs.
- For unequal-mass mesons (e.g., $D, B$ ), higher-twist DAs are more concentrated in the lower longitudinal momentum fraction ( $x$ ) region.
- For the $2S$ state, the nodal structure of the wave function manifests as a dip in the DAs, with the higher-twist DAs showing enhanced peaks near $x=0$ and suppressed peaks near $x=0.5$ .
- The study clarifies the definition of the twist-3 DA $\psi_{3;P}(x)$ from the pseudotensor current, showing that including the $(1-\rho_+)$ factor yields a narrower, higher peak compared to definitions omitting this factor.
Moments: The $\xi$ -moments ( $\langle \xi^n \rangle$ , where $\xi = 2x-1$ ) up to $n=6$ are computed for various heavy pseudoscalar mesons ( $D, D_s, \eta_c, B, B_s, B_c, \eta_b$ ) in both $1S$ and $2S$ states. The results indicate that odd moments grow in magnitude for higher twists, reflecting a shift of the DA away from the center.

Significance
The paper claims that its primary significance lies in establishing a self-consistent LFQM framework capable of simultaneously describing mass spectroscopy and wave-function-related observables (decay constants and DAs) without violating fundamental symmetries. By enforcing the BT construction and the $M \to M_0$ replacement, the authors resolve the long-standing issue of frame and component dependence in LFQM calculations. The explicit demonstration of the equality of decay constants from three different operators validates the model's robustness. Furthermore, the provision of twist-4 DAs and their moments for both ground and radially excited heavy mesons offers essential non-perturbative inputs for future QCD factorization calculations, particularly in low- $Q^2$ regions where higher-twist effects become relevant. The authors note that these results are intended to support the calculation of hard exclusive processes and suggest that the methodology could be extended to vector, scalar, and tensor mesons.

More like this