Search

We present our lattice QCD result for the long-distance part of the hadronic vacuum polarization contribution, $(a_\mu^{\rm hvp})^{\rm LD}$, to the muon $g-2$ in the time-momentum representation. This is the numerically dominant, and at the same time the most challenging part regarding statistical precision. Our calculation is based on ensembles with dynamical up, down and strange quarks, employing the O($a$)-improved Wilson fermion action with lattice spacings ranging from $0.035-0.099$ fm. In order to reduce statistical noise in the long-distance part of the correlator to the per-mille level, we apply low-mode averaging and combine it with an explicit spectral reconstruction. Our result is $(a_\mu^{\rm hvp})^{\rm LD} = 423.2(4.2)_{\rm stat}(3.4)_{\rm syst}\times 10^{-10}$ in isospin-symmetric QCD, where the pion decay constant is used to set the energy scale. When combined with our previous results for the short- and intermediate-distance window observables and after including all sub-dominant contributions as well as isospin-breaking corrections, we obtain the total leading-order hadronic vacuum polarization contribution as $a_\mu^{\rm hvp} = 724.9(5.0)_{\rm stat}(4.9)_{\rm syst}\times 10^{-10}$. Our result displays a tension of 3.9 standard deviations with the data-driven estimate published in the 2020 White Paper, but leads to a SM prediction for the total muon anomalous magnetic moment that agrees with the current experimental average., Comment: 55 pages, 14 figures, 8 tables

We present an update of our determination of the isovector charges $g_A^{u-d}$, $g_S^{u-d}$ and $g_T^{u-d}$, and the isovector twist-2 forward matrix elements $\langle x\rangle_{u-d}$, $\langle x\rangle_{\Delta u-\Delta d}$ and $\langle x\rangle_{\delta u-\delta d}$ on the $N_\mathrm{f}=2+1$ gauge ensembles generated by the Coordinated Lattice Simulations (CLS) effort. We have significantly extended our coverage of the parameter space by adding ensembles at the physical pion mass and fine lattice spacing, at nearly-physical pion masses and very fine lattice spacings, and at very large physical lattice volumes, enabling a well-controlled extrapolation to the physical point. Another major improvement is achieved owing to the extended range of source-sink separations, which allows us to perform two-state fits to summed correlator ratios, leading to a much higher level of control over excited-state effects. Systematic uncertainties from the chiral, continuum and infinite-volume extrapolations are incorporated via model averages based on the Akaike Information Criterion. Our final results at the physical point are $g_A^{u-d} = 1.254(19)_\mathrm{stat}(15)_\mathrm{sys}[24]_\mathrm{total}$, $g_S^{u-d} = 1.203(77)_\mathrm{stat}(81)_\mathrm{sys}[112]_\mathrm{total}$, $g_T^{u-d} = 0.993(15)_\mathrm{stat}(05)_\mathrm{sys}[16]_\mathrm{total}$, $\langle x\rangle_{u-d} = 0.153(15)_\mathrm{stat}(10)_\mathrm{sys}[17]_\mathrm{total}$, $\langle x\rangle_{\Delta u - \Delta d} = 0.207(15)_\mathrm{stat}(06)_\mathrm{sys}[16]_\mathrm{total}$, and $\langle x\rangle_{\delta u - \delta d} = 0.195(17)_\mathrm{stat}(15)_\mathrm{sys}[23]_\mathrm{total}$. While our results for the isovector charges are in excellent agreement with the FLAG\,21 averages, we note that our error for the tensor charge $g_T^{u-d}$ is considerably smaller., Comment: 20 pages, 11 figures

We present results for the short-distance window observable of the hadronic vacuum polarization contribution to the muon $g-2$, computed via the time-momentum representation (TMR) in lattice QCD. A key novelty of our calculation is the reduction of discretization effects by a suitable subtraction applied to the TMR kernel function, which cancels the leading $x_0^4$-behaviour at short distances. To compensate for the subtraction, one must substitute a term that can be reliably computed in perturbative QCD. We apply this strategy to our data for the vector current collected on ensembles generated with $2+1$ flavours of O($a$)-improved Wilson quarks at six values of the lattice spacing and pion masses in the range $130-420\,$MeV. Our estimate at the physical point contains a full error budget and reads $(a_\mu^{\rm hvp})^{\rm SD}=68.85(14)_{\rm stat}\,(42)_{\rm syst}\cdot10^{-10}$, which corresponds to a relative precision of 0.7\%. We discuss the implications of our result for the observed tensions between lattice and data-driven evaluations of the hadronic vacuum polarization., Comment: 37 pages, 10 figures, 9 tables

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. 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 chiral perturbation theory including vector mesons amended by models for lattice artefacts. From these fits, we determine the electric and magnetic radii and the magnetic moments of the proton and neutron, as well as the Zemach radius of the proton. 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 results for the magnetic moments of the proton and neutron are in good agreement with the experimental values and have a relative precision of about $2.4\,\%$ and $3.7\,\%$, respectively. For the electromagnetic radii of the proton, we achieve a precision at the $1.5\,\%$ level., Comment: 10 pages, 4 figures, contribution to the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab

