⚛️ nuclear theory

Light cone distributions of $P$ -wave quarkonia

Using the basis light-front quantization formalism, this study reveals that relativistic effects induce novel "W"-shaped structures in the leading-twist distribution amplitudes of $P$ -wave axial vector quarkonia, highlighting the necessity of relativistic frameworks for accurately modeling heavy quarkonium production and structure.

Original authors: Yang Li, Tianyang Hu, Xianghui Cao, Meijian Li

Published 2026-09-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Yang Li, Tianyang Hu, Xianghui Cao, Meijian Li

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 held together by a force so powerful it defies everyday intuition. This force, carried by particles called gluons, binds together even smaller entities known as quarks. When quarks pair up, they form a family of particles called mesons. Some of these mesons are simple, sitting in their lowest energy state, while others are excited, vibrating with extra energy like a plucked guitar string. Physicists have long been interested in a specific group of these excited particles, known as P-wave quarkonia. These are heavy mesons made of a quark and its antimatter partner, spinning in a way that gives them a distinct shape and behavior. Understanding exactly how the quarks move inside these particles is crucial. It is the missing piece of the puzzle needed to interpret the massive amounts of data pouring in from modern particle colliders, where scientists smash particles together at near-light speeds to recreate the conditions of the early universe.

For decades, scientists have relied on a simplified way of thinking about these heavy particles, treating them almost like tiny, slow-moving planets orbiting a sun. This approach, known as the non-relativistic model, works well enough when the particles are heavy and moving slowly. However, as we look closer, we realize that inside these tiny particles, the quarks are moving incredibly fast, close to the speed of light. At such speeds, the rules of Einstein's relativity take over, and the simple planetary model begins to break down. The question has been: how much does this high-speed motion actually change the internal structure of these particles? Does it just tweak the numbers, or does it fundamentally reshape the way the particles are built?

A team of researchers has now answered this question by peering into the internal architecture of these P-wave quarkonia using a sophisticated mathematical framework called basis light-front quantization. Instead of treating the particles as slow-moving objects, they used a method that accounts for the fact that the quarks are moving at relativistic speeds. They focused on three main types of these particles: scalar mesons, which have no spin; axial vector mesons, which have a specific type of spin alignment; and tensor mesons, which have a more complex spin structure. By running detailed simulations, they mapped out the "light-cone distribution amplitudes" of these particles. In simple terms, this is a map that shows how the momentum of the particle is shared between the quark and the antiquark as they zip around inside.

The researchers found that while the overall shape of these particles looks somewhat familiar, resembling the predictions of the older, slower models, the details tell a very different story. The most striking discovery appeared in the axial vector mesons. In the old, non-relativistic view, these particles were expected to have a smooth, predictable shape. However, the new simulations revealed a surprising and complex pattern: a distinct "W" shape in the way the momentum is distributed. This shape was not predicted by the traditional models. It arises because the fast-moving quarks inside the particle mix different types of motion in a way that is only possible when relativity is taken into account. Specifically, the rapid motion causes the quarks to adopt a combination of orbital movements that would be impossible if they were moving slowly. This mixing creates a structural feature that is purely a product of high-speed dynamics.

This finding is significant because it proves that even in heavy particles, which were once thought to be simple enough to ignore relativistic effects, the speed of the quarks plays a critical role in defining their identity. The researchers observed that as the particles get heavier, moving from the lighter charmonium family to the heavier bottomonium family, the strange "W" shape begins to fade, and the particles start to look more like the simple, non-relativistic models predicted. This confirms that the effect is indeed driven by the speed of the quarks relative to the speed of light. The study also showed that these complex internal structures are not just mathematical curiosities; they are essential for accurately describing how these particles are created and how they behave in high-energy collisions.

The implications of this work extend to the very frontiers of experimental physics. When scientists at facilities like the Large Hadron Collider or the Belle II experiment observe these particles, they are looking for clues about the fundamental forces of nature. If the theoretical models used to interpret these observations ignore the relativistic "W" shape, they risk misreading the data. By providing a more accurate map of how these particles are built, this research offers a new, more precise tool for physicists. It allows them to better understand the hard exclusive processes—specific types of collisions where particles are produced without breaking apart—where these P-wave quarkonia play a starring role. The study does not just refine an old theory; it reveals a new layer of complexity in the subatomic world, showing that even in the heaviest particles, the dance of the quarks is far more intricate than previously imagined.

Technical Summary: Light Cone Distributions of P-Wave Quarkonia

Problem Statement
Light-cone distributions encode the non-perturbative structural information of hadrons and are essential for describing high-energy scattering experiments, particularly inclusive and exclusive processes at large momentum transfer. While Parton Distribution Functions (PDFs) and Light-Cone Distribution Amplitudes (LCDAs) are well-studied for ground states, there is renewed interest in the light-cone distributions of P-wave quarkonia (scalar, axial vector, and tensor states). Theoretical access to these quantities is challenging, especially for Euclidean formulations like Lattice QCD, due to the intrinsically Minkowskian nature of light-cone distributions. Furthermore, in a relativistic framework, spin ( $S$ ) and orbital angular momentum ( $L$ ) are not strictly conserved quantum numbers. This necessitates a rigorous examination of how angular excitation interplays with relativistic dynamics to dictate the shape of P-wave mesons, particularly in the transition region between relativistic and non-relativistic regimes where charmonia and bottomonia reside.

