Abstract The Local Unitarity (LU) formalism provides a constructive, integrand-level realisation of the Kinoshita-Lee-Nauenberg (KLN) theorem, by combining loop and phase-space integrals appearing in scattering cross-sections in such a way that their final-state infrared singularities cancel before integration. Supplemented with localised ultraviolet renormalisation, it enables the direct Monte Carlo integration of cross sections at arbitrary perturbative order in four-dimensional spacetime. In this paper, we present its application to the next-to-next-to-leading order (NNLO) QCD total cross sections for heavy-quark pair production in direct photon fusion, involving the contribution from 138 distinct forward-scattering diagrams where external photons couple only to heavy quarks. By combining NNLO QCD with next-to-leading order (NLO) electroweak (EW) corrections and next-to-leading power (NLP) Coulomb resummation, we obtain state-of-the-art predictions for top-, bottom-, and charm-quark production in ultraperipheral hadron collisions and at e + e − colliders.

We provide an overview of the status of Monte-Carlo event generators for high-energy particle physics. Guided by the experimental needs and requirements, we highlight areas of active development, and opportunities for future improvements. Particular emphasis is given to physics models and algorithms that are employed across a variety of experiments. These common themes in event generator development lead to a more comprehensive understanding of physics at the highest energies and intensities, and allow models to be tested against a wealth of data that have been accumulated over the past decades. A cohesive approach to event generator development will allow these models to be further improved and systematic uncertainties to be reduced, directly contributing to future experimental success. Event generators are part of a much larger ecosystem of computational tools. They typically involve a number of unknown model parameters that must be tuned to experimental data, while maintaining the integrity of the underlying physics models. Making both these data, and the analyses with which they have been obtained accessible to future users is an essential aspect of open science and data preservation. It ensures the consistency of physics models across a variety of experiments.

The recent experimental observation of Light-by-Light (LbL) scattering at the Large Hadron Collider has revived interest in this fundamental process, and especially of the accurate prediction of its cross-section, which we present here for the first time at Next-to-Leading Order (NLO) in both QCD and QED. We compare two radically different computational approaches, both exact in the fermion mass dependence, thus offering a strong cross-check of our results. The first approach is a fully analytic method to calculate compact and well-organized two-loop helicity amplitudes. The second one is entirely numerical and leverages the Local Unitarity construction. Our two calculations agree with each other and conclude that including the exact fermion mass contribution typically increases the size of the NLO corrections. Moreover, we find that the exact result converges slowly to the massless limit of the high-energy regime, thus emphasizing the importance of including the full mass dependence at NLO. We also compare our results with the ATLAS measurement of LbL in ultra-peripheral lead-lead collisions, and find that the inclusion of exact NLO corrections reduces, but does not eliminate, the existing tension with theoretical predictions.

We present an update of the Universal FeynRules Output model format, commonly known as the UFO format, that is used by several automated matrix-element generators and high-energy physics software. We detail different features that have been proposed as extensions of the initial format during the last ten years, and collect them in the current second version of the model format that we coin the Universal Feynman Output format. Following the initial philosophy of the UFO, they consist of flexible and modular additions to address particle decays, custom propagators, form factors, the renormalisation group running of parameters and masses, and higher-order quantum corrections.

We compute the next-to-leading-order (NLO) QCD corrections to the Higgs pT distribution in Higgs production in association with a jet via gluon fusion at the LHC, with exact dependence on the mass of the quark circulating in the heavy-quark loops. The NLO corrections are presented including the top-quark mass, and for the first time, the bottom-quark mass as well. Further, besides the on-shell mass scheme, we consider for the first time a running mass renormalisation scheme. The computation is based on amplitudes which are valid for arbitrary heavy-quark masses.

A bstract The Local Unitarity (LU) representation of differential cross-sections locally realises the cancellations of infrared singularities predicted by the Kinoshita-Lee-Nauenberg theorem. In this work we solve the two remaining challenges to enable practical higher-loop computations within the LU formalism. The first concerns the generalisation of the LU representation to graphs with raised propagators. The solution to this problem results in a generalisation of distributional Cutkosky rules. The second concerns the regularisation of ultraviolet and spurious soft singularities, solved using a fully automated and local renormalisation procedure based on Bogoliubov’s R -operation. We detail an all-order construction for the hybrid $$ \overline{\textrm{MS}} $$ MS ¯ and On-Shell scheme whose only analytic input is single-scale vacuum diagrams. We validate this novel technology by providing (semi-)inclusive results for two multi-leg processes at NLO, study limits of individual supergraphs up to N3LO and present the first physical NNLO cross-sections computed fully numerically in momentum-space, namely for the processes γ ∗ → jj and $$ {\gamma}^{\ast}\to t\overline{t} $$ γ ∗ → t t ¯ .

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}