In this report we present the status of the Mainz group's lattice QCD calculation of the pion transition form factor, which describes the interaction of an on-shell pion with two off-shell photons. This form factor is the main ingredient in the calculation of the pion-pole contribution to hadronic light-by-light scattering in the muon $g-2$. We use the $N_f = 2 + 1$ CLS gauge ensembles, and we update our previous work by including a physical pion mass ensemble (E250). We compute the transition form factor in a moving frame as well as in the pion rest frame in order to have access to a wider range of photon virtualities. In addition to the quark-line connected correlator we also compute the quark-line disconnected diagrams that contribute to the form factor. At the final stage of the analysis, the result on E250 will be combined with the previous work published in 2019 to extrapolate the form factor to the continuum and to physical quark masses., Comment: Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

We present the first lattice-QCD result for the Zemach and Friar radii of the proton and neutron. Our calculation includes both quark-connected and -disconnected diagrams and assesses all sources of systematic uncertainties arising from excited-state contributions, finite-volume effects and the continuum extrapolation. At the physical point, we obtain for the proton $r_Z^p = ( 1.013 \pm 0.010\ (\mathrm{stat}) \pm 0.012\ (\mathrm{syst}) )~\mathrm{fm}$ and $r_F^p = ( 1.301 \pm 0.012\ (\mathrm{stat}) \pm 0.014\ (\mathrm{syst}) )~\mathrm{fm}$. These numbers suggest small values of the Zemach and Friar radii of the proton, but are compatible with most of the experimental studies., Comment: 13 pages, 2 figures; v2: significantly extended version including, in addition to the Zemach radius of the proton, results for the Friar radius of the proton as well as for the Zemach and Friar radii of the neutron; v3: typos fixed, matches published version

We present lattice-QCD results for the electromagnetic form factors of the proton and neutron including both quark-connected and -disconnected contributions. The parametrization of the $Q^2$-dependence of the form factors is combined with the extrapolation to the physical point. In this way, we determine the electric and magnetic radii and the magnetic moments of the proton and neutron. For the proton, we obtain at the physical pion mass and in the continuum and infinite-volume limit $\sqrt{\langle r_E^2 \rangle^p} = 0.820(14)$ fm, $\sqrt{\langle r_M^2 \rangle^p} = 0.8111(89)$ fm, and $\mu_M^p = 2.739(66)$, where the errors include all systematics., Comment: 8 pages, 3 figures; v2: matches published version; for the accompanying paper, see arXiv:2309.06590 [hep-lat]

There is a long-standing discrepancy between different measurements of the electric and magnetic radii of the proton. Lattice QCD calculations are a well-suited tool for theoretical investigations of the structure of the nucleon from first principles. However, all previous lattice studies of the proton's electromagnetic radii have either neglected quark-disconnected contributions or were not extrapolated to the continuum and infinite-volume limit. Here, 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. From simultaneous fits to the $Q^2$-, pion-mass, lattice-spacing, and finite-volume dependence of the form factors, we determine the electric and magnetic radii and the magnetic moments of the proton and neutron. For the proton, we obtain as our final values $\langle r_E^2 \rangle^p = (0.672 \pm 0.014$ (stat)${} \pm 0.018$ (syst)$)$ fm${}^2$, $\langle r_M^2 \rangle^p = (0.658 \pm 0.012$ (stat)${} \pm 0.008$ (syst)$)$ fm${}^2$, and $\mu_M^p = 2.739 \pm 0.063$ (stat)${} \pm 0.018$ (syst). The magnetic moment is in good agreement with the experimental value, as is the one of the neutron. On the one hand, our result for the electric (charge) radius of the proton clearly points towards a small value, as favored by muonic hydrogen spectroscopy and the recent $ep$-scattering experiment by PRad. Our estimate for the magnetic radius, on the other hand, is well compatible with that inferred from the A1 $ep$-scattering experiment., Comment: 58 pages, 19 figures; v2: significantly extended crosschecks of systematic effects, matches published version; for the accompanying letter, see arXiv:2309.07491 [hep-lat]

The problem of having to reconstruct the decay rates and corresponding amplitudes of the single-exponential components of a noisy multi-exponential signal is common in many other areas of physics and engineering besides lattice field theory, and it can be helpful to study the methods devised and used for that purpose in those contexts in order to get a better handle on the problem of extracting masses and matrix elements from lattice correlators. Here we consider the use of Pad\'e and Pad\'e-Laplace methods, which have found wide use in laser fluorescence spectroscopy and beyond, emphasizing the importance of using robust Pad\'e approximants to avoid spurious poles. To facilitate the accurate evaluation of the Laplace transform required for the Pad\'e-Laplace method, we also present a novel approach to the numerical quadrature of multi-exponential functions., Comment: 6 pages, 4 PDF figures; poster presented at 39th International Symposium on Lattice Field Theory (Lattice2022), 8-13 August, 2022, Bonn, Germany

We present a lattice-QCD based analysis of the nucleon sigma terms using gauge ensembles with $N_f = 2 + 1$ flavors of ${\cal O}(a)$-improved Wilson fermions, with a complete error budget concerning excited-state contaminations, the chiral interpolation as well as finite-size and lattice spacing effects. We compute the sigma terms determined directly from the matrix elements of the scalar currents. The chiral interpolation is based on SU(3) baryon chiral perturbation theory using the extended on-mass shell renormalization scheme. For the pion nucleon sigma term, we obtain $\sigma_{\pi N} = (43.7\pm3.6)$ MeV, where the error includes our estimate of the aforementioned systematics. The tension with extractions based on dispersion theory persists at the 2.4-$\sigma$ level. For the strange sigma term, we obtain a non-zero value, $\sigma_s=(28.6\pm9.3)$ MeV., Comment: Matches published version

