Skip to content

Rapidity Renormalization and Two-Scale Evolution

Dimensional regularization separates modes with different invariant masses, but it cannot distinguish modes that have the same virtuality and are separated only by a large boost. Their separate soft and collinear matrix elements can therefore contain divergent integrals over rapidity even though the complete amplitude is finite. A rapidity regulator converts those divergences into counterterms and introduces a scale ν\nu whose evolution resums logarithms of rapidity ratios.

This page derives the rapidity counterterms of a one-loop massive Sudakov form factor in the η\eta regulator, obtains the coupled μ\mu and ν\nu anomalous dimensions, and verifies path independence explicitly. The regulator and individual sector functions are scheme dependent; cancellation in the sector sum and the integrability of two-scale evolution are the physical tests.

Required background. Matching onto Factorized Operator Structures supplies the soft and collinear sectors and their overlap subtractions. Evolution Kernels, Consistency Relations, and Resummation Architecture supplies the single-scale kernel logic. Rapidity Divergences, Glauber Exchange, and Factorization Limits explains why rapidity consistency does not by itself prove cancellation of Glauber exchange.

Choose null vectors n2=nˉ2=0n^2=\bar n^2=0 with n ⁣ ⁣nˉ=2n\!\cdot\!\bar n=2, and write momentum components in the order (n ⁣p,nˉ ⁣p,p)(n\!\cdot p,\bar n\!\cdot p,p_\perp). For a hard scale QQ, an infrared mass MM, and λ=M/Q1\lambda=M/Q\ll1, the relevant scalings are

pnQ(λ2,1,λ),pn2M2,pnˉQ(1,λ2,λ),pnˉ2M2,psQ(λ,λ,λ),ps2M2.\begin{aligned} p_n&\sim Q(\lambda^2,1,\lambda), &p_n^2&\sim M^2, \\ p_{\bar n}&\sim Q(1,\lambda^2,\lambda), &p_{\bar n}^2&\sim M^2, \\ p_s&\sim Q(\lambda,\lambda,\lambda), &p_s^2&\sim M^2. \end{aligned}

All three modes have the same virtuality. Their rapidities,

y(p)=12lnnˉ ⁣pn ⁣p,y(p)=\frac12\ln\left|\frac{\bar n\!\cdot p}{n\!\cdot p}\right|,

are instead centered near +ln(1/λ)+\ln(1/\lambda), ln(1/λ)-\ln(1/\lambda), and 00. On a fixed-invariant-mass hyperbola one may set k+=κeyk^+=\kappa e^y and k=κeyk^-=\kappa e^{-y}. Dimensional regularization controls the transverse and invariant-mass integrations but leaves

0dk+k+=dy\int_0^\infty\frac{dk^+}{k^+}=\int_{-\infty}^{\infty}dy

unregulated. The endpoints are the soft–collinear overlap viewed as an infinite boost, not an additional ultraviolet limit in k2k^2.

DivergenceDirection in loop momentumRegulator and scaleRenormalized evolution
UltravioletIncreasing invariant mass or transverse momentumd=42ϵd=4-2\epsilon and μ\muμd/dμ\mu\,d/d\mu
RapidityIncreasing boost at fixed invariant massη\eta and ν\nu in the scheme belowνd/dν\nu\,d/d\nu

Mixed poles such as 1/(ηϵ)1/(\eta\epsilon) are expected: a rapidity counterterm can itself have ultraviolet scale dependence. They do not identify the two divergences.

Consider the spacelike color-singlet quark current in a gauge theory whose gauge boson has mass MM, used here as an infrared regulator. Work in Feynman gauge and subtract both ϵ\epsilon and η\eta poles minimally. For an emission of momentum kk, the η\eta prescription modifies collinear and soft Wilson-line vertices schematically as

nˉμnˉ ⁣knˉμnˉ ⁣kw2(νnˉ ⁣k)η,\frac{\bar n^\mu}{\bar n\!\cdot k} \longrightarrow \frac{\bar n^\mu}{\bar n\!\cdot k}\, w^2\left(\frac{\nu}{|\bar n\!\cdot k|}\right)^\eta,

and

nμn ⁣knμn ⁣kw(ν2k3)η/2.\frac{n^\mu}{n\!\cdot k} \longrightarrow \frac{n^\mu}{n\!\cdot k}\, w\left(\frac{\nu}{|2k^3|}\right)^{\eta/2}.

The auxiliary ww keeps track of rapidity renormalization and is set to one after differentiation. The relative power of η\eta follows from the soft line having two eikonal attachments where the corresponding collinear graph has one. The order of limits is part of the scheme:

η0beforeϵ0,ηϵm0(m>0).\eta\to0\quad\text{before}\quad\epsilon\to0, \qquad \frac{\eta}{\epsilon^m}\to0\quad(m>0).

