Source author record

Claude Duhr

Claude Duhr appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

64works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

64 published item(s)

preprint2026arXiv

NNLOCAL: Completely Local Subtractions

The computation of higher-order corrections to cross sections relevant at LHC involves the evaluation of phase-space integrals that exhibit soft and collinear divergences. The subtraction of these divergences is a key ingredient to obtain fully-differential predictions for physical observables. We discuss a subtraction method to handle these divergences based on the construction of universal local counterterms. The integration of the counterterms is carried out analytically, giving a strong control on the numerical stability of our predictions. We implement our method in a numerical program, that we dub NNLOCAL, and validate it by computing the fully-differential NNLO cross-section for Higgs boson production in gluon-gluon fusion.

preprint2023arXiv

The SAGEX Review on Scattering Amplitudes, Chapter 3: Mathematical structures in Feynman integrals

Dimensionally-regulated Feynman integrals are a cornerstone of all perturbative computations in quantum field theory. They are known to exhibit a rich mathematical structure, which has led to the development of powerful new techniques for their computation. We review some of the most recent advances in our understanding of the analytic structure of multiloop Feynman integrals in dimensional regularisation. In particular, we give an overview of modern approaches to computing Feynman integrals using differential equations, and we discuss some of the properties of the functions that appear in the solutions. We then review how dimensional regularisation has a natural mathematical interpretation in terms of the theory of twisted cohomology groups, and how many of the well-known ideas about Feynman integrals arise naturally in this context. This is Chapter 3 of a series of review articles on scattering amplitudes, of which Chapter 0 [arXiv:2203.13011] presents an overview and Chapter 4 [arXiv:2203.13015] contains closely related topics.

preprint2022arXiv

Functions Beyond Multiple Polylogarithms for Precision Collider Physics

Feynman diagrams constitute one of the essential ingredients for making precision predictions for collider experiments. Yet, while the simplest Feynman diagrams can be evaluated in terms of multiple polylogarithms -- whose properties as special functions are well understood -- more complex diagrams often involve integrals over complicated algebraic manifolds. Such diagrams already contribute at NNLO to the self-energy of the electron, $t \bar{t}$ production, $γγ$ production, and Higgs decay, and appear at two loops in the planar limit of maximally supersymmetric Yang-Mills theory. This makes the study of these more complicated types of integrals of phenomenological as well as conceptual importance. In this white paper contribution to the Snowmass community planning exercise, we provide an overview of the state of research on Feynman diagrams that involve special functions beyond multiple polylogarithms, and highlight a number of research directions that constitute essential avenues for future investigation.

preprint2022arXiv

The Diagrammatic Coaction

The diagrammatic coaction underpins the analytic structure of Feynman integrals, their cuts and the differential equations they admit. The coaction maps any diagram into a tensor product of its pinches and cuts. These correspond respectively to differential forms defining master integrals, and integration contours which place a subset of the propagators on shell. In a canonical basis these forms and contours are dual to each other. In this talk I review our present understanding of this algebraic structure and its manifestation for dimensionally-regularized Feynman integrals that are expandable to polylogarithms around integer dimensions. Using one- and two-loop integral examples, I will explain the duality between forms and contours, and the correspondence between the local coaction acting on the Laurent coefficients in the dimensional regulator and the global coaction acting on generalised hypergeometric functions.

preprint2022arXiv

The Four-Loop Rapidity Anomalous Dimension and Event Shapes to Fourth Logarithmic Order

We obtain the quark and gluon rapidity anomalous dimension to fourth order in QCD. We calculate the N$^3$LO rapidity anomalous dimensions to higher order in the dimensional regulator and make use of the soft/rapidity anomalous dimension correspondence in conjunction with the recent determination of the N$^4$LO threshold anomalous dimensions to achieve our result. We show that the results for the quark and gluon rapidity anomalous dimensions at four loops are related by generalized Casimir scaling. Using the N$^4$LO rapidity anomalous dimension, we perform the resummation of the Energy-Energy Correlation in the back-to-back limit at N$^4$LL, achieving for the first time the resummation of an event shape at this logarithmic order. We present numerical results and observe a reduction of perturbative uncertainties on the resummed cross section to below 1%.

preprint2022arXiv

The Path forward to N$^3$LO

The LHC experiments will achieve percent level precision measurements of processes key to some of the most pressing questions of contemporary particle physics: What is the nature of the Higgs boson? Can we successfully describe the interaction of fundamental particles at high energies? Is there physics beyond the Standard Model at the LHC? The capability to predict and describe such observables at next-to-next-to-next-to-leading order (N$^3$LO) in QCD perturbation theory is paramount to fully exploit these experimental measurements. We describe the current status of N$^3$LO predictions and highlight their importance in the upcoming precision phase of the LHC. Furthermore, we identify key conceptual and mathematical developments necessary to see wide-spread N$^3$LO phenomenology come to fruition.

preprint2021arXiv

Lepton-pair production at hadron colliders at N$^3$LO in QCD

We compute for the first time the complete corrections at N$^3$LO in the strong coupling constant to the inclusive neutral-current Drell-Yan process including contributions from both photon and $Z$-boson exchange. Our main result is the computation of the QCD corrections to the inclusive production cross section of an axial-vector boson to third order in the strong coupling in a variant of QCD with five massless quark flavours. Since the axial anomaly does not cancel for an odd number of flavours, we also consistently include non-decoupling effects in the top-quark mass through three loops. We perform a phenomenological study of our results, and we present for the first time predictions for the inclusive Drell-Yan process at the LHC at this order in QCD perturbation theory.

