Determining the gluonic gravitational form factors of the proton

The proton is one of the main building blocks of all visible matter in the Universe¹. Among its intrinsic properties are its electric charge, mass and spin². These properties emerge from the complex dynamics of its fundamental constituents-quarks and gluons-described by the theory of quantum chromodynamics^3-5. The electric charge and spin of protons, which are shared among the quarks, have been investigated previously using electron scattering². An example is the highly precise measurement of the electric charge radius of the proton⁶. By contrast, little is known about the inner mass density of the proton, which is dominated by the energy carried by gluons. Gluons are hard to access using electron scattering because they do not carry an electromagnetic charge. Here we investigated the gravitational density of gluons using a small colour dipole, through the threshold photoproduction of the J/ψ particle. We determined the gluonic gravitational form factors of the proton^7,8 from our measurement. We used a variety of models^9-11 and determined, in all cases, a mass radius that is notably smaller than the electric charge radius. In some, but not all cases, depending on the model, the determined radius agrees well with first-principle predictions from lattice quantum chromodynamics¹². This work paves the way for a deeper understanding of the salient role of gluons in providing gravitational mass to visible matter.

The Higgs boson turns ten

The discovery of the Higgs boson, ten years ago, was a milestone that opened the door to the study of a new sector of fundamental physical interactions. We review the role of the Higgs field in the Standard Model of particle physics and explain its impact on the world around us. We summarize the insights into Higgs physics revealed so far by ten years of work, discuss what remains to be determined and outline potential connections of the Higgs sector with unsolved mysteries of particle physics.

Operator-Algebraic Renormalization and Wavelets

We report on a rigorous operator-algebraic renormalization group scheme and construct the free field with a continuous action of translations as the scaling limit of Hamiltonian lattice systems using wavelet theory. A renormalization group step is determined by the scaling equation identifying lattice observables with the continuum field smeared by compactly supported wavelets. Causality follows from Lieb-Robinson bounds for harmonic lattice systems. The scheme is related with the multiscale entanglement renormalization ansatz and augments the semicontinuum limit of quantum systems.

Fresh Extraction of the Proton Charge Radius from Electron Scattering

We present a novel method for extracting the proton radius from elastic electron-proton (ep) scattering data. The approach is based on interpolation via continued fractions augmented by statistical sampling and avoids any assumptions on the form of function used for the representation of data and subsequent extrapolation onto Q^{2}≃0. Applying the method to extant modern ep datasets, we find that all results are mutually consistent and, combining them, we arrive at r_{p}=0.847(8)  fm. This result compares favorably with values obtained from contemporary measurements of the Lamb shift in muonic hydrogen, transitions in electronic hydrogen, and muonic deuterium spectroscopy.

Resonance Electroproduction and the Origin of Mass

One of the greatest challenges within the Standard Model is to discover the source of visible mass. Indeed, this is the focus of a “Millennium Problem”, posed by the Clay Mathematics Institute. The answer is hidden within quantum chromodynamics (QCD); and it is probable that revealing the origin of mass will also explain the nature of confinement. In connection with these issues, this perspective will describe insights that have recently been drawn using contemporary methods for solving the continuum bound-state problem in relativistic quantum field theory and how they have been informed and enabled by modern experiments on nucleon-resonance electroproduction.

Empirical Consequences of Emergent Mass

The Lagrangian that defines quantum chromodynamics (QCD), the strong interaction piece of the Standard Model, appears very simple. Nevertheless, it is responsible for an astonishing array of high-level phenomena with enormous apparent complexity, e.g., the existence, number and structure of atomic nuclei. The source of all these things can be traced to emergent mass, which might itself be QCD’s self-stabilising mechanism. A background to this perspective is provided, presenting, inter alia, a discussion of the gluon mass and QCD’s process-independent effective charge and highlighting an array of observable expressions of emergent mass, ranging from its manifestations in pion parton distributions to those in nucleon electromagnetic form factors.

Quantum Gravity on the Computer: Impressions of a Workshop

