Search

We present a computation of the one-loop QCD corrections to top-quark pair production in association with a $W$ boson, including terms up to order $\varepsilon^2$ in dimensional regularization. Providing a first glimpse into the complexity of the corresponding two-loop amplitude, this result is a first step towards a description of this process at next-to-next-to-leading order (NNLO) in QCD. We perform a tensor decomposition and express the corresponding form factors in terms of a basis of independent special functions with compact rational coefficients, providing a structured framework for future developments. In addition, we derive an explicit analytic representation of the form factors, valid up to order $\varepsilon^0$, expressed in terms of logarithms and dilogarithms. For the complete set of special functions required, we obtain a semi-numerical solution based on generalized power series expansion., Comment: 37 pages, 1 figure, 2 tables. The paper is complemented by a Zenodo repository found at https://zenodo.org/records/14810498

We address the analytic computation of the two-loop scattering amplitudes for the production of two photons in parton-parton scattering, mediated by loops of heavy quarks. Due to the presence of integrals of elliptic type, both partonic channels have been previously computed using semi-numerical methods. In this paper, leveraging new advances in the theory of differential equations for elliptic Feynman integrals, we derive a canonical basis for all integrals involved and compute them in terms of independent iterated integrals over elliptic and polylogarithmic differential forms. We use this representation to showcase interesting cancellations in the physical expressions for the scattering amplitudes. Furthermore, we address their numerical evaluation by producing series expansion representations for the whole amplitudes, which we demonstrate to be fast and numerically reliable across a large region of the phase space., Comment: 48 pages, 6 figures, accompanied by a Zenodo respository submission found at https://doi.org/10.5281/zenodo.14733100

We present the first benchmark evaluation of the two-loop finite remainders for the production of a top-quark pair in association with a jet at hadron colliders in the gluon channel. We work in the leading colour approximation, and perform the numerical evaluation in the physical phase space. To achieve this result, we develop a new method for expressing the master integrals in terms of a (over-complete) basis of special functions that enables the infrared and ultraviolet poles to be cancelled analytically despite the presence of elliptic Feynman integrals. The special function basis makes it manifest that the elliptic functions appear solely in the finite remainder, and can be evaluated numerically through generalised series expansions. The helicity amplitudes are constructed using four dimensional projectors combined with finite-field techniques to perform integration-by-parts reduction, mapping to special functions and Laurent expansion in the dimensional regularisation parameter., Comment: 22 pages, 6 figures, 3 tables, ancillary files at https://zenodo.org/records/14473301

We compute the differential equations for the two remaining integral topologies contributing to the leading colour two-loop amplitudes for $pp \rightarrow t\bar{t}j$. We derive differential equations for the master integrals by solving the integration-by-parts identities over finite fields. Of the two systems of differential equations, one is presented in canonical '${\rm d} \log$' form, while the other is found to have an elliptic sector. For the elliptic topology we identify the relevant elliptic curve, and present the differential equations in a more general form which depends quadratically on $\epsilon$ and contains non-logarithmic one-forms in addition to the canonical ${\rm d} \log$'s. We solve the systems of differential equations numerically using generalised series expansions with the boundary terms obtained using the auxiliary mass flow method. A summary of all one-loop and two-loop planar topologies is presented including the list of alphabet letters for the '${\rm d} \log$' form systems and high-precision boundary values., Comment: 22 pages, 7 figures, 3 tables, 3 appendices

We report on recent progress for the QCD corrections to top quark pair plus jet production. In particular, we discuss a recent computation for the two-loop master integrals associated to a two-loop five-point pentagon-box integral configuration with one internal massive propagator, that contributes to top quark pair production in association with a jet in the QCD planar limit., Comment: Proceedings for Radcor 2023, 8 pages, 1 figure

We present the computation of the two-loop form factors for diphoton production in the quark annihilation channel. These quantities are relevant for the NNLO QCD corrections to diphoton production at LHC recently presented in arXiv:2308.10885. The computation is performed retaining full dependence on the mass of the heavy quark in the loops. The master integrals are evaluated by means of differential equations which are solved exploiting the generalised power series technique., Comment: 23 pages, 5 figures