We present an analysis of the pion-nucleon sigma term on the CLS ensembles with $N_f=2+1$ flavors of $\mathcal{O}(a)$-improved Wilson fermions. We perform a chiral interpolation based on ensembles with pion masses ranging from 130 MeV to roughly 350 MeV. The analysis covers four lattice spacings between $a\approx [0.05 \, \rm{fm}\ldots 0.09\, \rm{fm}]$, allowing for an estimate of systematics associated with lattice artefacts., Comment: Talk presented at the 39th International Symposium on Lattice Field Theory (Lattice 2022), 8-13 August, 2022, Bonn, Germany

We present the current status of our analysis of nucleon structure observables including isovector charges and twist-2 matrix elements as well as the nucleon mass. Results are computed on a large set of CLS $N_f=2+1$ gauge ensembles with $M_\pi\approx0.130\mathrm{MeV}\ldots 350\mathrm{MeV}$, four values of the lattice spacing $a=0.05\mathrm{fm}\ldots 0.09\mathrm{fm}$ and covering a large range of physical volumes. Compared to the results presented at last year's conference we have added data on a very fine and large box at small light quark mass ($T\times L^3=192 \times 96^3$, $M_\pi=172\mathrm{MeV}$, $a=0.05\mathrm{fm}$). Additional (intermediate) source-sink separations have been computed on the coarser ensembles, further increasing effective statistics and allowing for a more fine-grained control in the treatment of the excited-state contamination. Excited states in the nucleon matrix elements are tamed by a simultaneous, two-state fit ansatz using the summation method., Comment: 10 pages, 4 figures, Talk presented at the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

We present results for the hadronic running of the electromagnetic coupling and the weak mixing angle from simulations of lattice QCD with $N_f=2+1$ flavours of $O(a)$-improved Wilson fermions. Using two different discretisations of the vector current, we compute the quark-connected and -disconnected contributions to the hadronic vacuum polarisation (HVP) functions $\bar{\Pi}^{\gamma\gamma}$ and $\bar{\Pi}^{Z\gamma}$ for spacelike squared momenta $Q^2\leq 7$ $\mathrm{GeV}^2$. Our results are extrapolated to the physical point using ensembles at four lattice spacings, with pion masses ranging from 130 to 420 MeV. We observe a tension of up to 3.5 standard deviations between our lattice results for $\Delta\alpha_{\rm had}^{(5)}(-Q^2)$ and estimates based on the $\textit{R}$-ratio for space-like momenta in the range $Q^2=3-7\,\rm GeV^2$. To obtain an estimate for $\Delta\alpha_\mathrm{had}^{(5)}(M_Z^2)$, we employ the Euclidean split technique. The implications for comparison with global electroweak fits are assessed., Comment: 16 pages, 7 figures, proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

The isovector axial form factor of the nucleon plays a key role in interpreting data from long-baseline neutrino oscillation experiments. We present a lattice QCD calculation of this form factor, introducing a new method to directly extract its z-expansion from lattice correlators. Our final parameterization of the form factor, which extends up to spacelike virtualities of 0.7 GeV^2 with fully quantified uncertainties, agrees with previous lattice calculations but is significantly less steep than neutrino-deuterium scattering data suggests., Comment: Proceeding for the 39th International Symposium on Lattice Field Theory (Lattice2022), 8th-13th August 2022, Bonn, Germany

Following the publication of the new measurement of the anomalous magnetic moment of the muon, the discrepancy between experiment and the theory prediction from the $g-2$ theory initiative has increased to $4.2\,\sigma$. Recent lattice QCD calculations predict values for the hadronic vacuum polarization contribution that are larger than the data-driven estimates, bringing the Standard Model prediction closer to the experimental measurement. Euclidean time windows in the time-momentum representation of the hadronic vacuum polarization contribution to the muon $g-2$ can help clarify the discrepancy between the phenomenological and lattice predictions. We present our calculation of the intermediate distance window contribution using $N_\mathrm{f}=2+1$ flavors of O$(a)$ improved Wilson quarks. We employ ensembles at six lattice spacings below $0.1\,$fm and pion masses down to the physical value. We present a detailed study of the continuum limit, using two discretizations of the vector current and two independent sets of improvement coefficients. Our result at the physical point displays a tension of $3.9\,\sigma$ with a recent evaluation of the intermediate window based on the data-driven method., Comment: 10 pages, 5 figures, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

We present results for the electromagnetic form factors of the proton and neutron computed on the Coordinated Lattice Simulations (CLS) ensembles with $N_f = 2 + 1$ flavors of $\mathcal{O}(a)$-improved Wilson fermions and an $\mathcal{O}(a)$-improved conserved vector current. In order to estimate the excited-state contamination, we employ several source-sink separations and apply the summation method. The quark-disconnected diagrams entering the isoscalar quantities are computed explicitly. For this purpose, a stochastic estimation based on the one-end trick is performed, in combination with a frequency-splitting technique and the hopping-parameter expansion. By these means, we obtain a clear signal for the form factors including the quark-disconnected contributions, which have a statistically significant effect on our results. From the $Q^2$-dependence of the form factors, we determine the electric and magnetic charge radii and the magnetic moments of the proton and neutron. The chiral interpolation is carried out by simultaneously fitting the pion mass and $Q^2$-dependence of our form factor data directly to the expressions resulting from covariant chiral perturbation theory including vector mesons. 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., Comment: 10 pages, 6 figures, contribution to the 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany; v2: fixed a bug in the fit code, slightly different results for the neutron, results for the proton unchanged. arXiv admin note: substantial text overlap with arXiv:2110.10626

