Hybrid Renormalization with Gradient Flow for Baryon Quasi-Distribution Amplitudes

Basic Information

• Paper ID: 2510.14325
• Title: Hybrid Renormalization with Gradient Flow for Baryon Quasi-Distribution Amplitudes
• Authors: Jia-Lu Zhang, Mu-Hua Zhang
• Classification: hep-ph (High Energy Physics - Phenomenology)
• Publication Date: October 16, 2025
• Paper Link: https://arxiv.org/abs/2510.14325

Abstract

This paper establishes factorization relations between baryon quasi-distribution amplitudes (quasi-DAs) defined via gradient flow and their renormalized counterparts in the $\overline{MS}$ scheme. Beyond the small flow-time limit, the authors perform complete one-loop calculations, obtaining the full matching kernel and associated linear divergences of the Wilson lines in the flowed quasi-distribution amplitudes. Based on this relation, the authors formulate a hybrid renormalization scheme combining ratio schemes with gradient flow matching. This framework enables systematic connection between lattice QCD matrix elements and continuum baryon quasi-distribution amplitudes, reducing renormalization uncertainties and providing explicit guidance for practical lattice analyses.

Research Background and Motivation

Core Problems

1. Non-perturbative nature of baryon distribution amplitudes: Baryon light-cone distribution amplitudes (LCDAs) are key non-perturbative quantities describing the internal structure of hadrons, encoding how valence quarks share longitudinal momentum within the baryon. However, since their definition involves non-local light-cone correlations, they cannot be directly calculated in Euclidean lattice QCD.
2. Technical challenges in lattice calculations:
   • Signal-to-noise ratio problems in non-local correlation functions, particularly severe at large hadron momentum and long Wilson line separations
   • Renormalization difficulties: continuum schemes such as $\overline{MS}$ cannot be directly implemented on the lattice, requiring intermediate non-perturbative schemes
   • Baryon quasi-distribution amplitudes depend on two independent momentum fractions, making them two-dimensional quantities that significantly amplify signal-to-noise ratio problems

Research Significance

Baryon LCDAs play crucial roles in QCD factorization theorems, particularly in baryon form factors at large momentum transfer and exclusive meson production processes. Understanding the parton structure of baryons is essential for comprehending QCD from first principles.

Limitations of Existing Methods

1. Traditional hybrid renormalization schemes: While theoretically consistent, they are numerically demanding, particularly for the two-dimensional nature of baryon distribution amplitudes
2. Self-renormalization procedures: Require parameter fitting in equations, introducing additional systematic uncertainties
3. RI/MOM scheme: Requires precise gauge fixing to control systematic uncertainties

Core Contributions

1. Establishes factorization relations for gradient flow baryon quasi-distribution amplitudes: Connecting flowed quasi-distribution amplitudes with their counterparts in the $\overline{MS}$ scheme
2. Completes full one-loop calculations beyond the small flow-time limit: Obtaining the complete matching kernel and Wilson line linear divergences
3. Proposes an improved hybrid renormalization scheme: Replacing self-renormalization with gradient flow matching, improving signal-to-noise ratio and achieving theoretically controlled matching
4. Provides a systematic framework for lattice QCD implementation: Reducing renormalization uncertainties and providing practical analysis guidance

Methodology Details

Task Definition

Establishing systematic connection from baryon quasi-distribution amplitudes calculated in lattice QCD to continuum baryon light-cone distribution amplitudes, with input being non-local matrix elements calculated on the lattice and output being physical baryon distribution amplitudes.

Theoretical Framework

1. Definitions of Baryon LCDAs and Quasi-Distribution Amplitudes

Baryon light-cone distribution amplitudes are defined through non-local hadron-vacuum matrix elements: $H(z_1, z_2, z_3)_{\alpha\beta\gamma} = \epsilon_{ijk}\langle 0|f^{i'}_\alpha(z_1n)W^{i'i}(z_1n, z_0n) \times g^{j'}_\beta(z_2n)W^{j'j}(z_2n, z_0n) \times h^{k'}_\gamma(z_3n)W^{k'k}(z_3n, z_0n)|B(P)\rangle$