Computer simulations allow us to explore non-perturbative phenomena in physics. This has the potential to help us understand quantum gravity. Finding a theory of quantum gravity is a hard problem, but, in the last several decades, many promising and intriguing approaches that utilize or might benefit from using numerical methods were developed. These approaches are based on very different ideas and assumptions, yet they face the common challenge to derive predictions and compare them to data. In March 2018, we held a workshop at the Nordic Institute for Theoretical Physics (NORDITA) in Stockholm gathering experts in many different approaches to quantum gravity for a workshop on “Quantum gravity on the computer”. In this article, we try to encapsulate some of the discussions held and talks given during this workshop and combine them with our own thoughts on why and how numerical approaches will play an important role in pushing quantum gravity forward. The last section of the article is a road map providing an outlook of the field and some intentions and goalposts that were debated in the closing session of the workshop. We hope that it will help to build a strong numerical community reaching beyond single approaches to combine our efforts in the search for quantum gravity.

Decoding the phase structure of QCD via particle production at high energy

Recent studies based on lattice Monte Carlo simulations of quantum chromodynamics (QCD)-the theory of strong interactions-have demonstrated that at high temperature there is a phase change from confined hadronic matter to a deconfined quark-gluon plasma in which quarks and gluons can travel distances that greatly exceed the size of hadrons. Here we show that the phase structure of such strongly interacting matter can be decoded by analysing particle production in high-energy nuclear collisions within the framework of statistical hadronization, which accounts for the thermal distribution of particle species. Our results represent a phenomenological determination of the location of the phase boundary of strongly interacting matter, and imply quark-hadron duality at this boundary.

A per-cent-level determination of the nucleon axial coupling from quantum chromodynamics

The axial coupling of the nucleon, g_A, is the strength of its coupling to the weak axial current of the standard model of particle physics, in much the same way as the electric charge is the strength of the coupling to the electromagnetic current. This axial coupling dictates the rate at which neutrons decay to protons, the strength of the attractive long-range force between nucleons and other features of nuclear physics. Precision tests of the standard model in nuclear environments require a quantitative understanding of nuclear physics that is rooted in quantum chromodynamics, a pillar of the standard model. The importance of g_A makes it a benchmark quantity to determine theoretically-a difficult task because quantum chromodynamics is non-perturbative, precluding known analytical methods. Lattice quantum chromodynamics provides a rigorous, non-perturbative definition of quantum chromodynamics that can be implemented numerically. It has been estimated that a precision of two per cent would be possible by 2020 if two challenges are overcome^1,2: contamination of g_A from excited states must be controlled in the calculations and statistical precision must be improved markedly^2-10. Here we use an unconventional method ¹¹ inspired by the Feynman-Hellmann theorem that overcomes these challenges. We calculate a g_A value of 1.271 ± 0.013, which has a precision of about one per cent.

Fluctuations of conserved charges from imaginary chemical potential

When comparing lattice calculation to experimental data from heavy ion collision experiments, the higher order fluctuations of conserved charges are important observables. An efficient way to study these fluctuations is to determine them from simulations at imaginary chemical potential. In this talk we present results up to the six order derivative in μB (with up to eighth order included in the fit), calculated on a 483 × 12 lattice with staggered fermions using different values of μB while μS = μQ = 0.

Single flavour filtering for RHMC in BQCD

Filtering algorithms for two degenerate quark flavours have advanced to the point that, in 2+1 flavour simulations, the cost of the strange quark is significant compared with the light quarks. This makes efficient filtering algorithms for single flavour actions highly desirable, in particular when considering 1+1+1 flavour simulations for QED+QCD. Here we discuss methods for filtering the RHMC algorithm that are implemented within BQCD, an open-source Fortran program for Hybrid Monte Carlo simulations.

Hadron thermodynamics from imaginary chemical potentials

We use 4stout improved staggered lattice data at imaginary chemical potentials to calculate fugacity expansion coefficients in finite temperature QCD. We discuss the phenomenological interpretation of our results within the hadron resonance gas (HRG) model, and the hints they give us about the hadron spectrum. We also discuss features of the higher order coefficients that are not captured by the HRG. This conference contribution is based on our recent papers [1, 2].,

Lattice QCD results for the HVP contribution to the anomalous magnetic moments of leptons