The energy dependency (running) of the strength of electromagnetic interactions $\alpha$ plays an important role in precision tests of the Standard Model. The running of the former to the $Z$ pole is an input quantity for global electroweak fits, while the running of the mixing angle is susceptible to the effects of Beyond Standard Model physics, particularly at low energies. We present a computation of the hadronic vacuum polarization (HVP) contribution to the running of these electroweak couplings at the non-perturbative level in lattice QCD, in the space-like regime up to $Q^2$ momentum transfers of $7\,\mathrm{GeV}^2$. This quantity is also closely related to the HVP contribution to the muon $g-2$. We observe a tension of up to $3.5$ standard deviation between our lattice results for $\Delta\alpha^{(5)}_{\mathrm{had}}(-Q^2)$ and estimates based on the $R$-ratio for $Q^2$ in the $3$ to $7\,\mathrm{GeV}^2$ range. The tension is, however, strongly diminished when translating our result to the $Z$ pole, by employing the Euclidean split technique and perturbative QCD, which yields $\Delta\alpha^{(5)}_{\mathrm{had}}(M_Z^2)=0.027\,73(15)$. This value agrees with results based on the $R$-ratio within the quoted uncertainties, and can be used as an alternative to the latter in global electroweak fits., Comment: 6 pages, 2 figures, v2: accepted version with minor changes, talk presented at the 41st International Conference on High Energy physics - ICHEP2022

Prompted by an observation about the integral of exponential functions of the form $f(x)=\lambda\mathrm{e}^{\alpha x}$, we investigate the possibility to exactly integrate families of functions generated from a given function by scaling or by affine transformations of the argument using nonlinear generalizations of quadrature formulae. The main result of this paper is that such formulae can be explicitly constructed for a wide class of functions, and have the same accuracy as Newton-Cotes formulae based on the same nodes. We also show how Newton-Cotes formulae emerge as the linear case of our general formalism, and demonstrate the usefulness of the nonlinear formulae in the context of the Pad\'e-Laplace method of exponential analysis., Comment: 14 pages, 3 figures (24 pdf files)

The isovector axial form factor of the nucleon plays a key role in interpreting data from long-baseline neutrino oscillation experiments. We perform a lattice-QCD based calculation of this form factor, introducing a new method to directly extract its $z$-expansion from lattice correlators. Our final parametrization of the form factor, which extends up to spacelike virtualities of $0.7\,{\rm GeV}^2$ with fully quantified uncertainties, agrees with previous lattice calculations but is significantly less steep than neutrino-deuterium scattering data suggests., Comment: 6 pages + 6 pages supplementary material; (3+6) figures; 2 tables in suppl.mat

Euclidean time windows in the integral representation of the hadronic vacuum polarization contribution to the muon $g-2$ serve to test the consistency of lattice calculations and may help in tracing the origins of a potential tension between lattice and data-driven evaluations. In this paper, we present results for the intermediate time window observable computed using O($a$) improved Wilson fermions at six values of the lattice spacings below 0.1\,fm and pion masses down to the physical value. Using two different sets of improvement coefficients in the definitions of the local and conserved vector currents, we perform a detailed scaling study which results in a fully controlled extrapolation to the continuum limit without any additional treatment of the data, except for the inclusion of finite-volume corrections. To determine the latter, we use a combination of the method of Hansen and Patella and the Meyer-Lellouch-L\"uscher procedure employing the Gounaris-Sakurai parameterization for the pion form factor. We correct our results for isospin-breaking effects via the perturbative expansion of QCD+QED around the isosymmetric theory. Our result at the physical point is $a_\mu^{\mathrm{win}}=(237.30\pm0.79_{\rm stat}\pm1.22_{\rm syst})\times10^{-10}$, where the systematic error includes an estimate of the uncertainty due to the quenched charm quark in our calculation. Our result displays a tension of 3.9$\sigma$ with a recent evaluation of $a_\mu^{\mathrm{win}}$ based on the data-driven method., Comment: 43 pages, 9 figures, 10 tables; version accepted for publication: extended discussion of finite-volume corrections. Results and conclusions unchanged

