Use of computer algebra in the calculation of Feynman diagrams

The increasing precision of current and future experiments in high-energy physics requires a likewise increase in the accuracy of the calculation of theoretical predictions, in order to find evidence for possible deviations of the generally accepted Standard Model of elementary particles and interactions. Calculating the experimentally measurable cross sections of scattering and decay processes to a higher accuracy directly translates into including higher order radiative corrections in the calculation. The large number of particles and interactions in the full Standard Model results in an exponentially growing number of Feynman diagrams contributing to any given process in higher orders. Additionally, the appearance of multiple independent mass scales makes even the calculation of single diagrams non-trivial. For over two decades now, the only way to cope with these issues has been to rely on the assistance of computers. The aim of the xloops project is to provide the necessary tools to automate the calculation procedures as far as possible, including the generation of the contributing diagrams and the evaluation of the resulting Feynman integrals. The latter is based on the techniques developed in Mainz for solving one- and two-loop diagrams in a general and systematic way using parallel/orthogonal space methods. These techniques involve a considerable amount of symbolic computations. During the development of xloops it was found that conventional computer algebra systems were not a suitable implementation environment. For this reason, a new system called GiNaC has been created, which allows the development of large-scale symbolic applications in an object-oriented fashion within the C++ programming language. This system, which is now also in use for other projects besides xloops, is the main focus of this thesis. The implementation of GiNaC as a C++ library sets it apart from other algebraic systems. Our results prove that a highly efficient symbolic manipulator can be designed in an object-oriented way, and that having a very fine granularity of objects is also feasible. The xloops-related parts of this work consist of a new implementation, based on GiNaC, of functions for calculating one-loop Feynman integrals that already existed in the original xloops program, as well as the addition of supplementary modules belonging to the interface between the library of integral functions and the diagram generator.

The massless two-loop two-point function and zeta functions in counterterms of Feynman diagrams

The main part of this thesis describes a method of calculating the massless two-loop two-point function which allows expanding the integral up to an arbitrary order in the dimensional regularization parameter epsilon by rewriting it as a double Mellin-Barnes integral. Closing the contour and collecting the residues then transforms this integral into a form that enables us to utilize S. Weinzierl's computer library nestedsums. We could show that multiple zeta values and rational numbers are sufficient for expanding the massless two-loop two-point function to all orders in epsilon. We then use the Hopf algebra of Feynman diagrams and its antipode, to investigate the appearance of Riemann's zeta function in counterterms of Feynman diagrams in massless Yukawa theory and massless QED. The class of Feynman diagrams we consider consists of graphs built from primitive one-loop diagrams and the non-planar vertex correction, where the vertex corrections only depend on one external momentum. We showed the absence of powers of pi in the counterterms of the non-planar vertex correction and diagrams built by shuffling it with the one-loop vertex correction. We also found the invariance of some coefficients of zeta functions under a change of momentum flow through these vertex corrections.

First principle study of ions in water and ionic liquids

Die beiden in dieser Arbeit betrachteten Systeme, wässrige Lösungen von Ionen und ionische Flüssigkeiten, zeigen vielfältige Eigenschaften und Anwendungsmöglichkeiten, im Gegensatz zu anderen Systemen. Man findet sie beinahe überall im normalen Leben (Wasser), oder ihre Bedeutung wächst (ioinische Flüssigkeiten). Der elektronische Anteil und der atomare Anteil wurden getrennt voneinander untersucht und im Zusammenhang analysiert. Mittels dieser Methode konnten die in dem jeweiligen System auftretenden Mechanismen genauer untersucht werden. Diese Methode wird "Multiscale Modeling" genannt, dabei werden die Untereinheiten eines Systems genauer betrachtet, wie in diesem Fall die elektronischen and atomaren Teilsystem. Die Ergebnisse, die aus den jeweiligen Betrachtungen hervorgehen, zeigen, dass, im Falle von hydratisierten Ionen die Wasser-Wasser Wechselwirkungen wesentlich stärker sind als die elektrostatischen Wechselwirkung zwischen Wasser und dem Ion. Anhand der Ergebnisse ergibt sich, dass normale nicht-polarisierbare Modelle ausreichen, um Ionen-Wasser Lösungen zu beschreiben. Im Falle der ionischen Flüssigkeiten betrachten wir die elektronische Ebene mittels sehr genauer post-Hartree-Fock Methoden und DFT, deren Ergebnisse dann mit denen auf molekularer Ebene (mithilfe von CPMD/klassischer MD) in Beziehung gesetzt werden. Die bisherigen Ergebnisse zeigen, dass die Wasserstoff-Brückenbindungen im Fall der ionischen Flüssigkeiten nicht vernachässigt werden können. Weiterhin hat diese Studie herausgefunden, dass die klassischen Kraftfelder die Elektrostatik (Dipol- und Quadrupolmomente) nicht genau genug beschreibt. Die Kombination des mikroskopischen Mechanismus und der molekularen Eigenschaften ist besonders sinnvoll um verschiedene Anhaltspunkte von Simualtionen (z.B. mit klassische Molekular-Dynamik) oder Experimenten zu liefern oder solche zu erklären.