preprint2021arXiv

Meromorphic modular forms and the three-loop equal-mass banana integral

We consider a class of differential equations for multi-loop Feynman integrals which can be solved to all orders in dimensional regularisation in terms of iterated integrals of meromorphic modular forms. We show that the subgroup under which the modular forms transform can naturally be identified with the monodromy group of a certain second-order differential operator. We provide an explicit decomposition of the spaces of modular forms into a direct sum of total derivatives and a basis of modular forms that cannot be written as derivatives of other functions, thereby generalising a result by one of the authors form the full modular group to arbitrary finite-index subgroups of genus zero. Finally, we apply our results to the two- and three-loop equal-mass banana integrals, and we obtain in particular for the first time complete analytic results for the higher orders in dimensional regularisation for the three-loop case, which involves iterated integrals of meromorphic modular forms.

preprint2020arXiv

A double integral of dlog forms which is not polylogarithmic

Feynman integrals are central to all calculations in perturbative Quantum Field Theory. They often give rise to iterated integrals of dlog-forms with algebraic arguments, which in many cases can be evaluated in terms of multiple polylogarithms. This has led to certain folklore beliefs in the community stating that all such integrals evaluate to polylogarithms. Here we discuss a concrete example of a double iterated integral of two dlog-forms that evaluates to a period of a cusp form. The motivic versions of these integrals are shown to be algebraically independent from all multiple polylogarithms evaluated at algebraic arguments. From a mathematical perspective, we study a mixed elliptic Hodge structure arising from a simple geometric configuration in $\mathbb{P}^2$, consisting of a modular plane elliptic curve and a set of lines which meet it at torsion points, which may provide an interesting worked example from the point of view of periods, extensions of motives, and L-functions.

preprint2020arXiv

Higgs production in bottom quark fusion: Matching the 4- and 5-flavour schemes to third order in the strong coupling

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.

preprint2020arXiv

Tree-level splitting amplitudes for a quark into four collinear partons

We compute in conventional dimensional regularisation the tree-level splitting amplitudes for a quark parent in the limit where four partons become collinear to each other. This is part of the universal infrared behaviour of the QCD scattering amplitudes at next-to-next-to-next-to-leading order (${\rm N^3LO}$) in the strong coupling constant. Further, we consider the iterated limit when $m'$ massless partons become collinear to each other within a bigger set of $m$ collinear partons, as well as the limits when one gluon or a $q\bar{q}$ pair or two gluons become soft within a set of $m$ collinear partons.

preprint2019arXiv

Algorithms and tools for iterated Eisenstein integrals

We present algorithms to work with iterated Eisenstein integrals that have recently appeared in the computation of multi-loop Feynman integrals. These algorithms allow one to analytically continue these integrals to all regions of the parameter space, and to obtain fast converging series representations in each region. We illustrate our approach on the examples of hypergeometric functions that evaluate to iterated Eisenstein integrals as well as the well-known sunrise graph.

preprint2019arXiv

From positive geometries to a coaction on hypergeometric functions

It is well known that Feynman integrals in dimensional regularization often evaluate to functions of hypergeometric type. Inspired by a recent proposal for a coaction on one-loop Feynman integrals in dimensional regularization, we use intersection numbers and twisted homology theory to define a coaction on certain hypergeometric functions. The functions we consider admit an integral representation where both the integrand and the contour of integration are associated with positive geometries. As in dimensionally-regularized Feynman integrals, endpoint singularities are regularized by means of exponents controlled by a small parameter $ε$. We show that the coaction defined on this class of integral is consistent, upon expansion in $ε$, with the well-known coaction on multiple polylogarithms. We illustrate the validity of our construction by explicitly determining the coaction on various types of hypergeometric ${}_{p+1}F_p$ and Appell functions.

preprint2019arXiv

Three-loop contributions to the $ρ$ parameter and iterated integrals of modular forms

We compute fully analytic results for the three-loop diagrams involving two different massive quark flavours contributing to the $ρ$ 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.

preprint2016arXiv

CP-even scalar boson production via gluon fusion at the LHC

In view of the searches at the LHC for scalar particle resonances in addition to the 125 GeV Higgs boson, we present the cross section for a CP-even scalar produced via gluon fusion at N3LO in perturbative QCD assuming that it couples directly to gluons in an effective theory approach. We refine our prediction by taking into account the possibility that the scalar couples to the top-quark and computing the corresponding contributions through NLO in perturbative QCD. We assess the theoretical uncertainties of the cross section due to missing higher-order QCD effects and we provide the necessary information for obtaining the cross section value and uncertainty from our results in specific scenarios beyond the Standard Model. We also give detailed results for the case of a 750 GeV scalar, which will be the subject of intense experimental studies.

preprint2016arXiv

High precision determination of the gluon fusion Higgs boson cross-section at the LHC

