by Lee Siu-keung. / Parallel title in Chinese. / Thesis (M.Phil.)--Chinese University of Hong Kong, 1992. / Includes bibliographical references (leaves 162-164). / Chapter Chapter 1 --- Introduction --- p.5 / Chapter Chapter 2 --- Path Integrals / Chapter 2.1 --- Path´ؤIntegral Approach to Quantum Mechanics --- p.8 / Chapter 2.2 --- Path´ؤIntegral Approach to Statistical Mechanics --- p.14 / Chapter 2.3 --- Variational Principle --- p.18 / Chapter 2.4 --- "Variational Method Proposed by Giachetti and Tognetti, and by Feynman and Kleinert" / Chapter 2.4.1 --- Effective Classical Partition Function --- p.24 / Chapter 2.4.2 --- Particle Distribution Function From Effective Classical Potential --- p.34 / Chapter Chapter 3 --- Systematic Perturbation Corrections to the Variational Approximation Proposed in Section2.4 / Chapter 3.1 --- Formalism / Chapter 3.1.1 --- Free Energy --- p.38 / Chapter 3.1.2 --- Particle Distribution Function --- p.49 / Chapter 3.2 --- Second Order Correction to Free Energy --- p.53 / Chapter 3.3 --- First Order Correction to Particle Distribution Function --- p.60 / Chapter Chapter 4 --- Examples and Results / Chapter 4.1 --- Quartic Anharmonic Oscillator / Chapter 4.1.1 --- "Free Energy, Internal Energy and Specific Heat" --- p.69 / Chapter 4.1.2 --- Particle Distribution Function --- p.87 / Chapter 4.2 --- Symmetric Double-well Potential / Chapter 4.2.1 --- "Free Energy, Internal Energy and Specific Heat" --- p.88 / Chapter 4.2.2 --- Particle Distribution Function --- p.106 / Chapter 4.3 --- Quartic-cubic Anharmonic Potential / Chapter 4.3.1 --- Free Energy --- p.108 / Chapter 4.3.2 --- Particle Distribution Function --- p.115 / Chapter Chapter 5 --- Application to the One-dimensional Ginzburg-Landau Model / Chapter 5.1 --- Introduction --- p.120 / Chapter 5.2 --- Exact Partition Function and Free Energy Per Unit Length --- p.123 / Chapter 5.3 --- Zeroth Order Approximation to Free Energy Per Unit Length --- p.126 / Chapter 5.4 --- Exact Specific Heat --- p.133 / Chapter 5.5 --- Zeroth Order Approximation to Specific Heat --- p.139 / Chapter Chapter 6 --- Conclusion --- p.141 / Chapter Appendix I --- Functional Calculus - Differentiation --- p.145 / Chapter Appendix II --- Evaluation of Feynman Propagator Δf(τ) --- p.147 / Chapter Appendix III --- Vanishing of the First Order Correction-βf1 --- p.150 / Chapter Appendix IV --- Numerical Method for the Energy Eigenvalues and Eigenfunctions of the One-dimensional Schroedinger Equation with ax2 + bx4 Potential --- p.153 / Chapter Appendix V --- Numerical Integrations with imaginary Ω --- p.158 / References --- p.162 / Figures --- p.165

The Feynman path integral formulation of quantum mechanics is a path integral representation for a propagator or probability amplitude in going between two points in space-time. The wave function is expressed in terms of an integral equation from which the Schrodinger equation can be derived. On taking the limit h — 0, the method of stationary phase can be applied and Newton's second law of motion is obtained. Also, the condition the phase vanishes leads to the Hamilton - Jacobi equation. The secondary objective of this paper is to study ways of relating quantum mechanics and classical mechanics. The Ehrenfest theorem is applied to a particle in an electromagnetic field. Expressions are found which are the hermitian Lorentz force operator, the hermitian torque operator, and the hermitian power operator.

