Two-Loop DGLAP Splitting Functions from Light Cone Perturbation Theory

We perform a two-loop calculation in Light Cone Perturbation Theory (LCPT) to evaluate the next-to-leading order nonsinglet splitting function. Our calculation demonstrates the methodology and feasibility of performing higher order calculations in LCPT. Since in Hamiltonian perturbation theory the longitudinal $k^+$ momentum is always positive, poles in $1/k^+$ can be regularized by a simple cutoff which cancels in physical results, without any associated ambiguities. For transverse momentum integrals we use dimensional regularization. Developing methods for loop calculations in LCPT paves the way for a systematical, automatizable procedure for precision calculations in this framework with a transparent physical partonic interpretation. This can provide a standard framework in higher order calculations in the gluon saturation regime of QCD.