We present the most precise value for the Higgs boson cross-section in the gluon-fusion production mode at the LHC. Our result is based on a perturbative expansion through N$^3$LO in QCD, in an effective theory where the top-quark is assumed to be infinitely heavy, while all other Standard Model quarks are massless. We combine this result with QCD corrections to the cross-section where all finite quark-mass effects are included exactly through NLO. In addition, electroweak corrections and the first corrections in the inverse mass of the top-quark are incorporated at three loops. We also investigate the effects of threshold resummation, both in the traditional QCD framework and following a SCET approach, which resums a class of $π^2$ contributions to all orders. We assess the uncertainty of the cross-section from missing higher-order corrections due to both perturbative QCD effects beyond N$^3$LO and unknown mixed QCD-electroweak effects. In addition, we determine the sensitivity of the cross-section to the choice of parton distribution function (PDF) sets and to the parametric uncertainty in the strong coupling constant and quark masses. For a Higgs mass of $m_H = 125~{\rm GeV}$ and an LHC center-of-mass energy of $13~{\rm TeV}$, our best prediction for the gluon fusion cross-section is \[ σ= 48.58\,{\rm pb} {}^{+2.22\, {\rm pb}\, (+4.56\%)}_{-3.27\, {\rm pb}\, (-6.72\%)} \mbox{ (theory)} \pm 1.56 \,{\rm pb}\, (3.20\%) \mbox{ (PDF+$α_s$)} \]

preprint2016arXiv

Jet production in the CoLoRFulNNLO method: event shapes in electron-positron collisions

We present the CoLoRFulNNLO method to compute higher order radiative corrections to jet cross sections in perturbative QCD. We apply our method to the computation of event shape observables in electron-positron collisions at NNLO accuracy and validate our code by comparing our predictions to previous results in the literature. We also calculate for the first time jet cone energy fraction at NNLO.

preprint2016arXiv

Long-distance singularities in multi-leg scattering amplitudes

We report on the recent completion of the three-loop calculation of the soft anomalous dimension in massless gauge-theory scattering amplitudes. This brings the state-of-the-art knowledge of long-distance singularities in multi-leg QCD amplitudes with any number of massless particles to three loops. The result displays some novel features: this is the first time non-dipole corrections appear, which directly correlate the colour and kinematic degrees of freedom of four coloured partons. We find that non-dipole corrections appear at three loops also for three coloured partons, but these are independent of the kinematics. The final result is remarkably simple when expressed in terms of single-valued harmonic polylogarithms, and it satisfies several non-trivial constraints. In particular, it is consistent with the high-energy limit behaviour and it satisfies the expected factorization properties in two-particle collinear limits.

preprint2016arXiv

Multi-Regge kinematics and the moduli space of Riemann spheres with marked points

We show that scattering amplitudes in planar N = 4 Super Yang-Mills in multi-Regge kinematics can naturally be expressed in terms of single-valued iterated integrals on the moduli space of Riemann spheres with marked points. As a consequence, scattering amplitudes in this limit can be expressed as convolutions that can easily be computed using Stokes' theorem. We apply this framework to MHV amplitudes to leading-logarithmic accuracy (LLA), and we prove that at L loops all MHV amplitudes are determined by amplitudes with up to L + 4 external legs. We also investigate non-MHV amplitudes, and we show that they can be obtained by convoluting the MHV results with a certain helicity flip kernel. We classify all leading singularities that appear at LLA in the Regge limit for arbitrary helicity configurations and any number of external legs. Finally, we use our new framework to obtain explicit analytic results at LLA for all MHV amplitudes up to five loops and all non-MHV amplitudes with up to eight external legs and four loops.

preprint2016arXiv

Three-jet production in electron-positron collisions using the CoLoRFulNNLO method

We introduce a subtraction method for jet cross sections at next-to-next-to-leading order (NNLO) accuracy in the strong coupling and use it to compute event shapes in three-jet production in electron-positron collisions. We validate our method on two event shapes, thrust and C-parameter, which are already known in the literature at NNLO accuracy and compute for the first time oblateness and the energy-energy correlation at the same accuracy.

preprint2016arXiv

Three-loop corrections to the soft anomalous dimension in multi-leg scattering

We present the three-loop result for the soft anomalous dimension governing long-distance singularities of multi-leg gauge-theory scattering amplitudes of massless partons. We compute all contributing webs involving semi-infinite Wilson lines at three loops and obtain the complete three-loop correction to the dipole formula. We find that non-dipole corrections appear already for three coloured partons, where the correction is a constant without kinematic dependence. Kinematic dependence appears only through conformally-invariant cross ratios for four coloured partons or more, and the result can be expressed in terms of single-valued harmonic polylogarithms of weight five. While the non-dipole three-loop term does not vanish in two-particle collinear limits, its contribution to the splitting amplitude anomalous dimension reduces to a constant, and it only depends on the colour charges of the collinear pair, thereby preserving strict collinear factorization properties. Finally we verify that our result is consistent with expectations from the Regge limit.

preprint2015arXiv

Computing decay rates for new physics theories with FeynRules and MadGraph5/aMC@NLO

We present new features of the FeynRules and MadGraph5/aMC@NLO programs for the automatic computation of decay widths that consistently include channels of arbitrary final-state multiplicity. The implementations are generic enough so that they can be used in the framework of any quantum field theory, possibly including higher-dimensional operators. We extend at the same time the conventions of the Universal FeynRules Output (or UFO) format to include decay tables and information on the total widths. We finally provide a set of representative examples of the usage of the new functions of the different codes in the framework of the Standard Model, the Higgs Effective Field Theory, the Strongly Interacting Light Higgs model and the Minimal Supersymmetric Standard Model and compare the results to available literature and programs for validation purposes.

preprint2015arXiv

Higgs boson decay into b-quarks at NNLO accuracy

