A C++ code to solve the DGLAP equations applied to Ultra High Energy Cosmic Rays

Abstract
We solve numerically the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations for the evolution of fragmentation functions using the Laguerre method. We extend this method to include supersymmetric evolution. The solution to the DGLAP equations is particularly interesting to calculate the expected spectra of Ultra High Energy Cosmic Rays in models where they are produced by the decay of a massive particle X, M_X> 10^(12)GeV.

Title of program: evolve
Catalogue Id: ADPQ_v1_0

Nature of problem
In order to predict the spectra of UHECR produced by the decay of a dark matter superheavy particle with mass MX one needs to calculate fragmentation functions at the energy scale MX. These can be calculated from measured low energy fragmentation functions using the DGLAP equations.

Versions of this program held in the CPC repository in Mendeley Data
ADPQ_v1_0; evolve; 10.1016/S0010-4655(01)00458-1

This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)