Taking the limits in the opposite order leaves the fixed-virtuality hyperbola before removing the rapidity cutoff and mixes the intended sector boundaries.

The leading-power current factorizes as

Jμ=H(Q2,μ)Jn(M;μ,ν/Q)γμJnˉ(M;μ,ν/Q)S(M;μ,ν/M).J^\mu =H(Q^2,\mu)\, J_n(M;\mu,\nu/Q)\, \gamma_\perp^\mu\, J_{\bar n}(M;\mu,\nu/Q)\, S(M;\mu,\nu/M).

Define, only for this example,

aαsCFπ=g2CF4π2,RϵeγEϵΓ(ϵ)(μM)2ϵ.a\equiv\frac{\alpha_s C_F}{\pi} =\frac{g^2C_F}{4\pi^2}, \qquad R_\epsilon\equiv e^{\gamma_E\epsilon}\Gamma(\epsilon) \left(\frac{\mu}{M}\right)^{2\epsilon}.

The rapidity-pole parts of the one-loop renormalization factors are

Zn(η)=1+a2ηRϵ,Znˉ(η)=1+a2ηRϵ,ZS(η)=1aηRϵ.\begin{aligned} Z_n^{(\eta)}&=1+\frac{a}{2\eta}R_\epsilon, \\ Z_{\bar n}^{(\eta)}&=1+\frac{a}{2\eta}R_\epsilon, \\ Z_S^{(\eta)}&=1-\frac{a}{\eta}R_\epsilon. \end{aligned}

Pure 1/ϵ1/\epsilon terms are omitted from this display. Since

Rϵ=1ϵ+lnμ2M2+O(ϵ),R_\epsilon =\frac1\epsilon+\ln\frac{\mu^2}{M^2} +\mathcal O(\epsilon),

these factors contain both the mixed pole and the finite logarithm needed for rapidity evolution. At order aa,

δZn(η)+δZnˉ(η)+δZS(η)=0.\delta Z_n^{(\eta)} +\delta Z_{\bar n}^{(\eta)} +\delta Z_S^{(\eta)}=0.

Thus the 1/η1/\eta poles and all ν\nu dependence cancel in the sector sum, which is boost invariant. The η\eta-regulated soft-bin integrals are scaleless, but they must still be defined: their role is to place the soft–collinear boundary consistently before the scaleless result is set to zero. Chiu, Jain, Neill, and Rothstein derive the regulated integrals, limit order, counterterms, and cancellation in Chiu et al. 2012, §§ 4–4.2, preprint pp. 7–14, Open PDF.

For each renormalized low-energy factor Fi{Jn,Jnˉ,S}F_i\in\{J_n,J_{\bar n},S\}, define

dlnFidlnμ=γμi,dlnFidlnν=γνi.\frac{d\ln F_i}{d\ln\mu}=\gamma_\mu^i, \qquad \frac{d\ln F_i}{d\ln\nu}=\gamma_\nu^i.

For symmetric external momenta nˉ ⁣p1=n ⁣p2=Q\bar n\!\cdot p_1=n\!\cdot p_2=Q, the one-loop anomalous dimensions are

Factorγμi/a\gamma_\mu^i/aγνi/a\gamma_\nu^i/a
JnJ_n34+lnνQ\dfrac34+\ln\dfrac{\nu}{Q}12lnμ2M2\dfrac12\ln\dfrac{\mu^2}{M^2}
JnˉJ_{\bar n}34+lnνQ\dfrac34+\ln\dfrac{\nu}{Q}12lnμ2M2\dfrac12\ln\dfrac{\mu^2}{M^2}
SSlnμ2ν2\ln\dfrac{\mu^2}{\nu^2}lnμ2M2-\ln\dfrac{\mu^2}{M^2}

Rapidity independence of the current gives

γνn+γνnˉ+γνS=0.\gamma_\nu^n+\gamma_\nu^{\bar n}+\gamma_\nu^S=0.

The hard factor has no rapidity anomalous dimension. Its virtuality anomalous dimension cancels the low-energy sum,

γμH=a(lnμ2Q2+32),γμH+γμn+γμnˉ+γμS=0.\gamma_\mu^H =-a\left(\ln\frac{\mu^2}{Q^2}+\frac32\right), \qquad \gamma_\mu^H+\gamma_\mu^n+ \gamma_\mu^{\bar n}+\gamma_\mu^S=0.

Independence of the order of renormalization requires a vanishing RG curl,

(lnμ+β(αs)αs)γνi=γμilnν.\left(\frac{\partial}{\partial\ln\mu} +\beta(\alpha_s)\frac{\partial}{\partial\alpha_s}\right) \gamma_\nu^i =\frac{\partial\gamma_\mu^i}{\partial\ln\nu}.

