Nonrelativistic QCD Factorization Framework

• Nonrelativistic QCD factorization is a framework that separates perturbative short-distance processes from universal long-distance matrix elements in heavy-quarkonium production and decay.
• It utilizes an effective field theory expansion in the heavy quark’s relative velocity, enabling systematic inclusion of higher-order QCD corrections.
• The approach underpins precision phenomenology through global fits and detailed predictions for cross sections, polarization, and endpoint production challenges.

Nonrelativistic QCD (NRQCD) factorization is a theoretical framework developed to describe the production and decay of heavy-quarkonium states in Quantum Chromodynamics (QCD) where the heavy-quark mass $m_Q$ generates a hierarchy of energy scales. The factorization rigorously separates the perturbative short-distance dynamics, calculable in QCD, from universal long-distance matrix elements (LDMEs) that encode nonperturbative hadronization, with an effective field theory (EFT) expansion in the relative velocity $v$ of the heavy quark in the bound state. NRQCD factorization underlies the modern approach to $J/\psi$, $\Upsilon$, and other quarkonium production and decay, and has been the subject of extensive analytic and phenomenological studies across hadroproduction, photoproduction, lepto/photoproduction, and heavy-flavor decays.

1. Core Structure of NRQCD Factorization

The NRQCD factorization framework formulates the inclusive production (or decay) cross section of a heavy-quarkonium state $H$ in a high-energy process as a sum over intermediate $Q\bar Q$ Fock states $n$, labeled by their spin, angular momentum, and color:

$$\sigma[H] = \sum_n \hat\sigma[Q\bar Q(n) + X]\, \langle O^H(n)\rangle$$

where

• $\hat\sigma[Q\bar Q(n) + X]$ is the perturbatively calculable short-distance coefficient (SDC) for producing a $Q\bar Q$ pair in state $n$ at a scale $\mu \sim m_Q$,
• $\langle O^H(n)\rangle$ is the nonperturbative LDME governing the hadronization of the $Q\bar Q$ pair in state $n$ into the physical quarkonium $H$,
• the sum over $n$ includes both color-singlet (CS) and color-octet (CO) quantum numbers and runs up to a stated order in the $v$-expansion.

This structure is valid for production (as in $pp\to H+X$, $\gamma p\to H+X$, $e^+e^-\to H+X$, etc.) and for quarkonium decays ($H\to X$). The LDMEs are organized according to NRQCD velocity scaling rules and, in principle, are universal—independent of the process—up to corrections higher order in $v$ or $\alpha_s$ (Brambilla et al., 2024, Butenschoen et al., 2012, Brambilla et al., 2022, Butenschoen et al., 2011, Sun et al., 2017, Butenschoen et al., 2012, Zheng et al., 2021, He et al., 2019).

2. Short-Distance Coefficients and Operator Expansion

Short-distance coefficients are determined by matching full QCD onto NRQCD at the scale $\mu\sim m_Q$, integrating out high-momentum modes. The leading SDCs for each Fock state $n$ are calculated perturbatively in $\alpha_s$, with higher-order QCD corrections (NLO, NNLO) improving accuracy and allowing for systematic uncertainty estimates:

• At leading order (LO), SDCs correspond to Born-level (tree) diagrams for $Q\bar Q(n)$ production.
• Next-to-leading order (NLO) includes virtual (one-loop) corrections and real-emission processes; both ultraviolet (UV) and infrared (IR) singularities must be handled by a combination of renormalization (typically $\overline{\rm MS}$) and cancellation between real and virtual diagrams (Brambilla et al., 2024, Butenschoen et al., 2011, Butenschoen et al., 2012, Butenschoen et al., 2022).
• P-wave (and higher) channels require operator renormalization to absorb IR poles into the appropriate LDMEs (Butenschoen et al., 2014).

The operator expansion organizes four-fermion NRQCD operators according to their velocity and color structure:

