Two-loop Amplitudes Research Articles

Caravel: A C++ framework for the computation of multi-loop amplitudes with numerical unitarity

We present the first public version of Caravel, a C++17 framework for the computation of multi-loop scattering amplitudes in quantum field theory, based on the numerical unitarity method. Caravel is composed of modules for the D-dimensional decomposition of integrands of scattering amplitudes into master and surface terms, the computation of tree-level amplitudes in floating point or finite-field arithmetic, the numerical computation of one- and two-loop amplitudes in QCD and Einstein gravity, and functional reconstruction tools. We provide programs that showcase Caravel's main functionalities and allow to compute selected one- and two-loop amplitudes. Program summaryProgram Title:CaravelCPC Library link to program files:https://doi.org/10.17632/rfjrxrb3rk.1Developer's repository link:https://gitlab.com/caravel-public/caravel.gitLicensing provisions: GPLv3Programming language: C++External dependencies:• Required:Python3 [1], meson [2]• Optional:Doxygen [3], Eigen [4], GiNaC [5], GMP [6], Lapack [7], MPFR [8], MPI [9], PentagonLibrary [10,11], QD [12]Nature of problem: The computation of multi-loop multi-particle scattering amplitudes in quantum field theorySolution method: The multi-loop numerical unitarity method, functional reconstruction algorithmsAdditional comments including restrictions and unusual features: Current version includes tools employed in previous calculations, with the aim of showcasing details of the algorithms employed. Computations are organized by provided data files.

Two-loop leading-color helicity amplitudes for three-photon production at the LHC

We calculate all planar contributions to the two-loop massless helicity amplitudes for the process qoverline{q} → γγγ. The results are presented in fully analytic form in terms of the functional basis proposed recently by Chicherin and Sotnikov. With this publication we provide the two-loop contributions already used by us in the NNLO QCD calculation of the LHC process pp → γγγ [Chawdhry et al. (2019)]. Our results agree with a recent calculation of the same amplitude [Abreu et al. (2020)] which was performed using different techniques. We combine several modern computational techniques, notably, analytic solutions for the IBP identities, finite-field reconstruction techniques as well as the recent approach [Chen (2019)] for efficiently projecting helicity amplitudes. Our framework appears well-suited for the calculation of two-loop multileg amplitudes for which complete sets of master integrals exist.

Top quark contribution to two-loop helicity amplitudes for Z boson pair production in gluon fusion

We compute the top quark contribution to the two-loop amplitude for on-shell Z boson pair production in gluon fusion, gg → ZZ. Exact dependence on the top quark mass is retained. For each phase space point the integral reduction is performed numerically and the master integrals are evaluated using the auxiliary mass flow method, allowing fast computation of the amplitude with very high precision.

Two-loop helicity amplitudes for gg \u2192 ZZ with full top-quark mass effects

We calculate the two-loop QCD corrections to gg → ZZ involving a closed top-quark loop. We present a new method to systematically construct linear combinations of Feynman integrals with a convergent parametric representation, where we also allow for irreducible numerators, higher powers of propagators, dimensionally shifted integrals, and subsector integrals. The amplitude is expressed in terms of such finite integrals by employing syzygies derived with linear algebra and finite field techniques. Evaluating the amplitude using numerical integration, we find agreement with previous expansions in asymptotic limits and provide ab initio results also for intermediate partonic energies and non-central scattering at higher energies.

Mixed EW-QCD two-loop amplitudes for qoverline{q}to {mathrm{ell}}^{+}{mathrm{ell}}^{-} and \u03b35 scheme independence of multi-loop corrections

We perform a dedicated study of the qoverline{q} -initiated two-loop electroweak-QCD Drell-Yan scattering amplitude in dimensional regularization schemes for vanishing light quark and lepton masses. For the relative order α and αs one-loop Standard Model corrections, details of our comparison to the original literature are given. The infrared pole terms of the mixed two-loop amplitude are governed by a known generalization of the dipole formula and we show explicitly that exactly the same two-loop polarized hard scattering functions are obtained in both the standard ’t Hooft-Veltman-Breitenlohner-Maison γ5 scheme and Kreimer’s anticommuting γ5 scheme.

Locally finite two-loop amplitudes for off-shell multi-photon production in electron-positron annihilation