We compute the fully differential decay rate of the Standard Model Higgs boson into b-quarks at next-to-next-to-leading order (NNLO) accuracy in alpha_S. We employ a general subtraction scheme developed for the calculation of higher order perturbative corrections to QCD jet cross sections, which is based on the universal infrared factorization properties of QCD squared matrix elements. We show that the subtractions render the various contributions to the NNLO correction finite. In particular, we demonstrate analytically that the sum of integrated subtraction terms correctly reproduces the infrared poles of the two-loop double virtual contribution to this process. We present illustrative differential distributions obtained by implementing the method in a parton level Monte Carlo program. The basic ingredients of our subtraction scheme, used here for the first time to compute a physical observable, are universal and can be employed for the computation of more involved processes.

preprint2015arXiv

Higgs boson gluon-fusion production in N3LO QCD

We present the cross-section for the production of a Higgs boson at hadron-colliders at next-to- next-to-next-to-leading order (N3LO) in perturbative QCD. The calculation is based on a method to perform a series expansion of the partonic cross-section around the threshold limit to an arbitrary order. We perform this expansion to sufficiently high order to obtain the value of the hadronic cross at N3LO in the large top-mass limit. For renormalisation and factorisation scales equal to half the Higgs mass, the N3LO corrections are of the order of +2.2%. The total scale variation at N3LO is 3%, reducing the uncertainty due to missing higher order QCD corrections by a factor of three.

preprint2015arXiv

Soft Expansion of Double-Real-Virtual Corrections to Higgs Production at N$^3$LO

We present methods to compute higher orders in the threshold expansion for the one-loop production of a Higgs boson in association with two partons at hadron colliders. This process contributes to the N$^3$LO Higgs production cross section beyond the soft-virtual approximation. We use reverse unitarity to expand the phase-space integrals in the small kinematic parameters and to reduce the coefficients of the expansion to a small set of master integrals. We describe two methods for the calculation of the master integrals. The first was introduced for the calculation of the soft triple-real radiation relevant to N$^3$LO Higgs production. The second uses a particular factorization of the three body phase-space measure and the knowledge of the scaling properties of the integral itself. Our result is presented as a Laurent expansion in the dimensional regulator, although some of the master integrals are computed to all orders in this parameter.

preprint2014arXiv

Bootstrapping six-gluon scattering in planar ${\cal N}=4$ super-Yang-Mills theory

We describe the hexagon function bootstrap for solving for six-gluon scattering amplitudes in the large $N_c$ limit of ${\cal N}=4$ super-Yang-Mills theory. In this method, an ansatz for the finite part of these amplitudes is constrained at the level of amplitudes, not integrands, using boundary information. In the near-collinear limit, the dual picture of the amplitudes as Wilson loops leads to an operator product expansion which has been solved using integrability by Basso, Sever and Vieira. Factorization of the amplitudes in the multi-Regge limit provides additional boundary data. This bootstrap has been applied successfully through four loops for the maximally helicity violating (MHV) configuration of gluon helicities, and through three loops for the non-MHV case.

preprint2014arXiv

FeynRules 2.0 - A complete toolbox for tree-level phenomenology

FeynRules is a Mathematica-based package which addresses the implementation of particle physics models, which are given in the form of a list of fields, parameters and a Lagrangian, into high-energy physics tools. It calculates the underlying Feynman rules and outputs them to a form appropriate for various programs such as CalcHEP, FeynArts, MadGraph, Sherpa and Whizard. Since the original version, many new features have been added: support for two-component fermions, spin-3/2 and spin-2 fields, superspace notation and calculations, automatic mass diagonalization, completely general FeynArts output, a new universal FeynRules output interface, a new Whizard interface, automatic 1 to 2 decay width calculation, improved speed and efficiency, new guidelines for validation and a new web-based validation package. With this feature set, FeynRules enables models to go from theory to simulation and comparison with experiment quickly, efficiently and accurately.

preprint2014arXiv

From multiple unitarity cuts to the coproduct of Feynman integrals

We develop techniques for computing and analyzing multiple unitarity cuts of Feynman integrals, and reconstructing the integral from these cuts. We study the relations among unitarity cuts of a Feynman integral computed via diagrammatic cutting rules, the discontinuity across the corresponding branch cut, and the coproduct of the integral. For single unitarity cuts, these relations are familiar. Here we show that they can be generalized to sequences of unitarity cuts in different channels. Using concrete one- and two-loop scalar integral examples we demonstrate that it is possible to reconstruct a Feynman integral from either single or double unitarity cuts. Our results offer insight into the analytic structure of Feynman integrals as well as a new approach to computing them.

preprint2014arXiv

Higgs boson gluon-fusion production at threshold in N3LO QCD

We present the cross-section for the threshold production of the Higgs boson at hadron-colliders at next-to-next-to-next-to-leading order (N3LO) in perturbative QCD. We present an analytic expression for the partonic cross-section at threshold and the impact of these corrections on the numerical estimates for the hadronic cross-section at the LHC. With this result we achieve a major milestone towards a complete evaluation of the cross-section at N3LO which will reduce the theoretical uncertainty in the determination of the strengths of the Higgs boson interactions.

preprint2014arXiv

Higgs boson gluon-fusion production beyond threshold in N3LO QCD

In this article, we compute the gluon fusion Higgs boson cross-section at N3LO through the second term in the threshold expansion. This calculation constitutes a major milestone towards the full N3LO cross section. Our result has the best formal accuracy in the threshold expansion currently available, and includes contributions from collinear regions besides subleading corrections from soft and hard regions, as well as certain logarithmically enhanced contributions for general kinematics. We use our results to perform a critical appraisal of the validity of the threshold approximation at N3LO in perturbative QCD.

preprint2014arXiv