We present lattice QCD results by the Budapest-Marseille-Wuppertal (BMW) Collaboration for the leading-order contribution of the hadron vacuum polarization (LOHVP) to the anomalous magnetic moments of all charged leptons. Calculations are performed with u, d, s and c quarks at their physical masses, in volumes of linear extent larger than 6 fm, and at six values of the lattice spacing, allowing for controlled continuum extrapolations. All connected and disconnected contributions are calculated for not only the muon but also the electron and tau anomalous magnetic moments. Systematic uncertainties are thoroughly discussed and comparisons with other calculations and phenomenological estimates are made.

Fundamental properties of resonances

All resonances, from hydrogen nuclei excited by the high-energy gamma rays in deep space to newly discovered particles produced in Large Hadron Collider, should be described by the same fundamental physical quantities. However, two distinct sets of properties are used to describe resonances: the pole parameters (complex pole position and residue) and the Breit-Wigner parameters (mass, width, and branching fractions). There is an ongoing decades-old debate on which one of them should be abandoned. In this study of nucleon resonances appearing in the elastic pion-nucleon scattering we discover an intricate interplay of the parameters from both sets, and realize that neither set is completely independent or fundamental on its own.

Up and Down Quark Masses and Corrections to Dashen's Theorem from Lattice QCD and Quenched QED

In a previous Letter [Borsanyi et al., Phys. Rev. Lett. 111, 252001 (2013)] we determined the isospin mass splittings of the baryon octet from a lattice calculation based on N_{f}=2+1 QCD simulations to which QED effects have been added in a partially quenched setup. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. Our ensembles include 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm, and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashen's theorem, we obtain ϵ=0.73(2)(5)(17), where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, m_{u}=2.27(6)(5)(4) and m_{d}=4.67(6)(5)(4)  MeV in the modified minimal subtraction scheme at 2  GeV and the isospin breaking ratios m_{u}/m_{d}=0.485(11)(8)(14), R=38.2(1.1)(0.8)(1.4), and Q=23.4(0.4)(0.3)(0.4). Our results exclude the m_{u}=0 solution to the strong CP problem by more than 24 standard deviations.

Lattice Computation of the Nucleon Scalar Quark Contents at the Physical Point

We present a QCD calculation of the u, d, and s scalar quark contents of nucleons based on 47 lattice ensembles with N_{f}=2+1 dynamical sea quarks, 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm, and pion masses down to 120 MeV. Using the Feynman-Hellmann theorem, we obtain f_{ud}^{N}=0.0405(40)(35) and f_{s}^{N}=0.113(45)(40), which translates into σ_{πN}=38(3)(3)  MeV, σ_{sN}=105(41)(37)  MeV, and y_{N}=0.20(8)(8) for the sigma terms and the related ratio, where the first errors are statistical and the second errors are systematic. Using isospin relations, we also compute the individual up and down quark contents of the proton and neutron (results in the main text).

Non-Abelian SU(2) Lattice Gauge Theories in Superconducting Circuits

We propose a digital quantum simulator of non-Abelian pure-gauge models with a superconducting circuit setup. Within the framework of quantum link models, we build a minimal instance of a pure SU(2) gauge theory, using triangular plaquettes involving geometric frustration. This realization is the least demanding, in terms of quantum simulation resources, of a non-Abelian gauge dynamics. We present two superconducting architectures that can host the quantum simulation, estimating the requirements needed to run possible experiments. The proposal establishes a path to the experimental simulation of non-Abelian physics with solid-state quantum platforms.

QCD and strongly coupled gauge theories: challenges and perspectives

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, [Formula: see text]- and [Formula: see text]-meson physics with the aim of making them easily accessible to the particle-physics community. More specifically, we report on the determination of the light-quark masses, the form factor [Formula: see text], arising in semileptonic [Formula: see text] transition at zero momentum transfer, as well as the decay-constant ratio [Formula: see text] of decay constants and its consequences for the CKM matrix elements [Formula: see text] and [Formula: see text]. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of [Formula: see text] and [Formula: see text] Chiral Perturbation Theory and review the determination of the [Formula: see text] parameter of neutral kaon mixing. The inclusion of heavy-quark quantities significantly expands the FLAG scope with respect to the previous review. Therefore, we focus here on [Formula: see text]- and [Formula: see text]-meson decay constants, form factors, and mixing parameters, since these are most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. In addition we review the status of lattice determinations of the strong coupling constant [Formula: see text].

