We compute the two-loop massless QCD corrections to the four-point amplitude $g+g rightarrow H+H$ resulting from effective operator insertions that describe the interaction of a Higgs boson with gluons in the infinite top quark mass limit. This amplitude is an essential ingredient to the third-order QCD corrections to Higgs boson pair production. We have implemented our results in a numerical code that can be used for further phenomenological studies.
We report consistent results for $Gamma(h rightarrow gamma gamma)$, $sigma(mathcal{G} ,mathcal{G}rightarrow h)$ and $Gamma(h rightarrow mathcal{G} ,mathcal{G})$ in the Standard Model Effective Field Theory (SMEFT) perturbing the SM by corrections $ma
thcal{O}(bar{v}_T^2/16 pi^2 Lambda^2)$ in the Background Field Method (BFM) approach to gauge fixing, and to $mathcal{O}(bar{v}_T^4/Lambda^4)$ using the geometric formulation of the SMEFT. We combine and modify recent results in the literature into a complete set of consistent results, uniforming conventions, and simultaneously complete the one loop results for these processes in the BFM. We emphasise calculational scheme dependence present across these processes, and how the operator and loop expansions are not independent beyond leading order. We illustrate several cross checks of consistency in the results.
A large set of techniques needed to compute decay rates at the two-loop level are derived and systematized. The main emphasis of the paper is on the two Standard Model decays H -> gamma gamma and H -> g g. The techniques, however, have a much wider r
ange of application: they give practical examples of general rules for two-loop renormalization; they introduce simple recipes for handling internal unstable particles in two-loop processes; they illustrate simple procedures for the extraction of collinear logarithms from the amplitude. The latter is particularly relevant to show cancellations, e.g. cancellation of collinear divergencies. Furthermore, the paper deals with the proper treatment of non-enhanced two-loop QCD and electroweak contributions to different physical (pseudo-)observables, showing how they can be transformed in a way that allows for a stable numerical integration. Numerical results for the two-loop percentage corrections to H -> gamma gamma, g g are presented and discussed. When applied to the process pp -> gg + X -> H + X, the results show that the electroweak scaling factor for the cross section is between -4 % and + 6 % in the range 100 GeV < Mh < 500 GeV, without incongruent large effects around the physical electroweak thresholds, thereby showing that only a complete implementation of the computational scheme keeps two-loop corrections under control.
We present the analytic evaluation of the two-loop corrections to the amplitude for the scattering of four fermions in Quantum Electrodynamics, $f^- + f^+ + F^- + F^+ to 0$, with $f$ and $F$ representing a massless and a massive lepton, respectively.
Dimensional regularization is employed to evaluate the loop integrals. Ultraviolet divergences are removed by renormalizing the coupling constant in the ${overline{text{MS}}}$-scheme, and the lepton mass as well as the external fields in the on-shell scheme. The analytic result for the renormalized amplitude is expressed as Laurent series around $d=4$ space-time dimensions, and contains Generalized Polylogarithms with up to weight four. The structure of the residual infrared divergences of the virtual amplitude is in agreement with the prediction of the Soft Collinear Effective Theory. Our analytic results are an essential ingredient for the computation of the scattering cross section for massive fermion-pair production in massless fermion-pair annihilation, i.e. $f^- f^+ to F^- F^+$, and crossing related processes such as the elastic scattering $f F to f F$, with up to Next-to-Next to Leading Order accuracy.
Two-loop electroweak corrections to the muon anomalous magnetic moment are automatically calculated by using GRACE-FORM system, as a trial to extend our system for two-loop calculation. We adopt the non-linear gauge (NLG) to check the reliability of
our calculation. In total 1780 two-loop diagrams consisting of 14 different topological types and 70 one-loop diagrams composed of counter terms are calculated. We check UV- and IR-divergences cancellation and the independence of the results from NLG parameters. As for the numerical calculation, we adopt trapezoidal rule with Double Exponential method (DE). Linear extrapolation method (LE) is introduced to regularize UV- and IR- divergence and to get finite values.
We extend useful properties of the $Htogammagamma$ unintegrated dual amplitudes from one- to two-loop level, using the Loop-Tree Duality formalism. In particular, we show that the universality of the functional form -- regardless of the nature of the
internal particle -- still holds at this order. We also present an algorithmic way to renormalise two-loop amplitudes, by locally cancelling the ultraviolet singularities at integrand level, thus allowing a full four-dimensional numerical implementation of the method. Our results are compared with analytic expressions already available in the literature, finding a perfect numerical agreement. The success of this computation plays a crucial role for the development of a fully local four-dimensional framework to compute physical observables at Next-to-Next-to Leading order and beyond.