At one loop the running of aa contributes only at higher order. The two nontrivial checks are therefore

γνnlnμ=a=γμnlnν,γνSlnμ=2a=γμSlnν.\begin{aligned} \frac{\partial\gamma_\nu^n}{\partial\ln\mu} &=a =\frac{\partial\gamma_\mu^n}{\partial\ln\nu}, \\ \frac{\partial\gamma_\nu^S}{\partial\ln\mu} &=-2a =\frac{\partial\gamma_\mu^S}{\partial\ln\nu}. \end{aligned}

The cusp anomalous dimension fixes these slopes beyond one loop; the noncusp integration constant requires a fixed-order calculation.

First application: an explicit path-independence check

Section titled “First application: an explicit path-independence check”

Evolve JnJ_n from (μi,νi)(\mu_i,\nu_i) to (μf,νf)(\mu_f,\nu_f) at fixed one-loop aa, and abbreviate

Lμ=lnμfμi,Lν=lnνfνi.L_\mu=\ln\frac{\mu_f}{\mu_i}, \qquad L_\nu=\ln\frac{\nu_f}{\nu_i}.

Path A first changes μ\mu at fixed νi\nu_i and then changes ν\nu at fixed μf\mu_f:

lnUnA=a[(34+lnνiQ)Lμ+12lnμf2M2Lν].\ln U_n^{A} =a\left[ \left(\frac34+\ln\frac{\nu_i}{Q}\right)L_\mu +\frac12\ln\frac{\mu_f^2}{M^2}L_\nu \right].

Path B reverses the order:

lnUnB=a[12lnμi2M2Lν+(34+lnνfQ)Lμ].\ln U_n^{B} =a\left[ \frac12\ln\frac{\mu_i^2}{M^2}L_\nu +\left(\frac34+\ln\frac{\nu_f}{Q}\right)L_\mu \right].

Their difference is the integral of the RG curl over the rectangle:

lnUnAUnB=a[LνLμ+LμLν]=0+O(a2).\ln\frac{U_n^A}{U_n^B} =a\left[-L_\nu L_\mu+L_\mu L_\nu\right] =0+\mathcal O(a^2).

The soft factor gives the same result with the corresponding 2a-2a curl, and the full rapidity kernel cancels among the three low-energy sectors. With a running coupling, one must integrate the consistency equation rather than insert the fixed-order γν\gamma_\nu at a scale where ln(μ/M)\ln(\mu/M) is large.

The natural boundary points are

μHQ,μnμnˉμSM,νnνnˉQ,νSM.\mu_H\sim Q, \qquad \mu_n\sim\mu_{\bar n}\sim\mu_S\sim M, \qquad \nu_n\sim\nu_{\bar n}\sim Q, \qquad \nu_S\sim M.

Although the low modes share μM\mu\sim M, their natural rapidity scales differ. In this particular one-loop example γνi\gamma_\nu^i vanishes at μ=M\mu=M, so a path that performs the ν\nu evolution there makes that segment trivial; a path that first raises μ\mu has a nonzero rapidity kernel. The equality above shows how the logarithm is redistributed without changing the result.

Evidence status, checked through 3 August 2026. The explicit counterterms above belong to the η\eta regulator and minimal subtraction. Other valid schemes can move finite terms and even hide the auxiliary rapidity scale, so the presence of an explicit ν\nu is not itself physical.

RepresentationHow the rapidity boundary is imposedWhat changes with the scheme
η\eta rapidity RGAnalytic powers of longitudinal momentum and an explicit ν\nuSector counterterms, boundary constants, and explicit γνi\gamma_\nu^i
Analytic regulator or collinear anomalyRegulator powers can combine sectors into an all-order expression with no visible ν\nuRapidity logarithms appear in an anomaly exponent or correlated sector product
Exponential regulatorA Laplace-space energy constraint damps the rapidity endpointsThe subtraction and boundary definition; the physical rapidity kernel is fixed once the hard scheme and symmetric sector treatment are fixed
Off-light-cone or δ\delta regulatorTilted Wilson lines or shifted eikonal denominatorsNonzero overlap subtractions and the allocation of soft contributions

An all-order analytic-regulator formula can be formally rapidity-scale independent while its finite-order truncation still carries an independent rapidity-scale uncertainty. Jaiswal and Okui derive this reconciliation in Jaiswal and Okui 2015, abstract and §§ I, IV.A, preprint pp. 1–7, Open PDF. Li, Neill, and Zhu show, under symmetric beam-sector treatment and correct subtractions, how regulator dependence is confined by the hard-function scheme while boundary constants remain scheme dependent in Li, Neill, and Zhu 2020, §§ II–III, preprint pp. 5–11, Open PDF.