Mathematical aspects of scattering amplitudes

In these lectures we discuss some of the mathematical structures that appear when computing multi-loop Feynman integrals. We focus on a specific class of special functions, the so-called multiple polylogarithms, and discuss introduce their Hopf algebra structure. We show how these mathematical concepts are useful in physics by illustrating on several examples how these algebraic structures are useful to perform analytic computations of loop integrals, in particular to derive functional equations among polylogarithms.

preprint2014arXiv

NNLO QCD corrections to $pp \to γ^* γ^*$ in the large $N_F$ limit

We compute the NNLO QCD corrections for the hadroproduction of a pair of off-shell photons in the limit of a large number of quark flavors. We perform a reduction of the two-loop amplitude to master integrals and calculate the latter analytically as a Laurent series in the dimensional regulator using modern integration methods. Real radiation corrections are evaluated numerically with a direct subtraction of infrared limits which we cast in a simple factorized form. The results presented here constitute a gauge invariant part of the full NNLO corrections but are not necessarily dominant. We view this calculation as a step towards a complete computation. Our partial corrections to the total cross-section are about $1\%-3\%$ and vary with the virtuality of the two off-shell photons.

preprint2014arXiv

The BFKL equation, Mueller-Navelet jets and single-valued harmonic polylogarithms

We introduce a generating function for the coefficients of the leading logarithmic BFKL Green's function in transverse-momentum space, order by order in alpha_s, in terms of single-valued harmonic polylogarithms. As an application, we exhibit fully analytic azimuthal-angle and transverse-momentum distributions for Mueller-Navelet jet cross sections at each order in alpha_s. We also provide a generating function for the total cross section valid to any number of loops.

preprint2014arXiv

The four-loop remainder function and multi-Regge behavior at NNLLA in planar N=4 super-Yang-Mills theory

We present the four-loop remainder function for six-gluon scattering with maximal helicity violation in planar N=4 super-Yang-Mills theory, as an analytic function of three dual-conformal cross ratios. The function is constructed entirely from its analytic properties, without ever inspecting any multi-loop integrand. We employ the same approach used at three loops, writing an ansatz in terms of hexagon functions, and fixing coefficients in the ansatz using the multi-Regge limit and the operator product expansion in the near-collinear limit. We express the result in terms of multiple polylogarithms, and in terms of the coproduct for the associated Hopf algebra. From the remainder function, we extract the BFKL eigenvalue at next-to-next-to-leading logarithmic accuracy (NNLLA), and the impact factor at NNNLLA. We plot the remainder function along various lines and on one surface, studying ratios of successive loop orders. As seen previously through three loops, these ratios are surprisingly constant over large regions in the space of cross ratios, and they are not far from the value expected at asymptotically large orders of perturbation theory.

preprint2014arXiv

Two-loop splitting amplitudes and the single-real contribution to inclusive Higgs production at N3LO

The factorisation of QCD matrix elements in the limit of two external partons becoming collinear is described by process-independent splitting amplitudes, which can be expanded systematically in perturbation theory. Working in conventional dimensional regularisation, we compute the two-loop splitting amplitudes for all simple collinear splitting processes, including subleading terms in the regularisation parameter. Our results are then applied to derive an analytical expression for the two-loop single-real contribution to inclusive Higgs boson production in gluon fusion to fourth order (N3LO) in perturbative QCD.

preprint2013arXiv

Leading singularities and off-shell conformal integrals

The three-loop four-point function of stress-tensor multiplets in N=4 super Yang-Mills theory contains two so far unknown, off-shell, conformal integrals, in addition to the known, ladder-type integrals. In this paper we evaluate the unknown integrals, thus obtaining the three-loop correlation function analytically. The integrals have the generic structure of rational functions multiplied by (multiple) polylogarithms. We use the idea of leading singularities to obtain the rational coefficients, the symbol - with an appropriate ansatz for its structure - as a means of characterising multiple polylogarithms, and the technique of asymptotic expansion of Feynman integrals to obtain the integrals in certain limits. The limiting behaviour uniquely fixes the symbols of the integrals, which we then lift to find the corresponding polylogarithmic functions. The final formulae are numerically confirmed. The techniques we develop can be applied more generally, and we illustrate this by analytically evaluating one of the integrals contributing to the same four-point function at four loops. This example shows a connection between the leading singularities and the entries of the symbol.

preprint2013arXiv

New developments in FeynRules

The program FeynRules is a Mathematica package developed to facilitate the implementation of new physics theories into high-energy physics tools. Starting from a minimal set of information such as the model gauge symmetries, its particle content, parameters and Lagrangian, FeynRules provides all necessary routines to extract automatically from the Lagrangian (that can also be computed semi-automatically for supersymmetric theories) the associated Feynman rules. These can be further exported to several Monte Carlo event generators through dedicated interfaces, as well as translated into a Python library, under the so-called UFO model format, agnostic of the model complexity, especially in terms of Lorentz and/or color structures appearing in the vertices or of number of external legs. In this work, we briefly report on the most recent new features that have been added to FeynRules, including full support for spin-3/2 fermions, a new module allowing for the automated diagonalization of the particle spectrum and a new set of routines dedicated to decay width calculations.

preprint2013arXiv

Real-Virtual contributions to the inclusive Higgs cross-section at N3LO