We present the key features relevant to the automated computation of all the leading- and next-to-leading order contributions to short-distance cross sections in a mixed-coupling expansion, with special emphasis on the first subleading NLO term in the QCD+EW scenario, commonly referred to as NLO EW corrections. We discuss, in particular, the FKS subtraction in the context of a mixed-coupling expansion; the extension of the FKS subtraction to processes that include final-state tagged particles, defined by means of fragmentation functions; and some properties of the complex mass scheme. We combine the present paper with the release of a new version of MadGraph5_aMC@NLO, capable of dealing with mixed-coupling expansions. We use the code to obtain illustrative inclusive and differential results for the 13-TeV LHC.

We propose a novel representation of differential scattering cross-sections that locally realises the direct cancellation of infrared singularities exhibited by its so-called real-emission and virtual degrees of freedom. We take advantage of the Loop-Tree Duality representation of each individual forward-scattering diagram and we prove that the ensuing expression is locally free of infrared divergences, applies at any perturbative order and for any process without initial-state collinear singularities. Divergences for loop momenta with large magnitudes are regulated using local ultraviolet counterterms that reproduce the usual Lagrangian renormalisation procedure of quantum field theories. Our representation is especially suited for a numerical implementation and we demonstrate its practical potential by computing fully numerically and without any IR counterterm the next-to-leading order accurate differential cross-section for the process e+e− → $$ d\overline{d} $$ . We also show first results beyond next-to-leading order by computing interference terms part of the N4LO-accurate inclusive cross-section of a 1 → 2 + X scalar scattering process.

We present analytic results for the partonic cross sections for the production of a Higgs boson via the fusion of two bottom quarks at N$^3$LO in QCD perturbation theory in the five-flavour scheme. We combine this perturbative result with NLO accurate predictions in the four-flavour scheme that include the full bottom quark mass dependence by appropriately removing any double-counting stemming from contributions included in both predictions. We thereby obtain state-of-the-art predictions for the inclusive production probability of a Higgs boson via bottom quark fusion at hadron colliders.

We introduce a novel construction of a contour deformation within the framework of Loop-Tree Duality for the numerical computation of loop integrals featuring threshold singularities in momentum space. The functional form of our contour deformation automatically satisfies all constraints without the need for fine-tuning. We demonstrate that our construction is systematic and efficient by applying it to more than 100 examples of finite scalar integrals featuring up to six loops. We also showcase a first step towards handling non-integrable singularities by applying our work to one-loop infrared divergent scalar integrals and to the one-loop amplitude for the ordered production of two and three photons. This requires the combination of our contour deformation with local counterterms that regulate soft, collinear and ultraviolet divergences. This work is an important step towards computing higher-order corrections to relevant scattering cross-sections in a fully numerical fashion.

Loop-Tree Duality (LTD) is a framework in which the energy components of all loop momenta of a Feynman integral are integrated out using residue theorem, resulting in a sum over tree-like structures. Originally, the LTD expression exhibits cancellations of non-causal thresholds between summands, also known as dual cancellations. As a result, the expression exhibits numerical instabilities in the vicinity of non-causal thresholds and for large loop momenta. In this work we derive a novel, generically applicable, Manifestly Causal LTD (cLTD) representation whose only thresholds are causal thresholds, i.e. it manifestly realizes dual cancellations. Consequently, this result also serves as a general proof for dual cancellations. We show that LTD, cLTD, and the expression stemming from Time Ordered Perturbation Theory (TOPT) are locally equivalent. TOPT and cLTD both feature only causal threshold singularities, however LTD features better scaling with the number of propagators. On top of the new theoretical perspectives offered by our representation, it has the useful property that the ultraviolet (UV) behaviour of the original 4D integrand is maintained for every summand. We show that the resulting LTD integrand expression is completely stable in the UV region which is key for practical applications of LTD to the computation of amplitudes and cross sections. We present explicit examples of the LTD expression for a variety of up to four-loop integrals and show that its increased computational complexity can be efficiently mitigated by optimising its numerical implementation. Finally, we provide computer code that automatically generates the LTD expression for an arbitrary topology.

We propose a procedure to cross-validate Monte Carlo implementations of the standard model effective field theory. It is based on the numerical comparison of squared amplitudes computed at specific phase-space and parameter points in pairs of implementations. Interactions are fully linearised in the effective field theory expansion. The squares of linear effective field theory amplitudes and their interference with standard-model contributions are compared separately. Such pairwise comparisons are primarily performed at tree level and a possible extension to the one-loop level is also briefly considered. We list the current standard model effective field theory implementations and the comparisons performed to date.