In this paper we consider the diphoton production in hadronic collisions at the next-to-next-to-leading order (NNLO) in perturbative QCD, taking into account for the first time the full top quark mass dependence up to two loops (full NNLO). We show selected numerical distributions, highlighting the kinematic regions where the massive corrections are more significant. We make use of the recently computed two-loop massive amplitudes for diphoton production in the quark annihilation channel. The remaining massive contributions at NNLO are also considered, and we comment on the weight of the different types of contributions to the full and complete result., Comment: 14 pages and 5 figures

We present the analytic expressions for the two-loop form factors for the production or decay of pseudo-scalar quarkonia, in a scheme where the quarks are produced at threshold. We consider the two-loop amplitude for the process $\gamma \gamma \leftrightarrow {^1S_0^{[1]}}$, that was previously known only numerically, as well as for the processes $gg \leftrightarrow {^1S_0^{[1]}}$, $\gamma g \leftrightarrow {^1S_0^{[8]}}$ and $gg \leftrightarrow {^1S_0^{[8]}}$, which have not been computed before. The two-loop corrections to $gg \leftrightarrow {^1S_0^{[1]}}$ are the last missing ingredients for a full NNLO calculation of $\eta_Q$ hadro-production. We discuss how the singularity structure of the amplitudes is affected by the threshold kinematics, which in particular introduces Coulomb singularities. In this context, we first show how the usual structure of the infrared singularities degenerates at threshold kinematics, and then extract the anomalous dimensions governing the Coulomb singularities for colour-singlet and octet channels, the latter being presented here for the first time. We give high-precision numerical results for the hard functions, which can be used for phenomenological studies of $\eta_Q$ production and decay at NNLO., Comment: 37 pages, 1 figure, v2: version published in JHEP

We consider the case of a two-loop five-point pentagon-box integral configuration with one internal massive propagator that contributes to top-quark pair production in association with a jet at hadron colliders. We construct the system of differential equations for all the master integrals in a canonical form where the analytic form is reconstructed from numerical evaluations over finite fields. We find that the system can be represented as a sum of d-logarithmic forms using an alphabet of 71 letters. Using high precision boundary values obtained via the auxiliary mass flow method, a numerical solution to the master integrals is provided using generalised power series expansions., Comment: 31 pages, 47 figures, ancillary material attached to the submission. Version v3 as published in JHEP

We compute the master integrals relevant for the two-loop corrections to pseudo-scalar quarkonium and leptonium production and decay. We present both analytic and high-precision numerical results. The analytic expressions are given in terms of multiple polylogarithms (MPLs), elliptic multiple polylogarithms (eMPLs) and iterated integrals of Eisenstein series. As an application of our results, we obtain for the first time an analytic expression for the two-loop amplitude for para-positronium decay to two photons at two loops., Comment: 34 pages, 6 figures, v2: references updated, version published in JHEP

We compute helicity amplitudes for the one-loop QCD corrections to top-quark pair production analytically in terms of a set of uniformly transcendental master integrals. We provide corrections up to $O(\epsilon^2)$ in the dimensional regulator for the first time which are relevant at NNLO. Four independent pentagon integral topologies appear in the complete description of the colour structure for which we provide numerical solutions using canonical form differential equations and the method of generalised power series expansions. Analytic forms of the boundary values are obtained in all cases except one where we find a one-dimensional integral representation., Comment: 37 pages. ancillary files containing analytic amplitudes and integral evaluation routines. v2 references added. v3 fixing minor typos in the text and in the ancillary files

We report on the two-loop amplitude computation for the mixed QCD-electroweak corrections to the process $gg \to Hg$, with exact dependence on the electroweak boson masses. This amplitude has been employed in the computation of next-to-leading order (NLO) mixed QCD-electroweak corrections to the Higgs-boson production rate in arXiv:2010.09451. The master integrals that appear in the amplitude are evaluated by means of generalized power series expansions, which allows for fast and high-precision numerical evaluation of the amplitude in the physical phase-space, proving to be a powerful tool for phenomenological applications., Comment: 23 pages, 2 figures

We review recent progress on operator mixing in the light of the theory of canonical forms for linear systems of differential equations and, in particular, of the Poincar\'e-Dulac theorem. We show that the matrix $A(g) = -\frac{\gamma(g)}{\beta(g)} =\frac{\gamma_0}{\beta_0}\frac{1}{g} + \cdots $ determines which different cases of operator mixing can occur, and we review their classification. We derive a sufficient condition for $A(g)$ to be set in the one-loop exact form $A(g) = \frac{\gamma_0}{\beta_0}\frac{1}{g}$. Finally, we discuss the consequences of the unitarity requirement in massless QCD-like theories, and we demonstrate that $\gamma_0$ is always diagonalizable if the theory is conformal invariant and unitary in its free limit at $g =0$., Comment: 8 pages, fixed typo in the references. Contribution to the 15th International Symposium on Radiative Corrections (RADCOR) and the XIX Workshop on Radiative Corrections for the LHC and Future Colliders (LoopFest). arXiv admin note: substantial text overlap with arXiv:2103.16220

