Black hole scattering from Worldline Quantum Field Theory and the Classical Double Copy of Spinning particles

This PhD thesis focuses on studying the classical scattering of massive/massless particles toward black holes, and investigating double copy relations between classical observables in gauge theories and gravity. This is done in the Post-Minkowskian approximation i.e. a perturbative expansion of observables controlled by the gravitational coupling constant κ = 32πGN, with GN being the Newtonian coupling constant. The investigation is performed by using the Worldline Quantum Field Theory (WQFT), displaying a worldline path integral describing the scattering objects and a QFT path integral in the Born approximation, describing the intermediate bosons exchanged in the scattering event by the massive/massless particles. We introduce the WQFT, by deriving a relation between the Kosower- Maybee-O’Connell (KMOC) limit of amplitudes and worldline path integrals, then, we use that to study the classical Compton amplitude and higher point amplitudes. We also present a nice application of our formulation to the case of Hard Thermal Loops (HTL), by explicitly evaluating hard thermal currents in gauge theory and gravity. Next we move to the investigation of the classical double copy (CDC), which is a powerful tool to generate integrands for classical observables related to the binary inspiralling problem in General Relativity. In order to use a Bern-Carrasco-Johansson (BCJ) like prescription, straight at the classical level, one has to identify a double copy (DC) kernel, encoding the locality structure of the classical amplitude. Such kernel is evaluated by using a theory where scalar particles interacts through bi-adjoint scalars. We show here how to push forward the classical double copy so to account for spinning particles, in the framework of the WQFT. Here the quantization procedure on the worldline allows us to fully reconstruct the quantum theory on the gravitational side. Next we investigate how to describe the scattering of massless particles off black holes in the WQFT.

On the path integral formulation and the evaluation of quantum statistical averages