The corresponding quasi-distribution amplitudes are defined through equal-time correlators: $\tilde{\Phi}_A(z_1, z_2, P^z, \mu) f_A P^\nu u_\Lambda(P) = \langle 0|u^T(z_1n_z)(C\gamma_5\gamma_\nu)d(z_2n_z)s(0)|\Lambda(P)\rangle$

2. Gradient Flow Formalism

Gradient flow evolves gauge and fermionic fields by solving diffusion equations: $\partial_t B_\mu = D_\nu G_{\nu\mu} + \kappa D_\mu \partial_\nu B_\nu$$\partial_t \chi = D_\mu D_\mu \chi - \kappa(\partial_\mu B_\mu)\chi$

The flow time $t$ introduces a physical scale $\sqrt{8t}$, systematically suppressing ultraviolet fluctuations while preserving long-distance physics.

3. Factorization Relations

The core factorization relation is: $\tilde{\Phi}^R(z_1, z_2, P^z, t) = e^{\delta m \tilde{z}} C_q(t, z_1, z_2, \mu) \tilde{\Phi}^{\overline{MS}}(z_1, z_2, P^z, \mu)$

where $e^{\delta m \tilde{z}}$ represents linear divergences from Wilson line self-energy corrections: $\delta m = -\frac{\alpha_s}{4\pi}C_F\frac{\sqrt{2\pi}}{\sqrt{t}} + O(\alpha_s^2)$

One-Loop Calculation Details

Matching Kernel Calculation

The complete finite flow-time matching kernel is: $C_q(z_1, z_2, t, \mu) = 1 - \frac{C_F\alpha_s}{16\pi}[18\log(2\mu^2 te^{\gamma_E}) + 6\log(432) + 12 + 7\text{Ei}(-\bar{z}_1^2) + 7\text{Ei}(-\bar{z}_2^2) + 6\text{Ei}(-\bar{z}_{12}^2) + \text{flow-time correction terms}] + O(\alpha_s^2)$

where $\bar{z}_i \equiv z_i/\sqrt{8t}$ and Ei denotes the exponential integral function.

Wilson Line Linear Divergences

The total linear divergence from self-energy diagrams is: $\delta m \tilde{z} = -\frac{C_F\alpha_s}{2\sqrt{\pi}}(\bar{z}_1 \text{erf}(\bar{z}_1) + \bar{z}_2 \text{erf}(\bar{z}_2) + \bar{z}_{12} \text{erf}(\bar{z}_{12}))$

Hybrid Renormalization Scheme

Region Partition Strategy

The calculation domain is partitioned based on spatial separations $(z_1, z_2)$ into:

• Fully perturbative region: $|z_1| < z_s, |z_2| < z_s$
• $z_1$ perturbative region: $|z_1| < z_s, |z_2| > 2z_s$
• $z_2$ perturbative region: $|z_1| > 2z_s, |z_2| < z_s$
• $z_1-z_2$ perturbative region: $|z_1-z_2| < z_s$
• Fully non-perturbative region: Other cases

Renormalization Prescription

Appropriate ratio constructions are applied to different regions, for example in the fully perturbative region: $\hat{\Phi}^H(z_1, z_2, P^z, \mu) = \frac{\hat{\Phi}^{\overline{MS}}(z_1, z_2, P^z, \mu)}{\hat{\Phi}^{\overline{MS}}(z_1, z_2, 0, \mu)}$

This eliminates short-distance logarithms $\ln z_1^2, \ln z_2^2, \ln(z_1-z_2)^2$.

Experimental Setup

Theoretical Calculation Verification

1. Small flow-time limit verification: The matching kernel simplifies to results from the auxiliary field method as $t \to 0$
2. Matching coefficients for single $\bar{h}_v\psi$ operators are exactly one-third of the three-operator case
3. Linear divergence structure recovers the usual $|z|$ dependence in the small flow-time limit

Numerical Stability Analysis

Figure 2 demonstrates comparison between the complete matching kernel and the small flow-time limit:

• Variation with flow time at fixed $\bar{z}_1 = 1, \bar{z}_2 = 2$
• The ratio $C_q(z_1,z_2,\mu,t)/\lim_{t\to 0}C_q(z_1,z_2,\mu,t)$ shows significant deviations at $z_1 = z_2$ and $z_2 = 0$