$$\mathcal{O}^H(n) = \chi^\dagger \mathcal{K}_n \psi\, a_H^\dagger a_H\, \psi^\dagger \mathcal{K}_n' \chi$$

Here, $\psi, \chi$ are Pauli spinors for the heavy quark and antiquark, $\mathcal{K}_n$ project onto spin–angular-momentum–color state $n$. LDMEs are expectation values of these operators in the QCD vacuum. Fock-state decomposition up to a given power of $v$ (typically $v^4$ relative to the CS) is standard (Brambilla et al., 2022, Sun et al., 2017).

3. Long-Distance Matrix Elements and Velocity Scaling

LDMEs encode the soft, nonperturbative transition probability for a $Q\bar Q(n)$ pair to hadronize into $H$. The velocity-scaling rules provide the parametric suppression associated with each channel:

• CS $S$-wave: $\langle O^{H}(^3S_1^{[1]})\rangle \sim v^3$
• CO $S$- and $P$-waves: $$\langle O^{H}(^1S_0^{[8]}), O^{H}(^3S_1^{[8]}), O^{H}(^3P_J^{[8]})\rangle \sim v^7$$

The precise numerical values of LDMEs are extracted from global fits to data, with heavy-quark spin symmetry (HQSS) providing further relations, e.g.:

$$\langle O^{\eta_c}(^3S_1^{[8]}) \rangle = \langle O^{J/\psi}(^1S_0^{[8]}) \rangle$$

$$\langle O^{\eta_c}(^1S_0^{[1]}) \rangle = \frac{1}{3} \langle O^{J/\psi}(^3S_1^{[1]}) \rangle$$

This systematics enables determination of LDMEs for various quarkonia by mapping from well-constrained states such as $J/\psi$ (Butenschoen et al., 2014, Butenschoen et al., 2022, Zheng et al., 2021). In potential NRQCD (pNRQCD), the LDMEs can be reduced to expressions involving the wavefunction at the origin and three universal gluonic correlators, further enhancing predictive power (Brambilla et al., 2022).

4. Applications: Inclusive Production, Polarization, and Decay

NRQCD factorization underlies the modern computation of cross sections and polarization observables for a variety of processes, including:

• Inclusive hadroproduction (e.g., $pp\to J/\psi+X$) (Brambilla et al., 2024, Butenschoen et al., 2022)
• Photoproduction and lepto/photoproduction ($\gamma p \to J/\psi + X$, $ep\to eJ/\psi X$) (Butenschoen et al., 2011, Sun et al., 2017, Butenschoen et al., 2012)
• Inclusive production in $Z$ and $\Upsilon$ decays, where the formalism is applied to total and differential rates (Zheng et al., 2021, He et al., 2019)
• Polarization measurements, with theory predictions for parameters $\lambda_\theta$, $\lambda_\phi$, and $\lambda_{\theta\phi}$, obtained from spin-density matrix elements $d\sigma_{ij}$ (helicity/Collins-Soper/target frames), are computed for comparison to experimental lepton angular distributions (Butenschoen et al., 2011, Brambilla et al., 2024, Butenschoen et al., 2012, Butenschoen et al., 2012)

Key features include:

• Marked improvement in the fit to global cross-section and polarization data upon inclusion of NLO corrections and all $O(v^4)$ CO contributions (Butenschoen et al., 2011, Butenschoen et al., 2012, Butenschoen et al., 2022).
• NLO NRQCD can achieve simultaneous description of high-$p_T$ hadroproduction, photoproduction (away from endpoint regions), and some polarization observables, although challenges remain—most notably in reconciling Tevatron $J/\psi$ polarization data (Butenschoen et al., 2012, Butenschoen et al., 2012, Brambilla et al., 2024).

5. Universality, Precision Fits, and Theoretical Developments

A central conjecture of NRQCD factorization is that LDMEs are universal, i.e., process-independent. Recent works have systematically tested this by:

