Accelerating Berends-Giele recursion for gluons in arbitrary dimensions over finite fields

This work provides a proof of concept for the computation of pure gluonic amplitudes in quantum chromodynamics (QCD) on graphics processing units (GPUs). The implementation relies on the Berends-Giele recursion algorithm and, for the first time on a GPU, enables the numerical computation of amplitudes in an arbitrary number of space-time dimensions and over finite fields. This demonstrates the advantages of hardware acceleration, not only for the computation of tree-level amplitudes for real-radiation processes in four dimensions over complex numbers but also for generating loop integrands for virtual corrections in $d$ dimensions over finite fields. The associated computer program is publicly available.