Die Auswertung von Mehrschleifen-Feynman-Integralen ist eine der größten Herausforderungen bei der Berechnung präziser theoretischer Vorhersagen für die am LHC gemessenen Wirkungsquerschnitte. In den vergangenen Jahren hat sich die Nutzung von Differentialgleichungen bei der Berechnung von Feynman-Integralen als sehr erfolgreich erwiesen. Es wurde dabei beobachtet, dass die von den Feynman-Integralen erfüllte Differentialgleichung oftmals in eine sogenannte kanonische Form transformiert werden kann, welche die Integration der Differentialgleichung mittels iterierter Integrale wesentlich vereinfacht. Das zentrale Ergebnis der vorliegenden Arbeit ist ein Algorithmus zur Berechnung rationaler Transformationen von Differentialgleichungen von Feynman-Integralen in eine kanonische Form. Neben der Existenz einer solchen rationalen Transformation stellt der Algorithmus keinerlei weitere Bedingungen an die Differentialgleichung. Insbesondere ist der Algorithmus auf Mehrskalenprobleme anwendbar und erlaubt eine rationale Abhängigkeit der Differentialgleichung vom dimensionalen Regulator. Bei der Anwendung des Algorithmus wird zunächst das Transformationsgesetz im dimensionalen Regulator entwickelt, um Differentialgleichungen für die Koeffizienten in der Entwicklung der Transformation herzuleiten. Diese Differentialgleichungen werden dann mit einem rationalen Ansatz für die gesuchte Transformation gelöst. Es wird zudem eine Implementation des Algorithmus in dem Mathematica Paket CANONICA vorgestellt, welches das erste veröffentlichte Programm dieser Art ist, das auf Mehrskalenprobleme anwendbar ist. CANONICAs Potential für moderne Mehrschleifenrechnungen wird anhand mehrerer nicht trivialer Mehrschleifen-Integraltopologien demonstriert. Die gezeigten Topologien hängen von bis zu drei Variablen ab und umfassen auch vormals ungelöste Topologien, die zu Korrekturen höherer Ordnung zum Wirkungsquerschnitt der Produktion einzelner Top-Quarks am LHC beitragen. / The evaluation of multi-loop Feynman integrals is one of the main challenges in the computation of precise theoretical predictions for the cross sections measured at the LHC. In recent years, the method of differential equations has proven to be a powerful tool for the computation of Feynman integrals. It has been observed that the differential equation of Feynman integrals can in many instances be transformed into a so-called canonical form, which significantly simplifies its integration in terms of iterated integrals. The main result of this thesis is an algorithm to compute rational transformations of differential equations of Feynman integrals into a canonical form. Apart from requiring the existence of such a rational transformation, the algorithm needs no further assumptions about the differential equation. In particular, it is applicable to problems depending on multiple kinematic variables and also allows for a rational dependence on the dimensional regulator. First, the transformation law is expanded in the dimensional regulator to derive differential equations for the coefficients of the transformation. Using an ansatz in terms of rational functions, these differential equations are then solved to determine the transformation. This thesis also presents an implementation of the algorithm in the Mathematica package CANONICA, which is the first publicly available program to compute transformations to a canonical form for differential equations depending on multiple variables. The main functionality and its usage are illustrated with some simple examples. Furthermore, the package is applied to state-of-the-art integral topologies appearing in recent multi-loop calculations. These topologies depend on up to three variables and include previously unknown topologies contributing to higher-order corrections to the cross section of single top-quark production at the LHC.