• Global fits to $\sim 1000$ data points for $J/\psi$, $\psi(2S)$ production, fitting the three main CO LDMEs at NLO (Butenschoen et al., 2022, Butenschoen et al., 2012).
• Applying fit-and-predict strategies that combine scale variations with covariance tracking in uncertainty bands, enhancing robustness of predictions across different processes (Brambilla et al., 2024).
• Using pNRQCD to reduce the LDME set to a minimal number of universal gluonic correlators plus wavefunction-at-origin, which in principle could be computed on the lattice (Brambilla et al., 2022).

The numerical extraction of LDMEs, especially for the $^3P_J^{[8]}$ channel, is delicate, requiring high-quality data at high $p_T$ to avoid negative cross sections or large cancellations (Brambilla et al., 2024, Butenschoen et al., 2022). In summary, the inclusion of all S and P-wave CO states at NLO, a fit-and-predict approach, and the pNRQCD minimal parametrization, together enable high-precision phenomenology and highlight the need for further theoretical work in endpoint (small $p_T$, $z\to1$) regimes.

6. Limitations, Extensions, and Open Challenges

The NRQCD factorization framework, while broadly successful, exhibits well-defined limitations and areas of ongoing research:

• In low transverse momentum ($p_T$) regimes, multiple soft scatterings and soft-gluon emissions violate simple factorization, necessitating transverse-momentum-dependent (TMD) factorization and the introduction of TMD shape functions (TMDShFs) as the appropriate nonperturbative input in the small-$q_T$ limit. The TMDShFs generalize LDMEs to functions of $q_T$, encode soft radiation, and reunite with standard LDMEs in the high-$p_T$ limit via OPE matching (Echevarria et al., 2024).
• Fixed-order NLO NRQCD fails to accurately describe production in endpoint regions ($z\to1$, $p_T\ll m_Q$) due to large logarithms and nonperturbative effects not captured by the OPE. LP/threshhold resummation and the inclusion of nonperturbative shape functions are required for accurate phenomenology (Brambilla et al., 2024, Butenschoen et al., 2012).
• Tensions remain in polarization observables, in particular strong discrepancies between NLO NRQCD predictions and Tevatron CDF polarization data for $J/\psi$, suggesting either a breakdown of universality or missing higher-order or nonperturbative mechanisms (Butenschoen et al., 2012, Butenschoen et al., 2012, Butenschoen et al., 2011).
• In certain processes, such as $\eta_c$ production at LHCb, NLO NRQCD with LDMEs determined from $J/\psi$ phenomenology drastically overshoots the data, while the color-singlet model succeeds, challenging the universality hypothesis (Butenschoen et al., 2014).

7. Outlook and Future Directions

Several directions are currently being pursued to address the challenges and extend the predictive scope of NRQCD factorization:

• Calculation of NNLO SDCs and systematic inclusion of large logarithms via LP/threshhold resummation to improve theoretical control at large $p_T$ and near kinematic endpoints (Brambilla et al., 2024).
• Development of TMD factorization and the calculation of TMDShFs at higher order, including their evolution and matching to LDMEs (Echevarria et al., 2024).
• Lattice calculations of pNRQCD universal correlators, aiming for first-principles determinations of the minimal set of nonperturbative parameters (Brambilla et al., 2022).
• Comprehensive testing of universality across diverse quarkonium states, decay and production modes, and in new associated production channels to further constrain the structure of LDMEs (Zheng et al., 2021, He et al., 2019, He et al., 2016).
• Precision measurements at emerging facilities such as the Electron-Ion Collider (EIC) and high-luminosity LHC will provide decisive data to refine, falsify, or extend NRQCD factorization.

NRQCD factorization remains the foundational framework for heavy-quarkonium production in high-energy QCD, supporting a broad phenomenology while motivating significant ongoing theoretical and experimental inquiry.