Spelling suggestions: "subject:"lattice"" "subject:"lattices""
41 |
The step scaling function of QCD at negative flavor numberGehrmann, Bernd 05 June 2002 (has links)
Wir untersuchen eine SU(3) Yang-Mills-Theorie mit einer Kopplung an ein bosonisches Spinorfeld. Diese als Bermion-Modell bekannte Theorie entspricht formal QCD mit minus zwei Quark-Flavors. Gegenüber der vollen QCD erfordert sie wesentlich weniger Computerzeit und ist deshalb als relativ kostengünstiges Testmodell geeignet. Im Mittelpunkt unseres Interesses steht die Step-Scaling-Funktion, die die Skalenabhängigkeit der laufenden Kopplung im Schrödinger-Funktional-Renormierungsschema beschreibt. Mit Hilfe einer nicht-perturbativen Finite-Size-Technik kann sie benutzt werden, um den Lambda-Parameter, der die Kopplung bei hohen Energien charakterisiert, aus experimentellen Daten bei niedrigen Energien zu bestimmen. Wir studieren im Detail die Gitterartefakte und die Kontinuumsextrapolation der aus Gittersimulationen bestimmten Step-Scaling-Funktion, wenn O(a)-Verbesserung nach Symanzik verwendet wird. Unsere Resultate stellen wir dem Fall von unverbesserten Bermionen und dynamischen Fermionen gegenüber, und vergleichen im Kontinuumslimes mit renormierter Störungstheorie. Weiterhin betrachten wir im Bermion-Modell die Step-Scaling-Funktion mit massiven Quarks. Nach dem Appelquist-Carazzone-Theorem erwartet man, daß Beiträge von Materiefeldern mit ansteigender Masse verschwinden, so daß die Step-Scaling-Funktion gegen den Fall reiner Eichtheorie konvergieren sollte. Wenn man nicht-perturbativ verschiedene effektive Theorien mit verschiedener Anzahl von aktiven Quarks über Massenschwellen hinweg verbinden will, sollten Gitterartefakte klein sein. Um die Durchführbarkeit einer solchen Methode zu testen, untersuchen wir die Step-Scaling-Funktion und ihre Gitterartefakte für verschiedene Massen. Für die Monte-Carlo-Simulation von verbesserten Bermionen entwickeln wir einen geeigneten Algorithmus und vergleichen seine Effizienz mit unbesserten Bermionen und mit voller QCD. Als vorbereitende Studie vergleichen wir die Effizienz verschiedener Algorithmen in reiner Eichtheorie. / As a computationally less costly test case for full QCD, we investigate an SU(3) Yang-Mills theory coupled to a bosonic spinor field. This theory corresponds to QCD with minus two quark flavors and is known as the bermion model. Our central object of interest is the step scaling function which describes the scale evolution of the running coupling in the Schrödinger functional scheme. With the help of a non-perturbative recursive finite size technique, it can be used to determine the Lambda parameter, which characterizes the coupling at high energy, from experimental input at low energies. We study in detail the lattice artefacts and the continuum extrapolation of the step scaling function from lattice simulations when O(a) improvement according to the Symanzik programme is used. Our results are compared to the unimproved bermion and dynamical fermion cases, and to renormalized perturbation theory in the continuum limit. For the bermion model, we also examine the step scaling function with massive quarks. According to the Appelquist-Carazzone theorem the contributions from matter fields are expected to vanish for large masses, such that the step scaling function converges to the pure gauge theory case. If one wants to connect non-perturbatively different effective theories with different numbers of active quarks over flavor thresholds, lattice artefacts should be reasonably small. In order to test the feasibility of such a method, we investigate the step scaling function and its lattice artefacts for several values of the mass. For the Monte Carlo simulation of improved bermions, we develop a suitable algorithm and compare its performance with unimproved bermions and full QCD. As a preparative study, we compare the efficiency of algorithms in pure gauge theory.
|
42 |
Chiral fermions on the latticeChiarappa, Thomas 09 December 2004 (has links)
Das Ziel dieser Doktorarbeit ist die Berechnung, aus ersten Prinzipien, von Niedrigenenergiekonstanten (NEK), welche die chirale Störungsheorie (ChST) parametrisieren, durch Simulationen auf dem Gitter. Diese Arbeit ist eine Pilotstudie und will deshalb nicht zu definitiven und präzisen quantitativen Vorhersagen führen, sondern wir möchten qualitative Leitlinien für zukünftige genauere Erforschungen der epsilon-Entwicklung der ChST geben. Letztere ist ein Gebiet in dem die chirale Symmetrie wiederhergestellt ist und wo die Compton-Wellenlänge der leichtesten Mesonen grösser als die lineare Dimension des endlichen Volumens ist. Die epsilon-Region ist weiterhin durch die wichtige Rolle der topologischen Ladung, nu, charakterisiert, so dass die Untersuchung der mesonischen Zweipunkt-Korrelationsfunktionen in definierten topologischen Sektoren vorzunehmen ist. Aus diesem Grunde haben wir chirale Fermionen auf dem Gitter simuliert, wobei der overlap Formalismus für den Dirac Operator, mit dem gewöhnliches Wilson Dirac Operator als Kern, in der quenched Approximation genommen wurde. Wir demonstrieren dass der Sektor mit neutraler topologischer Ladung, nu = 0, numerisch sehr schwierig zu untersuchen ist, weil die Daten von erheblichen Spitzen beeinflusst werden, die aufgrund sehr kleiner, nicht verschwindender Eigenwerte herrühren. Diese Beobachtung findet eine Übereinstimmung in einer Studie der Random Matrix Theory (RMT), die besagt, dass man eine Statistik von mehr als 10000 Konfigurationen benötigt, falls man auf kleine Eigenwerte empfindlich reagierende physikalische Grössen, im triviale topologische Sektoren simulieren möchte. Weniger Probleme ergeben sich in dem nicht-trivialen Sektoren |nu| > 0: wir präsentieren deshalb unsere Ergebnisse aus dem topologisches Sektor nu = 1 und den Daten aus das Sektor mit nu = 2 nur als Gegenprobe. Wir zeigen die Existenz einer unteren Schranke für das physikalisches Volumen bei V > 1 fm zur vierten Potenz für die Gültigkeit der ChST, übereinstimmend mit einer vorherigen Beobachtung aus der RMT. Wir beschränken deshalb unsere Untersuchung auf ein grösseres Gitter, wo wir die Daten des axialen Korrelators mit den Vorhersagen der ChST in der Valenzquark-Approximation fitten und eine stabile Bestimmung der gequenchten Pionen-Zerfallskonstante, F, erhalten. Die skalaren und pseudoskalaren Korrelationsfunktionen sind durch mehrere NEK parametrisiert. Dadurch ist der Vergleich mit den Vorhersagen der ChST schwieriger, und deshalb präsentieren wir nur einige Abschätzungen, die mit der entsprechenden Literatur verglichen werden. Zum Schluss präsentieren wir mögliche Implementierungen von verbesserten Algorithmen für die Inversion des overlap Operators, die die Simulationkosten reduzieren. / The aim of this thesis is the computation of Low Energy Constants (LEC) which parameterise Chiral Perturbation Theory (ChPT) from a first principles analysis via lattice simulations. The thesis provides a pilot study and will not give definitive and precise quantitative predictions, but rather our aim is to provide qualitative hints for future accurate investigations of the epsilon-expansion of ChPT, where chiral symmetry is restored and the Compton wavelength of the lightest meson is larger than the linear size of the finite volume. One of the property characterising the epsilon-regime is the important role played by the topological charge, nu, leading to the investigation of the two-point meson correlation functions in distinct topological sectors. To this end, we simulate chiral fermions on the lattice adopting the overlap formalism for the Dirac operator in the quenched approximation, with the kernel provided by the usual Wilson Dirac operator. We demonstrate that the neutral topological sector, nu = 0, is very difficult to explore numerically, as the data are affected by large spikes due to the presence of very small, non-zero eigenvalues. This observation is in agreement with a study of Random Matrix Theory (RMT), which indicates that a statistic of more than 10000 configurations is required when physical quantities sensitive to small eigenvalues are investigated in the neutral topological sector. Therefore, we present our results corresponding to the topological sector nu = 1. Due to the modest statistic, we only use the nu = 2 data as a crosscheck. We find a lower bound on the physical volume V > 1 fm to the four for ChPT to be used, in agreement with a previous observation using RMT. Restricting our attention to a larger lattice, we fit the data of the axial correlation function with the predictions of quenched ChPT, obtaining a stable determination of the quenched pion decay constant, F. The scalar and pseudoscalar correlation functions are parameterised by a larger number of LEC, rendering the comparison with ChPT predictions much more difficult and hence we present only some estimates that are compared with other determinations in the literature. Finally, we present possible implementations of improved algorithms used in the ``inversion'''' of the overlap operator, whose aim is to reduce the computational cost of the simulations.
|
43 |
Investigating new lattice approaches to the momentum and spin structure of the nucleonWiese, Christian 03 June 2016 (has links)
Diese Arbeit beschäftigt sich mit der Berechnung von für die Struktur des Nukleons relevanten Observablen, die experimentell durch inklusive und semi-inklusive Streuexperimente bestimmt werden können. Es werden zwei Pilotstudien erörtert, welche die Spin- und Impulsstruktur des Nukleons mithilfe von Gitter-QCD untersuchen. Hierfür wird der Twisted-Mass-Formalismus mit dynamischen Fermionen verwendet, um sicherzustellen, dass die untersuchten Größen einen verbesserten Kontinuumslimes aufweisen. Der erste Teil dieser Arbeit untersucht die Umsetzbarkeit einer Rechnung, die sich mit dem durchschnittlichen Impulsanteil der Gluonen im Nukleon auseinandersetzt. Diese Größe wurde bisher kaum im Rahmen der Gitter-QCD behandelt. In diesem Zusammenhang werden zwei verschiedene Gittermethoden untersucht: das Feynman-Hellman-Theorem, sowie die direkte Berechnung der relevanten Formfaktoren. Mithilfe der zweiten Methode und mehreren Iterationen des Schmierens der Eichlinks ist es möglich, statistisch aussagekräftige Resultate zu erhalten. Die zweite Studie beschäftigt sich mit der direkten Berechnung der vollständigen Impuls- und Spinverteilung von Quarks und Antiquarks im Nukleon. Hierfür wird untersucht, ob eine kürzlich publizierte Methode praktikabel ist, nach der eine räumliche Quasiverteilung zu berechnen und aus dieser die physikalische Verteilung abzuleiten ist. In diesem Zusammenhang wird der Einfluß des Schmierens der Eichlinks und unterschiedlicher Impulsboosts des Nukleons erprobt. Die anschließend berechneten Isovektor-Quarkverteilungen (unpolarisiert und polarisiert) weisen eine gute qualitative Übereinstimmung mit Verteilungen auf, die mithilfe von phänomenologischen Analysen bestimmt wurden. Zentrale Erkenntnis dieser Arbeit ist der Nachweis, dass es auf dem Gitter prinzipiell möglich ist, beide Observablen zu berechnen. Trotzdem muss noch erheblich mehr Arbeit aufgewendet werden, um verlässliche Resultate für diese Größen zu erhalten. / This thesis deals with the theoretical computation of nucleon structure observables as they can be experimentally obtained from inclusive and semi-inclusive scattering experiments. I present two exploratory studies on spin and momentum structure observables of the nucleon in the framework of lattice QCD. Throughout this work, I use the twisted mass formalism with dynamical fermions at maximal twist, which ensures an improved continuum limit scaling for the relevant quantities. In the first part, I investigate the feasibility of a lattice calculation of the gluons’ average momentum fraction in the nucleon, a quantity that is rarely studied in lattice QCD. For this purpose, I study two different methods, namely the Feynman-Hellman theorem and the direct computation of the relevant form factor. Applying the latter method and combining it with several steps of stout gauge link smearing, I obtain a statistically significant results for the gluon content. The second study is concerned with the direct computation of the full momentum and spin distribution of quarks and antiquarks within the nucleon. I investigate the feasibility of a recently published approach proposing the computation of a purely spatial quasi-distribution that can be related to the physical distribution. I test the influence of gauge link smearing and different nucleon momentum boosts on the lattice data. Ultimately, I obtain iso-vector quark distributions for the unpolarized and polarize case that featuring a decent qualitative agreement to quark distributions acquired from phenomenological fits. As a key result of this work, I demonstrate that the demanding calculation of gluon content and the novel approach of computing quark distributions directly within lattice QCD are feasible in principle, although significantly more effort has to be invested into obtaining accurate results with reliable uncertainties.
|
44 |
Form factors for semileptonic Bs → Kℓνdecays in lattice QCDBahr, Felix Tobias 02 December 2015 (has links)
Wir präsentieren eine Forschungsstudie zur Berechnung des Formfaktors f_+(q^2) für den semileptonischen Zerfall B_s -> K l nu in Gitter-QCD-Simulationen im großen Volumen mit zwei dynamischen Seequark-Flavours mit O(a)-verbesserten Wilson-Fermionen. Wir diskutieren die Berechnung relevanter Zwei- und Dreipunkt-Funktionen und betrachten komplementäre Methoden diese zu kombinieren, um den Formfaktor zu erhalten. Insbesondere stellen wir die Strategie eines kombinierten Fits vor, in den Datenpunkte aller Korrelatoren eingehen und der als Fitparameter Energien, Amplituden und den Formfaktor hat. Das b-Quark wird in HQET behandelt; unsere momentane Analyse konzentriert sich auf den statischen Grenzfall. Indes haben wir den Code für alle nötigen O(1/m_h)-Korrekturen entwickelt und die Messungen dieser durchgeführt; sie werden verwendet werden, sobald ihre Koeffizienten von der ALPHA-Kollaboration bestimmt worden sein werden. Um den Formfaktor auf allen Ensembles bei dem gleichen Wert des Impulsübertrags q^2 bestimmen zu können, führen wir getwistete Randbedingungen für das s- und das b-Quark ein, die ein freies Einstellen der Quarkimpulse und damit von q^2 ermöglichen. Wir führen Messungen auf einer Untermenge von N_f=2 CLS Eichkonfigurationen durch und erhalten den Formfaktor bei drei verschiedenen Gitterabständen und etwa gleicher Pionmasse von ungefähr 330 MeV. Wir benutzen diese, um eine Kontinuumsextrapolation durchzuführen, und beobachten, dass diese relativ flach in a^2 ist. Eine Messung bei einer unterschiedlichen Pionmasse deutet an, dass Quarkmassen-Effekte klein sind. Wir vergleichen unseren Kontinuumswert des Formfaktors mit kürzlich veröffentlichten Ergebnissen anderer Kollaborationen und stellen eine gute Übereinstimmung fest. / We present an exploratory study of the calculation of the form factor f_+(q^2) for the semileptonic decay B_s -> K l nu in large-volume lattice QCD simulations with two dynamical sea quark flavours using O(a) improved Wilson fermions. We discuss the computation of relevant two- and three-point functions and consider complementary methods how these can be combined to obtain the form factor. In particular, we put forward the strategy of a combined fit in which data of all correlators enter and which has as fit parameters energies and amplitudes of the correlators and the form factor. The b quark is treated in HQET; our present analysis focuses on the static limit. Meanwhile, we have developed the code and performed the measurements of all needed O(1/m_h) corrections which will be used as soon as their coefficients will have been computed by the ALPHA collaboration. In order to be able to measure the form factor at the same value of the momentum transfer q^2 on all ensembles, we impose twisted boundary conditions on the s and b quarks that allow for a free tuning of the quark momenta and thus of q^2. We perform measurements on a subset of N_f=2 CLS gauge configurations, obtaining the form factor at three different lattice spacings and roughly the same pion mass of about 330 MeV. Using these, we carry out a continuum extrapolation and observe that it is relatively flat in a^2. A measurement at a different pion mass indicates that quark mass effects are small. We compare our continuum value of the form factor with recently published results of other collaborations and observe a good agreement.
|
45 |
The renormalised quark mass in the Schrödinger functional of lattice QCDKurth, Stefan 04 September 2002 (has links)
Diese Arbeit befasst sich mit störungstheoretischen Rechnungen zur renormierten Quarkmasse im Schrödinger-Funktional mit nicht verschwindendemHintergrundfeld. Als Grundlage der Rechnungen werden das Schroedinger-Funktional und seinegrundlegenden Eigenschaften erläutert. Auch die O(a)-Verbesserung, die zu einem schnelleren Erreichen des Kontinuumslimes fuehren soll, wird in diesem Zusammenhang dargestellt.Des weiteren wird erklärt, aufwelche Weise das Schrödinger-Funktional dazu dient, das Skalenverhaltenrenormierter Größen ueber einen grossen Energiebereich zuuntersuchen. Das Skalenverhalten sowohl der renormierten Kopplung als auchder renormierten Quarkmassen wird in diesem Schema durch Step-Scaling-Funktionenbeschrieben. Die Definition der renormierten Kopplung wird dargestellt,ebenso die Definition der renormierten Masse, die mit Hilfe derPCAC-Relation ueber den Axialvektorstrom und die Pseudoskalardichte erfolgt. Die Skalenabhängigkeit der renormierten Massewird auf die Skalenabhängigkeit der Renormierungskonstanten derPseudoskalardichte zurueckgefuehrt. Breiten Raum nimmt die Berechnung verschiedenerKorrelationsfunktionen bis zur Ein-Loop-Ordnung in Stoerungstheorie ein. Mit Hilfe der soermittelten Koeffizienten wird die kritische Quarkmasse, bei der die renormierte Masse verschwindet, in Ein-Loop-Naeherung berechnet,ebenso der Ein-Loop-Koeffizent der Renormierungskonstanten der Pseudoskalardichte. Mit Hilfe dieses Koeffizienten wird aus der bekanntenanomalen Dimension in Zwei-Loop-Ordnung im MS-bar-Schemadie anomale Dimension im Schrödinger-Funktional berechnet. Als weitere Anwendung der Störungstheorie werden verschiedene Diskretisierungsfehler bestimmt. Die kritische Quarkmasse in Ein-Loop-Ordnunggeht in den Zwei-Loop-Koeffizienten des Diskretisierungfehlers der Step-Scaling-Funktion der renormierten Kopplung ein, der durchdie Abweichung dieser Funktion von ihrem Kontinuumslimes definiert ist.Verschiedene Diskretisierungsfehler der Strommasse, die durch die PCAC-Relationmit unrenormiertem Axialvektorstrom und Pseudoskalardichte definiert ist, werdenin Ein-Loop-Ordnung berechnet. Ein wichtiger Diskretisierungsfehler derrenormierten Quarkmasse ist die Abweichung ihrer Step-Scaling-Funktion vomKontinuumslimes. Dieser Fehler ist in Ein-Loop-Ordnung bislang nur mitverschwindendem Hintergrundfeld bekannt und wird in dieser Arbeit mitnicht verschwindendem Hintergrundfeld berechnet. / The renormalised quark mass in the Schroedinger functional is studied perturbatively with a non-vanishing background field. The framework in which the calculations are done is the Schroedinger functional. Its definition and basic properties are reviewed and it is shown how to make the theory converge faster towards its continuum limit by O(a) improvement. It is explained how the Schroedinger functional scheme avoids the implications of treating a large energy range on a single lattice in order to determine the scale dependence of renormalised quantities. The description of the scale dependence by the step scaling function is introduced both for the renormalised coupling and the renormalised quark masses. The definition of the renormalised coupling in the Schroedinger functional is reviewed, and the concept of the renormalised mass being defined by the axial current and density via the PCAC-relation is explained. The running of the renormalised mass described by its step scaling function is presented as a consequence of the fact that the renormalisation constant of the axial density is scale dependent. The central part of the thesis is the expansion of several correlation functions up to 1-loop order. The expansion coefficients are used to compute the critical quark mass at which the renormalised mass vanishes, as well as the 1-loop coefficient of the renormalisation constant of the axial density. Using the result for this renormalisation constant, the 2-loop anomalous dimension is obtained by conversion from the MS-bar-scheme. Another important application of perturbation theory carried out in this thesis is the determination of discretisation errors. The critical quark mass at 1-loop order is used to compute the deviation of the coupling's step scaling function from its continuum limit at 2-loop order. Several lattice artefacts of the current quark mass, defined by the PCAC relation with the unrenormalised axial current and density, are computed at 1-loop order. An essential property of the renormalised quark mass being computed in this thesis at 1-loop order is the deviation of its step scaling function from the continuum limit, which was so far only known for the zero background field case.
|
46 |
On the chirally rotated Schrödinger functional with Wilson fermionsLópez, Jénifer González 13 July 2011 (has links)
Viele Phaenomene in der Natur sind eng verknuepft mit dem Niederenergieverhalten der QCD und damit von nicht-perturbative Natur. Viele Groeßen benoetigen auch eine nicht-perturbative Renormierung. Als nicht-perturbative Renormierungsschema schlagen wir das chiral gedrehte Schroedingerfunktional, χSF, in einer Gitterregularisierung vor. Auf dem Baumgraphenniveau wird eine analytische Rechnung im Kontinuum und auf dem Gitter durchgefuehrt. Weitere Untersuchungen werden dann in der Valenzquark-Approximation der Gitter QCD durchgefuehrt. Eines der Hauptziele ist es dabei, die im χSF benoetigten Koeffizienten nicht-perturbativ so einzustellen, dass ein wohl-definierter Kontinuumlimes durchgefuehrt werden kann. Es wird gezeigt, dass solch eine Feineinstellung der Parameter des χSF durchfuehrbar ist und dass physikalische Groeßen nicht sensitiv auf die spezielle Wahl der Bedingung zur Einstellung der Parameter sind. Es wird gezeigt, dass das Skalierungsverhalten physikalischer Groeßen konsistent mit fuehrenden O(a2) Diskretisierungseffekten ist. Das Hauptergebnis dieser Arbeit ist der Nachweis, dass das χSF mit den hier berechneten Verbesserungskoeffizienten, zu einem korrekten Kontinuumlimes fuehrt. Dazu wurden drei unterschiedliche Werte der Renormierungsskala verwendet und mehrere uns interessierende physikalische Groeßen berechnet. Wir koennen deshalb den Schluss ziehen, dass das χSF ein viel versprechendes Renormierungsschema darstellt, um eine nicht-perturbative Renormierung vorzunehmen und dabei gleichzeitig die automatische O(a)-Verbesserung aufrecht erhalten. Dies eroeffnet den sehr wichtigen Ausblick, dass das χSF in zukuenftigen nicht-perturbativen Berechnungen von Renormierungskonstanten auch ueber die Valenzquark-Approximation hinaus eingesetzt werden kann. / There are many phenomena in nature which are closely linked to the low energy regime of QCD. Theoretically, these can be dealt with only by means of non-perturbative methods. Often, a non-perturbative renormalization of QCD is required. We employ a 4-dimensional lattice as a regulator of QCD. As a non-perturbative renormalization scheme, we propose the chirally rotated Schrödinger functional, χSF. We perform analytical studies at tree-level of perturbation theory, in the continuum and on the lattice. Beyond tree-level, all studies are performed in the quenched approximation of QCD. One of the main targets has been to perform the non-perturbative tuning of the two required coefficients of the χSF scheme, such that a well defined continuum limit can be reached. We demonstrate that the tuning is feasible and physical quantities are insensitive to the tuning condition. There are also a couple of improvement counterterms at the boundaries. However, besides these boundary O(a) effects, the χSF is expected to be compatible with bulk automatic O(a)-improvement. We show that the scaling behavior of physical quantities is consistent with automatic O(a)-improvement. The other most important achievement has been to demonstrate that the χSF, with the here computed tuning coefficients, leads to the correct continuum limit. For this, we have performed universality tests of the continuum limit, at three different values of the renormalization scale and through the computation of several physical quantities of interest. The conclusion of these results is that the χSF is a promising scheme to perform non-perturbative renormalizations while maintaining bulk automatic O(a)-improvement. This opens the most relevant prospect that the χSF can be safely used in future non-perturbative computations of renormalization factors also beyond the quenched approximation.
|
47 |
Critical slowing down and error analysis of lattice QCD simulationsVirotta, Francesco 07 May 2012 (has links)
In dieser Arbeit untersuchen wir das Critical Slowing down der Gitter-QCD Simulationen. Wir führen eine Vorstudie in der quenched Approximation durch, in der wir feststellen, dass unsere Schätzung der exponentiellen Autokorrelation wie $\tauexp(a) \sim a^{-5} $ skaliert, wobei $a$ der Gitterabstand ist. In unquenched Simulationen mit O(a)-verbesserten Wilson-Fermionen finden wir ein ähnliches Skalierungsgesetz. Die Diskussion wird von einem gro\ss{}en Satz an Ensembles sowohl in reiner Eichtheorie als auch in der Theorie mit zwei entarteten Seequarks unterstützt. Wir haben darüber hinaus die Wirkung von langsamen algorithmischen Modi in der Fehleranalyse des Erwartungswertes von typischen Gitter-QCD-Observablen (hadronische Matrixelemente und Massen) untersucht. Im Kontext der Simulationen, die durch langsame Modi betroffen sind, schlagen wir vor und testen eine Methode, um zuverlässige Schätzungen der statistischen Fehler zu bekommen. Diese Methode soll in dem typischen Simulationsbereich der Gitter-QCD helfen, nämlich dann, wenn die gesamte erfasste Statistik O(10)\tauexp ist. Dies ist der typische Fall bei Simulationen in der Nähe des Kontinuumslimes, wo der Rechenaufwand für die Erzeugung von zwei unabhängigen Datenpunkten sehr gro\ss{} sein kann. Schlie\ss{}lich diskutieren wir die Skalenbestimmung in N_f=2-Simulationen mit der Kaon Zerfallskonstante f_K als experimentellem Input. Die Methode wird zusammen mit einer gründlichen Diskussion der angewandten Fehleranalyse erklärt. Eine Beschreibung der öffentlich zugänglichen Software, die für die Fehleranalyse genutzt wurde, ist eingeschlossen. / In this work we investigate the critical slowing down of lattice QCD simulations. We perform a preliminary study in the quenched approximation where we find that our estimate of the exponential auto-correlation time scales as $\tauexp(a)\sim a^{-5}$, where $a$ is the lattice spacing. In unquenched simulations with O(a) improved Wilson fermions we do not obtain a scaling law but find results compatible with the behavior that we find in the pure gauge theory. The discussion is supported by a large set of ensembles both in pure gauge and in the theory with two degenerate sea quarks. We have moreover investigated the effect of slow algorithmic modes in the error analysis of the expectation value of typical lattice QCD observables (hadronic matrix elements and masses). In the context of simulations affected by slow modes we propose and test a method to obtain reliable estimates of statistical errors. The method is supposed to help in the typical algorithmic setup of lattice QCD, namely when the total statistics collected is of O(10)\tauexp. This is the typical case when simulating close to the continuum limit where the computational costs for producing two independent data points can be extremely large. We finally discuss the scale setting in Nf=2 simulations using the Kaon decay constant f_K as physical input. The method is explained together with a thorough discussion of the error analysis employed. A description of the publicly available code used for the error analysis is included.
|
48 |
Determination of csw in Nf = 3 + 1 Lattice QCD with massive Wilson fermionsStollenwerk, Felix 22 February 2017 (has links)
Um aussagekräftige, mit dem Experiment vergleichbare Resultate aus Berechnungen der Gitter-QCD zu erhalten, ist die Extrapolation zum Kontinuum unabdingbar. Das bewährte Symanzik-Verbesserungsprogramm führt zu einer systematischen Reduzierung der Ordnung von Cutoff-Effekten, die eine bessere Kontrolle über die genannten Fehler sowie größere und damit erschwinglichere Gitterabstände ermöglicht. Auf die Wilson-Fermionenwirkung bezogen bedarf es nur des Hinzufügens des Sheikholeslami-Wohlert-Terms mit dem O(a)-Verbesserungskoeffizienten csw. In der vorliegenden Arbeit wird eine Strategie zur nicht-perturbativen Bestimmung dieses Koeffizienten in der Theorie mit Nf=3+1 massiven Seequarks entwickelt. Diese ist in ein allgemeines, massenabhängiges Renormierungs- und Verbesserungsschema eingebettet, dessen Grundlagen dargelegt werden. Die Auferlegung der Verbesserungsbedingung, bei der die PCAC-Relation im Schrödinger-Funktional Verwendung findet, geschieht entlang einer Linie konstanter Physik, welche dem Charm-Quark näherungsweise seine physikalische Masse zuordnet. Dieser vergleichsweise aufwendige Ansatz hat zum Ziel, große, massenabhängige O(a^2)-Effekte in zukünftigen Simulationen im großen Volumen mit vier dynamischen Quarkspezies zu vermeiden. Die numerischen Resultate dieser Arbeit werden unter Verwendung der tree-level-verbesserten Lüscher-Weisz-Eichwirkung gewonnen. Da die sogenannte Gradient-Flow-Kopplung bei der Definition der Linie konstanter Physik Verwendung findet, wird in einer zusätzlichen Untersuchung die Wechselbeziehung dieser Kopplung mit der Topologischen Ladung beleuchtet, insbesondere im Bezug auf die unter den Namen Critical Slowing Down und Topology Freezing bekannten Phänomene. / In order to obtain sensible results from Lattice QCD that may be compared with experiment, extrapolation to the continuum is crucial. The well-established Symanzik improvement program systematically reduces the order of cutoff effects, allowing for better control of the aforementioned errors, as well as larger and thus more affordable lattice spacings. Applied to the Wilson fermion action, it entails the addition of the Sheikholeslami–Wohlert term with the O(a) improvement coefficient csw. In this work, a strategy is developed for the non-perturbative determination of csw in the theory with Nf=3+1 massive sea quarks. It is embedded in a general, mass-dependent renormalization and improvement scheme, for which we lay the foundations. The improvement condition, formulated by means of the PCAC relation in the Schrödinger Functional, is imposed along a line of constant physics that is designed to be close to the physical mass of the charm quark. The aim of this rather elaborate approach is to avoid large, mass-dependent O(a^2) effects in future large volume simulations with four dynamical quark species. The numerical results are worked out using the tree-level improved Lüscher–Weisz gauge action. Since the gradient flow coupling is employed in the definition of the line of constant physics, its interdependence with the topological charge in regard to critical slowing down and topology freezing is investigated in a supplemental study.
|
49 |
Gluon and ghost propagator studies in lattice QCD at finite temperatureAouane, Rafik 14 May 2013 (has links)
Die im infraroten Impulsbereich der Quantenchromodynamik (QCD) berechneten Gluon- und Ghost-Propagatoren spielen eine große Rolle für das sogenannte Confinement der Quarks und Gluonen. Sie sind Gegenstand intensiver Foschungen dank nicht-perturbativer Methoden basierend auf Dyson-Schwinger- (DS) und funktionalen Renormierungsgruppen-Gleichungen (FRG). Darüber hinaus sollte es deren Verhalten bei endlichen Temperaturen erlauben, den chiralen und Deconfinement-Phasenübergang bzw. das Crossover in der QCD besser aufzuklären. Unser Zugang beruht auf der gitter-diskretisierten QCD (LQCD), die es als ab-initio-Methode gestattet, verschiedenste störungstheoretisch nicht zugängliche QCD-Observablen der hadronischen Welt zu berechnen. Wir untersuchen das Temperaturverhalten der Gluon- und Ghost-Propagatoren in der Landau-Eichung für die reine Gluodynamik und die volle QCD. Für den Gluon-Propagator berechnen wir deren longitudinale (DL) sowie transversale (DT) Komponenten. Ziel ist es, Datensätze in Form von Fit-Formeln zu liefern, welche als Input für die DS- (oder FRG-) Gleichungen verwendet werden können. Wir beschäftigen uns mit der vollen (Nf=2) LQCD unter Verwendung der sogenannten twisted mass Fermiondiskretisierung. Von der tmfT-Kollaboration wurden uns dafür Eichfeldkonfigurationen für Temperaturen im Crossover-Bereich sowie jeweils für drei fixierte Pion-Massenwerte im Intervall [300, 500] MeV bereitgestellt. Schließlich berechnen wir innerhalb der reinen SU(3) Eichtheorie (bei T=0) den Landau Gluon-Propagator unter Verwendung verschiedener Eichfixierungskriterien. Unser Ziel ist es, den Einfluss von Eich-Kopien mit minimalen (nicht-trivialen) Eigenwerten des Faddeev-Popov-Operators zu verstehen. Eine solche Studie soll klären, wie Gribov-Kopien das Verhalten der Gluon- und Ghost-Propagatoren im infraroten Bereich prinzipiell beeinflussen. / Gluon and ghost propagators in quantum chromodynamics (QCD) computed in the infrared momentum region play an important role to understand quark and gluon confinement. They are the subject of intensive research thanks to non-perturbative methods based on Dyson-Schwinger (DS) and functional renormalization group (FRG) equations. Moreover, their temperature behavior might also help to explore the chiral and deconfinement phase transition or crossover within QCD at non-zero temperature. Our prime tool is the lattice discretized QCD (LQCD) providing a unique ab-initio non-perturbative approach to deal with the computation of various observables of the hadronic world. We investigate the temperature dependence of Landau gauge gluon and ghost propagators in pure gluodynamics and in full QCD. Regarding the gluon propagator, we compute its longitudinal DL as well its transversal DT components. The aim is to provide a data set in terms of fitting formulae which can be used as input for DS (or FRG) equations. We deal with full (Nf=2) LQCD with the twisted mass fermion discretization. We employ gauge field configurations provided by the tmfT collaboration for temperatures in the crossover region and for three fixed pion mass values in the range [300,500] MeV. Finally, within SU(3) pure gauge theory (at T=0) we compute the Landau gauge gluon propagator according to different gauge fixing criteria. Our goal is to understand the influence of gauge copies with minimal (non-trivial) eigenvalues of the Faddeev-Popov operator.
|
50 |
Investigating the large N limit of SU(N) Yang-Mills gauge theories on the latticeGarcía Vera, Miguel Francisco 02 August 2017 (has links)
In dieser Arbeit praesentieren wir Resultate der topologischen Suszeptibilitaet “chi” und untersuchen die Faktorisierung der reinen SU(N) Yang-Mills Eichtheorie im 't Hooft'schen Grenzwert grosser N. Ein entscheidender Teil der Berechnung von chi in der Gittereichtheorie ist die Abschaetzung des topologischen Ladungsdichtekorrelators, die durch ein schlechtes Signal-Rausch- Verhaeltnis beeintraechtigt ist. Um dieses Problem abzuschwaechen, fuehren wir einen neuen, auf einem mehrstufigen Vorgehen beruhenden Algorithmus ein, um die Korrelationsfunktion von Observablen zu berechnen, die mit dem Yang-Mills Gradientenfluss geglaettet wurden. Angewandt auf unsere Observablen, erhalten wir Ergebnisse, deren Fehlerskalierung besser ist, als die von herkoemmlichen Monte-Carlo Simulationen.
Wir bestimmen die topologische Suszeptibilitaet in der reinen Yang-Mills Eichtheorie fuer Eichgruppen mit N = 4,5,6 und drei verschiedenen Gitterabstaenden. Um das Einfrieren der Topologie zu umgehen, wenden wir offene Randbedingungen an. Zusaetzlich wenden wir die korrekte Definition der topologischen Ladungsdichte durch den Gradientenfluss an. Unser Endresultat im des Grenzfalls von grossen N repraesentiert eine neue Qualitaet in der Verifikation der Witten-Veneziano Formel.
Schliesslich benutzen wir die Gitterformulierung, um die Erwartungswertfaktorisierung des Produkts eichinvarianter Operatoren im Grenzwert grosser N zu verifizieren. Wir arbeiten mit durch den Yang-Mills Grandientenfluss geglaetteten Wilsonschleifen und Simulationen bis zur Eichgruppe SU(8). Die Extrapolationen zu grossen N sind in Ueberstimmung mit der Faktorisierung sowohl fuer endlichen Gitterabstand als auch in Kontinnumslimes. Unsere Daten erlauben uns nicht nur die Verifizierung der Faktorisierung, sondern auch einen hochpraezisen Test des 1/N Skalierungsverhaltens. Hier konnten wir das quadratische Skalierungsverhalten in 1/N finden, welches von 't Hooft vorhergesagt wurde. / In this thesis we present results for the topological susceptibility “chi”, and investigate the property of factorization in the 't Hooft large N limit of SU(N) pure Yang-Mills gauge theory. A key component in the lattice gauge theory computation of chi is the estimation of the topological charge density correlator, which is affected by a severe signal to noise problem. To alleviate this problem, we introduce a novel algorithm that uses a multilevel type approach to compute the correlation function of observables smoothed with the Yang-Mills gradient flow. When applied to our observables, the results show an scaling of the error which is better than the one of standard Monte-Carlo simulations.
We compute the topological susceptibility in the pure Yang-Mills gauge theory for the gauge groups with N = 4, 5, 6 and three different lattice spacings. In order to deal with the freezing of topology, we use open boundary conditions. In addition, we employ the theoretically sound definition of the topological charge density through the gradient flow. Our final result in the limit N to infinity, represents a new quality in the verification of the Witten-Veneziano formula.
Lastly, we use the lattice formulation to verify the factorization of the expectation value of the product of gauge invariant operators in the large N limit. We work with Wilson loops smoothed with the Yang-Mills gradient flow and simulations up to the gauge group SU(8). The large N extrapolations at finite lattice spacing and in the continuum are compatible with factorization. Our data allow us not only to verify factorization, but also to test the 1/N scaling up to very high precision, where we find it to agree very well with a quadratic series in 1/N as predicted originally by 't Hooft for the case of the pure Yang-Mills gauge theory.
|
Page generated in 0.0551 seconds