We study the singularity structure of two-loop QED amplitudes for the production of multiple off-shell photons in massless electron-positron annihilation and develop counterterms that remove their infrared and ultraviolet divergences point by point in the loop integrand. The remainders of the subtraction are integrable in four dimensions and can be computed in the future with numerical integration. The counterterms capture the divergences of the amplitudes and factorize in terms of the Born amplitude and the finite remainder of the one-loop amplitude. They consist of simple one- and two-loop integrals with at most three external momenta and can be integrated analytically in a simple manner with established methods. We uncover novel aspects of fully local IR factorization, where vertex and self energy subdiagrams must be modified by new symmetrizations over loop momenta, in order to expose their tree-like tensor structures and hence factorization of IR singularities prior to loop integration. This work is a first step towards isolating locally the hard contributions of generic gauge theory amplitudes and rendering them integrable in exactly four dimensions with numerical methods.

Non-planar form factors of generic local operators via on-shell unitarity and color-kinematics duality

Form factors, as quantities involving both local operators and asymptotic particle states, contain information of both the spectrum of operators and the on-shell amplitudes. So far the studies of form factors have been mostly focused on the large Nc planar limit, with a few exceptions of Sudakov form factors. In this paper, we discuss the systematical construction of full color dependent form factors with generic local operators. We study the color decomposition for form factors and discuss the general strategy of using on-shell unitarity cut method. As concrete applications, we compute the full two-loop non-planar minimal form factors for both half-BPS operators and non-BPS operators in the SU(2) sector in mathcal{N} = 4 SYM. Another important aspect is to investigate the color-kinematics (CK) duality for form factors of high-length operators. Explicit CK dual representation is found for the two-loop half-BPS minimal form factors with arbitrary number of external legs. The full-color two-loop form factor result provides an independent check of the infrared dipole formula for two-loop n-point amplitudes. By extracting the UV divergences, we also reproduce the known non-planar SU(2) dilatation operator at two loops. As for the finite remainder function, interestingly, the non-planar part is found to contain a new maximally transcendental part beyond the known planar result.

Two-loop leading colour QCD corrections to qoverline{q} \u2192 \u03b3\u03b3g and qg \u2192 \u03b3\u03b3q

We present the leading colour and light fermionic planar two-loop corrections for the production of two photons and a jet in the quark-antiquark and quark-gluon channels. In particular, we compute the interference of the two-loop amplitudes with the corresponding tree level ones, summed over colours and polarisations. Our calculation uses the latest advancements in the algorithms for integration-by-parts reduction and multivariate partial fraction decomposition to produce compact and easy-to-use results. We have implemented our results in an efficient C++ numerical code. We also provide their analytic expressions in Mathematica format.

Two-Loop Analysis of the Pion Mass Dependence of the ρ Meson.

Analyzing the pion mass dependence of ππ scattering phase shifts beyond the low-energy region requires the unitarization of the amplitudes from chiral perturbation theory. In the two-flavor theory, unitarization via the inverse-amplitude method (IAM) can be justified from dispersion relations, which is therefore expected to provide reliable predictions for the pion mass dependence of results from lattice QCD calculations. In this work, we provide compact analytic expression for the two-loop partial-wave amplitudes for J=0, 1, 2 required for the IAM at subleading order. To analyze the pion mass dependence of recent lattice QCD results for the P wave, we develop a fit strategy that for the first time allows us to perform stable two-loop IAM fits and assess the chiral convergence of the IAM approach. While the comparison of subsequent orders suggests a breakdown scale not much below the ρ mass, a detailed understanding of the systematic uncertainties of lattice QCD data is critical to obtain acceptable fits, especially at larger pion masses.

ZH production in gluon fusion: two-loop amplitudes with full top quark mass dependence

We present results for the two-loop helicity amplitudes entering the NLO QCD corrections to the production of a Higgs boson in association with a Z -boson in gluon fusion. The two-loop integrals, involving massive top quarks, are calculated numerically. Results for the interference of the finite part of the two-loop amplitudes with the Born amplitude are shown as a function of the two kinematic invariants on which the amplitudes depend.

The symbol and alphabet of two-loop NMHV amplitudes from overline{Q} equations