Anhydrous proton conducting polymer electrolytes based on polymeric ionic liquids

Imidazolium types of ionic liquids were immobilized by tethering it to acrylate backbone. These imidazolium salt containing acrylate monomers were polymerize at 70oC by free radical polymerization to give polymers poly(AcIm-n) with n being the side chain lenght. The chemical structure of the polymer electrolytes obtained by the described synthetic routes was investigated by NMR-spectroscopy. The polymers were doped with various amounts of H3PO4 and LiN(SO2CF3)2, to obtain poly(AcIm-n) x H3PO4 and poly(AcIm-2-Li) x LiN(SO2CF3)2. The TG curves show that the polymer electrolytes are thermally stable up to about 200◦C. DSC results indicates the softening effect of the length of the spacers (n) as well as phosphoric acid. The proton conductivity of the samples increase with x and reaches to 10-2 Scm-1 at 120oC for both poly(AcIm-2)2H3PO4 and poly(AcIm-6)2H3PO4. It was observed that the lithium ion conductivity of the poly(AcIm-2-Li) x LiN(SO2CF3)2 increases with blends (x) up to certain composition and then leveled off independently from blend content. The conductivity reaches to about 10-5 S cm-1 at 30oC and 10-3 at 100oC for poly(AcIm-2-Li) x LiN(SO2CF3)2 where x is 10. The phosphate and phosphoric acid functionality in the resulting polymers, poly(AcIm-n) x H3PO4, undergoes condensation leading to the formation of cross-linked materials at elevated temperature which may improve the mechanical properties to be used as membrane materials in fuel cells. High resolution nuclear magnetic resonance (NMR) spectroscopy was used to obtain information about hydrogen bonding in solids. The low Tg enhances molecular mobility and this leads to better resolved resonances in both the backbone region and side chain region. The mobile and immobile protons can be distinguished by comparing 1H MAS and 1H-DQF NMR spectra. The interaction of the protons which may contribute to the conductivity is observed from the 2D double quantum correlation (DQC) spectra.

Numerische und analytische Untersuchungen stark korrelierter fermionischer Mehrbandsysteme