In dieser Dissertation geht es um Schwinger-parametrische Feynmanintegrale in der Quantenelektrodynamik. Mittels einer Vielzahl von Methoden aus der Kombinatorik und Graphentheorie wird eine signifikante Vereinfachung des Integranden erreicht. Nach einer größtenteils in sich geschlossenen Einführung zu Feynmangraphen und -integralen wird die Herleitung der Schwinger-parametrischen Darstellung aus den klassischen Impulsraumintegralen ausführlich erläutert, sowohl für skalare Theorien als auch Quantenelektrodynamik. Es stellt sich heraus, dass die Ableitungen, die benötigt werden um Integrale aus der Quantenelektrodynamik in ihrer parametrischen Version zu formulieren, neue Graphpolynome enthalten, die auf Zykeln und minimalen Schnitten (engl. "bonds") basieren. Danach wird die Tensorstruktur der Quantenelektrodynamik, bestehend aus Dirac-Matrizen und ihren Spuren, durch eine diagrammatische Interpretation ihrer Kontraktion zu ganzzahligen Faktoren reduziert. Dabei werden insbesondere gefärbte Sehnendiagramme benutzt. Dies liefert einen parametrischen Integranden, der über bestimmte Teilmengen solcher Diagramme summierte Produkte von Zykel- und Bondpolynomen enthält. Weitere Untersuchungen der im Integranden auftauchenden Polynome decken Verbindungen zu Dodgson- und Spannwaldpolynomen auf. Dies wird benutzt um eine Identität zu beweisen, mit der sehr große Summen von Sehnendiagrammen in einer kurzen Form ausgedrückt werden können. Insbesondere führt dies zu Aufhebungen, die den Integranden massiv vereinfachen. / This thesis is concerned with the study of Schwinger parametric Feynman integrals in quantum electrodynamics. Using a variety of tools from combinatorics and graph theory, significant simplification of the integrand is achieved. After a largely self-contained introduction to Feynman graphs and integrals, the derivation of the Schwinger parametric representation from the standard momentum space integrals is reviewed in full detail for both scalar theories and quantum electrodynamics. The derivatives needed to express Feynman integrals in quantum electrodynamics in their parametric version are found to contain new types of graph polynomials based on cycle and bond subgraphs. Then the tensor structure of quantum electrodynamics, products of Dirac matrices and their traces, is reduced to integer factors with a diagrammatic interpretation of their contraction. Specifically, chord diagrams with a particular colouring are used. This results in a parametric integrand that contains sums of products of cycle and bond polynomials over certain subsets of such chord diagrams. Further study of the polynomials occurring in the integrand reveals connections to other well-known graph polynomials, the Dodgson and spanning forest polynomials. This is used to prove an identity that expresses some of the very large sums over chord diagrams in a very concise form. In particular, this leads to cancellations that massively simplify the integrand.

Wir geben eine rigorose Konstruktion von analytisch-regularisierten Feynman-Integralen im D-dimensionalen Minkowski-Raum als meromorphe Distributionen in den externen Impulsen, sowohl in der Impuls- als auch in der parametrischen Darstellung. Wir zeigen, dass ihre Pole durch die üblichen Power-counting Formeln gegeben sind, und dass ihr singulärer Träger in mikrolokalen Verallgemeinerungen der (+alpha)-Landauflächen enthalten ist. Als weitere Anwendungen geben wir eine Konstruktion von dimensional regularisierten Integralen im Minkowski-Raum und beweisen Diskontinuitätsformeln für parametrische Amplituden. / We give a rigorous construction of analytically regularized Feynman integrals in D-dimensional Minkowski space as meromorphic distributions in the external momenta, both in the momentum and parametric representation. We show that their pole structure is given by the usual power-counting formula and that their singular support is contained in a microlocal generalization of the alpha-Landau surfaces. As further applications, we give a construction of dimensionally regularized integrals in Minkowski space and prove discontinuity formula for parametric amplitudes.

