Electromagnetic form factors and radii of the nucleon from Lattice QCD and the proton radius puzzles

Abstract

In this thesis, we investigate the electromagnetic form factors of the proton and neutron in the framework of Quantum Chromodynamics (QCD). To perform a calculation of such low-energy quantities based on first principles, we employ the lattice regularization of QCD. In Lattice-QCD simulations of nucleon-structure observables, systematic errors are inherent due to the finite lattice spacing and volume as well as due to contamination by excited states. These can be controlled and removed in a systematic fashion: for the removal of excited-state contributions, a variety of dedicated methods exists, while discretization and finite-volume effects also need to be taken into account by performing a continuum and infinite-volume extrapolation. However, all previous lattice studies of the electromagnetic form factors of the nucleon have either neglected the numerically challenging quark-disconnected contributions or were not extrapolated to the continuum and infinite-volume limit.

We present results for the electromagnetic form factors of the proton and neutron computed on the (2 + 1)-flavor Coordinated Lattice Simulations (CLS) ensembles including both quark-connected and -disconnected contributions while, at the same time, controlling all sources of systematic uncertainties. For the excited-state analysis, we explore three complementary methods based on two-state fits to the effective form factors and on two different truncations of the summation method, respectively. The Q^2-, pion-mass, lattice-spacing and finite-volume dependence of our form factor data is fitted simultaneously to the expressions resulting from covariant baryon chiral perturbation theory including vector mesons amended by models for lattice artefacts. From these fits, we determine the electric, magnetic, Zemach and Friar radii as well as the magnetic moments of the proton and neutron. To assess the influence of systematic effects, we average over various cuts in the pion mass and the momentum transfer, as well as over different models for the lattice-spacing and finite-volume dependence, using weights derived from the Akaike Information Criterion (AIC).

Our ab-initio QCD results for the electromagnetic radii of the proton are of particular relevance in light of the so-called proton radius puzzle, i.e., the observation of a large tension between different experimental measurements of the proton's electric radius, which is after more than a decade of vigorous research still not completely explained. Also for the magnetic radius, analyses based on different data sets find discrepant results. In this situation, a firm theoretical prediction of the proton radii can contribute towards the clarification of the origins of the discrepancies.

Our results for the magnetic moments of the proton and neutron are in good agreement with the experimentally very precisely known values. For the radii of the proton, we achieve, including all systematic errors, a precision which enables a meaningful comparison to the various experiments and data-driven evaluations. On the one hand, our result for the electric radius of the proton clearly points towards a small value, as favored by muonic hydrogen spectroscopy, the recent ep-scattering experiment by PRad and data-driven dispersive analyses. Our estimate for the magnetic radius, on the other hand, is well compatible with that inferred from the A1 ep-scattering experiment by a z-expansion analysis and in tension with z-expansion on the remaining world data as well as with dispersive approaches.

In dieser Doktorarbeit werden die elektromagnetischen Formfaktoren des Protons und Neutrons im Rahmen der Quantenchromodynamik (QCD) untersucht. Um solche Niederenergiegrößen basierend auf ersten Prinzipien zu berechnen, wird die Gitterregularisierung der QCD verwendet. In Gitter-QCD-Simulationen von Nukleonstruktur-Observablen sind systematische Fehler inhärent aufgrund des endlichen Gitterabstands und -volumens sowie aufgrund der Kontamination durch angeregte Zustände. Diese können kontrolliert und auf systematische Weise entfernt werden: Für die Entfernung der Beiträge angeregter Zustände existieren eine Reihe dedizierter Methoden, während Diskretisierungs- und Endlich-Volumen-Effekte ebenfalls berücksichtigt werden müssen, indem man eine Extrapolation ins Kontinuum und zu unendlichem Volumen durchführt. Jedoch haben alle bisherigen Gitterstudien der elektromagnetischen Formfaktoren des Nukleons entweder die numerisch herausfordernden quark-unverbundenen Beiträge vernachlässigt oder waren nicht zum Kontinuums- und Unendlich-Volumen-Limes extrapoliert.

Wir präsentieren Resultate für die elektromagnetischen Formfaktoren des Protons und Neutrons, die auf den (2 + 1)-flavour Coordinated Lattice Simulations (CLS) Ensembles berechnet wurden. Sie berücksichtigen sowohl quark-verbundene als auch -unverbundene Beiträge und kontrollieren gleichzeitig alle Quellen systematischer Unsicherheiten. Für die Analyse angeregter Zustände werden drei komplementäre Methoden erforscht, die auf Zweizustandsfits an die effektiven Formfaktoren bzw. auf zwei verschiedenen Trunkierungen der Summationsmethode basieren. Die Q^2-, Pionmassen-, Gitterabstands- und Volumen-Abhängigkeit unserer Formfaktor-Daten wird simultan an die Ausdrücke gefittet, die aus kovarianter baryonischer chiraler Störungstheorie inklusive Vektormesonen resultieren, und die um Modelle für Gitterartefakte ergänzt werden. Mithilfe dieser Fits werden die elektrischen, magnetischen, Zemach- und Friar-Radien sowie die magnetischen Momente des Protons und Neutrons bestimmt. Um die Einflüsse systematischer Effekte zu beurteilen, wird über diverse Cuts in der Pionmasse und im Impulsübertrag sowie über verschiedene Modelle für die Gitterabstands- und Volumen-Abhängigkeit gemittelt. Hierzu werden Gewichte verwendet, die aus dem Akaike-Informationskriterium (AIC) hergeleitet sind.

Unsere ab-initio-Resultate für die elektromagnetischen Radien des Protons sind von besonderer Bedeutung im Lichte des sogenannten Proton-Radius-Puzzles, d.h. der großen Spannung zwischen verschiedenen experimentellen Messungen des elektrischen Radius des Protons, die nach mehr als einer Dekade intensiver Forschung immer noch nicht vollständig erklärt ist. Auch für den magnetischen Radius finden Analysen basierend auf verschiedenen Datensätzen diskrepante Ergebnisse. In dieser Situation kann eine solide theoretische Vorhersage der Protonradien dazu beitragen, die Ursprünge der Diskrepanzen zu klären.

Unsere Ergebnisse für die magnetischen Momente des Protons und Neutrons befinden sich in guter Übereinstimmung mit den experimentell sehr genau bekannten Werten. Für die Radien des Protons erreichen wir inklusive aller systematischer Fehler eine Präzision, die einen aussagekräftigen Vergleich mit den diversen Experimenten und datengetriebenen Bestimmungen erlaubt. Einerseits deutet unser Resultat für den elektrischen Radius des Protons klar auf einen kleinen Wert hin, wie er von der Spektroskopie myonischen Wasserstoffs, dem kürzlich erfolgten ep-Streuexperiment von PRad und datengetriebenen dispersiven Analysen bevorzugt wird. Andererseits ist unser Ergebnis für den magnetischen Radius gut kompatibel mit dem mithilfe einer z-expansion-Analyse aus dem A1 ep-Streuexperiment abgeleiteten Wert und in Spannung mit z-expansion auf den verbleibenden Welt-Daten sowie mit dispersiven Ansätzen.