Matrix product states for gauge field theories

The matrix product state formalism is used to simulate Hamiltonian lattice gauge theories. To this end, we define matrix product state manifolds which are manifestly gauge invariant. As an application, we study (1+1)-dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model, and are able to determine very accurately the ground-state properties and elementary one-particle excitations in the continuum limit. In particular, a novel particle excitation in the form of a heavy vector boson is uncovered, compatible with the strong coupling expansion in the continuum. We also study full quantum nonequilibrium dynamics by simulating the real-time evolution of the system induced by a quench in the form of a uniform background electric field.

Astronomical reach of fundamental physics

Using basic physical arguments, we derive by dimensional and physical analysis the characteristic masses and sizes of important objects in the universe in terms of just a few fundamental constants. This exercise illustrates the unifying power of physics and the profound connections between the small and the large in the cosmos we inhabit. We focus on the minimum and maximum masses of normal stars, the corresponding quantities for neutron stars, the maximum mass of a rocky planet, the maximum mass of a white dwarf, and the mass of a typical galaxy. To zeroth order, we show that all these masses can be expressed in terms of either the Planck mass or the Chandrasekar mass, in combination with various dimensionless quantities. With these examples, we expose the deep interrelationships imposed by nature between disparate realms of the universe and the amazing consequences of the unifying character of physical law.

Measurement of the 1H(γ, p)π0 reaction using a novel nucleon spin polarimeter

We report the first large-acceptance measurement of polarization transfer from a polarized photon beam to a recoiling nucleon. The measurement pioneers a novel polarimetry technique, which can be applied to many other nuclear and hadron physics experiments. The commissioning reaction of 1H(γ, p)π0 in the range 0.4<Eγ<1.4  GeV validates the technique and provides essential new data to constrain the excitation spectrum of the nucleon.

Isospin splittings in the light-baryon octet from lattice QCD and QED

While electromagnetic and up-down quark mass difference effects on octet baryon masses are very small, they have important consequences. The stability of the hydrogen atom against beta decay is a prominent example. Here, we include these effects by adding them to valence quarks in a lattice QCD calculation based on Nf=2+1 simulations with five lattice spacings down to 0.054 fm, lattice sizes up to 6 fm, and average up-down quark masses all the way down to their physical value. This allows us to gain control over all systematic errors, except for the one associated with neglecting electromagnetism in the sea. We compute the octet baryon isomultiplet mass splittings, as well as the individual contributions from electromagnetism and the up-down quark mass difference. Our results for the total splittings are in good agreement with experiment.

Model-independent extraction of the pole and Breit-Wigner resonance parameters

We show that a slightly modified Breit-Wigner formula can successfully describe the total cross section even for the broad resonances, from the light ρ(770) to the heavy Z boson. In addition to the mass, width, and branching fraction, we include another resonance parameter that turns out to be directly related to the pole residue phase. The new formula has two mathematically equivalent forms: one with the pole and the other with the Breit-Wigner parameters.

Freeze-out parameters: lattice meets experiment

We present our results for ratios of higher order fluctuations of electric charge as functions of the temperature. These results are obtained in a system of 2+1 quark flavors at physical quark masses and continuum extrapolated. We compare them to preliminary data on higher order moments of the net electric charge distribution from the STAR collaboration. This allows us to determine the freeze-out temperature and chemical potential from first principles. We also show continuum-extrapolated results for ratios of higher order fluctuations of baryon number. These will allow us to test the consistency of the approach, by comparing them to the corresponding experimental data (once they become available) and thus, extracting the freeze-out parameters in an independent way.

Discrete differential geometry: the nonplanar quadrilateral mesh

We consider the problem of constructing a discrete differential geometry defined on nonplanar quadrilateral meshes. Physical models on discrete nonflat spaces are of inherent interest, as well as being used in applications such as computation for electromagnetism, fluid mechanics, and image analysis. However, the majority of analysis has focused on triangulated meshes. We consider two approaches: discretizing the tensor calculus, and a discrete mesh version of differential forms. While these two approaches are equivalent in the continuum, we show that this is not true in the discrete case. Nevertheless, we show that it is possible to construct mesh versions of the Levi-Civita connection (and hence the tensorial covariant derivative and the associated covariant exterior derivative), the torsion, and the curvature. We show how discrete analogs of the usual vector integral theorems are constructed in such a way that the appropriate conservation laws hold exactly on the mesh, rather than only as approximations to the continuum limit. We demonstrate the success of our method by constructing a mesh version of classical electromagnetism and discuss how our formalism could be used to deal with other physical models, such as fluids.

