Vacuum correlators at short distances from lattice QCD

TL;DR

This paper introduces a novel lattice QCD method using high-temperature correlators to accurately compute vacuum polarization at short distances, reducing discretization errors for electroweak precision tests.

Contribution

It proposes a new approach leveraging static screening correlators at high temperature to control discretization errors in short-distance vacuum polarization calculations.

Findings

Identified a specific lattice artifact in vacuum polarization calculations.

Demonstrated the method's effectiveness with non-perturbative lattice QCD data.

Achieved simulations with lattice spacings down to 0.033 fm.

Abstract

Non-perturbatively computing the hadronic vacuum polarization at large photon virtualities and making contact with perturbation theory enables a precision determination of the electromagnetic coupling at the $Z$ pole, which enters global electroweak fits. In order to achieve this goal ab initio using lattice QCD, one faces the challenge that, at the short distances which dominate the observable, discretization errors are hard to control. Here we address challenges of this type with the help of static screening correlators in the high-temperature phase of QCD, yet without incurring any bias. The idea is motivated by the observations that (a) the cost of high-temperature simulations is typically much lower than their vacuum counterpart, and (b) at distances $x_3$ far below the inverse temperature $1/T$, the operator-product expansion guarantees the thermal correlator of two local currents to deviate from the vacuum correlator by a relative amount that is power-suppressed in $(x_3\:T)$. The method is first investigated in lattice perturbation theory, where we point out the appearance of an O$(a^2 \log(1/a))$ lattice artifact in the vacuum polarization with a prefactor that we calculate. It is then applied to non-perturbative lattice QCD data with two dynamical flavors of quarks. Our lattice spacings range down to 0.049 fm for the vacuum simulations and down to 0.033 fm for the simulations performed at a temperature of 250 MeV.