We compute the hadronic running of the electromagnetic and weak couplings in lattice QCD with $N_{\mathrm{f}}=2+1$ flavors of $\mathcal{O}(a)$ improved Wilson fermions. Using two different discretizations of the vector current, we compute the quark-connected and -disconnected contributions to the hadronic vacuum polarization (HVP) functions $\bar{\Pi}^{\gamma\gamma}$ and $\bar{\Pi}^{\gamma Z}$ for Euclidean squared momenta $Q^2\leq 7\,\mathrm{GeV}^2$. Gauge field ensembles at four values of the lattice spacing and several values of the pion mass, including its physical value, are used to extrapolate the results to the physical point. The ability to perform an exact flavor decomposition allows us to present the most precise determination to date of the $\mathrm{SU}(3)$-flavor-suppressed HVP function $\bar{\Pi}^{08}$ that enters the running of $\sin^2\theta_{\mathrm{W}}$. Our results for $\bar{\Pi}^{\gamma\gamma}$, $\bar{\Pi}^{\gamma Z}$ and $\bar{\Pi}^{08}$ are presented in terms of rational functions for continuous values of $Q^2$ below $7 \,\mathrm{GeV}^2$. We observe a tension of up to $3.5$ standard deviation between our lattice results for $\Delta\alpha^{(5)}_{\mathrm{had}}(-Q^2)$ and estimates based on the $R$-ratio for space-like momenta in the range $3$--$7\,\mathrm{GeV}^2$. The tension is, however, strongly diminished when translating our result to the $Z$ pole, by employing the Euclidean split technique and perturbative QCD, which yields $\Delta\alpha^{(5)}_{\mathrm{had}}(M_Z^2)=0.027\,73(15)$ and agrees with results based on the $R$-ratio within the quoted uncertainties., Comment: 63 pages, 17 figures, 11 tables. v2: version accepted for publication in JHEP

We present an update of the scale setting for $N_f=2+1$ flavor QCD using gradient flow scales and pseudo-scalar decay constants. We analyze the latest ensembles with $2+1$ flavors of non-perturbatively improved Wilson fermions generated by CLS for improved precision. Special care is taken to correct for mistuning by measuring directly the mass derivatives of the various observables. We determine $t_0$ with input taken from a combination of leptonic decay rates of the Pion and the Kaon., Comment: 10 pages, 6 figures, The 38th International Symposium on Lattice Field Theory, LATTICE2021

We present the analysis of isovector axial vector nucleon form factors on a set of $N_f=2+1$ CLS ensembles with $\mathcal O(a)$-improved Wilson fermions and L\"uscher-Weisz gauge action. The set of ensembles covers a pion mass range of $130-353\,$MeV with lattice spacings between $0.05\,$fm and $0.09\,$fm. In particular, the set includes a $L/a=96$ ensemble at the physical pion mass. For the purpose of the form factor extraction, we employ both the summed operator insertion method (summation method) and explicit two-state fits in order to account for excited-state contributions to the nucleon correlation functions. To describe the $Q^{2}$-behavior of the form factors, we perform $z$-expansion fits. Finally, we present HBChPT-inspired chiral and continuum extrapolations of the axial charge and radius., Comment: 10 pages, 3 figures, 1 table. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present results for the isoscalar electromagnetic form factors of the nucleon computed on the Coordinated Lattice Simulations (CLS) ensembles with $N_\mathrm{f} = 2 + 1$ flavors of $\mathcal{O}(a)$-improved Wilson fermions and an $\mathcal{O}(a)$-improved conserved vector current. In order to estimate the excited-state contamination, we employ several source-sink separations and apply the summation method. For the computation of the quark-disconnected diagrams, a stochastic estimation based on the one-end trick is performed, in combination with a frequency-splitting technique and the hopping parameter expansion. By these means, we obtain a clear signal for the form factors including the quark-disconnected contributions, which have a statistically significant effect on our results., Comment: 9 pages, 5 figures, contribution to the 38th International Symposium on Lattice Field Theory (LATTICE21), 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology; enlarged font size in figures

Preliminary results are presented for nucleon isovector charges and twist-2 matrix elements which have been obtained employing an improved analysis strategy to deal with excited-state contamination. The set of CLS $N_f=2+1$ gauge ensembles in this study has been extended compared to our 2018 calculation, including an ensemble at physical quark masses. Besides the addition of new ensembles, the number of gauge configurations and measurements has been increased on several of the existing ensembles and the analysis has been extended to include additional source-sink separations. The ensembles cover a range of the light quark mass corresponding to $M_\pi\approx 0.130\,\mathrm{MeV} \ldots 350\,\mathrm{MeV}$, four values of the lattice spacing $a\approx0.05\,\mathrm{fm}\ldots0.09\,\mathrm{fm}$ and a large range of volumes. Results at the physical point are computed for each observable from a combined chiral, continuum and finite-volume extrapolation., Comment: 9 pages, 2 figures, 2 tables. Talk presented at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

As present and future experiments, on both the energy and precision frontiers, look to identify new physics beyond the Standard Model, we require more precise determinations of fundamental quantities, like the QED and electroweak couplings at various momenta. These can be obtained either entirely from experimental measurements, or from one such measurement at a particular virtuality combined with the couplings' virtuality dependence computed within the SM. Thus, a precise, entirely theoretical determination of the running couplings is highly desirable, even more since the preliminary results of the E989 experiment in Fermilab were published. We give results for the hadronic contribution to the QED running coupling $\alpha(Q^2)$ and weak mixing angle $\sin^2\theta_W(Q^2)$ in the space-like energy region $(0, 7]~\text{GeV}^2$ with a total relative uncertainty of $2\%$ at energies $Q^2 \ll 1~\text{GeV}^2$, and $1\%$ at $Q^2 > 1~\text{GeV}^2$., Comment: 9 pages, 3 figures, talk presented at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