In dieser Arbeit werden vier unterschiedliche, stark korrelierte, fermionische Mehrbandsysteme untersucht. Es handelt sich dabei um ein Mehrstörstellen-Anderson-Modell, zwei Hubbard-Modelle sowie ein Mehrbandsystem, wie es sich aus einer ab initio-Beschreibung für ein korreliertes Halbmetall ergibt.rnrnDie Betrachtung des Mehrstörstellen-Anderson-Modells konzentriert sich auf die Untersuchung des Einflusses der Austauschwechselwirkung und der nicht-lokalen Korrelationen zwischen zwei Störstellen in einem einfach-kubischen Gitter. Das zentrale Resultat ist die Abstandsabhängigkeit der Korrelationen der Störstellenelektronen, welche stark von der Gitterdimension und der relativen Position der Störstellen abhängen. Bemerkenswert ist hier die lange Reichweite der Korrelationen in der Diagonalrichtung des Gitters. Außerdem ergibt sich, dass eine antiferromagnetische Austauschwechselwirkung ein Singulett zwischen den Störstellenelektronen gegenüber den Kondo-Singuletts der einzelnen Störstellen favorisiert und so den Kondo-Effekt der einzelnen Störstellen behindert.rnrnEin Zweiband-Hubbard-Modell, das Jz-Modell, wird im Hinblick auf seine Mott-Phasen in Abhängigkeit von Dotierung und Kristallfeldaufspaltung auf dem Bethe-Gitter untersucht. Die Entartung der Bänder ist durch eine unterschiedliche Bandbreite aufgehoben. Wichtigstes Ergebnis sind die Phasendiagramme in Bezug auf Wechselwirkung, Gesamtfüllung und Kristallfeldparameter. Im Vergleich zu Einbandmodellen kommen im Jz-Modell sogenannte orbital-selektive Mott-Phasen hinzu, die, abhängig von Wechselwirkung, Gesamtfüllung und Kristallfeldparameter, einerseits metallischen und andererseits isolierenden Charakter haben. Ein neuer Aspekt ergibt sich durch den Kristallfeldparameter, der die ionischen Einteilchenniveaus relativ zueinander verschiebt, und für bestimmte Werte eine orbital-selektive Mott-Phase des breiten Bands ermöglicht. Im Vergleich mit analytischen Näherungslösungen und Einbandmodellen lassen sich generische Vielteilchen- und Korrelationseffekte von typischen Mehrband- und Einteilcheneffekten differenzieren.rnrnDas zweite untersuchte Hubbard-Modell beschreibt eine magneto-optische Falle mit einer endlichen Anzahl Gitterplätze, in welcher fermionische Atome platziert sind. Es wird eine z-antiferromagnetische Phase unter Berücksichtigung nicht-lokaler Vielteilchenkorrelationen erhalten, und dabei werden bekannte Ergebnisse einer effektiven Einteilchenbeschreibung verbessert.rnrnDas korrelierte Halbmetall wird im Rahmen einer Mehrbandrechnung im Hinblick auf Korrelationseffekte untersucht. Ausgangspunkt ist eine ab initio-Beschreibung durch die Dichtefunktionaltheorie (DFT), welche dann durch die Hinzunahme lokaler Korrelationen ergänzt wird. Die Vielteilcheneffekte werden an Hand einer einfachen Wechselwirkungsnäherung verdeutlicht, und für ein Wechselwirkungsmodell in sphärischer Symmetrie präzisiert. Es ergibt sich nur eine schwache Quasiteilchenrenormierung. Besonders für röntgenspektroskopische Experimente wird eine gute Übereinstimmung erzielt.rnrnDie numerischen Ergebnisse für das Jz-Modell basieren auf Quanten-Monte-Carlo-Simulationen im Rahmen der dynamischen Molekularfeldtheorie (DMFT). Für alle anderen Systeme wird ein Mehrband-Algorithmus entwickelt und implementiert, welcher explizit nicht-diagonale Mehrbandprozesse berücksichtigt.rnrn

Ionic liquids in a multiscale study

Trotz des hohen Interesse an Ionischen Flüssigkeiten wird das zielgerichtete Design und die Anwendung Ionischer Flüssigkeiten durch fehlendes grundlegendes Verständnis erschwert. Deshalb wurde die Balance der molekularen Wechselwirkungen in Ionischen Flüssigkeiten studiert, um die Eigenschaften dieser zu verstehen und die Kraftfeldentwicklung im Rahmen des Multiskalenansatzes zu systematisieren. Es wurden reine Imidazolium-basierte Ionische Flüssigkeiten, Mischungen mit kleinen Molekülen und eine protische Ionische Flüssigkeit mit ab-initio-Methoden, hauptsächlich Car-Parrinello-Molekulardynamik, untersucht. Weiterhin wurden Eigenschaften der Flüssigphase mit denen von Ionenpaaren verglichen.rnIm Fokus standen die molekularen elektrostatischen Eigenschaften und es wurde gezeigt, dass Coulomb-Wechselwirkungen zu einzigartigen Charakteristika führten. So waren die Ionen-Nettoladungen stets reduziert, die molekularen Dipolmomentverteilungen sehr breit, elektronische Polarisation war entscheidend. Die elektrostatischen Eigenschaften waren allgemein lokal auf molekularen Größen- und Zeitskalen und hingen stark von Phasenzustand und Zusammensetzung ab. Für andere molekulare Eigenschaften, wie der Neigung zu dispersiven Kontakten oder Wasserstoffbrücken, wurde gezeigt, dass sie einen entscheidenden Einfluss auf die Feinstruktur Ionischer Flüssigkeiten hatten. Das Gleichgewicht der Wechselwirkungen zeigte sich auch in Leistungsspektren, die sich aus den ab-initio-Molekulardynamiksimulationen ergaben. Diese boten einen neuen Weg für den Vergleich zum Experiment und für einen Einblick in die schnelle Dynamik Ionischer Flüssigkeiten.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.