We study the symbol and the alphabet for two-loop NMHV amplitudes in planar mathcal{N} = 4 super-Yang-Mills from the overline{Q} equations, which provide a first-principle method for computing multi-loop amplitudes. Starting from one-loop N2MHV ratio functions, we explain in detail how to use overline{Q} equations to obtain the total differential of two-loop n-point NMHV amplitudes, whose symbol contains letters that are algebraic functions of kinematics for n ≥ 8. We present explicit formula with nice patterns for the part of the symbol involving algebraic letters for all multiplicities, and we find 17 − 2m multiplicative-independent letters for a given square root of Gram determinant, with 0 ≤ m ≤ 4 depending on the number of particles involved in the square root. We also observe that these algebraic letters can be found as poles of one-loop four-mass leading singularities with MHV or NMHV trees. As a byproduct of our algebraic results, we find a large class of components of two-loop NMHV, which can be written as differences of two double-pentagon integrals, particularly simple and free of square roots. As an example, we present the complete symbol for n = 9 whose alphabet contains 59 × 9 rational letters, in addition to the 11 × 9 independent algebraic ones. We also give all-loop NMHV last-entry conditions for all multiplicities.

A note on letters of Yangian invariants

Motivated by reformulating Yangian invariants in planar mathcal{N} = 4 SYM directly as d log forms on momentum-twistor space, we propose a purely algebraic problem of determining the arguments of the d log’s, which we call “letters”, for any Yangian invariant. These are functions of momentum twistors Z ’s, given by the positive coordinates α’s of parametrizations of the matrix C(α), evaluated on the support of polynomial equations C(α) · Z = 0. We provide evidence that the letters of Yangian invariants are related to the cluster algebra of Grassmannian G(4, n), which is relevant for the symbol alphabet of n-point scattering amplitudes. For n = 6, 7, the collection of letters for all Yangian invariants contains the cluster mathcal{A} coordinates of G(4, n). We determine algebraic letters of Yangian invariant associated with any “four-mass” box, which for n = 8 reproduce the 18 multiplicative-independent, algebraic symbol letters discovered recently for two-loop amplitudes.

Leading-color two-loop QCD corrections for three-photon production at hadron colliders

We compute the two-loop helicity amplitudes for the production of three photons at hadron colliders in QCD at leading-color. Using the two-loop numerical unitarity method coupled with analytic reconstruction techniques, we obtain the decomposition of the two-loop amplitudes in terms of master integrals in analytic form. These expressions are valid to all orders in the dimensional regulator. We use them to compute the two-loop finite remainders, which are given in a form that can be efficiently evaluated across the whole physical phase space. We further package these results in a public code which assembles the helicity-summed squared two-loop remainders, whose numerical stability across phase-space is demonstrated. This is the first time that a five-point two-loop process is publicly available for immediate phenomenological applications.

Contribution of third generation quarks to two-loop helicity amplitudes for W boson pair production in gluon fusion

We compute the contribution of third generation quarks (t, b) to the two-loop amplitude for on-shell W boson pair production in gluon fusion gg → WW. We present plots for the amplitude across partonic phase space as well as reference values for two kinematic points. The master integrals are efficiently evaluated by numerically solving a system of ordinary differential equations.

Triphoton production at hadron colliders in NNLO QCD

We present next-to-next-to-leading-order (NNLO) QCD corrections to the production of three isolated photons in hadronic collisions at the fully differential level. We employ qT subtraction within Matrix and an efficient implementation of analytic two-loop amplitudes in the leading-colour approximation to achieve the first on-the-fly calculation for this process at NNLO accuracy. Numerical results are presented for proton–proton collisions at energies ranging from 7 TeV to 100 TeV. We find full agreement with the 8 TeV results of Ref. [1] and confirm that NNLO corrections are indispensable to describe ATLAS 8 TeV data. In addition, we demonstrate the significance of NNLO corrections for future precision studies of triphoton production at higher collision energies.

Multi-Regge limit of the two-loop five-point amplitudes in mathcal{N} = 4 super Yang-Mills and mathcal{N} = 8 supergravity