We present a lattice calculation of the leading hadronic contribution to the anomalous magnetic moment of the muon. This work is based on a subset of the CLS ensembles with $N_f = 2+1$ dynamical quarks and a quenched charm quark. Noise reduction techniques are used to improve significantly the statistical precision of the dominant light quark contribution. The main source of systematic error comes from finite size effects which are estimated using the formalism described in Ref. [7] and based on our knowledge of the timelike pion form factor. The strange and charm quark contributions are under control and an estimate of the quark-disconnected contribution is included. Isospin breaking effects will be studied in a future publication but are included in the systematic error using an estimate based on published lattice results. Our final result, $a_\mu^{\rm hvp} = (720.0\pm 12.4 \pm 6.8)\times 10^{-10}$, has a precision of about 2%., Comment: 7 pages. Contribution to the 37th International Symposium on Lattice Field Theory (LATTICE2019), 16-22 June 2019, Wuhan, China

We report on our calculation of the strange contribution to the vector and axial vector form factors. The strange charge radii, magnetic moment, and axial charge are extracted by model independent $z$-expansion fits to the $Q^2$-dependence of the respective form factors. Furthermore, the isoscalar contribution to the axial and tensor charge is investigated by combining the calculation of connected and disconnected diagrams. The required renormalization is performed with the Rome-Southampton method. We make use of the CLS $N_f=2+1$ $\mathcal{O}(a)$-improved Wilson fermion ensembles. Results are reported for pion masses in the range $m_\pi=200-360\,\text{MeV}$ and lattice spacings $a=0.05-0.086\,\text{fm}$., Comment: 7 pages, 11 figures, 1 table, Proceedings for the 37th International Symposium on Lattice Field Theory, 16-22 June 2019, Wuhan, China

We present results from a lattice QCD study of nucleon matrix elements at vanishing momentum transfer for local and twist-2 isovector operator insertions. Computations are performed on gauge ensembles with non-perturbatively improved $N_f=2+1$ Wilson fermions, covering four values of the lattice spacing and pion masses down to $M_\pi\approx200$MeV. Several source-sink separations (typically ~1.0fm to ~1.5fm) allow us to assess excited-state contamination. Results on individual ensembles are obtained from simultaneous two-state fits across all observables and all available source-sink separations with the energy gap as a common fit parameter. Renormalization has been performed non-perturbatively using the Rome-Southampton method for all but the finest lattice spacing for which an extrapolation has been used. Physical results are quoted in the $\overline{MS}$ scheme at a scale of $\mu=2$GeV and are obtained from a combined chiral, continuum and finite-size extrapolation. For the nucleon isovector axial, scalar and tensor charges we find physical values of $g_A^{u-d}=1.242(25)_\text{stat}(\genfrac{}{}{0pt}{2}{+00}{-31})_\text{sys}$, $g_S^{u-d}=1.13(11)_\text{stat}(\genfrac{}{}{0pt}{2}{+07}{-06})_\text{sys}$ and $g_T^{u-d}=0.965(38)_\text{stat}(\genfrac{}{}{0pt}{2}{+13}{-41})_\text{sys}$, respectively, where individual systematic errors in each direction from the chiral, continuum and finite-size extrapolation have been added in quadrature. Our final results for the isovector average quark momentum fraction and the isovector helicity and transversity moments are given by $\langle x\rangle_{u-d}=0.180(25)_\text{stat}(\genfrac{}{}{0pt}{2}{+14}{-06})_\text{sys}$, $\langle x\rangle_{\Delta u-\Delta d}=0.221(25)_\text{stat}(\genfrac{}{}{0pt}{2}{+10}{-00})_\text{sys}$ and $\langle x\rangle_{\delta u-\delta d}=0.212(32)_\text{stat}(\genfrac{}{}{0pt}{2}{+20}{-10})_\text{sys}$, respectively., Comment: 35 pages, 13 figures; matching version accepted for publication in PRD

The comparison of the theoretical and experimental determinations of the anomalous magnetic moment of the muon $(g-2)_\mu$ constitutes one of the strongest tests of the Standard Model at low energies. In this article, we compute the leading hadronic contribution to $(g-2)_\mu$ using lattice QCD simulations employing Wilson quarks. Gauge field ensembles at four different lattice spacings and several values of the pion mass down to its physical value are used. We apply the O($a$) improvement programme with two discretizations of the vector current to better constrain the approach to the continuum limit. The electromagnetic current correlators are computed in the time-momentum representation. In addition, we perform auxiliary calculations of the pion form factor at timelike momenta in order to better constrain the tail of the isovector correlator and to correct its dominant finite-size effect. For the numerically dominant light-quark contribution, we have rescaled the lepton mass by the pion decay constant computed on each lattice ensemble. We perform a combined chiral and continuum extrapolation to the physical point, and our final result is $ a_\mu^{\rm hvp}=(720.0\pm12.4_{\rm stat}\,\pm9.9_{\rm syst})\cdot10^{-10}$. It contains the contributions of quark-disconnected diagrams, and the systematic error has been enlarged to account for the missing isospin-breaking effects., Comment: 30 pages, 10 figures

We improve the non-relativistic QCD (NRQCD) action by comparing the dispersion relation to that of the continuum through $\mathcal{O}(p^6)$ in perturbation theory. The one-loop matching coefficients of the $\mathcal{O}(p^4)$ kinetic operators are determined, as well as the scale at which to evaluate $\alpha_s$ in the $V$-scheme for each quantity. We utilise automated lattice perturbation theory using twisted boundary conditions as an infrared regulator. The one-loop radiative corrections to the mass renormalisation, zero-point energy and overall energy-shift of an NRQCD $b$-quark are also found. We also explore how a Fat$3$-smeared NRQCD action and changes of the stability parameter $n$ affect the coefficients. Finally, we use gluon field ensembles at multiple lattice spacing values, all of which include $u$, $d$, $s$ and $c$ quark vacuum polarisation, to test how the improvements affect the non-perturbatively determined $\Upsilon(1S)$ and $\eta_b(1S)$ kinetic masses, and the tuning of the $b$ quark mass., Comment: 19 pages, 14 figures, 9 tables