Experimental Results

Main Theoretical Results

1. Complete one-loop matching kernel: Provides exact expressions beyond the small flow-time approximation
2. Wilson line linear divergences: Gives complete treatment within the gradient flow framework
3. Optimized region partition: Eliminates singularities in the matching kernel through ratio schemes

Method Verification

1. Auxiliary field method consistency: Small flow-time limit completely agrees with known results
2. Renormalization group properties: Matching kernel satisfies expected renormalization group equations
3. Physical scale separation: The hierarchy $\sqrt{8t} \ll |z_1|, |z_2|, |z_1-z_2| \ll \Lambda_{QCD}^{-1}$ is maintained

Numerical Improvement Effects

Advantages of the gradient flow method compared to traditional approaches:

• Improved signal-to-noise ratio for non-local operators
• Eliminated need for additional renormalization
• Provides gauge-invariant ultraviolet regularization

Related Work

Development of Quasi-Distribution Methods

1. LaMET framework: Proposed by Ji et al., connecting lattice and continuum through large-momentum effective theory
2. Pseudo-distribution method: Alternative approach proposed by Radyushkin
3. Gradient flow applications in quasi-PDFs: Factorization relations for quasi-PDFs established by Brambilla et al.

Baryon Distribution Amplitude Research

1. Theoretical development: Systematic studies by Braun et al. establishing theoretical foundations of baryon LCDAs
2. Lattice calculations: Recent significant progress in lattice calculations of baryon LCDAs
3. Renormalization schemes: Development of hybrid and self-renormalization schemes

Conclusions and Discussion

Main Conclusions

1. Establishes complete factorization relations between gradient flow baryon quasi-distribution amplitudes and the $\overline{MS}$ scheme
2. Provides exact one-loop matching kernels beyond the small flow-time limit
3. Develops numerically stable and theoretically controlled hybrid renormalization scheme
4. Provides systematic framework for lattice QCD implementation

Limitations

1. One-loop approximation: Current calculations limited to one-loop level; higher-order corrections may be important
2. Flow-time selection: Requires balance between $t\Lambda_{QCD}^2 \ll 1$ and $|z_i| \gg \sqrt{8t}$
3. Linear renormalization group coefficients: Wilson line linear divergences affected by linear renormalization group coefficients, requiring leading-order renormalization group summation

Future Directions

1. Lattice numerical implementation: Applying theoretical framework to actual lattice QCD calculations
2. Higher-order calculations: Extending to two-loop level for improved precision
3. Other hadrons: Generalizing method to other baryon and meson distribution amplitudes

In-Depth Evaluation

Strengths

1. Theoretical rigor: Provides complete one-loop calculations beyond small flow-time approximation
2. Methodological innovation: Cleverly combines gradient flow and ratio schemes, addressing unique challenges of baryon quasi-distribution amplitudes
3. Practical value: Provides operational theoretical framework for lattice QCD calculations
4. Computational completeness: Detailed contributions from all Feynman diagrams and final results

Weaknesses

1. Complexity: Two-dimensional renormalization scheme remains quite complex; practical implementation may face challenges
2. Parameter tuning: Selection of hybrid cutoff $z_s$ requires careful balance of different systematic effects
3. Higher-order uncertainties: Magnitude of higher-order corrections to one-loop results unknown

Impact

1. Theoretical contribution: Provides important theoretical tools for lattice QCD research on baryon structure
2. Methodological value: Application of gradient flow in renormalization of non-local operators has universal significance
3. Experimental relevance: Will ultimately contribute to understanding baryon internal structure, relevant to experimental measurements

Applicable Scenarios

This method is particularly suitable for:

• Lattice QCD calculations of baryon distribution amplitudes
• Research on non-local operators requiring high-precision renormalization
• Applications of gradient flow methods in hadron physics

References

The paper cites 101 related references, covering important works in quasi-distribution methods, gradient flow theory, lattice QCD calculations and related fields, demonstrating comprehensive understanding and accurate grasp of the field's development.