We compute the contributions to the N3LO inclusive Higgs boson cross-section from the square of one-loop amplitudes with a Higgs boson and three QCD partons as external states. Our result is a Taylor expansion in the dimensional regulator epsilon, where the coefficients of the expansion are analytic functions of the ratio of the Higgs boson mass and the partonic center of mass energy and they are valid for arbitrary values of this ratio. We also perform a threshold expansion around the limit of soft-parton radiation in the final state. The expressions for the coefficients of the threshold expansion are valid for arbitrary values of the dimension. As a by-product of the threshold expansion calculation, we have developed a soft expansion method at the integrand level by identifying the relevant soft and collinear regions for the loop-momentum.

preprint2013arXiv

Soft triple-real radiation for Higgs production at N3LO

We present the two first terms in the threshold expansion of Higgs production partonic cross-sections at hadron colliders for processes with three partons in the final state. These are contributions to the inclusive Higgs cross-section in gluon fusion at N3LO. We have developed a new technique for the expansion of the squared matrix-elements around the soft limit and for the reduction of the required phase-space integrals to only ten single-scale master integrals. We compute the master integrals building upon modern techniques for the integration of multidimensional integrals in dimensional regularization. Our results constitute an important step towards a systematic computation of the Higgs boson cross-section as an expansion around the threshold limit.

preprint2013arXiv

The two-loop soft current in dimensional regularization

The soft current describes the factorization behavior of quantum chromodynamics (QCD) scattering amplitudes in the limit of vanishing energy of one of the external partons. It is process-independent and can be expanded in a perturbative series in the coupling constant. To all orders in the dimensional regularization parameter, we compute the two-loop correction to the soft current for processes involving two hard partons.

preprint2012arXiv

From Lagrangians to Events: Computer Tutorial at the MC4BSM-2012 Workshop

This is a written account of the computer tutorial offered at the Sixth MC4BSM workshop at Cornell University, March 22-24, 2012. The tools covered during the tutorial include: FeynRules, LanHEP, MadGraph, CalcHEP, Pythia 8, Herwig++, and Sherpa. In the tutorial, we specify a simple extension of the Standard Model, at the level of a Lagrangian. The software tools are then used to automatically generate a set of Feynman rules, compute the invariant matrix element for a sample process, and generate both parton-level and fully hadronized/showered Monte Carlo event samples. The tutorial is designed to be self-paced, and detailed instructions for all steps are included in this write-up. Installation instructions for each tool on a variety of popular platforms are also provided.

preprint2012arXiv

Hopf algebras, coproducts and symbols: an application to Higgs boson amplitudes

We show how the Hopf algebra structure of multiple polylogarithms can be used to simplify complicated expressions for multi-loop amplitudes in perturbative quantum field theory and we argue that, unlike the recently popularized symbol-based approach, the coproduct incorporates information about the zeta values. We illustrate our approach by rewriting the two-loop helicity amplitudes for a Higgs boson plus three gluons in a simplified and compact form involving only classical polylogarithms.

preprint2012arXiv

Infrared singularities and the high-energy limit

We review recent results on the high-energy limit of gauge amplitudes, which can be derived from the universal properties of their infrared singularities. Using the dipole formula, a compact ansatz for infrared singularities of massless gauge amplitudes, and taking the high-energy limit, we provide a simple expression for the soft factor of a generic high-energy amplitude, valid to leading power in $t/s$ and to all logarithmic orders. This gives a direct and general proof of leading-logarithmic Reggeization for infrared divergent contributions to the amplitude, and it shows how Reggeization breaks down at NNLL level. We further show how the dipole formula constrains the high-energy limit of multi-particle amplitudes in multi-Regge kinematics, and how, on the other hand, Regge theory constrains possible corrections to the dipole formula.

preprint2012arXiv

Infrared singularities in the high-energy limit

We use our current understanding of the all-order singularity structure of gauge theory amplitudes to probe their high-energy limit. Our starting point is the dipole formula, a compact ansatz for the soft anomalous dimension matrix of massless multi-particle amplitudes. In the high-energy limit, we find a simple and general expression for the infrared factor generating all soft and collinear singularities of the amplitude, which is valid to leading power in $|t|/s$ and to all logarithmic orders. This leads to a direct and general proof of leading-logarithmic Reggeization for infrared divergent contributions to the amplitude. Furthermore, we can prove explicitly that the simplest form of Reggeization, based on the absence of Regge cuts in the complex angular momentum plane, breaks down at the NNLL level. Finally, we note that the known features of the high-energy limit can be used to constrain possible corrections to the dipole formula, starting at the three-loop order.

preprint2012arXiv

Introducing an interface between WHIZARD and FeynRules

While Monte Carlo event generators like WHIZARD have become indispensable tools in studying the impact of new physics on collider observables over the last decades, the implementation of new models in such packages has remained a rather awkward and error-prone process. Recently, the FeynRules package was introduced which greatly simplifies this process by providing a single unified model format from which model implementations for many different Monte Carlo codes can be derived automatically. In this note, we present an interface which extends FeynRules to provide this functionality also for the WHIZARD package, thus making WHIZARD's strengths and performance easily available to model builders.

preprint2012arXiv

NNLO phase space master integrals for two-to-one inclusive cross sections in dimensional regularization

We evaluate all phase space master integrals which are required for the total cross section of generic 2 -> 1 processes at NNLO as a series expansion in the dimensional regulator epsilon. Away from the limit of threshold production, our expansion includes one order higher than what has been available in the literature. At threshold, we provide expressions which are valid to all orders in terms of Gamma functions and hypergeometric functions. These results are a necessary ingredient for the renormalization and mass factorization of singularities in 2 -> 1 inclusive cross sections at NNNLO in QCD.