Nous étudions des fluctuations qui sont omniprésentes dans des collisions entre particules aux hautes énergies. Ces fluctuations peuvent être de nature classique ou quantique et nous allons considérer ces deux cas. D'abord, nous étudions les fluctuations quantiques qui sont présentes dans des collisions entre protons. Celles-ci sont calculables en théorie quantique des champs, et nous allons nous concentrer sur une certaine classe de diagrammes dans ce cadre. Dans un second temps nous allons étudier des fluctuations qui sont présentes dans des collisions entre particules plus lourdes que le proton. Celles-ci sont décrites par les lois quantiques de la nature qui donnent les positions des nucléons dans le noyau, ou bien des fluctuations classiques, d'origine thermique, qui affectent l'évolution hydrodynamique du milieu produit dans une collision. Les fluctuations dans des collisions entre protons peuvent être calculées analytiquement jusqu'à un certain ordre en théorie quantique des champs. Nous allons nous concentrer sur des diagrammes à une boucle, d'une topologie donnée. Ces diagrammes aux boucles donnent des intégrales, qui typiquement sont difficiles à calculer. Nous allons démontrer comment des outils des mathématiques modernes peuvent être utilisés pour faciliter leur évaluation. En particulier, nous allons étudier des relations entre des coupures d'un diagramme, la discontinuité à travers d'un branchement et le coproduit. Nous allons démontrer comment l'intégrale originale peut être reconstruit à partir de l'information contenue dans le coproduit. Nous nous attendons à ce que ces méthodes seront utiles pour le calcul des diagrammes avec des topologies plus difficiles et ainsi aident au calcul des nouvelles amplitudes de diffusion. A la fin, nous étudions les deux types de fluctuations qui ont lieu dans des collisions entre ions lourds. Celles-ci sont liées soit à l'état initial de la matière, soit à l'état intermédiaire produit dans une telle collision. Les fluctuations de l'état initial ont été mesurées expérimentalement, et on voit qu'elles donnent lieu à des non-Gaussianités dans le spectre final de particules. Nous allons démontrer comment ces non-Gaussianités peuvent être comprises comme des positions et des énergies d'interaction aléatoires des 'sources' dans les noyaux entrant en collision. En plus, nous étudions le bruit hydrodynamique dans le milieu produit juste après une collision. Le comportement de ce milieu est celui d'un fluide à basse viscosité. / We study fluctuations that are omnipresent in high-energy particle collisions. These fluctuations can be either of either classical or quantum origin and we will study both. Firstly, we consider the type of quantum fluctuations that arise in proton-proton collisions. These are computable perturbatively in quantum field theory and we will focus on a specific class of diagrams in this set-up. Secondly, we will consider the fluctuations that are present in collisions between nuclei that can be heavier than protons. These are the quantum laws of nature that describe the positions of nucleons within a nucleus, but also the hydrodynamic fluctuations of classical, thermal origin that affect the evolution of the medium produced in heavy-ion collisions. The fluctuations arising in proton-proton collisions can be computed analytically up to a certain order in perturbative quantum field theory. We will focus on one-loop diagrams of a fixed topology. Loop diagrams give rise to integrals that typically are hard to evaluate. We show how modern mathematical methods can be used to ease their computation. We will study the relations among unitarity cuts of a diagram, the discontinuity across the corresponding branch cut and the coproduct. We show how the original integral corresponding to a given diagram can be reconstructed from the information contained in the coproduct. We expect that these methods can be applied to solve more complicated topologies and help in the computation of new amplitudes in the future. Finally, we study the two types of fluctuations arising in heavy-ion collisions. These are related either to the initial state or the intermediate state of matter produced in such collisions. The initial state fluctuations are experimentally observed to give rise to non-Gaussianities in the final-state spectra. We show how these non-Gaussianities can be explained by the random position and interaction energy of `sources' in the colliding nuclei. Furthermore, we investigate the effect of hydrodynamical noise in the evolution of the medium produced just after a collision. This medium behaves like a fluid with a very low viscosity, and so the corresponding evolution is hydrodynamical.