Is the proton radius a player in the redefinition of the International System of Units?

It is now recognized that the International System of Units (SI units) will be redefined in terms of fundamental constants, even if the date when this will occur is still under debate. Actually, the best estimate of fundamental constant values is given by a least-squares adjustment, carried out under the auspices of the Committee on Data for Science and Technology (CODATA) Task Group on Fundamental Constants. This adjustment provides a significant measure of the correctness and overall consistency of the basic theories and experimental methods of physics using the values of the constants obtained from widely differing experiments. The physical theories that underlie this adjustment are assumed to be valid, such as quantum electrodynamics (QED). Testing QED, one of the most precise theories is the aim of many accurate experiments. The calculations and the corresponding experiments can be carried out either on a boundless system, such as the electron magnetic moment anomaly, or on a bound system, such as atomic hydrogen. The value of fundamental constants can be deduced from the comparison of theory and experiment. For example, using QED calculations, the value of the fine structure constant given by the CODATA is mainly inferred from the measurement of the electron magnetic moment anomaly carried out by Gabrielse's group. (Hanneke et al. 2008 Phys. Rev. Lett. 100, 120801) The value of the Rydberg constant is known from two-photon spectroscopy of hydrogen combined with accurate theoretical quantities. The Rydberg constant, determined by the comparison of theory and experiment using atomic hydrogen, is known with a relative uncertainty of 6.6×10(-12). It is one of the most accurate fundamental constants to date. A careful analysis shows that knowledge of the electrical size of the proton is nowadays a limitation in this comparison. The aim of muonic hydrogen spectroscopy was to obtain an accurate value of the proton charge radius. However, the value deduced from this experiment contradicts other less accurate determinations. This problem is known as the proton radius puzzle. This new determination of the proton radius may affect the value of the Rydberg constant . This constant is related to many fundamental constants; in particular, links the two possible ways proposed for the redefinition of the kilogram, the Avogadro constant N(A) and the Planck constant h. However, the current relative uncertainty on the experimental determinations of N(A) or h is three orders of magnitude larger than the 'possible' shift of the Rydberg constant, which may be shown by the new value of the size of the proton radius determined from muonic hydrogen. The proton radius puzzle will not interfere in the redefinition of the kilogram. After a short introduction to the properties of the proton, we will describe the muonic hydrogen experiment. There is intense theoretical activity as a result of our observation. A brief summary of possible theoretical explanations at the date of writing of the paper will be given. The contribution of the proton radius puzzle to the redefinition of SI-based units will then be examined.

Δ(1232) axial charge and form factors from lattice QCD

We present the first calculation on the Δ axial vector and pseudoscalar form factors using lattice QCD. Two Goldberger-Treiman relations are derived and examined. A combined chiral fit is performed to the nucleon axial charge, N to Δ axial transition coupling constant and Δ axial charge.

Nuclear physics from lattice QCD at strong coupling

We study numerically the strong coupling limit of lattice QCD with one flavor of massless staggered quarks. We determine the complete phase diagram as a function of temperature and chemical potential, including a tricritical point. We clarify the nature of the low temperature dense phase, which is strongly bound "nuclear" matter. This strong binding is explained by the nuclear potential, which we measure. Finally, we determine, from this first-principles limiting case of QCD, the masses of "atomic nuclei" up to A=12 "carbon".

Disentangling the dynamical origin of P11 nucleon resonances

We show that two almost degenerate poles near the piDelta threshold and the next higher mass pole in the P11 partial wave of piN scattering evolve from a single bare state through its coupling with piN, etaN, and pipiN reaction channels. This finding provides new information on understanding the dynamical origins of the Roper N{*}(1440) and N{*}(1710) resonances listed by Particle Data Group. Our results for the resonance poles in other piN partial waves are also presented.