Order-$v^2$ relativistic corrections to heavy-quark fragmentation into $P$-wave quarkonium states

Within the framework of nonrelativistic QCD (NRQCD) factorization,and based on the Collins–Soper operator definition of fragmentation functions, we present a systematic calculation of the fragmentation functions for a heavy quark fragmenting into color-singlet $P$-wave quarkonium states. After reproducing and confirming the known leading-order results, we further compute the relativistic corrections up to order $\mathcal{O}(v^{2})$. Our analysis applies both to quarkonium systems composed of heavy quarks with the same flavor and to $B_c$-type mesons formed by heavy quarks of different flavors. Numerical results show that, for all color-singlet $P$-wave channels, the $\mathcal{O}(v^{2})$ relativistic corrections give sizable negative contributions over most of the momentum-fraction $z$ region. We further compute inclusive cross sections for $P$-wave quarkonium plus charmed hadrons in $e^+e^-$ annihilation via the single photon process up to $\mathcal{O}(v^{2})$ by applying our obtained fragmentation functions, and the resulting predictions are consistent with the full fixed-order results in the high-energy region.

💡 Research Summary

**
In this work the authors present a comprehensive calculation of the fragmentation functions describing a heavy quark turning into a color‑singlet (P)-wave quarkonium state, within the framework of non‑relativistic QCD (NRQCD) factorization and using the Collins–Soper operator definition of fragmentation functions. After reproducing the known leading‑order (LO) results, they systematically derive the relativistic corrections of order (v^{2}) (where (v) is the relative velocity of the heavy quark–antiquark pair) for all relevant channels: the spin‑singlet (^{1}P_{1}) and the spin‑triplet (^{3}P_{J}) with (J=0,1,2).

The calculation is performed in light‑cone coordinates and in the axial gauge ((\hat n!\cdot!A=0)), which eliminates the Wilson line in the Collins–Soper definition and reduces the contributing diagrams to a single gluon‑exchange topology. The heavy‑quark pair’s momenta are parametrized by the total momentum (P) and the relative momentum (q); the latter is expanded in powers of (q^{2}\sim v^{2}). The authors construct the appropriate spin‑projectors for both equal‑mass and unequal‑mass systems (the latter covering the (B_{c}) family) and expand the amplitude up to (O(q^{2})). The squared amplitude is then matched to the NRQCD operator expansion, which contains the leading‑order four‑fermion operators (O^{(0)}) and the relativistic correction operators (P^{(0)}).

Through this matching they extract the short‑distance coefficients (SDCs) (F_{n}(z)) (LO) and (G_{n}(z)) (relativistic correction) for each channel. The SDCs are presented analytically in the appendices for both equal‑mass and unequal‑mass cases. Numerical evaluation shows that the (G_{n}(z)) functions are negative over most of the momentum‑fraction (z) range and can reach 20–30 % of the magnitude of the LO term, especially around (z\sim0.4)–(0.6). Consequently, the relativistic corrections lead to sizable reductions of the fragmentation probabilities compared with the pure LO predictions.

To demonstrate phenomenological relevance, the authors apply the newly obtained fragmentation functions to the process (e^{+}e^{-}\to\gamma^{*}\to H+X_{c}), where (H) is a (P)-wave quarkonium state and (X_{c}) denotes a charmed hadron. They compute the inclusive cross sections by integrating the fragmentation functions over (z) and compare the results with full fixed‑order (NLO) calculations. In the high‑energy regime ((\sqrt{s}) much larger than the heavy‑quark mass) the fragmentation approximation reproduces the fixed‑order results very well, confirming that fragmentation dominates the production of high‑(p_T) quarkonia.

Overall, the paper fills a notable gap in the theoretical description of heavy‑quark fragmentation into (P)-wave quarkonia by providing the first complete (O(v^{2})) relativistic corrections for both equal‑mass and unequal‑mass systems. The results are essential for improving the precision of quarkonium production predictions at current and future colliders (LHC, RHIC, Belle II) and lay the groundwork for higher‑order studies, such as (O(v^{4})) relativistic effects or NNLO QCD corrections.