Two-Loop Renormalization-Group Evolution for the Nucleon Distribution Amplitude

Technical Summary: Two-Loop Renormalization-Group Evolution for the Nucleon Distribution Amplitude

Problem Statement
The nucleon light-cone distribution amplitude (LCDA), $\Phi_N$, is a fundamental non-perturbative object required for the systematic description of hard exclusive reactions involving nucleons, such as electromagnetic form factors and semileptonic decays (e.g., $\Lambda_b \to p \ell \bar{\nu}_\ell$). While the one-loop (leading-order) renormalization-group (RG) evolution kernel for the leading-twist nucleon distribution amplitude was established over forty years ago, the next-to-leading-order (NLO) QCD correction to this three-particle RG kernel has remained elusive. This gap prevents the completion of next-to-leading-logarithmic (NLL) corrections for nucleon form factors within the hard-collinear factorization framework. The primary challenge lies in the technical complexity of performing two-loop ultraviolet (UV) renormalization for non-local baryonic operators, particularly due to the presence of evanescent operators that vanish in four dimensions but are necessary for consistent dimensional regularization.

Methodology
The authors employ a modern effective field theory approach to compute the NLO evolution kernel for the leading-twist nucleon distribution amplitude. The methodology proceeds through the following steps:

1. Operator Basis and Renormalization: The study utilizes a renormalized three-particle light-ray operator matrix element. To handle UV divergences, the basis of collinear operators is enlarged to include evanescent operators ($O_2, O_3$) alongside the physical operator ($O_1$). The renormalized physical operator is expressed as a linear combination of bare operators, involving a matrix of renormalization constants $Z_{ij}$.
2. Two-Loop Calculation: The authors compute the QCD matrix elements of these operators at the two-loop order ($O(\alpha_s^2)$) using dimensional regularization to capture UV divergences and a non-vanishing mass for internal particles to regulate infrared singularities.
   • A total of 70 Feynman diagrams were generated for the physical operator $\Pi_1$.
   • Vector and tensor integrals were decomposed using the Passarino-Veltman method.
   • Dirac and color algebra were reduced using QCD equations of motion and on-shell conditions.
   • The resulting scalar integrals were reduced to a set of 20 two-loop master integrals using integration-by-parts relations and the Laporta algorithm (implemented via the FIRE package).
3. Derivation of the Evolution Kernel: The NLO evolution kernel $H^{(1)}$ is derived using a "master formula" that relates the kernel to the renormalization constants $Z_{ij}$. Crucially, the calculation accounts for the mixing between evanescent and physical operators, specifically the finite renormalization constant $Z^{(1,0)}_{21}$, which is essential for determining the correct two-loop anomalous dimension despite the vanishing of evanescent operators in $D=4$.
4. Analytic Solution: To solve the resulting integro-differential evolution equation, the authors apply a conformal partial wave expansion. The nucleon distribution amplitude is expanded in terms of orthogonal polynomials $P_{Mm}$, which are eigenfunctions of the leading-order kernel $H^{(0)}$. This transforms the RG equation into a system of ordinary differential equations for the local moments $\Psi_{Mm}$.
5. Scheme Conversion: The authors derive the two-loop matching relation between their renormalization prescription (the "EO scheme," containing evanescent operators) and the established Krankl-Manashov (KM) scheme.

Key Contributions and Results

• First Determination of the Two-Loop Kernel: The paper presents the first explicit calculation of the two-loop RG evolution kernel ($H^{(1)}$) for the leading-twist nucleon distribution amplitude. The kernel is expressed in terms of primitive kernels ($V^{(1), n}_{LC}, V^{(1), n}_{2P}, V^{(1), n}_{3P}$) involving color factors $C_F, C_A, \beta_0$.
• Analytic Evolution Solution: An analytic solution for the scale dependence of the normalization coefficients and shape parameters is constructed at the NLL accuracy. The evolution matrix is derived explicitly, including the necessary anomalous dimension matrices $L^{(0)}$ and $L^{(1)}$.
• Scheme Independence: The work provides the explicit two-loop conversion factor between the EO and KM schemes, confirming the scheme independence of physical observables like the nucleon electromagnetic form factor.
• Numerical Implications:
  • The authors analyze the RG evolution using three sample models for the initial distribution amplitude (COZ, LAT25, and ABO1) at a reference scale $\mu_0 = 1.0$ GeV.
  • They find that the inclusion of NLL corrections leads to noticeable impacts (approximately 20%) on the normalized shape parameters at intermediate renormalization scales ($\mu \in [3.0, 10.0]$ GeV).
  • These effects are significantly more pronounced for the nucleon distribution amplitude than for the pion or B-meson distribution amplitudes.
  • The NLL corrections substantially affect the theoretical predictions for the Dirac nucleon form factors ($F_1^p$ and $F_1^n$), particularly at higher momentum transfers ($Q^2$). The neutron form factor is observed to be more sensitive to these resummation effects than the proton form factor.

Significance
The paper claims that determining the two-loop RG evolution of the nucleon distribution amplitude constitutes the "last missing ingredient" for achieving complete NLL corrections to nucleon form factors in the hard-collinear factorization framework. By providing the analytic solution and the explicit kernel, the authors enable more precise theoretical predictions for flagship hadron observables. The work demonstrates that neglecting these two-loop effects can lead to significant deviations in phenomenological predictions, thereby justifying the necessity of full NLO computations for the nucleon distribution amplitude. The authors suggest that extending this RG analysis to the full baryon octet and decuplet will be beneficial for exploring the partonic structure of composite hadrons and improving the precision of QCD descriptions for various hard exclusive reactions.