7 resultados para Numerical calculations
em ArchiMeD - Elektronische Publikationen der Universität Mainz - Alemanha
Resumo:
The goal of this thesis is the acceleration of numerical calculations of QCD observables, both at leading order and next–to–leading order in the coupling constant. In particular, the optimization of helicity and spin summation in the context of VEGAS Monte Carlo algorithms is investigated. In the literature, two such methods are mentioned but without detailed analyses. Only one of these methods can be used at next–to–leading order. This work presents a total of five different methods that replace the helicity sums with a Monte Carlo integration. This integration can be combined with the existing phase space integral, in the hope that this causes less overhead than the complete summation. For three of these methods, an extension to existing subtraction terms is developed which is required to enable next–to–leading order calculations. All methods are analyzed with respect to efficiency, accuracy, and ease of implementation before they are compared with each other. In this process, one method shows clear advantages in relation to all others.
Resumo:
This thesis deals with three different physical models, where each model involves a random component which is linked to a cubic lattice. First, a model is studied, which is used in numerical calculations of Quantum Chromodynamics.In these calculations random gauge-fields are distributed on the bonds of the lattice. The formulation of the model is fitted into the mathematical framework of ergodic operator families. We prove, that for small coupling constants, the ergodicity of the underlying probability measure is indeed ensured and that the integrated density of states of the Wilson-Dirac operator exists. The physical situations treated in the next two chapters are more similar to one another. In both cases the principle idea is to study a fermion system in a cubic crystal with impurities, that are modeled by a random potential located at the lattice sites. In the second model we apply the Hartree-Fock approximation to such a system. For the case of reduced Hartree-Fock theory at positive temperatures and a fixed chemical potential we consider the limit of an infinite system. In that case we show the existence and uniqueness of minimizers of the Hartree-Fock functional. In the third model we formulate the fermion system algebraically via C*-algebras. The question imposed here is to calculate the heat production of the system under the influence of an outer electromagnetic field. We show that the heat production corresponds exactly to what is empirically predicted by Joule's law in the regime of linear response.
Resumo:
Plasmonen stellen elektromagnetische Moden in metallischen Strukturen dar, in denen die quasifreien Elektronen im Metall kollektiv oszillieren. Während des letzten Jahrzehnts erfuhr das Gebiet der Plasmonik eine rasante Entwicklung, basierend auf zunehmenden Fortschritten der Nanostrukturierungsmethoden und spektroskopischen Untersuchungsmethoden, die zu der Möglichkeit von systematischen Einzelobjektuntersuchungen wohldefinierter Nanostrukturen führte. Die Anregung von Plasmonen resultiert neben einer radiativen Verstärkung der optischen Streuintensität im Fernfeld in einer nicht-radiativen Überhöhung der Feldstärke in unmittelbarer Umgebung der metallischen Struktur (Nahfeld), die durch die kohärente Ladungsansammlung an der metallischen Oberfläche hervorgerufen wird. Das optische Nahfeld stellt folglich eine bedeutende Größe für das fundamentale Verständnis der Wirkung und Wechselwirkung von Plasmonen sowie für die Optimierung plasmonbasierter Applikationen dar. Die große Herausforderung liegt in der Kompliziertheit des experimentellen Zugangs zum Nahfeld, der die Entwicklung eines grundlegenden Verständisses des Nahfeldes verhinderte.rnIm Rahmen dieser Arbeit wurde Photoemissionselektronenmikroskopie (PEEM) bzw. -mikrospektroskopie genutzt, um ortsaufgelöst die Eigenschaften nahfeld-induzierter Elektronenemission zu bestimmen. Die elektrodynamischen Eigenschaften der untersuchten Systeme wurden zudem mit numerischen, auf der Finiten Integrationsmethode basierenden Berechnungen bestimmt und mit den experimentellen Resultaten verglichen.rnAg-Scheiben mit einem Durchmesser von 1µm und einer Höhe von 50nm wurden mit fs-Laserstrahlung der Wellenlänge 400nm unter verschiedenen Polarisationszuständen angeregt. Die laterale Verteilung der infolge eines 2PPE-Prozesses emittierten Elektronen wurde mit dem PEEM aufgenommen. Aus dem Vergleich mit den numerischen Berechnungen lässt sich folgern, dass sich das Nahfeld an unterschiedlichen Stellen der metallischen Struktur verschiedenartig ausbildet. Insbesondere wird am Rand der Scheibe bei s-polarisierter Anregung (verschwindende Vertikalkomponente des elektrischen Felds) ein Nahfeld mit endlicher z-Komponente induziert, während im Zentrum der Scheibe das Nahfeld stets proportional zum einfallenden elektrischen Feld ist.rnWeiterhin wurde erstmalig das Nahfeld optisch angeregter, stark gekoppelter Plasmonen spektral (750-850nm) untersucht und für identische Nanoobjekte mit den entsprechenden Fernfeldspektren verglichen. Dies erfolgte durch Messung der spektralen Streucharakteristik der Einzelobjekte mit einem Dunkelfeldkonfokalmikroskop. Als Modellsystem stark gekoppelter Plasmonen dienten Au Nanopartikel in sub-Nanometerabstand zu einem Au Film (nanoparticle on plane, NPOP). Mit Hilfe dieser Kombination aus komplementären Untersuchungsmethoden konnte erstmalig die spektrale Trennung von radiativen und nicht-radiativen Moden stark gekoppelter Plasmonen nachgewiesen werden. Dies ist insbesondere für Anwendungen von großer Relevanz, da reine Nahfeldmoden durch den unterdrückten radiativen Zerfall eine große Lebensdauer besitzen, so dass deren Verstärkungswirkung besonders lange nutzbar ist. Ursachen für die Unterschiede im spektralen Verhalten von Fern- und Nahfeld konnten durch numerische Berechnungen identifiziert werden. Sie zeigten, dass das Nahfeld nicht-spärischer NPOPs durch die komplexe Oszillationsbewegung der Elektronen innerhalb des Spaltes zwischen Partikel und Film stark ortsabhängig ist. Zudem reagiert das Nahfeld stark gekoppelter Plasmonen deutlich empfindlicher auf strukturelle Störstellen des Resonators als die Fernfeld-Response. Ferner wurde der Elektronenemissionsmechanismus als optischer Feldemissionsprozess identifiziert. Um den Vorgang beschreiben zu können, wurde die Fowler-Nordheim Theorie der statischen Feldemission für den Fall harmonisch oszillierender Felder modifiziert.
Resumo:
Gegenstand dieser Arbeit ist die nummerische Berechnung von Schleifenintegralen welche in höheren Ordnungen der Störungstheorie auftreten.rnAnalog zur reellen Emission kann man auch in den virtuellen Beiträgen Subtraktionsterme einführen, welche die kollinearen und soften Divergenzen des Schleifenintegrals entfernen. Die Phasenraumintegration und die Schleifenintegration können dann in einer einzigen Monte Carlo Integration durchgeführt werden. In dieser Arbeit zeigen wir wie eine solche numerische Integration unter zu Hilfenahme einer Kontourdeformation durchgeführt werden kann. Ausserdem zeigen wir wie man die benötigeten Integranden mit Rekursionsformeln berechnen kann.
Resumo:
In dieser Arbeit stelle ich Aspekte zu QCD Berechnungen vor, welche eng verknüpft sind mit der numerischen Auswertung von NLO QCD Amplituden, speziell der entsprechenden Einschleifenbeiträge, und der effizienten Berechnung von damit verbundenen Beschleunigerobservablen. Zwei Themen haben sich in der vorliegenden Arbeit dabei herauskristallisiert, welche den Hauptteil der Arbeit konstituieren. Ein großer Teil konzentriert sich dabei auf das gruppentheoretische Verhalten von Einschleifenamplituden in QCD, um einen Weg zu finden die assoziierten Farbfreiheitsgrade korrekt und effizient zu behandeln. Zu diesem Zweck wird eine neue Herangehensweise eingeführt welche benutzt werden kann, um farbgeordnete Einschleifenpartialamplituden mit mehreren Quark-Antiquark Paaren durch Shufflesummation über zyklisch geordnete primitive Einschleifenamplituden auszudrücken. Ein zweiter großer Teil konzentriert sich auf die lokale Subtraktion von zu Divergenzen führenden Poltermen in primitiven Einschleifenamplituden. Hierbei wurde im Speziellen eine Methode entwickelt, um die primitiven Einchleifenamplituden lokal zu renormieren, welche lokale UV Counterterme und effiziente rekursive Routinen benutzt. Zusammen mit geeigneten lokalen soften und kollinearen Subtraktionstermen wird die Subtraktionsmethode dadurch auf den virtuellen Teil in der Berechnung von NLO Observablen erweitert, was die voll numerische Auswertung der Einschleifenintegrale in den virtuellen Beiträgen der NLO Observablen ermöglicht. Die Methode wurde schließlich erfolgreich auf die Berechnung von NLO Jetraten in Elektron-Positron Annihilation im farbführenden Limes angewandt.
Resumo:
This thesis is concerned with calculations in manifestly Lorentz-invariant baryon chiral perturbation theory beyond order D=4. We investigate two different methods. The first approach consists of the inclusion of additional particles besides pions and nucleons as explicit degrees of freedom. This results in the resummation of an infinite number of higher-order terms which contribute to higher-order low-energy constants in the standard formulation. In this thesis the nucleon axial, induced pseudoscalar, and pion-nucleon form factors are investigated. They are first calculated in the standard approach up to order D=4. Next, the inclusion of the axial-vector meson a_1(1260) is considered. We find three diagrams with an axial-vector meson which are relevant to the form factors. Due to the applied renormalization scheme, however, the contributions of the two loop diagrams vanish and only a tree diagram contributes explicitly. The appearing coupling constant is fitted to experimental data of the axial form factor. The inclusion of the axial-vector meson results in an improved description of the axial form factor for higher values of momentum transfer. The contributions to the induced pseudoscalar form factor, however, are negligible for the considered momentum transfer, and the axial-vector meson does not contribute to the pion-nucleon form factor. The second method consists in the explicit calculation of higher-order diagrams. This thesis describes the applied renormalization scheme and shows that all symmetries and the power counting are preserved. As an application we determine the nucleon mass up to order D=6 which includes the evaluation of two-loop diagrams. This is the first complete calculation in manifestly Lorentz-invariant baryon chiral perturbation theory at the two-loop level. The numerical contributions of the terms of order D=5 and D=6 are estimated, and we investigate their pion-mass dependence. Furthermore, the higher-order terms of the nucleon sigma term are determined with the help of the Feynman-Hellmann theorem.
Resumo:
In vielen Teilgebieten der Mathematik ist es w"{u}nschenswert, die Monodromiegruppe einer homogenen linearen Differenzialgleichung zu verstehen. Es sind nur wenige analytische Methoden zur Berechnung dieser Gruppe bekannt, daher entwickeln wir im ersten Teil dieser Arbeit eine numerische Methode zur Approximation ihrer Erzeuger.rnIm zweiten Abschnitt fassen wir die Grundlagen der Theorie der Uniformisierung Riemannscher Fl"achen und die der arithmetischen Fuchsschen Gruppen zusammen. Auss erdem erkl"aren wir, wie unsere numerische Methode bei der Bestimmung von uniformisierenden Differenzialgleichungen dienlich sein kann. F"ur arithmetische Fuchssche Gruppen mit zwei Erzeugern erhalten wir lokale Daten und freie Parameter von Lam'{e} Gleichungen, welche die zugeh"origen Riemannschen Fl"achen uniformisieren. rnIm dritten Teil geben wir einen kurzen Abriss zur homologischen Spiegelsymmetrie und f"uhren die $widehat{Gamma}$-Klasse ein. Wir erkl"aren wie diese genutzt werden kann, um eine Hodge-theoretische Version der Spiegelsymmetrie f"ur torische Varit"aten zu beweisen. Daraus gewinnen wir Vermutungen "uber die Monodromiegruppe $M$ von Picard-Fuchs Gleichungen von gewissen Familien $f:mathcal{X}rightarrow bbp^1$ von $n$-dimensionalen Calabi-Yau Variet"aten. Diese besagen erstens, dass bez"uglich einer nat"urlichen Basis die Monodromiematrizen in $M$ Eintr"age aus dem K"orper $bbq(zeta(2j+1)/(2 pi i)^{2j+1},j=1,ldots,lfloor (n-1)/2 rfloor)$ haben. Und zweitens, dass sich topologische Invarianten des Spiegelpartners einer generischen Faser von $f:mathcal{X}rightarrow bbp^1$ aus einem speziellen Element von $M$ rekonstruieren lassen. Schliess lich benutzen wir die im ersten Teil entwickelten Methoden zur Verifizierung dieser Vermutungen, vornehmlich in Hinblick auf Dimension drei. Dar"uber hinaus erstellen wir eine Liste von Kandidaten topologischer Invarianten von vermutlich existierenden dreidimensionalen Calabi-Yau Variet"aten mit $h^{1,1}=1$.