A Matrix Approach to Numerical Solution of the DGLAP Evolution Equations

Abstract

A matrix-based approach to numerical integration of the DGLAP evolution equations is presented. The method arises naturally on discretisation of the Bjorken x variable, a necessary procedure for numerical integration. Owing to peculiar properties of the matrices involved, the resulting equations take on a particularly simple form and may be solved in closed analytical form in the variable t=ln(alpha_0/alpha). Such an approach affords parametrisation via data x bins, rather than fixed functional forms. Thus, with the aid of the full correlation matrix, appraisal of the behaviour in different x regions is rendered more transparent and free of pollution from unphysical cross-correlations inherent to functional parametrisations. Computationally, the entire programme results in greater speed and stability; the matrix representation developed is extremely compact. Moreover, since the parameter dependence is linear, fitting is very stable and may be performed analytically in a single pass over the data values.

13 pages, no figures, typeset with revtex4 and uses packages: acromake, amssymb

Cited by in corpus (8)

• QCDNUM: Fast QCD Evolution and Convolution
• A Higher Order Perturbative Parton Evolution Toolkit (HOPPET)
• A posteriori inclusion of parton density functions in NLO QCD final-state calculations at hadron colliders: The APPLGRID Project
• Fast evolution of parton distributions
• ChiliPDF: Chebyshev Interpolation for Parton Distributions
• Nucleon and nuclear structure functions with non-perturbative and higher order perturbative QCD effects
• Comparative analysis of analytical solutions for $F_2^P$ at small $x$ in DGLAP approach
• Global Analyses of Generalized Parton Distributions with Diverse PDF Inputs