Domain decomposition and multilevel integration for fermions

The numerical computation of many hadronic correlation functions is exceedingly difficult due to the exponentially decreasing signal-to-noise ratio with the distance between source and sink. Multilevel integration methods, using independent updates of separate regions in space-time, are known to be able to solve such problems but have so far been available only for pure gauge theory. We present first steps into the direction of making such integration schemes amenable to theories with fermions, by factorizing a given observable via an approximated domain decomposition of the quark propagator. This allows for multilevel integration of the (large) factorized contribution to the observable, while its (small) correction can be computed in the standard way.Comment: 14 pages, 6 figures, v2: published version, talk presented at the 34th annual International Symposium on Lattice Field Theory, 24-30 July 2016, University of Southampton, U

Local multiboson factorization of the quark determinant

We discuss the recently proposed multiboson domain-decomposed factorization of the gauge-field dependence of the fermion determinant in lattice QCD. In particular, we focus on the case of a lattice divided in an arbitrary number of thick time slices. As a consequence, multiple space-time regions can be updated independently. This allows to address the exponential degradation of the signal-to-noise ration of correlation functions with multilevel Monte Carlo sampling. We show numerical evidence of the effectiveness of a two-level integration for pseudoscalar propagators with momentum and for vector propagators, in a two active regions setup. These results are relevant to lattice computation of the hadronic contributions to the anomalous magnetic moment of the muon and to heavy meson decay form factors.Comment: 8 pages, 4 figures, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spai

Stochastic locality and master-field simulations of very large lattices

In lattice QCD and other field theories with a mass gap, the field variables in distant regions of a physically large lattice are only weakly correlated. Accurate stochastic estimates of the expectation values of local observables may therefore be obtained from a single representative field. Such master-field simulations potentially allow very large lattices to be simulated, but require various conceptual and technical issues to be addressed. In this talk, an introduction to the subject is provided and some encouraging results of master-field simulations of the SU(3) gauge theory are reported.Comment: Talk given at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain; LaTeX source with 6 figure

Testing the strength of the UA(1)\text{U}_A(1) anomaly at the chiral phase transition in two-flavour QCD

We study the thermal transition of QCD with two degenerate light flavours by lattice simulations using $\mathcal{O}(a)$-improved Wilson quarks. Particular emphasis lies on the pattern of chiral symmetry restoration, which we probe via the static screening correlators. On $32^3$ volumes we observe that the screening masses in transverse iso-vector vector and axial-vector channels become degenerate at the transition temperature. The splitting between the screening masses in iso-vector scalar and pseudoscalar channels is strongly reduced compared to the splitting at zero temperature and is actually consistent with zero within uncertainties. In this proceedings article we extend our studies to matrix elements and iso-singlet correlation functions. Furthermore, we present results on larger volumes, including first results at the physical pion mass.Comment: 10 pages, 9 figures, invited contribution to the 9th International Workshop on Chiral Dynamics, Sept. 17-21, 2018, Duke University, Durham, NC, US

The leading hadronic contribution to (g−2)μ(g-2)_\mu from lattice QCD with Nf=2+1N_{\rm f}=2+1 flavours of O(aa) improved Wilson quarks

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 figure

Deep inelastic scattering on the quark-gluon plasma