Les calculs de précision dans les théories de jauge jouent un rôle très important pour l’étude de la physique du Modèle Standard et au-delà dans les super-collisionneurs de particules comme le LHC, TeVatron et ILC. Par conséquent, il est extrêmement important de fournir des outils du calcul d’amplitudes à une boucle stables, rapides, efficaces et hautement automatisés. Cette thèse a pour but de développer la bibliothèque d’intégrales Golem95. Cette bibliothèque est un programme écrit en Fortran95, qui contient tous les ingrédients nécessaires pour calculer une intégrale scalaire ou tensorielle à une boucle avec jusqu’à six pattes externes. Golem95 utilise une méthode traditionnelle de réduction (réduction à la Golem) qui réduit les facteurs de forme en des intégrales de base redondantes qui peuvent être scalaires (sans paramètres de Feynman au numérateur) ou tensorielles (avec des paramètres de Feynman au numérateur); ce formalisme permet d’éviter les problèmes de l’instabilité numérique engendrés par des singularités factices dues à l’annulation des déterminants de Gram. En plus, cette bibliothèque peut être interfacée avec des programmes du calcul automatique basés sur les méthodes d’unitarité comme GoSam par exemple. Les versions antérieures de Golem95 ont été conçues pour le calcul des amplitudes sans masses internes. Le but de ce travail de thèse est de généraliser cette bibliothèque pour les configurations les plus générales (les masses complexes sont incluses), et de fournir un calcul numériquement stable dans les régions problématique en donnant une représentation intégrale unidimensionnelle stable pour chaque intégrale de base de Golem95. / Higher order corrections in gauge theories play a crucial role in studying physics within the standard model and beyond at TeV colliders, like LHC, TeVatron and ILC. Therefore, it is of extreme importance to provide tools for next-to-leading order amplitude computation which are fast, stable, efficient and highly automatized. This thesis aims at developing the library of integrals Golem95. This library is a program written in Fortran95, it contains all the necessary ingredients to calculate any one-loop scalar or tensorial integral with up to six external legs. Golem95 uses the traditional reduction method (Golem reduction) to reduce the form factors into redundant basic integrals, which can be scalar (without Feynman parameters in the numerator) or tensorial (with Feynman parameter in the numerator); this formalism allows us to avoid the problems of numerical instabilities generated by the spurious singularities induced by the vanishing of the Gram determinants. In addition, this library can be interfaced with automatic programs of NLO calculation based on the unitarity inspired reduction methods as GoSam for example. Earlierversions of Golem95 were designed for the calculation of amplitudes without internal masses. The purpose of this thesis is to extend this library for more general configurations (complex masses are supported); and to provide numerically stable calculation in the problematic regions (det(G) → 0), by providing a stable one-dimensional integral representation for each Golem95 basic integral.

In dieser Dissertation erweitern wir die Ideen des Yangian-Bootstrap-Algorithmus auf Feynman-Diagramme mit massiven Teilchen. Ausgehend von der massiven dual-konformen Symmetrie der N = 4 Super-Yang-Mills Theorie auf dem Coulomb-Zweig konstruieren wir einen Satz von bilokalen Yangian Level-Eins Generatoren und zeigen, dass sie eine unendliche Anzahl von planaren ein- und zwei-Schleifen-Diagrammen vernichten. Wir beschreiben außerdem wie der dual-konforme Level-Eins Impuls-Operator auf eine massive Verallgemeinerung des gewöhnlichen spezial-konformen Generators im Impulsraum abgebildet wird. Als nächstes wenden wir den Yangian-Bootstrap-Algorithmus mit großem Erfolg auf eine Reihe von massiven Ein-Schleifen-Diagrammen mit verallgemeinerten Propagatorexponenten und in beliebiger Anzahl von Raumdimensionen an. Im Spezialfall der dual-konformen Integrale, deren Propagatorexponenten sich zur Raumdimension addieren, finden wir neue sehr einfache Darstellungen durch hypergeometrische Funktionen, die eine natürliche Verallgemeinerung für Diagramme mit beliebig vielen äußeren Punkten erlauben. Außerdem diskutieren wir Aspekte des Yangian-Bootstrap-Algorithmus in Minkowski-Raumzeit am Beispiel des masselosen Box-Integrals. Wir zeigen, dass dessen Yangian-Symmetrie gemeinsam mit seinen diskreten Permutationssymmetrien das Box-Integrals bis auf 12 unbestimmte Konstanten komplett festlegt. Schließlich schlagen wir vor, dass das Auftreten von Yangian-Symmetrie in massiven Fischnetz-Diagrammen mit deren Rolle als Ein-Spur-Streuamplituden in einer massiven Fischnetz-Theorie zusammenhängen könnte. In Analogie mit der masselosen Fischnetz-Theorie zeigen wir, wie diese Theorie als Deformation der N = 4 Super-Yang-Mills Theorie auf dem Coulomb-Zweig definiert werden kann. Wir diskutieren eine bestimmte Klasse von planaren Grenzfällen, in der die off-shell Streuamplituden der Theorie eine massive dual-konforme Symmetrie sowie Yangian-Symmetrie aufweisen. / In this dissertation, we extend the ideas of the Yangian bootstrap algorithm to massive Feynman diagrams. Based on the massive dual-conformal symmetry of Coulomb branch N = 4 super-Yang-Mills theory, we construct a set of bi-local Yangian level-one generators and show that they annihilate infinite classes of massive planar Feynman integrals at one and two loops. We also describe how the dual-conformal level-one momentum generator maps to a massive deformation of the ordinary momentum space special conformal generator. We then apply the Yangian bootstrap to a set of massive one-loop integrals with generalised propagator powers and in an arbitrary number of space dimensions to great success. In the special case of dual-conformal integrals, whose propagator powers sum to the space dimension, we find very simple novel hypergeometric structures, suggesting a natural generalisation to diagrams with an arbitrary number of external points. In the particular case of the massless box integral we also discuss elements of the Yangian bootstrap in Minkowski space. We show that its Yangian and discrete permutation symmetries constrain it up to 12 undetermined constants. We then derive the values of these constants via analytic continuation from the box integral in the Euclidean region. Finally, we provide evidence that the appearance of Yangian symmetry for massive fishnet diagrams is related to their role as colour-ordered scattering amplitudes in a massive fishnet theory. We show how to construct this theory from Coulomb branch N = 4 super-Yang-Mills theory, paralleling the original construction of the massless fishnet theory. We discuss how a particular class of planar limits leads to the emergence of massive dual-conformal symmetry as well as massive Yangian symmetry for the theory’s off-shell scattering amplitudes.

