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

Abstract

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.

Code at https://github.com/Amps-GPU/BG-Trees Matched published version

References in corpus (36)

• Array Programming with NumPy
• Scattering Amplitudes with Open Loops
• Event Generation with Sherpa 2.2
• Generalized Unitarity and One-Loop Amplitudes in N=4 Super-Yang-Mills
• An Automated Implementation of On-Shell Methods for One-Loop Amplitudes
• Full one-loop amplitudes from tree amplitudes
• A novel approach to integration by parts reduction
• Scattering amplitudes over finite fields and multivariate functional reconstruction
• A Numerical Unitarity Formalism for Evaluating One-Loop Amplitudes
• RECOLA: REcursive Computation of One-Loop Amplitudes
• Towards a Basis for Planar Two-Loop Integrals
• Two-loop Integrand Decomposition into Master Integrals and Surface Terms
• A Two-Loop Five-Gluon Helicity Amplitude in QCD
• To $d$, or not to $d$: Recent developments and comparisons of regularization schemes
• Hepta-Cuts of Two-Loop Scattering Amplitudes
• A GPU compatible quasi-Monte Carlo integrator interfaced to pySecDec
• Two-Loop Four-Gluon Amplitudes with the Numerical Unitarity Method
• A Complete Two-Loop, Five-Gluon Helicity Amplitude in Yang-Mills Theory
• Challenges in Monte Carlo event generator software for High-Luminosity LHC
• Simulation of vector boson plus many jet final states at the high luminosity LHC
• Multi-gluon one-loop amplitudes using tensor integrals
• Caravel: A C++ Framework for the Computation of Multi-Loop Amplitudes with Numerical Unitarity
• Numerical Scattering Amplitudes with pySecDec
• Les Houches 2021: Physics at TeV Colliders: Report on the Standard Model Precision Wishlist
• Efficient Numerical Evaluation of Feynman Integral
• Subleading Poles in the Numerical Unitarity Method at Two Loops
• VegasFlow: accelerating Monte Carlo simulation across multiple hardware platforms
• Design and engineering of a simplified workflow execution for the MG5aMC event generator on GPUs and vector CPUs
• Fast calculation of HELAS amplitudes using graphics processing unit (GPU)
• Calculation of HELAS amplitudes for QCD processes using graphics processing unit (GPU)
• Accelerating LHC event generation with simplified pilot runs and fast PDFs
• Developments in Performance and Portability for MadGraph5_aMC@NLO
• Thread-Scalable Evaluation of Multi-Jet Observables
• Numerical evaluation of two-loop QCD helicity amplitudes for $gg\to t\bar{t}g$ at leading colour
• Computational challenges for MC event generation
• MadFlow: automating Monte Carlo simulation on GPU for particle physics processes