preprint2012arXiv

Single-valued harmonic polylogarithms and the multi-Regge limit

We argue that the natural functions for describing the multi-Regge limit of six-gluon scattering in planar N=4 super Yang-Mills theory are the single-valued harmonic polylogarithmic functions introduced by Brown. These functions depend on a single complex variable and its conjugate, (w,w*). Using these functions, and formulas due to Fadin, Lipatov and Prygarin, we determine the six-gluon MHV remainder function in the leading-logarithmic approximation (LLA) in this limit through ten loops, and the next-to-LLA (NLLA) terms through nine loops. In separate work, we have determined the symbol of the four-loop remainder function for general kinematics, up to 113 constants. Taking its multi-Regge limit and matching to our four-loop LLA and NLLA results, we fix all but one of the constants that survive in this limit. The multi-Regge limit factorizes in the variables (ν,n) which are related to (w,w*) by a Fourier-Mellin transform. We can transform the single-valued harmonic polylogarithms to functions of (ν,n) that incorporate harmonic sums, systematically through transcendental weight six. Combining this information with the four-loop results, we determine the eigenvalues of the BFKL kernel in the adjoint representation to NNLLA accuracy, and the MHV product of impact factors to NNNLLA accuracy, up to constants representing beyond-the-symbol terms and the one symbol-level constant. Remarkably, only derivatives of the polygamma function enter these results. Finally, the LLA approximation to the six-gluon NMHV amplitude is evaluated through ten loops.

preprint2012arXiv

Three-mass triangle integrals and single-valued polylogarithms

We study one and two-loop triangle integrals with massless propagators and all external legs off shell. We show that there is a kinematic region where the results can be expressed in terms of a basis of single-valued polylogarithms in one complex variable. The relevant space of single-valued functions can be determined a priori and the results take strikingly a simple and compact form when written in terms of this basis. We study the properties of the basis functions and illustrate how one can easily analytically continue our results to all kinematic regions where the external masses have the same sign.

preprint2012arXiv

UFO - The Universal FeynRules Output

We present a new model format for automatized matrix-element generators, the so- called Universal FeynRules Output (UFO). The format is universal in the sense that it features compatibility with more than one single generator and is designed to be flexible, modular and agnostic of any assumption such as the number of particles or the color and Lorentz structures appearing in the interaction vertices. Unlike other model formats where text files need to be parsed, the information on the model is encoded into a Python module that can easily be linked to other computer codes. We then describe an interface for the Mathematica package FeynRules that allows for an automatic output of models in the UFO format.

preprint2011arXiv

A superspace module for the FeynRules package

We describe an additional module for the Mathematica package FeynRules that allows for an easy building of any N=1 supersymmetric quantum field theory, directly in superspace. After the superfield content of a specific model has been implemented, the user can study the properties of the model, such as the supersymmetric transformation laws of the associated Lagrangian, directly in Mathematica. While the model dependent parts of the latter, i.e., the soft supersymmetry-breaking Lagrangian and the superpotential, have to be provided by the user, the model independent pieces, such as the gauge interaction terms, are derived automatically. Using the strengths of the Feynrules program, it is then possible to derive all the Feynman rules associated to the model and implement them in all the Feynman diagram calculators interfaced to FeynRules in a straightforward way.

preprint2011arXiv

An infrared approach to Reggeization

We present a new approach to Reggeization of gauge amplitudes based on the universal properties of their infrared singularities. Using the "dipole formula", a compact ansatz for all infrared singularities of massless amplitudes, we study Reggeization of singular contributions to high-energy amplitudes for arbitrary color representations, and any logarithmic accuracy. We derive leading-logarithmic Reggeization for general cross-channel color exchanges, and we show that Reggeization breaks down for the imaginary part of the amplitude at next-to-leading logarithms and for the real part at next-to-next-to-leading logarithms. Our formalism applies to multiparticle amplitudes in multi-Regge kinematics, and constrains possible corrections to the dipole formula starting at three loops.

preprint2011arXiv

From polygons and symbols to polylogarithmic functions

We present a review of the symbol map, a mathematical tool that can be useful in simplifying expressions among multiple polylogarithms, and recall its main properties. A recipe is given for how to obtain the symbol of a multiple polylogarithm in terms of the combinatorial properties of an associated rooted decorated polygon. We also outline a systematic approach to constructing a function corresponding to a given symbol, and illustrate it in the particular case of harmonic polylogarithms up to weight four. Furthermore, part of the ambiguity of this process is highlighted by exhibiting a family of non-trivial elements in the kernel of the symbol map for arbitrary weight.

preprint2011arXiv

The infrared structure of gauge theory amplitudes in the high-energy limit

We develop an approach to the high-energy limit of gauge theories based on the universal properties of their infrared singularities. Our main tool is the dipole formula, a compact ansatz for the all-order infrared singularity structure of scattering amplitudes of massless partons. By taking the high-energy limit, we show that the dipole formula implies Reggeization of infrared-singular contributions to the amplitude, at leading logarithmic accuracy, for the exchange of arbitrary color representations in the cross channel. We observe that the real part of the amplitude Reggeizes also at next-to-leading logarithmic order, and we compute the singular part of the two-loop Regge trajectory, which is universally expressed in terms of the cusp anomalous dimension. Our approach provides tools to study the high-energy limit beyond the boundaries of Regge factorization: thus we show that Reggeization generically breaks down at next-to-next-to-leading logarithmic accuracy, and provide a general expression for the leading Reggeization-breaking operator. Our approach applies to multiparticle amplitudes in multi-Regge kinematics, and it also implies new constraints on possible corrections to the dipole formula, based on the Regge limit.

