BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Stefan Weinzierl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We study a two-loop four-point function with one internal mass. This Feynman integral is one of the simplest Feynman integrals depending on two elliptic curves. We transform the associated differential equation into an $\varepsilon$-form. We study the entries of the differential equation, and in particular the entries which depend on both elliptic curves.
In this talk we discuss the interplay of two elliptic curves, which occur in different sub-sectors of Feynman integrals. We analyse a particular Feynman integral depending on two elliptic curves and derive an associated differential equation in $\varepsilon$-form. We discuss the mixed entries of the differential equation, which depend on both elliptic curves.
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.
We show that methods developed in the context of perturbative calculations can be transferred to non-perturbative calculations. We demonstrate that correlation functions on the lattice can be computed with the method of differential equations, supplemented with techniques from twisted cohomology. We derive differential equations for the variation with the coupling or -- more generally -- with the parameters of the action. Already simple examples show that the differential equation with respect to the coupling has an essential singularity at zero coupling and a regular singularity at infinite coupling. The properties of the differential equation at zero coupling can be used to prove that the perturbative series is only an asymptotic series.
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.
We investigate the behaviour of elliptic Feynman integrals under modular transformations. This has a practical motivation: Through a suitable modular transformation we can achieve that the nome squared is a small quantity, leading to fast numerical evaluations. Contrary to the case of multiple polylogarithms, where it is sufficient to consider just variable transformations for the numerical evaluations of multiple polylogarithms, it is more natural in the elliptic case to consider a combination of a variable transformation (i.e. a modular transformation) together with a redefinition of the master integrals. Thus we combine a coordinate transformation on the base manifold with a basis transformation in the fibre. Only in the combination of the two transformations we stay within the same class of functions.
We study linear relations among correlation functions on a lattice obtained from integration-by-parts identities. We use the framework of twisted cocycles and determine for a scalar theory a basis of correlation functions, in which all other correlation functions may be expressed.
Using loop-tree duality, we relate a renormalised $n$-point $l$-loop amplitude in a quantum field theory to a phase-space integral of a regularised $l$-fold forward limit of a UV-subtracted $(n+2l)$-point tree-amplitude-like object. We show that up to three loops the latter object is easily computable from recurrence relations. This defines an integrand of the loop amplitude with a global definition of the loop momenta. Field and mass renormalisation are performed in the on-shell scheme.
The computation of Feynman integrals often involves square roots. One way to obtain a solution in terms of multiple polylogarithms is to rationalize these square roots by a suitable variable change. We present a program that can be used to find such transformations. After an introduction to the theoretical background, we explain in detail how to use the program in practice.
We solve the two-loop sunrise integral with unequal masses systematically to all orders in the dimensional regularisation parameter $\varepsilon$. In order to do so, we transform the system of differential equations for the master integrals to an $\varepsilon$-form. The sunrise integral with unequal masses depends on three kinematical variables. We perform a change of variables to standard coordinates on the moduli space ${\mathcal M}_{1,3}$ of a genus one Riemann surface with three marked points. This gives us the solution as iterated integrals on $\overline{\mathcal M}_{1,3}$. On the hypersurface $τ=\mbox{const}$ our result reduces to elliptic polylogarithms. In the equal mass case our result reduces to iterated integrals of modular forms.
A walk on sunset boulevard can teach us about transcendental functions associated to Feynman diagrams. On this guided tour we will see multiple polylogarithms, differential equations and elliptic curves. A highlight of the tour will be the generalisation of the polylogarithms to the elliptic setting and the all-order solution for the sunset integral in the equal mass case.
Inspired by the recent progress in the field of scattering amplitudes, we discuss hypothetical particles which can be characterised as the double copies of fermions -- in the same way gravitons can be viewed as double copies of gauge bosons. As the gravitons, these hypothetical particles interact only through gravitational interactions. We present two equivalent methods for the computation of the relevant scattering amplitudes. The hypothetical particles can be massive and non-relativistic.
In these lectures I discuss Hopf algebras and Dyson-Schwinger equations. The lectures start with an introduction to Hopf algebras, followed by a review where Hopf algebras occur in particles physics. The final part of these lectures is devoted to the relation between Hopf algebras and Dyson-Schwinger equations.
Numerical approaches to higher-order calculations often employ subtraction terms, both for the real emission and the virtual corrections. These subtraction terms have to be added back. In this paper we show that at NLO the real subtraction terms, the virtual subtraction terms, the integral representations of the field renormalisation constants and -- in the case of initial-state partons -- the integral representation for the collinear counterterm can be grouped together to give finite integrals, which can be evaluated numerically. This is useful for an extension towards NNLO.
In this talk we review relations and representations of primitive QCD tree amplitudes. Topics covered include the BCJ relations, the CHY representation, and the KLT relations. We will put a special emphasis on how these relations and representations generalise from pure Yang-Mills theory to QCD. The generalisation of the KLT relations from pure Yang-Mills to QCD includes the case of massive quarks. On the gravity side we then obtain hypothetical particles interacting with gravitational strength, which can be massive and non-relativistic.
We present a method to compute the Laurent expansion of the two-loop sunrise integral with equal non-zero masses to arbitrary order in the dimensional regularisation $\varepsilon$. This is done by introducing a class of functions (generalisations of multiple polylogarithms to include the elliptic case) and by showing that all integrations can be carried out within this class of functions.
We show that the Laurent series of the two-loop kite integral in $D=4-2\varepsilon$ space-time dimensions can be expressed in each order of the series expansion in terms of elliptic generalisations of (multiple) polylogarithms. Using differential equations we present an iterative method to compute any desired order. As an example, we give the first three orders explicitly.
We summarize recent computations with a class of elliptic generalizations of polylogarithms, arising from the massive sunrise integral. For the case of arbitrary masses we obtain results in two and four space-time dimensions. The iterated integral structure of our functions allows us to furthermore compute the equal mass case to arbitrary order.
This paper investigates how tree-level amplitudes with massless quarks, gluons and/or massless scalars transforming under a single copy of the gauge group can be expressed in the context of the scattering equations as a sum over the inequivalent solutions of the scattering equations. In the case where the amplitudes satisfy cyclic invariance, KK- and BCJ-relations the only modification is the generalisation of the permutation invariant function $E(z,p,\varepsilon)$. We present a method to compute the modified $\hat{E}(z,p,\varepsilon)$. The most important examples are tree amplitudes in ${\mathcal N}=4$ SYM and QCD amplitudes with one quark-antiquark pair and an arbitrary number of gluons. QCD amplitudes with two or more quark-antiquark pairs do not satisfy the BCJ-relations and require in addition a generalisation of the Parke-Taylor factors $C_σ(z)$. The simplest case of the QCD tree-level four-point amplitude with two quark-antiquark pairs is discussed explicitly.
In this talk I will focus on theoretical issues related to high precision determinations of the top mass. Several mass definitions are reviewed and their respective advantages and disadvantages are discussed. Precision determinations of the top mass will require a short-distance mass definition. I will summarise current work in this direction.
The fundamental BCJ-relation is a linear relation between primitive tree amplitudes with different cyclic orderings. The cyclic orderings differ by the insertion place of one gluon. The coefficients of the fundamental BCJ-relation are linear in the Lorentz invariants $2 p_i p_j$. The BCJ-relations are well established for pure gluonic amplitudes as well as for amplitudes in ${\mathcal N}=4$ super-Yang-Mills theory. Recently, it has been conjectured that the BCJ-relations hold also for QCD amplitudes. In this paper we give a proof of this conjecture. The proof is valid for massless and massive quarks.
In this paper we construct a CHY representation for all tree-level primitive QCD amplitudes. The quarks may be massless or massive. We define a generalised cyclic factor $\hat{C}(w,z)$ and a generalised permutation invariant function $\hat{E}(z,p,\varepsilon)$. The amplitude is then given as a contour integral encircling the solutions of the scattering equations with the product $\hat{C} \hat{E}$ as integrand. Equivalently, it is given as a sum over the inequivalent solutions of the scattering equations, where the summand consists of a Jacobian times the product $\hat{C} \hat{E}$. This representation separates information: The generalised cyclic factor does not depend on the helicities of the external particles, the generalised permutation invariant function does not depend on the ordering of the external particles.
In this talk, we discuss our recent computation of the two-loop sunrise integral with arbitrary non-zero particle masses. In two space-time dimensions, we arrive at a result in terms of elliptic dilogarithms. Near four space-time dimensions, we obtain a result which furthermore involves elliptic generalizations of Clausen and Glaisher functions.
In this talk we discuss the solution for the sunrise integral around two and four space-time dimensions in terms of a generalised elliptic version of the multiple polylogarithms. In two space-time dimensions we obtain a sum of three elliptic dilogarithms. The arguments of the elliptic dilogarithms have a nice geometric interpretation. In four space-time dimensions the sunrise integral can be expressed with the $ε^0$- and $ε^1$-solution around two dimensions, mass derivatives thereof and simpler terms.
We present the result for the finite part of the two-loop sunrise integral with unequal masses in four space-time dimensions in terms of the ${\mathcal O}(\varepsilon^0)$-part and the ${\mathcal O}(\varepsilon^1)$-part of the sunrise integral around two space-time dimensions. The latter two integrals are given in terms of elliptic generalisations of Clausen and Glaisher functions. Interesting aspects of the result for the ${\mathcal O}(\varepsilon^1)$-part of the sunrise integral around two space-time dimensions are the occurrence of depth two elliptic objects and the weights of the individual terms.
In this talk we discuss the color decomposition of tree-level and one-loop QCD amplitudes with arbitrary numbers of quarks and gluons. We present a method for the decomposition of partial amplitudes into primitive amplitudes, which is based on shuffle relations and is purely combinatorial. Closed formulae are derived, which do not require the inversion of a system of linear equations.
We present the decomposition of QCD partial amplitudes into primitive amplitudes at one-loop level and tree level for arbitrary numbers of quarks and gluons. Our method is based on shuffle relations. This method is purely combinatorial and does not require the inversion of a system of linear equations.
In this talk we report on first results from the NLO computation of Z production in association with five jets in hadron-hadron collisions. The results are obtained with the help of the numerical method, where apart from the phase space integration also the integration over the loop momentum is performed numerically. In addition we discuss several methods and techniques for the improvement of the Monte Carlo integration.
This paper addresses the question, whether the solutions of the scattering equations in four space-time dimensions can be expressed as rational functions of the momentum twistor variables. This is the case for $n\le5$ external particles. For general $n$ there are always two solutions, which are rational functions of the momentum twistor variables. However, the remaining solutions are in general not rational. In the case $n=6$ the remaining four solutions can be expressed as algebraic functions. These four solutions are constructed explicitly in this paper.
We present the two-loop sunrise integral with arbitrary non-zero masses in two space-time dimensions in terms of elliptic dilogarithms. We find that the structure of the result is as simple and elegant as in the equal mass case, only the arguments of the elliptic dilogarithms are modified. These arguments have a nice geometric interpretation.
In this pedagogical note I will discuss one-loop integrals, where (i) different regions of the integration region lead to divergences and (ii) where these divergences cancel in the sum over all regions. These integrals cannot be calculated without regularisation, in spite of the fact that they yield a finite result. A typical example where such integrals occur is the decay H --> gamma gamma.
The BCJ decomposition is a highly non-trivial property of gauge theories. In this paper we systematically construct an effective Lagrangian, whose Feynman rules automatically produce the BCJ numerators. The effective Lagrangian contains non-local terms. The difference between the standard Yang-Mills Lagrangian and the effective Lagrangian simplifies to zero.
In this talk, we discuss the speed-up of numerical calculations of jet observables by replacing the usual sum over all helicity amplitudes with an integral over a parametrisation for the parton polarisations called random polarisations. Random polarisations are a linear combination of helicity eigenstates multiplied by a phase factor depending on a so-called helicity angle. Instead of a summation over discrete helicities, random polarisations require an integration over the helicity angle. By combining this integral with the final-state phase space integral, we only have to evaluate one squared amplitude per phase space point instead of $2^n$ helicity amplitudes, where $n$ is the total number of particles in the process. While the technique itself has been known since 1998, so far there has been no way of using it with dipole subtraction, which is probably the most-used method for dealing with infrared divergences in NLO calculations. After giving detailed reasons for this statement, we propose a solution to this problem in terms of extending the existing subtraction method by a new term.
We present a method to construct a suitable contour deformation in loop momentum space for multi-loop integrals. This contour deformation can be used to perform the integration for multi-loop integrals numerically. The integration can be performed directly in loop momentum space without the introduction of Feynman or Schwinger parameters. The method can be applied to finite multi-loop integrals and to divergent multi-loop integrals with suitable subtraction terms. The algorithm extends techniques from the one-loop case to the multi-loop case. Examples at two and three loops are discussed explicitly.
In these lectures I discuss Feynman graphs and the associated Feynman integrals. Of particular interest are the classes functions, which appear in the evaluation of Feynman integrals. The most prominent class of functions is given by multiple polylogarithms. The algebraic properties of multiple polylogarithms are reviewed in the second part of these lectures. The final part of these lectures is devoted to Feynman integrals, which cannot be expressed in terms of multiple polylogarithms. Methods from algebraic geometry provide tools to tackle these integrals.
There is a fruitful interplay between algebraic geometry on the one side and perturbative quantum field theory on the other side. I review the main relevant mathematical concepts of periods, Hodge structures and Picard-Fuchs equations and discuss the connection with Feynman integrals.
We present a systematic method to derive an ordinary differential equation for any Feynman integral, where the differentiation is with respect to an external variable. The resulting differential equation is of Fuchsian type. The method can be used within fixed integer space-time dimensions as well as within dimensional regularisation. We show that finding the differential equation is equivalent to solving a linear system of equations. We observe interesting factorisation properties of the D-dimensional Picard-Fuchs operator when D is specialised to integer dimensions.
We discuss the analytical solution of the two-loop sunrise graph with arbitrary non-zero masses in two space-time dimensions. The analytical result is obtained by solving a second-order differential equation. The solution involves elliptic integrals and in particular the solutions of the corresponding homogeneous differential equation are given by periods of an elliptic curve.
We derive a second-order differential equation for the two-loop sunrise graph in two dimensions with arbitrary masses. The differential equation is obtained by viewing the Feynman integral as a period of a variation of a mixed Hodge structure, where the variation is with respect to the external momentum squared. The fibre is the complement of an elliptic curve. From the fact that the first cohomology group of this elliptic curve is two-dimensional we obtain a second-order differential equation. This is an improvement compared to the usual way of deriving differential equations: Integration-by-parts identities lead only to a coupled system of four first-order differential equations.
We present a method, which constructs a suitable deformation vector in loop momentum space, when the loop integration is done numerically with the help of the subtraction method. The method presented here extends previously discussed techniques from the massless case to the general case of arbitrary masses in the loop.
In this paper we discuss techniques, which lead to a significant improvement of the efficiency of the Monte Carlo integration, when one-loop QCD amplitudes are calculated numerically with the help of the subtraction method and contour deformation. The techniques discussed are: holomorphic and non-holomorphic division into sub-channels, optimisation of the integration contour, improvement of the ultraviolet subtraction terms, importance sampling and antithetic variates in loop momentum space, recurrence relations.
We extend the dipole formalism for massless and massive partons to random polarisations of the external partons. The dipole formalism was originally formulated for spin-summed matrix elements and later extended to individual helicity eigenstates. For efficiency reasons one wants to replace the spin sum by a smooth integration over additional variables. This requires the extension of the dipole formalism to random polarisations. In this paper we derive the modified subtraction terms. We only modify the real subtraction terms, the integrated subtraction terms do not require any modifications.
In this talk I will give an overview on theoretical aspects of top quark physics. The focus lies on top pair production and single top production.
This article discusses the occurences of one-loop amplitudes within a next-to-next-to-leading order calculation. In an NNLO calculation the one-loop amplitude enters squared and one would therefore naively expect that the O(epsilon)- and O(epsilon^2)-terms of the one-loop amplitudes are required. I show that the calculation of these terms can be avoided if a method is known, which computes the O(epsilon^0)-terms of the finite remainder function of the two-loop amplitude.
This article gives results on several jet algorithms in electron-positron annihilation: Considered are the exclusive sequential recombination algorithms Durham, Geneva, Jade-E0 and Cambridge, which are typically used in electron-positron annihilation. In addition also inclusive jet algorithms are studied. Results are provided for the inclusive sequential recombination algorithms Durham, Aachen and anti-kt, as well as the infrared-safe cone algorithm SISCone. The results are obtained in perturbative QCD and are NNNLO for the two-jet rates, NNLO for the three-jet rates, NLO for the four-jet rates and LO for the five-jet rates.
We present next-to-leading order corrections in the leading colour approximation for jet rates in electron-positron annihilation up to seven jets. The results for the two-, three-, and four-jet rates agree with known results. The NLO jet rates have been known previously only up to five jets. The results for the six- and seven-jet rate are new. The results are obtained by a new and efficient method based on subtraction and numerical integration.
The SISCone jet algorithm is a seedless infrared-safe cone jet algorithm. There exists an implementation which is highly optimised for a large number of final state particles. However, in fixed-order perturbative calculations with a small number of final state particles, it turns out that the computer time needed for the jet clustering of this implementation is comparable to the computer time of the matrix elements. This article reports on an implementation of the SISCone algorithm optimised for low particle multiplicities.
We show that a well-known simple formula for the explicit infrared poles of one-loop QCD amplitudes has a corresponding simple counterpart in unintegrated form. The unintegrated formula approximates the integrand of one-loop QCD amplitudes in all soft and collinear singular regions. It thus defines a local counter-term for the infrared singularities and can be used as an ingredient for the numerical calculation of one-loop amplitudes.
The integrand of any multi-loop integral is characterised after Feynman parametrisation by two polynomials. In this review we summarise the properties of these polynomials. Topics covered in this article include among others: Spanning trees and spanning forests, the all-minors matrix-tree theorem, recursion relations due to contraction and deletion of edges, Dodgson's identity and matroids.
In these lectures I will give an introduction to Feynman integrals. In the first part of the course I review the basics of the perturbative expansion in quantum field theories. In the second part of the course I will discuss more advanced topics: Mathematical aspects of loop integrals related to periods, shuffle algebras and multiple polylogarithms are covered as well as practical algorithms for evaluating Feynman integrals.
The strong coupling constant is a fundamental parameter of nature. It can be extracted from experiments measuring three-jet events in electron-positron annihilation. For this extraction precise theoretical calculations for jet rates and event shapes are needed. In this talk I will discuss the NNLO calculation for these observables.
We present an algorithm for the numerical calculation of one-loop QCD amplitudes. The algorithm consists of subtraction terms, approximating the soft, collinear and ultraviolet divergences of one-loop amplitudes and a method to deform the integration contour for the loop integration into the complex space. The algorithm is formulated at the amplitude level and does not rely on Feynman graphs. Therefore all required ingredients can be calculated efficiently using recurrence relations. The algorithm applies to massless partons as well as to massive partons.
Starting from the standard Lagrangian for a SU(2) x U(1) gauge theory plus a Higgs field we derive the corresponding "maximal helicity violating" (MHV) Lagrangian. From this MHV Lagrangian one deduces simple diagrammatic rules for the calculation of multi-particle scattering amplitudes. We arrive at the MHV Lagrangian by a canonical change of the field variables in the light-cone gauge. We comment on the modifications which occur in a spontaneously broken gauge theory as compared to a pure (unbroken) Yang-Mills theory.
We consider on-shell recursion relations for all Born QCD amplitudes. This includes amplitudes with several pairs of quarks and massive quarks. We give a detailed description on how to shift the external particles in spinor space and clarify the allowed helicities of the shifted legs. We proof that the corresponding meromorphic functions vanish at z --> infinity. As an application we obtain compact expressions for helicity amplitudes including a pair of massive quarks, one negative helicity gluon and an arbitrary number of positive helicity gluons.
We consider one-loop tensor and scalar integrals, which occur in a massless quantum field theory and we report on the implementation into a numerical program of an algorithm for the automated computation of these one-loop integrals. The number of external legs of the loop integrals is not restricted. All calculations are done within dimensional regularization.
We review recent developments for the calculation of Born amplitudes in QCD. This includes the computation of gluon helicity amplitudes from MHV vertices and an approach based on scalar propagators and a set of three- and four-valent vertices. The latter easily generalizes to amplitudes with any number of quark pairs. The quarks may be massless or massive.
We study the systematic inclusion of QED corrections in the evolution of parton distributions. O(alpha) corrections modify the evolution equation for parton distributions. They introduce additional parton distributions, like the photon distribution in the nucleon, and lead to additional mixing effects. We discuss the modifications for a realistic model of N_f,up up-type flavours and N_f,down down-type flavours. We have implemented these corrections into a numerical program and we quantify the size of these effects in a toy model. The corrections reach the order of 1%.
Perturbative calculations at next-to-next-to-leading order for multi-particle final states require a method to cancel infrared singularities. I discuss how to setup the subtraction method at NNLO.
In the recent years, Hopf algebras have been introduced to describe certain combinatorial properties of quantum field theories. I will give a basic introduction to these algebras and review some occurrences in particle physics.
Expansion of higher transcendental functions in a small parameter are needed in many areas of science. For certain classes of functions this can be achieved by algebraic means. These algebraic tools are based on nested sums and can be formulated as algorithms suitable for an implementation on a computer. Examples, such as expansions of generalized hypergeometric functions or Appell functions are discussed. As a further application, we give the general solution of a two-loop integral, the so-called C-topology, in terms of multiple nested sums. In addition, we discuss some important properties of nested sums, in particular we show that they satisfy a Hopf algebra.
We extend the dipole formalism of Catani and Seymour to QCD processes involving heavy fermions. We give the appropriate subtraction terms together with their integrated counterpart. All calculations are done within dimensional regularization. Our formulae can be used with three variants of dimensional regularization (conventional dimensional regularization, the 't Hooft-Veltman scheme or the four-dimensional scheme).
We suggest a new approach for calculating heavy-to-light form factors. The method is based on light cone sum rules (LCSR) and covers the whole kinematical range of momentum transfer. The derivation of the new sum rule uses a suitable combination of double and single dispersion integrals. As an example we give numerical results for the form factor f^+ for the B -> pion transition.
We consider a U(1) x SU(2) gauge theory on the four-dimensional manifold S^1 x S^3. If we make the assumption that only gauge transformations connected to the identity are allowed, the winding numbers of U(1) around S^1 and of SU(2) around S^3 become topological conserved quantities. We derive the effective theory for non-trivial winding numbers if all distances are small compared to the radii of the spheres. In the non-abelian case the gauge bosons become massive.