Recent low-xx work illustrates why regulator choice remains context dependent. Three regulators formulated in k+k^+, kk^-, or rapidity lead naturally to evolution in the corresponding variable and produce different finite-order collinear-log patterns. That result is specific to low-xx dipole factorization rather than a proof about every SCET observable; it supplies a current example of the same general separation between physical predictions and regulator-dependent organization. See Altinoluk, Beuf, and Penttala 2026, § 2.2 and § 6, preprint pp. 6–8 and 49–50, Open PDF.

The regulator-independent requirements are narrower and stronger: the regulated sector sum must reproduce the relevant full-theory region, all auxiliary poles and scales must cancel in a physical quantity, the two evolution directions must be integrable to the calculated order, and predicted symmetries must be restored when regulators are removed. No regulator can repair a missing leading mode or an uncanceled Glauber contribution.

The collinear and soft-II points on the line a+b=2a+b=2 share one virtuality but occupy different rapidities. This is precisely when the ν\nu direction in the final evolution box becomes necessary. The diagram also emphasizes the prior checks: homogeneous mode expansion, multipole expansion, and overlap subtraction must be correct before a rapidity counterterm has a physical interpretation.

Collinear and soft-II modes lie on one virtuality line at different rapidities, ultrasoft lies at lower virtuality, and the construction adds multipole and overlap tests before factorized evolution.

Mode locations are shown in the exponents aa and bb of (n ⁣p/Q,nˉ ⁣p/Q)(λa,λb)(n\!\cdot p/Q,\bar n\!\cdot p/Q)\sim(\lambda^a,\lambda^b), with the transverse exponent written in each label. The line a+b=2a+b=2 contains nn-collinear, nˉ\bar n-collinear, and soft-II scalings of virtuality Q2λ2Q^2\lambda^2; their separation along the line is a rapidity separation. Ultrasoft momentum has virtuality Q2λ4Q^2\lambda^4, while hard fluctuations are matched at Q2Q^2. The points are alternatives selected by a hierarchy and observable, not a universal simultaneous field list. A consistent construction requires homogeneous fields, multipole and overlap expansion, sector matching, and μ\mu and, when needed, ν\nu evolution, with explicit factorization checks. The diagram is schematic and not to scale.

Uncanceled 1/η1/\eta poles. Check that the mode list is complete, the soft and collinear regulator powers are compatible, and every zero-bin has been defined before being evaluated. Changing ν\nu cannot cancel a missing sector.

A nonzero μ–ν curl. Counterterms, anomalous dimensions, or perturbative orders have been mixed. Compare DμγνD_\mu\gamma_\nu with lnνγμ\partial_{\ln\nu}\gamma_\mu sector by sector before solving either evolution equation.

Residual boost or regulator dependence. The physical sector sum is incomplete or has inconsistent Wilson-line orientations and subtractions. Individual functions need not be boost invariant, but the stated observable must be.

Glauber sensitivity. Rapidity-pole cancellation is necessary but not sufficient for factorization. If a pinched Glauber region survives the measurement or cut sum, retain the corresponding operator or restrict the theorem’s scope rather than interpreting it as a regulator artifact.

  1. Why can ordinary μ\mu evolution not resum a soft–collinear rapidity logarithm when the modes have equal virtuality?

    Solution

    The modes have the same natural invariant-mass scale and hence the same natural μ\mu, but their ratios p+/pp^+/p^- differ parametrically. A second scale conjugate to that boost separation is required.

  2. What does a nonzero (μ,ν)(\mu,\nu) evolution commutator indicate?

    Solution

    It signals inconsistent anomalous dimensions, regulator counterterms, overlap subtraction, or perturbative truncation. Exact evolution between fixed endpoints cannot depend on the path in the (μ,ν)(\mu,\nu) plane.

  • Altinoluk, Tolga, Guillaume Beuf, and Jani Penttala. 2026. “Exploring Rapidity Regularization Schemes at Low xx with the DIS Longitudinal Structure Function.” Journal of High Energy Physics 2026 (4): 125. DOI. Open PDF.

  • Chiu, Jui-yu, Ambar Jain, Duff Neill, and Ira Z. Rothstein. 2012. “A Formalism for the Systematic Treatment of Rapidity Logarithms in Quantum Field Theory.” Journal of High Energy Physics 2012 (5): 084. DOI. Open PDF.

  • Jaiswal, Prerit, and Takemichi Okui. 2015. “Re-emergence of Rapidity Scale Uncertainty in SCET.” Physical Review D 92 (7): 074035. DOI. Open PDF.

  • Li, Ye, Duff Neill, and Hua Xing Zhu. 2020. “An Exponential Regulator for Rapidity Divergences.” Nuclear Physics B 960: 115193. DOI. Open PDF.