preprint2011arXiv

The one-loop six-dimensional hexagon integral with three massive corners

We compute the six-dimensional hexagon integral with three non-adjacent external masses analytically. After a simple rescaling, it is given by a function of six dual conformally invariant cross-ratios. The result can be expressed as a sum of 24 terms involving only one basic function, which is a simple linear combination of logarithms, dilogarithms, and trilogarithms of uniform degree three transcendentality. Our method uses differential equations to determine the symbol of the function, and an algorithm to reconstruct the latter from its symbol. It is known that six-dimensional hexagon integrals are closely related to scattering amplitudes in N=4 super Yang-Mills theory, and we therefore expect our result to be helpful for understanding the structure of scattering amplitudes in this theory, in particular at two loops.

preprint2010arXiv

A Two-Loop Octagon Wilson Loop in N = 4 SYM

In the planar N = 4 supersymmetric Yang-Mills theory at weak coupling, we perform the first analytic computation of a two-loop eight-edged Wilson loop embedded into the boundary of AdS3. Its remainder function is given as a function of uniform transcendental weight four in terms of a constant plus a product of four logarithms. We compare to the strong-coupling result, and test a conjecture on the universality of the remainder function proposed in the literature.

preprint2010arXiv

An Analytic Result for the Two-Loop Hexagon Wilson Loop in N = 4 SYM

In the planar N=4 supersymmetric Yang-Mills theory, the conformal symmetry constrains multi-loop n-edged Wilson loops to be basically given in terms of the one-loop n-edged Wilson loop, augmented, for n greater than 6, by a function of conformally invariant cross ratios. We identify a class of kinematics for which the Wilson loop exhibits exact Regge factorisation and which leave invariant the analytic form of the multi-loop n-edged Wilson loop. In those kinematics, the analytic result for the Wilson loop is the same as in general kinematics, although the computation is remarkably simplified with respect to general kinematics. Using the simplest of those kinematics, we have performed the first analytic computation of the two-loop six-edged Wilson loop in general kinematics.

preprint2010arXiv

The Two-Loop Hexagon Wilson Loop in N = 4 SYM

In the planar N=4 supersymmetric Yang-Mills theory, the conformal symmetry constrains multi-loop n-edged Wilson loops to be given in terms of the one-loop n-edged Wilson loop, augmented, for n greater than 6, by a function of conformally invariant cross ratios. That function is termed the remainder function. In a recent paper, we have displayed the first analytic computation of the two-loop six-edged Wilson loop, and thus of the corresponding remainder function. Although the calculation was performed in the quasi-multi-Regge kinematics of a pair along the ladder, the Regge exactness of the six-edged Wilson loop in those kinematics entails that the result is the same as in general kinematics. We show in detail how the most difficult of the integrals is computed, which contribute to the six-edged Wilson loop. Finally, the remainder function is given as a function of uniform transcendental weight four in terms of Goncharov polylogarithms. We consider also some asymptotic values of the remainder function, and the value when all the cross ratios are equal.

preprint2009arXiv

A comprehensive approach to new physics simulations

We describe a framework to develop, implement and validate any perturbative Lagrangian-based particle physics model for further theoretical, phenomenological and experimental studies. The starting point is FeynRules, a Mathematica package that allows to generate Feynman rules for any Lagrangian and then, through dedicated interfaces, automatically pass the corresponding relevant information to any supported Monte Carlo event generator. We prove the power, robustness and flexibility of this approach by presenting a few examples of new physics models (the Hidden Abelian Higgs Model, the general Two-Higgs-Doublet Model, the most general Minimal Supersymmetric Standard Model, the Minimal Higgsless Model, Universal and Large Extra Dimensions, and QCD-inspired effective Lagrangians) and their implementation/validation in FeynArts/FormCalc, CalcHep, MadGraph/MadEvent, and Sherpa.

preprint2009arXiv

The one-loop pentagon to higher orders in epsilon

We compute the one-loop scalar massless pentagon integral I_5^{6-2 eps} in D=6-2\eps dimensions in the limit of multi-Regge kinematics. This integral first contributes to the parity-odd part of the one-loop N=4 five-point MHV amplitude m_5^{(1)} at O(eps). In the high energy limit defined, the pentagon integral reduces to double sums or equivalently two-fold Mellin-Barnes integrals. By determining the O(eps) contribution to I_5^{6-2 eps}, one therefore gains knowledge of m_5^{(1)} through to O(eps^2) which is necessary for studies of the iterative structure of N=4 SYM amplitudes beyond one-loop. One immediate application is the extraction of the one-loop gluon-production vertex through to O(eps^2) and the iterative construction of the two-loop gluon-production vertex through to finite terms which is described in a companion paper. The analytic methods we have used for evaluating the pentagon integral in the high energy limit may also be applied to the hexagon integral and may ultimately give information on the form of the remainder function.

preprint2006arXiv

Color-dressed recursive relations for multi-parton amplitudes

Remarkable progress inspired by twistors has lead to very simple analytic expressions and to new recursive relations for multi-parton color-ordered amplitudes. We show how such relations can be extended to include color and present the corresponding color-dressed formulation for the Berends-Giele, BCF and a new kind of CSW recursive relations. A detailed comparison of the numerical efficiency of the different approaches to the calculation of multi-parton cross sections is performed.