We provide an interpretation of the structure functions of a thermal medium such as the quark-gluon plasma in terms of the scattering of an incoming electron on the medium via the exchange of a spacelike photon. We then focus on the deep-inelastic scattering (DIS) regime, and formulate the corresponding moment sum rules obeyed by the structure functions. Accordingly, these moments are given by the thermal expectation value of twist-two operators, which is computable from first principles in lattice QCD for the first few moments. We also show how lattice QCD calculations can be used to probe how large the photon virtuality needs to be in order for the Bjorken scaling of structure functions to set in. Finally, we provide the parton-model interpretation of the structure functions in the Bjorken limit and test its consistency. As in DIS on the proton, the kinematic variable $x$ is proportional to the longitudinal momentum carried by the partons, however $x$ ranges from zero to infinity. Choosing the parton momentum parametrization to be $x T u$ where $u$ is the fluid four-velocity and $T$ its temperature in the rest frame, the parton distribution function for a plasma of non-interacting quarks is proportional to $x \log(1+e^{-x/2})$.We provide an interpretation of the structure functions of a thermal medium such as the quark-gluon plasma in terms of the scattering of an incoming electron on the medium via the exchange of a spacelike photon. We then focus on the deep-inelastic scattering (DIS) regime, and formulate the corresponding moment sum rules obeyed by the structure functions. Accordingly, these moments are given by the thermal expectation value of twist-two operators, which is computable from first principles in lattice QCD for the first few moments. We also show how lattice QCD calculations can be used to probe how large the photon virtuality needs to be in order for the Bjorken scaling of structure functions to set in. Finally, we provide the parton-model interpretation of the structure functions in the Bjorken limit and test its consistency. As in DIS on the proton, the kinematic variable x is proportional to the longitudinal momentum carried by the partons, however x ranges from zero to infinity. Choosing the parton momentum parametrization to be xT u where u is the fluid four-velocity and T its temperature in the rest frame, the parton distribution function for a plasma of non-interacting quarks is proportional to x log(1 + e$^{−x/2}$)

The leading hadronic contribution to the running of the Weinberg angle using covariant coordinate-space methods

We present a preliminary study of the leading hadronic contribution to the running of the Weinberg angle $\theta_{\mathrm{W}}$. The running is extracted from the correlation function of the electromagnetic current with the vector part of the weak neutral current using both the standard time-momentum representation method and the Lorentz-covariant coordinate-space method recently introduced by Meyer. Both connected and disconnected contributions have been computed on $N_{\mathrm{f}}=2+1$ non-perturbatively $O(a)$-improved Wilson fermions configurations from the CLS initiative. Similar covariant coordinate-space methods can be used to compute the leading hadronic contribution to the anomalous magnetic moment of the muon $(g-2)_\mu$ and to the running of the QED coupling $\alpha$.Comment: 7 pages, 2 figures, talk presented at The 36th Annual International Symposium on Lattice Field Theory, July 22-28, 2018, East Lansing, MI, US

Photon emissivity of the quark-gluon plasma: a lattice QCD analysis of the transverse channel

We present results for the thermal photon emissivity of the quark-gluon plasma derived from spatially transverse vector correlators computed in lattice QCD at a temperature of 250 MeV. The analysis of the spectral functions, performed at fixed spatial momentum, is based on continuum-extrapolated correlators obtained with two flavours of dynamical Wilson fermions. We compare the next-to-leading order perturbative QCD correlators, as well as the ${\cal N}=4$ supersymmetric Yang-Mills correlators at infinite coupling, to the correlators from lattice QCD and find them to lie within $\sim10\%$ of each other. We then refine the comparison, performing it at the level of filtered spectral functions obtained model-independently via the Backus-Gilbert method. Motivated by these studies, for frequencies $\omega\lesssim2.5\,$GeV we use fit ans\"atze to the spectral functions that perform well when applied to mock data generated from the NLO QCD or from the strongly-coupled SYM spectral functions, while the high-frequency part, $\omega\gtrsim 2.5\,$GeV, is matched to NLO QCD. We compare our results for the photon emissivity to our previous analysis of a different vector channel at the same temperature. We obtain the most stringent constraint at photon momenta around $k\simeq0.8\,$GeV, for which we find a differential photon emission rate per unit volume of $d\Gamma_\gamma/d^3k = (\alpha_{\rm em}/(\exp(k/T)-1))\times (2.2 \pm 0.8 ) \times 10^{-3}\,{\rm GeV}$.Comment: 26 pages, 13 figures, 1 tabl