In dieser Dissertation untersuchen wir unterschiedliche Aspekte im Zusammenhang mit Korrelationsfunktionen in der Fischnetz-Theorie. Zunächst betrachten wir einen der einfachsten Korrelatoren der Fischnetz Theorie, das konforme Box-Integral, in Minkowski Signatur. Während dieses Integral in Euklidischer Signatur eine konforme Symmetrie aufweist, wird diese Symmetrie in Minkowski-Raumzeit subtil gebrochen. Wir beschreiben die Brechung der konformen Symmetrie quantitativ, indem wir die funktionale Form des Box-Integrals in allen kinematischen Regionen untersuchen. Ausserdem untersuchen wir das Ausmass zu dem das Box integral durch seine Yangian-Symmetrie festgelegt ist. Als nächstes widmen wir uns den Basso-Dixon-Graphen, die ebenfalls konforme Vier-Punkt-Integrale sind und Verallgemeinerungen des Box-Integrals zu höheren Schleifenordnungen darstellen. Wir leiten die Yangian-Ward-Identitäten ab, die diese Klasse von Integralen erfüllen. Die Ward-Identitäten sind einhomogene Erweiterungen der partiellen Differentialgleichungen, die im homogenen Fall durch Appell-Hypergeometrische Funktionen gelöst werden. Die Ward-Identitäten können natürlicherweise auf eine Ein-Parameter-Familie von D-dimensionalen Integralen erweitert werden, die Korrelatoren in der verallgemeinerten Fischnetz-Theorie von Kazakov und Olivucci darstellen. Schliesslich untersuchen wir den Dilatationsoperator in einem Drei-Skalar-Sektor der Fischnetztheorie, der auch als Eklektisches Modell bezeichnet wird. In diesem Sektor der Dilatationsoperator nimmt nicht--diagonalisierbare Form an. Das führt dazu, dass die Zwei-Punkt-Korrelationsfunktionen eine logarithmische Abhängigkeit von der Raumzeitseparierung der Operatoren annimmt. Unter Zuhilfenahme von kombinatorischen Argumenten führen wir eine generierende Funktion ein, die das Jordan-Block-Spektrum eines verwandten Modells, der hypereklektischen Spinkette, vollständig charakterisiert. / We study various aspects of correlation functions in fishnet theory. We begin with the study of the simplest correlator in theory theory, represented by the conformal box integral, in Minkowski space. While this integral is conformally invariant in Euclidean space, this symmetry is subtly broken in Minkowski space. We quantify the extent to which conformal symmetry is broken by analysing the functional form of the box in each kinematic region. We propose a new method to calculate the box integral directly in Minkowski space, by introducing a family of configurations with two points at infinity. Furthermore, we investigate the extent to which the box integral is constrained by Yangian symmetry. We constrain the functional form of the box integral in all kinematic regions up to twelve undetermined constants, which we fix by three separate analytic continuations from the Euclidean region. Next, we study the Basso-Dixon graphs, which represent higher-loop versions of the box integral. We derive and study Yangian Ward identities for this class of integrals. These take the form of inhomogeneous extensions of the partial differential equations defining the Appell hypergeometric functions. The Ward identities naturally generalise to a one-parameter family of D dimensional integrals representing correlators in a generalised fishnet theory. Finally, we study the dilatation operator in a particular three scalar sector of the fishnet theory, which has been dubbed the eclectic model. This dilatation operator is non-diagonalisable in this sector. This leads to logarithmic spacetime dependence in the corresponding two-point functions. Using combinatorial arguments, we introduce a generating function which fully characterises the Jordan block spectrum of a related model: the hypereclectic spin chain. This function is found by purely combinatorial means and can be expressed in terms of the q-binomial coefficient.