The MadGraph5_aMC@NLO framework aims to automate all types of leading- and next-to-leading-order-accurate simulations for any user-defined model that stems from a renormalisable Lagrangian. In this paper, we present all of the key ingredients of such models in the context of supersymmetric theories. In order to do so, we extend the FeynRules package by giving it the possibility of dealing with different renormalisation options that are relevant to supersymmetric models. We also show how to deal with the problem posed by the presence of narrow resonances, thus generalising the so-called on-shell subtraction approaches. We extensively compare our total rate results with those of both Prospino2 and Resummino, and present illustrative applications relevant to the 13 TeV LHC, both at the total-rate and differential levels. The computer programmes that we have used to obtain the predictions presented here are all publicly available.

We compute mixed QCD-weak corrections to inclusive Higgs production at the LHC from the partonic process $$ gg\to Hq\overline{q} $$ . We start from the UV- and IR-finite oneloop weak amplitude and consider its interference with the corresponding one-loop QCD amplitude. This contribution is a $$ \mathcal{O}\left({\alpha}_s\alpha \right) $$ correction to the leading-order gluon-fusion cross section, and was not numerically assessed in previous works. We also compute the cross section from the square of this weak amplitude, suppressed by $$ \mathcal{O}\left({\alpha}^2\right) $$ . Finally, we consider contributions from the partonic process gq → Hq, which are one order lower in αs, as a reference for the size of terms which are not enhanced by the large gluon luminosity. We find that, given the magnitude of the uncertainties on current state-of-the-art predictions for Higgs production, all contributions computed in this work can be safely ignored, both fully inclusively and in the boosted Higgs regime. This result supports the approximate factorisation of QCD and weak corrections to that process.

Loop-tree duality (LTD) offers a promising avenue to numerically integrate multiloop integrals directly in momentum space. It is well established at one loop, but there have been only sparse numerical results at two loops. We provide a formal derivation for a novel multiloop LTD expression and study its threshold singularity structure. We apply our findings numerically to a diverse set of up to four-loop finite topologies with kinematics for which no contour deformation is needed. We also lay down the ground work for constructing such a deformation. Our results serve as an important stepping stone towards a generalized and efficient numerical implementation of LTD, which is applicable to the computation of virtual corrections.

This note summarizes the state-of-the-art predictions for the cross sections expected for Higgs boson production in the 27 TeV proton-proton collisions of a high-energy LHC, including a full theoretical uncertainty analysis. It also provides projections for the progress that may be expected on the timescale of the high-luminosity LHC and an assessment of the main limiting factors to further reduction of the remaining theoretical uncertainties. https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCHXSWG

This report comprises the outcome of five working groups that have studied the physics potential of the high-luminosity phase of the LHC (HL-LHC) and the perspectives for a possible future high-energy LHC (HE-LHC).The working groups covered a broad range of topics: Standard Model measurements, studies of the properties ofthe Higgs boson, searches for phenomena beyond the Standard Model, flavor physics of heavy quarks and leptonsand studies of QCD matter at high density and temperature.The work is prepared as an input to the ongoing process of updating the European Strategy for Particle Physics,a process that will be concluded in May 2020.

The discovery of the Higgs boson in 2012, by the ATLAS and CMS experiments, was a success achieved with only a percent of the entire dataset foreseen for the LHC. It opened a landscape of possibilities in the study of Higgs boson properties, Electroweak Symmetry breaking and the Standard Model in general, as well as new avenues in probing new physics beyond the Standard Model. Six years after the discovery, with a conspicuously larger dataset collected during LHC Run 2 at a 13 TeV centre-of-mass energy, the theory and experimental particle physics communities have started a meticulous exploration of the potential for precision measurements of its properties. This includes studies of Higgs boson production and decays processes, the search for rare decays and production modes, high energy observables, and searches for an extended electroweak symmetry breaking sector. This report summarises the potential reach and opportunities in Higgs physics during the High Luminosity phase of the LHC, with an expected dataset of pp collisions at 14 TeV, corresponding to an integrated luminosity of 3 ab$^{-1}$. These studies are performed in light of the most recent analyses from LHC collaborations and the latest theoretical developments. The potential of an LHC upgrade, colliding protons at a centre-of-mass energy of 27 TeV and producing a dataset corresponding to an integrated luminosity of 15 ab$^{-1}$, is also discussed.