We report on our ongoing project to determine the leading-order hadronic vacuum polarisation contribution to the muon $g-2$, using ensembles with $N_f=2+1$ flavours of O($a$) improved Wilson quarks generated by the CLS effort, with pion masses down to the physical value. We employ O($a$) improved versions of the local and conserved vector currents to compute the contributions of the light, strange and charm quarks to $(g-2)_\mu$, using the time-momentum representation. We perform a detailed investigation of the systematic effects arising from constraining the long-distance regime of the vector correlator. To this end we make use of auxiliary calculations in the iso-vector channel using distillation and the L\"uscher formalism. Our results are corrected for finite-volume effects by computing the timelike pion form factor in finite and infinite volume. For certain parameter choices, the corrections computed in this way can also be confronted with results determined on different volumes. Currently, the overall precision of our results is limited by the uncertainties in the lattice scale., Comment: To appear in: Proceedings of the 36th International Symposium on Lattice Field Theory (Lattice 2018), Michigan State University, East Lansing, USA; 7 pages

We present preliminary results for strangeness form factors of the nucleon computed on the CLS ensembles with $N_f=2+1$ flavors of O($a$)-improved Wilson fermions. Our calculations are performed at two values of the lattice spacing ($a \in \{0.064, 0.086\} \,\mathrm{fm}$) and three values of the pion mass ($m_\pi \in \{ 200, 280, 340\} \,\mathrm{MeV}$). The determination of strangeness form factors proceeds by computing quark-disconnected diagrams, for which we employ hierarchical probing in four dimensions in order to deal with this most challenging part of the calculation. Furthermore, we investigate several source-sink separations to check for excited-state contamination., Comment: 7 pages, 8 figures, 1 table, Proceedings for the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

We present a nucleon structure analysis including local isovector charges as well as twist-2 operator insertions. Computations are performed on CLS ensembles with $N_f=2+1$ Wilson fermions, covering four values of the lattice spacing and pion masses down to $M_\pi \approx 200\,\mathrm{MeV}$. Several source-sink separations (typically $\sim 1.0\,\mathrm{fm}$ to $\sim 1.5\,\mathrm{fm}$) allow us to assess excited-state contaminations. Results on each ensemble are obtained from simultaneous two-state fits including all observables and all available source-sink separations with the mass gap as a free fit parameter. Finally, the chiral, continuum and finite size extrapolation is performed to extract physical results., Comment: 7 pages, 4 figures, 1 table. Talk presented at the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

The hadronic contribution to the eight forward amplitudes of light-by-light scattering ($\gamma^*\gamma^*\to \gamma^*\gamma^*$) is computed in lattice QCD. Via dispersive sum rules, the amplitudes are compared to a model of the $\gamma^*\gamma^*\to {\rm hadrons}$ cross sections in which the fusion process is described by hadronic resonances. Our results thus provide an important test for the model estimates of hadronic light-by-light scattering in the anomalous magnetic moment of the muon, $a_\mu^{\rm HLbL}$. Using simple parametrizations of the resonance $M\to \gamma^*\gamma^*$ transition form factors, we determine the corresponding monopole and dipole masses by performing a global fit to all eight amplitudes. Together with a previous dedicated calculation of the $\pi^0\to \gamma^*\gamma^*$ transition form factor, our calculation provides valuable information for phenomenological estimates of $a_\mu^{\rm HLbL}$. The presented calculations are performed in two-flavor QCD with pion masses extending down to 190\,MeV at two different lattice spacings. In addition to the fully connected Wick contractions, on two lattice ensembles we also compute the (2+2) disconnected class of diagrams, and find that their overall size is compatible with a parameter-free, large-$N$ inspired prediction, where $N$ is the number of colors. Motivated by this observation, we estimate in the same way the disconnected contribution to $a_\mu^{\rm HLbL}$., Comment: 41 pages, 10 figures with up to 8 panels

We present a lattice calculation of the nucleon iso-vector axial and induced pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical flavours of non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions and an $\mathcal{O}(a)$-improved axial current together with the pseudoscalar density. Excited-state effects in the extraction of the form factors are treated using a variety of methods, with a detailed discussion of their respective merits. The chiral and continuum extrapolation of the results is performed both using formulae inspired by Heavy Baryon Chiral Perturbation Theory (HBChPT) and a global approach to the form factors based on a chiral effective theory (EFT) including axial vector mesons. Our results indicate that careful treatment of excited-state effects is important in order to obtain reliable results for the axial form factors of the nucleon, and that the main remaining error stems from the systematic uncertainties of the chiral extrapolation. As final results, we quote $g_{\rm A} = 1.278 \pm 0.068\genfrac{}{}{0pt}{1}{+0.000}{-0.087}$, $\langle r_{\rm A}^2\rangle = 0.360 \pm 0.036\genfrac{}{}{0pt}{1}{+0.080}{-0.088}~\mathrm{fm}^2$, and $g_{\rm P} = 7.7 \pm 1.8 \genfrac{}{}{0pt}{1}{+0.8}{-2.0}$ for the axial charge, axial charge radius and induced pseudoscalar charge, respectively, where the first error is statistical and the second is systematic., Comment: 27 pages, 13 figures (23 PDF files); uses revtex4-1.cls; v4: corrections to affiliation addresses, matches published version