Wir untersuchen Feynman-Integrale in der Darstellung mit Schwinger-Parametern und leiten rekursive Integralgleichungen für masselose 3- und 4-Punkt-Funktionen her. Eigenschaften der analytischen (und dimensionalen) Regularisierung werden zusammengefasst und wir beweisen, dass in der Euklidischen Region jedes Feynman-Integral als eine Linearkombination konvergenter Feynman-Integrale geschrieben werden kann. Dies impliziert, dass man stets eine Basis aus konvergenten Masterintegralen wählen kann und somit divergente Integrale nicht selbst berechnet werden müssen. Weiterhin geben wir eine in sich geschlossene Darstellung der Theorie der Hyperlogarithmen und erklären detailliert die nötigen Algorithmen, um diese für die Berechnung mehrfacher Integrale anzuwenden. Wir definieren eine neue Methode um die Singularitäten solcher Integrale zu bestimmen und stellen ein Computerprogramm vor, welches die Integrationsalgorithmen implementiert. Unser Hauptresultat ist die Konstruktion unendlicher Familien masseloser 3- und 4-Punkt-Funktionen (diese umfassen unter anderem alle Leiter-Box-Graphen und deren Minoren), deren Feynman-Integrale zu allen Ordnungen in der epsilon-Entwicklung durch multiple Polylogarithmen dargestellt werden können. Diese Integrale können mit dem vorgestellten Programm explizit berechnet werden. Die Arbeit enthält interessante Beispiele von expliziten Ergebnissen für Feynman-Integrale mit bis zu 6 Schleifen. Insbesondere präsentieren wir den ersten exakt bestimmten Gegenterm in masseloser phi^4-Theorie, der kein multipler Zetawert ist sondern eine Linearkombination multipler Polylogarithmen, ausgewertet an primitiven sechsten Einheitswurzeln (und geteilt durch die Quadratwurzel aus 3). Zu diesem Zweck beweisen wir ein Paritätsresultat über die Zerlegbarkeit der Real- und Imaginärteile solcher Zahlen in Produkte und Beiträge geringerer Tiefe (depth). / We study Feynman integrals in the representation with Schwinger parameters and derive recursive integral formulas for massless 3- and 4-point functions. Properties of analytic (including dimensional) regularization are summarized and we prove that in the Euclidean region, each Feynman integral can be written as a linear combination of convergent Feynman integrals. This means that one can choose a basis of convergent master integrals and need not evaluate any divergent Feynman graph directly. Secondly we give a self-contained account of hyperlogarithms and explain in detail the algorithms needed for their application to the evaluation of multivariate integrals. We define a new method to track singularities of such integrals and present a computer program that implements the integration method. As our main result, we prove the existence of infinite families of massless 3- and 4-point graphs (including the ladder box graphs with arbitrary loop number and their minors) whose Feynman integrals can be expressed in terms of multiple polylogarithms, to all orders in the epsilon-expansion. These integrals can be computed effectively with the presented program. We include interesting examples of explicit results for Feynman integrals with up to 6 loops. In particular we present the first exactly computed counterterm in massless phi^4 theory which is not a multiple zeta value, but a linear combination of multiple polylogarithms at primitive sixth roots of unity (and divided by the square-root of 3). To this end we derive a parity result on the reducibility of the real- and imaginary parts of such numbers into products and terms of lower depth.