In previous work, the two-loop five-point amplitudes in mathcal{N} = 4 super Yang-Mills theory and mathcal{N} = 8 supergravity were computed at symbol level. In this paper, we compute the full functional form. The amplitudes are assembled and simplified using the analytic expressions of the two-loop pentagon integrals in the physical scattering region. We provide the explicit functional expressions, and a numerical reference point in the scattering region. We then calculate the multi-Regge limit of both amplitudes. The result is written in terms of an explicit transcendental function basis. For certain non-planar colour structures of the mathcal{N} = 4 super Yang-Mills amplitude, we perform an independent calculation based on the BFKL effective theory. We find perfect agreement. We comment on the analytic properties of the amplitudes.

Symbol alphabets from plabic graphs

Symbol alphabets of n-particle amplitudes in mathcal{N} = 4 super-Yang-Mills theory are known to contain certain cluster variables of Gr(4, n) as well as certain algebraic functions of cluster variables. In this paper we suggest an algorithm for computing these symbol alphabets from plabic graphs by solving matrix equations of the form C ∙ Z = 0 to associate functions on Gr(m, n) to parameterizations of certain cells of Gr(k, n) indexed by plabic graphs. For m = 4 and n = 8 we show that this association precisely reproduces the 18 algebraic symbol letters of the two-loop NMHV eight-particle amplitude from four plabic graphs.

Two-loop superstring five-point amplitudes. Part I. Construction via chiral splitting and pure spinors

The full two-loop amplitudes for five massless states in Type II and Heterotic superstrings are constructed in terms of convergent integrals over the genus-two moduli space of compact Riemann surfaces and integrals of Green functions and Abelian differentials on the surface. The construction combines elements from the BRST cohomology of the pure spinor formulation and from chiral splitting with the help of loop momenta and homology invariance. The α′ → 0 limit of the resulting superstring amplitude is shown to be in perfect agreement with the previously known amplitude computed in Type II supergravity. Investigations of the α′ expansion of the Type II amplitude and comparisons with predictions from S-duality are relegated to a first companion paper. A construction from first principles in the RNS formulation of the genus-two amplitude with five external NS states is relegated to a second companion paper.

Quark mass effects in two-loop Higgs amplitudes

We provide two two-loop amplitudes relevant for precision Higgs physics. The first is the two-loop amplitude for Higgs boson production through gluon fusion with exact dependence on the top quark mass up to squared order in the dimensional regulator ε. The second result we provide is the two-loop amplitude for the decay of a Higgs boson into a pair of massive bottom quarks through the Higgs-to-gluon coupling in the infinite top mass limit. Both amplitudes are computed by finding canonical bases of master integrals, which we evaluate explicitly in terms of harmonic polylogarithms. We obtain the bare, renormalized and IR-subtracted amplitude and provide the results in terms of building blocks suitable for changing renormalization schemes.

Polarized mathrm{q}overline{mathrm{q}} \u2192 Z+Higgs amplitudes at two loops in QCD: the interplay between vector and axial vector form factors and a pitfall in applying a non-anticommuting \u03b35

We consider QCD corrections to two loops for the polarized amplitudes of qoverline{q} → Z + Higgs boson. First we show how the polarized amplitudes of boverline{b} → Zh associated with a non-vanishing b-quark Yukawa coupling and a scalar or pseudoscalar Higgs boson h can be built up solely from vector form factors (FF) of properly grouped classes of diagrams, bypassing completely the need of explicitly manipulating γ5 in dimensional regularization (up to a few “anomalous”, i.e., triangle diagrams). We determine the contributions of the triangle diagrams in the heavy top limit. We present the analytic results of the vector FF and the triangle-diagram contributions to the axial vector FF, which are sufficient for deriving the two-loop QCD amplitudes for boverline{b} → Zh with a CP-even and CP-odd Higgs boson h. We derive the respective Ward identity for these amplitudes, which are subsequently verified to two-loop order in QCD using these FF. In addition, the FF of a class of corrections to qoverline{q} → ZH proportional to the top-Yukawa coupling are obtained analytically to two-loop order in QCD in the heavy-top limit using the Higgs-gluon effective Lagrangian where the top quark is integrated out. We address a pitfall that occurs when applying the non-anticommutating γ5 prescription to this class of contributions that has been overlooked so far in the literature. We attribute this issue to the fact that the absence of certain heavy-mass expanded diagrams in the infinite-mass limit of a scattering amplitude with an axial vector current depends on the particular γ5 prescription in use.