We complete the computation of the two-loop helicity amplitudes for the production of three photons at hadron colliders, including all contributions beyond the leading-color approximation. We reconstruct the analytic form of the amplitudes from numerical finite-field samples obtained with the numerical unitarity method. This method requires as input surface terms for all relevant five-point non-planar integral topologies, which we obtain by solving the associated syzygy problem in embedding space. The numerical samples are used to constrain compact spinor-helicity ans\"atze, which are optimized by taking advantage of the known one-loop analytic structure. We make our analytic results available in a public C++ library, which is suitable for immediate phenomenological applications. We estimate that the inclusion of the subleading-color contributions will decrease the size of the two-loop corrections by about 30% to 50%, and the NNLO cross sections by a few percent, compared to the results in the leading-color approximation.
@article{arxiv.2305.17056,
title = {Two-Loop QCD Corrections for Three-Photon Production at Hadron Colliders},
author = {Samuel Abreu and Giuseppe De Laurentis and Harald Ita and Maximillian Klinkert and Ben Page and Vasily Sotnikov},
journal= {arXiv preprint arXiv:2305.17056},
year = {2023}
}