Methodology
The authors investigate the leading-twist LCDAs of P-wave quarkonia using the Basis Light-Front Quantization (BLFQ) formalism. The study employs an effective light-cone Hamiltonian ( $H_{lc}$ ) that combines light-front holographic QCD confinement with a short-distance one-gluon exchange (OGE) interaction. This interaction is designed to generate proper hadronic spin structures.

The methodology proceeds as follows:

Covariant Structure Derivation: The authors derive the general covariant light-front structures for scalar ( $0^{++}$ ), axial vector ( $1^{++}$ and $1^{+-}$ ), and tensor ( $2^{++}$ ) quarkonia. They parameterize the Light-Front Wave Functions (LFWFs) using invariant functions ( $\phi_i$ ) that respect discrete symmetries, specifically mirror parity ( $mP$ ) and charge conjugation ( $C$ ). Unlike non-relativistic models which rely solely on Melosh rotations, this approach incorporates explicit light-front orientation-dependent terms to ensure strict Lorentz covariance.
Symmetry Constraints: The analysis highlights that relativistic spin-orbit coupling mixes states with different $L$ and $S$ . For instance, while the non-relativistic limit predicts pure $^3P_0$ , $^3P_1$ , $^1P_1$ , and $^3P_2$ configurations, the relativistic framework allows for partial-wave mixing (e.g., S/D-wave mixing in axial vectors) even within the valence $q\bar{q}$ sector, constrained by angular momentum addition and parity rules.
Numerical Implementation: The LFWFs are obtained as eigenfunctions of the light-front invariant mass squared operator. The wave functions are expanded in a basis of transverse (Laguerre polynomials) and longitudinal (Jacobi polynomials) modes. The basis size is truncated ( $N_{max}$ ) to impose an ultraviolet cutoff ( $\Lambda_{UV}$ ), which is tuned to be roughly equivalent to the mass of the quarkonium state ( $\Lambda_{UV} \approx M_{Q\bar{Q}}$ ).
Extraction of LCDAs: The leading-twist LCDAs are extracted by projecting the LFWFs onto specific spin configurations (specifically those with zero orbital angular momentum projection, $L_z=0$ ) and integrating over transverse momentum. The normalization conditions for these amplitudes are determined by the quantum numbers of the interpolating currents (vector, axial-vector, or tensor) and the charge conjugation parity of the meson.

Key Contributions and Results

Unified Formalism: The paper establishes a unified formalism where the leading-twist LCDAs of various P-wave mesons are directly related to specific $L_z=0$ components of their LFWFs.
Relativistic Partial-Wave Mixing: A central finding is the emergence of novel partial-wave components driven by relativistic dynamics. Specifically, the LFWFs of axial vector mesons ( $\chi_{c1}, \chi_{b1}$ and $h_c, h_b$ ) contain robust S/D-wave components (e.g., $\psi^{(m_j=0)}_{\uparrow\downarrow-\downarrow\uparrow}$ ) that are strictly absent in non-relativistic models.
"W"-Shaped Structures: The most notable result is the observation of distinct "W"-shaped structures in the distribution amplitudes of axial vector mesons. These shapes arise from the interference between relativistically induced S-wave and D-wave components. This feature cannot be explained by non-relativistic dynamics, which predict simpler profiles.
Convergence to Non-Relativistic Limits: As the constituent quark mass increases from charmonium to bottomonium, the computed LCDAs systematically converge toward the predictions of the non-relativistic quark model. In the heavy quark limit, the longitudinal and transverse distributions for tensor mesons converge to a single universal function, and the various P-wave LCDAs simplify into symmetric or antisymmetric sets based on longitudinal parity.
Normalization Nuances: The study clarifies normalization conditions for LCDAs based on the C-parity of the meson and the current. For C-even scalar and tensor mesons, the local matrix elements of the relevant currents vanish, requiring normalization to the first moment of the distribution amplitude. Conversely, for the C-odd axial vector mesons ( $h$ ), the tensor current matches the meson's C-parity, resulting in a non-zero local matrix element that allows the transverse LCDA to be normalized to its zeroth moment.

Significance
The paper claims that these findings provide essential non-perturbative inputs for analyzing hard exclusive processes involving P-wave quarkonia. The identification of "W"-shaped structures and other relativistic features highlights the importance of using relativistic frameworks, such as BLFQ, to understand heavy quarkonium production and structure at modern high-energy colliders (e.g., BESIII, Belle, LHC). The results serve as a quantitative signature of relativistic dynamics that can, in principle, be probed experimentally. The authors position this work as a step toward bridging fundamental dynamics with experimental observables, with future extensions planned for calculating transition form factors and higher-twist distributions.

More like this