We work out the interplay between lowest-order perturbative computations in the 't Hooft coupling, $g^2=g^2_{YM} N$, operator mixing, renormalization-group (RG) improved ultraviolet (UV) asymptotics of leading-order (LO) nonplanar/planar contributions to $2$-point correlators, and nonperturbative large-$N$ expansion of perturbatively massless QCD-like theories. As concrete examples, we compute to the lowest perturbative order in $SU(N)$ YM theory the ratios, $r_i$, of LO-nonplanar to planar contributions to the $2$-point correlators in the orthogonal basis in the coordinate representation of the gauge-invariant dimension-$8$ scalar operators and all the twist-$2$ operators. We demonstrate that -- if $\frac{\gamma_0}{\beta_0}$ has no LO-nonplanar contribution, with $\gamma_0$ and $\beta_0$ the one-loop coefficients of the anomalous-dimension matrix and beta function respectively -- $r_i$ actually coincides with the corresponding ratio in the large-$N$ expansion of the RG-improved UV asymptotics of the $2$-point correlators, provided that a certain canonical nonresonant diagonal renormalization scheme exists for the corresponding operators. Contrary to the aforementioned scalar operators, for the first $10^3$ twist-$2$ operators we actually verify the above conditions, and we get the universal value $r_i=-\frac{1}{N^2}$. Hence, nonperturbatively such $r_i$ must coincide with the UV asymptotics of the ratio of the glueball self-energy loop to the glueball tree contribution to the $2$-point correlators above. As a consequence, the universality of $r_i$ reflects the universality of the effective coupling in the nonperturbative large-$N$ YM theory for the twist-$2$ operators in the coordinate representation., Comment: 24 pages, no figures

Recently, a geometric approach to operator mixing in massless QCD-like theories -- that involves canonical forms based on the Poincare'-Dulac theorem for the linear system that defines the renormalized mixing matrix in the coordinate representation $Z(x,\mu)$ -- has been advocated in arXiv:2103.15527 . As a consequence, a classification of operator mixing in four cases -- depending on the canonical forms of $- \frac{\gamma(g)}{\beta(g)}$, with $\gamma(g)=\gamma_0 g^2+\cdots$ the matrix of the anomalous dimensions and $\beta(g)=-\beta_0 g^3 + \cdots$ the beta function -- has been proposed: (I) nonresonant $\frac{\gamma_0}{\beta_0}$ diagonalizable, (II) resonant $\frac{\gamma_0}{\beta_0}$ diagonalizable, (III) nonresonant $\frac{\gamma_0}{\beta_0}$ nondiagonalizable, (IV) resonant $\frac{\gamma_0}{\beta_0}$ nondiagonalizable. In particular, in arXiv:2103.15527 a detailed analysis of the case (I) -- where operator mixing reduces to all orders of perturbation theory to the multiplicatively renormalizable case -- has been provided. In the present paper, following the aforementioned approach, we work out in the remaining three cases the canonical forms for $- \frac{\gamma(g)}{\beta(g)}$ to all orders of perturbation theory, the corresponding UV asymptotics of $Z(x,\mu)$, and the physics interpretation. We also work out in detail physical realizations of the cases (I) and (II)., Comment: 35 pages, formulas unchanged, but some comments on the UV asymptotics corrected, physical realizations of the resonant case and new references added

We present for the first time the exact NLO QCD corrections to the light-quark part of the mixed QCD-EW contributions to Higgs production via gluon fusion at LHC13, with exact EW-boson mass dependence. The relevant two-loop real-emission matrix element is computed using a dynamic one-dimensional series expansion strategy whose stability and speed allows for a numerical phase-space integration using local IR subtraction counterterms. For $\mu_R=\mu_F=M_H$, we find: \begin{equation}\sigma^{(\alpha_s^2\alpha^2+\alpha_s^3\alpha^2)}_{g g\rightarrow H+X} = 1.467(2)^{\;+18.7\%}_{\;-14.6\%}\;(\mu_R\;\text{var.})\;\pm 2\%\;(\text{PDF}) \ \textrm{pb},\end{equation} which we use to provide the best result including an estimate of suppressed contributions: \begin{equation}\sigma^{(\text{EW},\textrm{best})}_{p p\rightarrow H+X} = 2.11 \pm 0.28 \ (\textrm{theory}) \ \mathrm{pb}.\end{equation}, Comment: 7 pages, 4 figures