Four problems of physical interest have been solved in this thesis using the path integral formalism. Using the trigonometric expansion method of Burton and de Borde (1955), we found the kernel for two interacting one dimensional oscillators• The result is the same as one would obtain using a normal coordinate transformation, We next introduced the method of Papadopolous (1969), which is a systematic perturbation type method specifically geared to finding the partition function Z, or equivalently, the Helmholtz free energy F, of a system of interacting oscillators. We applied this method to the next three problems considered• First, by summing the perturbation expansion, we found F for a system of N interacting Einstein oscillators^ The result obtained is the same as the usual result obtained by Shukla and Muller (1972) • Next, we found F to 0(Xi)f where A is the usual Tan Hove ordering parameter* The results obtained are the same as those of Shukla and Oowley (1971), who have used a diagrammatic procedure, and did the necessary sums in Fourier space* We performed the work in temperature space• Finally, slightly modifying the method of Papadopolous, we found the finite temperature expressions for the Debyecaller factor in Bravais lattices, to 0(AZ) and u(/K/ j,where K is the scattering vector* The high temperature limit of the expressions obtained here, are in complete agreement with the classical results of Maradudin and Flinn (1963) .

One-loop n-point equivalence among negative-dimensional, Mellin-Barnes and Feynman parametrization approaches to Feynman integrals

We show that at one-loop order, negative-dimensional, Mellin-Barnes (MB) and Feynman parametrization (FP) approaches to Feynman loop integral calculations are equivalent. Starting with a generating functional, for two and then for n-point scalar integrals, we show how to reobtain MB results, using negative-dimensional and FP techniques. The n-point result is valid for different masses, arbitrary exponents of propagators and dimension.

A systematization for one-loop 4D Feynman integrals

We present a strategy for the systematization of manipulations and calculations involving divergent (or not) Feynman integrals, typical of the one-loop perturbative solutions of QFT, where the use of an explicit regularization is avoided. Two types of systematization are adopted. The divergent parts are put in terms of a small number of standard objects, and a set of structure functions for the finite parts is also defined. Some important properties of the finite structures, specially useful in the verification of relations among Green's functions, are identified. We show that, in fundamental (renormalizable) theories, all the finite parts of two-, three- and four-point functions can be written in terms of only three basic functions while the divergent parts require (only) five objects. The final results obtained within the proposed strategy can be easily converted into those corresponding to any specific regularization technique providing an unified point of view for the treatment of divergent Feynman integrals. Examples of physical amplitudes evaluation and their corresponding symmetry relations verification are presented as well as generalizations of our results for the treatment of Green's functions having an arbitrary number of points are considered.

Non-planar double-box, massive and massless pentabox Feynman integrals in the negative-dimensional approach

The negative-dimensional integration method is a technique which can be applied, with success, in usual covariant gauge calculations. We consider three two-loop diagrams: the scalar massless non-planar double-box with six propagators and the scalar pentabox in two cases, where six virtual particles have the same mass, and in the case all of them are massless. Our results are given in terms of hypergeometric functions of Mandelstam variables and also for arbitrary exponents of propagators and dimension D.

NDIM achievements: massive, arbitrary tensor rank and N-loop insertions in Feynman integrals

One of the main difficulties in studying quantum field theory, in the perturbative regime, is the calculation of D-dimensional Feynman integrals. In general, one introduces the so-called Feynman parameters and, associated with them, the cumbersome parametric integrals. Solving these integrals beyond the one-loop level can be a difficult task. The negative-dimensional integration method (NDIM) is a technique whereby such a problem is dramatically reduced. We present the calculation of two-loop integrals in three different cases: scalar ones with three different masses, massless with arbitrary tensor rank, with and N insertions of a two-loop diagram.

Feynman integrals with tensorial structure in the negative dimensional integration scheme

The negative-dimensional integration method (NDIM) is revealing itself as a very useful technique for computing massless and/or massive Feynman integrals, covariant and noncovanant alike. Up until now however, the illustrative calculations done using such method have been mostly covariant scalar integrals/without numerator factors. We show here how those integrals with tensorial structures also can be handled straightforwardly and easily. However, contrary to the absence of significant features in the usual approach, here the NDIM also allows us to come across surprising unsuspected bonuses. Toward this end, we present two alternative ways of working out the integrals and illustrate them by taking the easiest Feynman integrals in this category that emerge in the computation of a standard one-loop self-energy diagram. One of the novel and heretofore unsuspected bonuses is that there are degeneracies in the way one can express the final result for the referred Feynman integral.

Loop Integrals in Three Outstanding Gauges: Feynman, Light-Cone, and Coulomb

We apply the negative dimensional integration method (NDIM) to three outstanding gauges: Feynman, light-cone, and Coulomb gauges. Our aim is to show that NDIM is a very suitable technique to deal with loop integrals, regardless of which gauge choice that originated them. In the Feynman gauge we perform scalar two-loop four-point massless integrals; in the light-cone gauge we calculate scalar two-loop integrals contributing to two-point functions without any kind of prescriptions, since NDIM can abandon such devices - this calculation is the first test of our prescriptionless method beyond one-loop order; and finally, for the Coulomb gauge we consider a four-propagator massless loop integral, in the split-dimensional regularization context. © 2001 Academic Press.

Mathematical aspects of Feynman integrals

In the present dissertation we consider Feynman integrals in the framework of dimensional regularization. As all such integrals can be expressed in terms of scalar integrals, we focus on this latter kind of integrals in their Feynman parametric representation and study their mathematical properties, partially applying graph theory, algebraic geometry and number theory. The three main topics are the graph theoretic properties of the Symanzik polynomials, the termination of the sector decomposition algorithm of Binoth and Heinrich and the arithmetic nature of the Laurent coefficients of Feynman integrals.rnrnThe integrand of an arbitrary dimensionally regularised, scalar Feynman integral can be expressed in terms of the two well-known Symanzik polynomials. We give a detailed review on the graph theoretic properties of these polynomials. Due to the matrix-tree-theorem the first of these polynomials can be constructed from the determinant of a minor of the generic Laplacian matrix of a graph. By use of a generalization of this theorem, the all-minors-matrix-tree theorem, we derive a new relation which furthermore relates the second Symanzik polynomial to the Laplacian matrix of a graph.rnrnStarting from the Feynman parametric parameterization, the sector decomposition algorithm of Binoth and Heinrich serves for the numerical evaluation of the Laurent coefficients of an arbitrary Feynman integral in the Euclidean momentum region. This widely used algorithm contains an iterated step, consisting of an appropriate decomposition of the domain of integration and the deformation of the resulting pieces. This procedure leads to a disentanglement of the overlapping singularities of the integral. By giving a counter-example we exhibit the problem, that this iterative step of the algorithm does not terminate for every possible case. We solve this problem by presenting an appropriate extension of the algorithm, which is guaranteed to terminate. This is achieved by mapping the iterative step to an abstract combinatorial problem, known as Hironaka's polyhedra game. We present a publicly available implementation of the improved algorithm. Furthermore we explain the relationship of the sector decomposition method with the resolution of singularities of a variety, given by a sequence of blow-ups, in algebraic geometry.rnrnMotivated by the connection between Feynman integrals and topics of algebraic geometry we consider the set of periods as defined by Kontsevich and Zagier. This special set of numbers contains the set of multiple zeta values and certain values of polylogarithms, which in turn are known to be present in results for Laurent coefficients of certain dimensionally regularized Feynman integrals. By use of the extended sector decomposition algorithm we prove a theorem which implies, that the Laurent coefficients of an arbitrary Feynman integral are periods if the masses and kinematical invariants take values in the Euclidean momentum region. The statement is formulated for an even more general class of integrals, allowing for an arbitrary number of polynomials in the integrand.

Picard-Fuchs equations of dimensionally regulated Feynman integrals

Zusammenfassung In der vorliegenden Arbeit besch¨aftige ich mich mit Differentialgleichungen von Feynman– Integralen. Ein Feynman–Integral h¨angt von einem Dimensionsparameter D ab und kann f¨ur ganzzahlige Dimension als projektives Integral dargestellt werden. Dies ist die sogenannte Feynman–Parameter Darstellung. In Abh¨angigkeit der Dimension kann ein solches Integral divergieren. Als Funktion in D erh¨alt man eine meromorphe Funktion auf ganz C. Ein divergentes Integral kann also durch eine Laurent–Reihe ersetzt werden und dessen Koeffizienten r¨ucken in das Zentrum des Interesses. Diese Vorgehensweise wird als dimensionale Regularisierung bezeichnet. Alle Terme einer solchen Laurent–Reihe eines Feynman–Integrals sind Perioden im Sinne von Kontsevich und Zagier. Ich beschreibe eine neue Methode zur Berechnung von Differentialgleichungen von Feynman– Integralen. ¨ Ublicherweise verwendet man hierzu die sogenannten ”integration by parts” (IBP)– Identit¨aten. Die neue Methode verwendet die Theorie der Picard–Fuchs–Differentialgleichungen. Im Falle projektiver oder quasi–projektiver Variet¨aten basiert die Berechnung einer solchen Differentialgleichung auf der sogenannten Griffiths–Dwork–Reduktion. Zun¨achst beschreibe ich die Methode f¨ur feste, ganzzahlige Dimension. Nach geeigneter Verschiebung der Dimension erh¨alt man direkt eine Periode und somit eine Picard–Fuchs–Differentialgleichung. Diese ist inhomogen, da das Integrationsgebiet einen Rand besitzt und daher nur einen relativen Zykel darstellt. Mit Hilfe von dimensionalen Rekurrenzrelationen, die auf Tarasov zur¨uckgehen, kann in einem zweiten Schritt die L¨osung in der urspr¨unglichen Dimension bestimmt werden. Ich beschreibe außerdem eine Methode, die auf der Griffiths–Dwork–Reduktion basiert, um die Differentialgleichung direkt f¨ur beliebige Dimension zu berechnen. Diese Methode ist allgemein g¨ultig und erspart Dimensionswechsel. Ein Erfolg der Methode h¨angt von der M¨oglichkeit ab, große Systeme von linearen Gleichungen zu l¨osen. Ich gebe Beispiele von Integralen von Graphen mit zwei und drei Schleifen. Tarasov gibt eine Basis von Integralen an, die Graphen mit zwei Schleifen und zwei externen Kanten bestimmen. Ich bestimme Differentialgleichungen der Integrale dieser Basis. Als wichtigstes Beispiel berechne ich die Differentialgleichung des sogenannten Sunrise–Graphen mit zwei Schleifen im allgemeinen Fall beliebiger Massen. Diese ist f¨ur spezielle Werte von D eine inhomogene Picard–Fuchs–Gleichung einer Familie elliptischer Kurven. Der Sunrise–Graph ist besonders interessant, weil eine analytische L¨osung erst mit dieser Methode gefunden werden konnte, und weil dies der einfachste Graph ist, dessen Master–Integrale nicht durch Polylogarithmen gegeben sind. Ich gebe außerdem ein Beispiel eines Graphen mit drei Schleifen. Hier taucht die Picard–Fuchs–Gleichung einer Familie von K3–Fl¨achen auf.

Rational algorithms for the decomposition of Feynman integrals via intersection theory

The decomposition of Feynman integrals into a basis of independent master integrals is an essential ingredient of high-precision theoretical predictions, that often represents a major bottleneck when processes with a high number of loops and legs are involved. In this thesis we present a new algorithm for the decomposition of Feynman integrals into master integrals with the formalism of intersection theory. Intersection theory is a novel approach that allows to decompose Feynman integrals into master integrals via projections, based on a scalar product between Feynman integrals called intersection number. We propose a new purely rational algorithm for the calculation of intersection numbers of differential $n-$forms that avoids the presence of algebraic extensions. We show how expansions around non-rational poles, which are a bottleneck of existing algorithms for intersection numbers, can be avoided by performing an expansion in series around a rational polynomial irreducible over $\mathbb{Q}$, that we refer to as $p(z)-$adic expansion. The algorithm we developed has been implemented and tested on several diagrams, both at one and two loops.

Negative dimensional approach for scalar two-loop three-point and three-loop two-point integrals

The well-known D-dimensional Feynman integrals were shown, by Halliday and Ricotta, to be capable of undergoing analytic continuation into the domain of negative values for the dimension of space-time. Furthermore, this could be identified with Grassmannian integration in positive dimensions. From this possibility follows the concept of negative-dimensional integration for loop integrals in field theories. Using this technique, we evaluate three two-loop three-point scalar integrals, with five and six massless propagators, with specific external kinematic configurations (two legs on-shell), and four three-loop two-point scalar integrals. These results are given for arbitrary exponents of propagators and dimension, in Euclidean space, and the particular cases compared to results published in the literature.

Probing negative-dimensional integration: Two-loop covariant vertex and one-loop light-cone integrals

The negative-dimensional integration method (NDIM) seems to be a very promising technique for evaluating massless and/or massive Feynman diagrams. It is unique in the sense that the method gives solutions in different regions of external momenta simultaneously. Moreover, it is a technique whereby the difficulties associated with performing parametric integrals in the standard approach are transferred to a simpler solving of a system of linear algebraic equations, thanks to the polynomial character of the relevant integrands. We employ this method to evaluate a scalar integral for a massless two-loop three-point vertex with all the external legs off-shell, and consider several special cases for it, yielding results, even for distinct simpler diagrams. We also consider the possibility of NDIM in non-covariant gauges such as the light-cone gauge and do some illustrative calculations, showing that for one-degree violation of covariance (i.e. one external, gauge-breaking, light-like vector n μ) the ensuing results are concordant with the ones obtained via either the usual dimensional regularization technique, or the use of the principal value prescription for the gauge-dependent pole, while for two-degree violation of covariance - i.e. two external, light-like vectors n μ, the gauge-breaking one, and (its dual) n * μ - the ensuing results are concordant with the ones obtained via causal constraints or the use of the so-called generalized Mandelstam-Leibbrandt prescription. © 1999 Elsevier Science B.V.

Prescriptionless light-cone integrals

Perturbative quantum gauge field theory as seen within the perspective of physical gauge choices such as the light-cone gauge entails the emergence of troublesome poles of the type (k · n)-α in the Feynman integrals. These come from the boson field propagator, where α = 1, 2, ⋯ and nμ is the external arbitrary four-vector that defines the gauge proper. This becomes an additional hurdle in the computation of Feynman diagrams, since any graph containing internal boson lines will inevitably produce integrands with denominators bearing the characteristic gauge-fixing factor. How one deals with them has been the subject of research over decades, and several prescriptions have been suggested and tried in the course of time, with failures and successes. However, a more recent development at this fronteer which applies the negative dimensional technique to compute light-cone Feynman integrals shows that we can altogether dispense with prescriptions to perform the calculations. An additional bonus comes to us attached to this new technique, in that not only it renders the light-cone prescriptionless but, by the very nature of it, it can also dispense with decomposition formulas or partial fractioning tricks used in the standard approach to separate pole products of the type (k · n)-α[(k - p) · n]-β (β = 1, 2, ⋯). In this work we demonstrate how all this can be done.