We present updated results on the nucleon electromagnetic form factors and axial coupling calculated using CLS ensembles with $N_\mathrm{f}=2+1$ dynamical flavours of Wilson fermions. The measurements are performed on large, fine lattices with a pseudoscalar mass reaching down to 200 MeV. The truncated-solver method is employed to reduce the variance of the measurements. Estimation of the matrix elements is challenging due to large contamination from excited states and further investigation is necessary to bring these effects under control., Comment: corrected reference; 7 pages, 4 figures, talk presented at 34th annual International Symposium on Lattice Field Theory, 24-30 July 2016, University of Southampton, UK

We study the performance of all-mode-averaging (AMA) when used in conjunction with a locally deflated SAP-preconditioned solver, determining how to optimize the local block sizes and number of deflation fields in order to minimize the computational cost for a given level of overall statistical accuracy. We find that AMA enables a reduction of the statistical error on nucleon charges by a factor of around two at the same cost when compared to the standard method. As a demonstration, we compute the axial, scalar and tensor charges of the nucleon in $N_f=2$ lattice QCD with non-perturbatively O(a)-improved Wilson quarks, using O(10,000) measurements to pursue the signal out to source-sink separations of $t_s\sim 1.5$ fm. Our results suggest that the axial charge is suffering from a significant amount (5-10%) of excited-state contamination at source-sink separations of up to $t_s\sim 1.2$ fm, whereas the excited-state contamination in the scalar and tensor charges seems to be small., Comment: 24 pages, 11 figures (25 pdf files); v3:version to appear in Nucl.Phys.B

We present preliminary results on the electromagnetic form factors and axial charge of the nucleon from ensembles generated by the CLS effort with $N_\mathrm{f}=2+1$ flavours of non-perturbatively $\mathrm{O}(a)$-improved Wilson fermions and open temporal boundary conditions. Systematic effects due to excited-state contamination are accounted for using both two-state fits and the method of summed operator insertions. This exploratory analysis demonstrates the viability of obtaining precision baryon observables with $N_\mathrm{f}=2+1$ flavours of Wilson fermions on fine lattices, aiming towards controlled chiral and continuum limits in the future., Comment: 7 pages, 4 figures; contribution to the proceedings of the 33rd International Symposium on Lattice Field Theory, July 14-18, 2015, Kobe, Japan

We report calculations of hadronic light-by-light scattering amplitudes via lattice QCD evaluation of Euclidean four-point functions of vector currents. These initial results include only the fully quark-connected contribution. Particular attention is given to the case of forward scattering, which can be related via dispersion relations to the $\gamma^* \gamma^* \to$ hadrons cross section, and thus allows lattice data to be compared with phenomenology. We also present a strategy for computing the hadronic light-by-light contribution to the muon anomalous magnetic moment., Comment: 7 pages, 4 figures. Talk presented at the 33rd International Symposium on Lattice Field Theory (Lattice 2015), 14-18 July 2015, Kobe International Conference Center, Kobe, Japan

We present a calculation of the hindered M$1$ $\Upsilon(2S) \to \eta_b(1S) \gamma$ decay rate using lattice non-relativistic QCD. The calculation includes spin-dependent relativistic corrections to the NRQCD action through $\mathcal{O}(v^6)$ in the quark's relative velocity, relativistic corrections to the leading order current which mediates the transition through the quark's magnetic moment, radiative corrections to the leading spin-magnetic coupling and for the first time a full error budget. We also use gluon field ensembles at multiple lattice spacing values, all of which include $u$, $d$, $s$ and $c$ quark vacuum polarisation. Our result for the branching fraction is $\mathcal{B}(\Upsilon(2S)\to\eta_b(1S)\gamma) = 5.4(1.8)\times 10^{-4} $, which agrees with the current experimental value., Comment: 26 pages, 11 figures. Updated to match version published in Phys.Rev.D

We extend our study of the pion scalar radius in two-flavour lattice QCD to include two additional lattice spacings as well as lighter pion masses, enabling us to perform a combined chiral and continuum extrapolation. We find discretisation artefacts to be small for the radius, and confirm the importance of the disconnected diagrams in reproducing the correct chiral behaviour. Our final result for the scalar radius of the pion at the physical point is $\left\langle r^2\right\rangle^\pi_{\rm S}=0.600\pm0.052$ fm$^2$, corresponding to a value of $\overline{\ell}_4=4.54\pm0.30$ for the low-energy constant $\overline{\ell}_4$ of NLO chiral perturbation theory., Comment: 4 pages, 4 figures, uses svjour.cls

We perform a lattice QCD calculation of the hadronic light-by-light scattering amplitude in a broad kinematical range. At forward kinematics, the results are compared to a phenomenological analysis based on dispersive sum rules for light-by-light scattering. The size of the pion pole contribution is investigated for momenta of typical hadronic size. The presented numerical methods can be used to compute the hadronic light-by-light contribution to the anomalous magnetic moment of the muon. Our calculations are carried out in two-flavor QCD with the pion mass in the range of 270 to 450MeV, and contain so far only the diagrams with fully connected quark lines., Comment: 5 pages, 5 figures