Analytic solution to leading order coupled DGLAP evolution equations: A new perturbative QCD tool

Journal Article · Tue Mar 01 00:00:00 EST 2011 · Physical Review. D, Particles Fields

DOI:https://doi.org/10.1103/PHYSREVD.83.054009· OSTI ID:21578185

Block, Martin M; Durand, Loyal; Ha, Phuoc; McKay, Douglas W ^[1]

---

1. Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208 (United States)

We have analytically solved the LO perturbative QCD singlet DGLAP equations [V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972)][G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977)][Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977)] using Laplace transform techniques. Newly developed, highly accurate, numerical inverse Laplace transform algorithms [M. M. Block, Eur. Phys. J. C 65, 1 (2010)][M. M. Block, Eur. Phys. J. C 68, 683 (2010)] allow us to write fully decoupled solutions for the singlet structure function F{sub s}(x,Q{sup 2}) and G(x,Q{sup 2}) as F{sub s}(x,Q{sup 2})=F{sub s}(F{sub s0}(x{sub 0}),G{sub 0}(x{sub 0})) and G(x,Q{sup 2})=G(F{sub s0}(x{sub 0}),G{sub 0}(x{sub 0})), where the x{sub 0} are the Bjorken x values at Q{sub 0}{sup 2}. Here F{sub s} and G are known functions--found using LO DGLAP splitting functions--of the initial boundary conditions F{sub s0}(x){identical_to}F{sub s}(x,Q{sub 0}{sup 2}) and G{sub 0}(x){identical_to}G(x,Q{sub 0}{sup 2}), i.e., the chosen starting functions at the virtuality Q{sub 0}{sup 2}. For both G(x) and F{sub s}(x), we are able to either devolve or evolve each separately and rapidly, with very high numerical accuracy--a computational fractional precision of O(10{sup -9}). Armed with this powerful new tool in the perturbative QCD arsenal, we compare our numerical results from the above equations with the published MSTW2008 and CTEQ6L LO gluon and singlet F{sub s} distributions [A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Eur. Phys. J. C 63, 189 (2009)], starting from their initial values at Q{sub 0}{sup 2}=1 GeV{sup 2} and 1.69 GeV{sup 2}, respectively, using their choice of {alpha}{sub s}(Q{sup 2}). This allows an important independent check on the accuracies of their evolution codes and, therefore, the computational accuracies of their published parton distributions. Our method completely decouples the two LO distributions, at the same time guaranteeing that both G and F{sub s} satisfy the singlet coupled DGLAP equations. It also allows one to easily obtain the effects of the starting functions on the evolved gluon and singlet structure functions, as functions of both Q{sup 2} and Q{sub 0}{sup 2}, being equally accurate in devolution (Q{sup 2}<Q{sub 0}{sup 2}) as in evolution (Q{sup 2}>Q{sub 0}{sup 2}). Further, it can also be used for nonsinglet distributions, thus giving LO analytic solutions for individual quark and gluon distributions at a given x and Q{sup 2}, rather than the numerical solutions of the coupled integral-differential equations on a large, but fixed, two-dimensional grid that are currently available.

• OSTI ID: 21578185
• Journal Information: Physical Review. D, Particles Fields, Vol. 83, Issue 5; Other Information: DOI: 10.1103/PhysRevD.83.054009; (c) 2011 American Institute of Physics; ISSN 0556-2821
• Country of Publication: United States
• Language: English

Related Subjects

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS
ACCURACY
ALGORITHMS
ANALYTICAL SOLUTION
BOUNDARY CONDITIONS
COMPARATIVE EVALUATIONS
DIFFERENTIAL EQUATIONS
DISTRIBUTION
EVOLUTION
GLUONS
LAPLACE TRANSFORMATION
NUMERICAL SOLUTION
QUANTUM CHROMODYNAMICS
STRUCTURE FUNCTIONS
TWO-DIMENSIONAL CALCULATIONS
BOSONS
EQUATIONS
EVALUATION
FIELD THEORIES
FUNCTIONS
INTEGRAL TRANSFORMATIONS
MATHEMATICAL LOGIC
MATHEMATICAL SOLUTIONS
QUANTUM FIELD THEORY
TRANSFORMATIONS