We present the analytic computation of the master integrals associated to certain two-loop non-planar topologies, which are needed to complete the evaluation of the last two color coefficients for the top-pair production in the quark-annihilation channel, which are not yet known analytically. The master integrals have been computed exploiting the differential equations method in canonical form. The solution is given as a series expansion in the dimensional regularization parameter through to weight four, the expansion coefficients are given in terms of multiple polylogarithms., Comment: 11 pages, 2 figures, talk given at RADCOR2019, based on arXiv:1904.10834

We compute fully analytic results for the three-loop diagrams involving two different massive quark flavours contributing to the $\rho$ parameter in the Standard Model. We find that the results involve exactly the same class of functions that appears in the well-known sunrise and banana graphs, namely elliptic polylogarithms and iterated integrals of modular forms. Using recent developments in the understanding of these functions, we analytically continue all the iterated integrals of modular forms to all regions of the parameter space, and in each region we obtain manifestly real and fast-converging series expansions for these functions., Comment: 37 pages, 2 figures

We present the analytic calculation of the Master Integrals for the two-loop, non-planar topologies that enter the calculation of the amplitude for top-quark pair hadroproduction in the quark-annihilation channel. Using the method of differential equations, we expand the integrals in powers of the dimensional regulator $\epsilon$ and determine the expansion coefficients in terms of generalized harmonic polylogarithms of two dimensionless variables through to weight four., Comment: 28 pages, 2 figures, ancillary files included

We verify, both perturbatively and nonperturbatively asymptotically in the ultraviolet (UV), a special case of a low-energy theorem of the NSVZ type in QCD-like theories, recently derived in arXiv:1701.07833, that relates the logarithmic derivative with respect to the gauge coupling, or the logarithmic derivative with respect to the renormalization-group (RG) invariant scale, of an $n$-point correlator of local operators in one side to an $n+1$-point correlator with the insertion of $Tr F^2$ at zero momentum in the other side. Our computation involves the operator product expansion (OPE) of the scalar glueball operator, $Tr F^2$, in massless QCD, worked out perturbatively in arXiv:1209.1516 -- and in its RG-improved form in the present paper -- by means of which we extract both the perturbative divergences and the nonperturbative UV asymptotics in both sides. We also discuss the role of the contact terms in the OPE, both finite and divergent, discovered some years ago in arXiv:1209.1516, in relation to the low-energy theorem. Besides, working the other way around by assuming the low-energy theorem for any 2-point correlator of a multiplicatively renormalizable gauge-invariant operator, we compute in a massless QCD-like theory the corresponding perturbative OPE to the order of $g^2$ and nonperturbative asymptotics. The low-energy theorem has a number of applications: to the renormalization in asymptotically free QCD-like theories, both perturbatively and nonperturbatively in the large-$N$ 't Hooft and Veneziano expansions, and to the way the open/closed string duality may or may not be realized in the would-be solution by canonical string theories for QCD-like theories, both perturbatively and in the 't Hooft large-$N$ expansion. Our computations will also enter further developments based on the low-energy theorem., Comment: Some arguments extended and minor typos corrected, paper as published in JHEP

We present the analytic calculation of the planar master integrals which contribute to compute the two-loop light-fermion electroweak corrections to the production of a Higgs boson in association with a jet in gluon-gluon fusion. The complete dependence on the electroweak-boson mass is retained. The master integrals are evaluated by means of the differential equations method and the analytic results are expressed in terms of multiple polylogarithms up to weight four., Comment: 21 pages, ancillary files

We present the analytic calculation of the Master Integrals necessary to compute the planar massive QCD corrections to Di-photon (and Di-jet) production at hadron colliders. The masters are evaluated by means of the differential equations method and expressed in terms of multiple polylogarithms and one- or two-fold integrals of polylogarithms and irrational functions, up to transcendentality four., Comment: 20 pages, ancillary files