Source author record

Harvey B. Meyer

Harvey B. Meyer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

67works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

67 published item(s)

preprint2023arXiv

Window observable for the hadronic vacuum polarization contribution to the muon $g-2$ from lattice QCD

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üscher 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_μ^{\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$σ$ with a recent evaluation of $a_μ^{\mathrm{win}}$ based on the data-driven method.

preprint2022arXiv

Electromagnetic form factors of the proton and neutron from $N_f = 2 + 1$ lattice QCD

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.

preprint2022arXiv

Intermediate window observable for the hadronic vacuum polarization contribution to the muon $g-2$ from O$(a)$ improved Wilson quarks

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\,σ$. 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\,σ$ with a recent evaluation of the intermediate window based on the data-driven method.

preprint2022arXiv

The charm-quark contribution to light-by-light scattering in the muon $(g-2)$ from lattice QCD

We compute the hadronic light-by-light scattering contribution to the muon $g-2$ from the charm quark using lattice QCD. The calculation is performed on ensembles generated with dynamical $(u,d,s)$ quarks at the SU(3)$_{\rm f}$ symmetric point with degenerate pion and kaon masses of around 415 MeV. It includes the connected charm contribution, as well as the leading disconnected Wick contraction, involving the correlation between a charm and a light-quark loop. Cutoff effects turn out to be sizeable, which leads us to use lighter-than-physical charm masses, to employ a broad range of lattice spacings reaching down to 0.039 fm and to perform a combined charm-mass and continuum extrapolation. We use the $η_c$ meson to define the physical charm-mass point and obtain a final value of $a_μ^{\rm HLbL,c} = (2.8\pm 0.5) \times 10^{-11}$, whose uncertainty is dominated by the systematics of the extrapolation. Our result is consistent with the estimate based on a simple charm-quark loop, whilst being free of any perturbative scheme dependence on the charm mass. The mixed charm-light disconnected contraction contributes a small negative amount to the final value.

preprint2022arXiv

The hadronic contribution to the running of the electromagnetic coupling and electroweak mixing angle

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 $α(Q^2)$ and weak mixing angle $\sin^2θ_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$.

preprint2022arXiv

The hadronic running of the electromagnetic coupling and the electroweak mixing angle from lattice QCD

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Π^{γγ}$ and $\barΠ^{γ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Π^{08}$ that enters the running of $\sin^2θ_{\mathrm{W}}$. Our results for $\barΠ^{γγ}$, $\barΠ^{γZ}$ and $\barΠ^{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 $Δα^{(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 $Δα^{(5)}_{\mathrm{had}}(M_Z^2)=0.027\,73(15)$ and agrees with results based on the $R$-ratio within the quoted uncertainties.

preprint2022arXiv

The isovector axial form factor of the nucleon from lattice QCD

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.

preprint2021arXiv

Vacuum correlators at short distances from lattice QCD

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.

preprint2019arXiv

The hadronic contribution to the running of the electromagnetic coupling and the electroweak mixing angle

The electromagnetic coupling $α$ and the electroweak mixing angle $θ_{\mathrm{W}}$ are parameters of the Standard Model (SM) that enter precision SM tests and play a fundamental rôle in beyond SM physics searches. Their values are energy dependent, and non-perturbative hadronic contributions are the main source of uncertainty to the theoretical knowledge of the running with energy. We present a lattice study of the leading hadronic contribution to the running of $α$ and $\sin^2θ_{\mathrm{W}}$. The former is related to the hadronic vacuum polarization (HVP) function of electromagnetic currents, and the latter to the HVP mixing of the electromagnetic current with the vector part of the weak neutral currents. We use the time-momentum representation (TMR) method to compute the HVP on the lattice, estimating both connected and disconnected contributions on $N_{\mathrm{f}}=2+1$ non-perturbatively $O(a)$-improved Wilson fermions ensembles from the Coordinated Lattice Simulations (CLS) initiative. The use of different lattice spacings and quark masses allows us to reliably extrapolate the results to the physical point.

preprint2016arXiv

A variational method for spectral functions

The Generalized Eigenvalue Problem (GEVP) has been used extensively in the past in order to reliably extract energy levels from time-dependent Euclidean correlators calculated in Lattice QCD. We propose a formulation of the GEVP in frequency space. Our approach consists of applying the model-independent Backus-Gilbert method to a set of Euclidean two-point functions with common quantum numbers. A GEVP analysis in frequency space is then applied to a matrix of estimators that allows us, among other things, to obtain particular linear combinations of the initial set of operators that optimally overlap to different local regions in frequency. We apply this method to lattice data from NRQCD. This approach can be interesting both for vacuum physics as well as for finite-temperature problems.

preprint2016arXiv

Lattice calculation of the pion transition form factor $π^0 \to γ^* γ^*$

We calculate the $π^0\to γ^*γ^*$ transition form factor ${\cal F}_{π^0γ^*γ^*}(q_1^2,q_2^2)$ in lattice QCD with two flavors of quarks. Our main motivation is to provide the input to calculate the $π^0$-pole contribution to hadronic light-by-light scattering in the muon $(g-2)$, $a_μ^{\rm HLbL;π^0}$. We therefore focus on the region where both photons are spacelike up to virtualities of about $1.5~$GeV$^2$, which has so far not been experimentally accessible. Results are obtained in the continuum at the physical pion mass by a combined extrapolation. We reproduce the prediction of the chiral anomaly for real photons with an accuracy of about $8-9\%$. We also compare to various recently proposed models and find reasonable agreement for the parameters of some of these models with their phenomenological values. Finally, we use the parametrization of our lattice data by these models to calculate $a_μ^{\rm HLbL;π^0}$.

preprint2016arXiv

Lattice calculation of the pion transition form factor $π^0 \to γ^* γ^*$

We calculate the pion transition form factor ${\cal F}_{π^0γ^*γ^*}(q_1^2,q_2^2)$, which describe the interaction of an on-shell pion with two off-shell photons, using lattice QCD simulations with two degenerate flavors of dynamical quarks. 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$, $a_μ^{\mathrm{HLbL}; π^0}$. We focus our study on the spacelike region with photon virtualities up to $1.5~\mathrm{GeV}^2$, not yet measured experimentally. Several lattice spacings and pion masses are used to extrapolate the results to the physical point and a comparison with different phenomenological models is performed. Finally, we use our extrapolated form factor to provide a lattice determinaiton of $a_μ^{\mathrm{HLbL}; π^0}$.

preprint2016arXiv

Position-space approach to hadronic light-by-light scattering in the muon $g-2$ on the lattice

The anomalous magnetic moment of the muon currently exhibits a discrepancy of about three standard deviations between the experimental value and recent Standard Model predictions. The theoretical uncertainty is dominated by the hadronic vacuum polarization and the hadronic light-by-light (HLbL) scattering contributions, where the latter has so far only been fully evaluated using different models. To pave the way for a lattice calculation of HLbL, we present an expression for the HLbL contribution to $g-2$ that involves a multidimensional integral over a position-space QED kernel function in the continuum and a lattice QCD four-point correlator. We describe our semi-analytic calculation of the kernel and test the approach by evaluating the $π^0$-pole contribution in the continuum.

preprint2016arXiv

Static and non-static vector screening masses

Thermal screening masses of the conserved vector current are calculated both in a weak-coupling approach and in lattice QCD. The inverse of a screening mass can be understood as the length scale over which an external electric field is screened in a QCD medium. The comparison of screening masses both in the zero and non-zero Matsubara frequency sectors shows good agreement of the perturbative and the lattice results. Moreover, at $T\approx 508\mathrm{MeV}$ the lightest screening mass lies above the free result ($2πT$), in agreement with the $\mathcal{O}(g^2)$ weak-coupling prediction.

preprint2015arXiv

Charge transport and vector meson dissociation across the thermal phase transition in lattice QCD with two light quark flavors

We compute and analyze correlation functions in the isovector vector channel at vanishing spatial momentum across the deconfinement phase transition in lattice QCD. The simulations are carried out at temperatures $T/T_c=0.156, 0.8, 1.0, 1.25$ and $1.67$ with $T_c\simeq203$MeV for two flavors of Wilson-Clover fermions with a zero-temperature pion mass of $\simeq270$MeV. Exploiting exact sum rules and applying a phenomenologically motivated ansatz allows us to determine the spectral function $ρ(ω,T)$ via a fit to the lattice correlation function data. From these results we estimate the electrical conductivity across the deconfinement phase transition via a Kubo formula and find evidence for the dissociation of the $ρ$ meson by resolving its spectral weight at the available temperatures. We also apply the Backus-Gilbert method as a model-independent approach to this problem. At any given frequency, it yields a local weighted average of the true spectral function. We use this method to compare kinetic theory predictions and previously published phenomenological spectral functions to our lattice study.

preprint2015arXiv

Direct calculation of hadronic light-by-light scattering

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 $γ^* γ^* \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.

preprint2015arXiv

Lattice QCD calculation of hadronic light-by-light scattering

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.

preprint2015arXiv

Non-perturbative improvement of the vector current in Wilson lattice QCD

Many observables of interest in lattice QCD are extracted from correlation functions involving the vector current. If Wilson fermions are used, it is therefore of practical importance that, besides the action, the current be O($a$) improved in order to remove the leading discretization errors from the observables. Here we introduce and apply a new method to determine the improvement coefficient for the two most widely used discretizations of the current.

preprint2015arXiv

Nucleon electromagnetic form factors and axial charge from CLS $N_\mathrm{f}=2+1$ ensembles

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.

preprint2015arXiv

Study of the hadronic contributions to the running of the QED coupling and the weak mixing angle

The electromagnetic coupling receives significant contributions to its running from non-perturbative QCD effects. We present an update of a lattice QCD study of the Adler function and of its application to the determination of the running of the QED coupling. We perform a high-statistics computation with two flavours of O$(a)$ improved Wilson fermions in a large range of momentum transfer $Q^2$. The running of the electromagnetic coupling, including contributions from $u$, $d$, $s$ and $c$ valence quarks, is compared to phenomenological determinations at intermediate $Q^2$ values. An extension of this study to the determination of the hadronic contributions to the running of the weak mixing angle is also described.

preprint2015arXiv

The pion quasiparticle in the low-temperature phase of QCD

We investigate the properties of the pion quasiparticle in the low-temperature phase of two-flavor QCD on the lattice with support from chiral effective theory. We find that the pion quasiparticle mass is significantly reduced compared to its value in the vacuum, in contrast to the static screening mass, which increases with temperature. By a simple argument, the two masses are expected to determine the quasiparticle dispersion relation near the chiral limit. Analyzing two-point functions of the axial charge density at non-vanishing spatial momentum, we find that the predicted dispersion relation and the residue of the pion pole are simultaneously consistent with the lattice data at low momentum. The test, based on fits to the correlation functions, is confirmed by a second analysis using the Backus-Gilbert method.

preprint2015arXiv

The pion quasiparticle in the low-temperature phase of QCD

We investigate the properties of the pion quasiparticle in the low-temperature phase of two-flavor QCD on the lattice with support from chiral effective theory. We find that the pion quasiparticle mass is significantly reduced compared to its value in the vacuum, by contrast with the static screening mass, which increases with temperature. By a simple argument, near the chiral limit the two masses are expected to determine the quasiparticle dispersion relation. Analyzing two-point functions of the axial charge density at non-vanishing spatial momentum, we find that the predicted dispersion relation and the residue of the pion pole are simultaneously consistent with the lattice data at low momentum. The test, based on fits to the correlation functions, is confirmed by a second analysis using the Backus-Gilbert method.

preprint2014arXiv

Chiral dynamics in the low-temperature phase of QCD

We investigate the low-temperature phase of QCD and the crossover region with two light flavors of quarks. The chiral expansion around the point $(T,m=0)$ in the temperature vs. quark-mass plane indicates that a sharp real-time excitation exists with the quantum numbers of the pion. An exact sum rule is derived for the thermal modification of the spectral function associated with the axial charge density; the (dominant) pion pole contribution obeys the sum rule. We determine the two parameters of the pion dispersion relation using lattice QCD simulations and test the applicability of the chiral expansion. The time-dependent correlators are also analyzed using the Maximum Entropy Method, yielding consistent results. Finally, we test the predictions of the chiral expansion around the point $(T=0,m=0)$ for the temperature dependence of static observables.

preprint2014arXiv

Chiral dynamics in the low-temperature phase of QCD

We investigate the low-temperature phase of QCD and the crossover region with two light flavors of quarks. The chiral expansion around the point $(T, m_q = 0)$ in the temperature vs. quark-mass plane indicates that a sharp real-time excitation exists with the quantum numbers of the pion. We determine its dispersion relation and test the applicability of the chiral expansion. The time-dependent correlators are also analyzed using the Maximum Entropy Method (MEM), yielding consistent results. Finally, we test the predictions of ordinary chiral perturbation theory around the point $(T = 0, m_q = 0)$ for the temperature dependence of static observables. Around the crossover temperature, we find that all quantities considered depend only mildly on the quark mass in the considered range 8MeV $\leq \bar{m}^{\bar{\text{MS}}} \leq$ 15MeV.

preprint2014arXiv

Hadronic contributions to the muon anomalous magnetic moment Workshop. $(g-2)_μ$: Quo vadis? Workshop. Mini proceedings

We present the mini-proceedings of the workshops Hadronic contributions to the muon anomalous magnetic moment: strategies for improvements of the accuracy of the theoretical prediction and $(g-2)_μ$: Quo vadis?, both held in Mainz from April 1$^{\rm rst}$ to 5$^{\rm th}$ and from April 7$^{\rm th}$ to 10$^{\rm th}$, 2014, respectively.

preprint2014arXiv

Lattice QCD Studies of the Leading Order Hadronic Contribution to the Muon $g-2$

The anomalous magnetic moment of the muon, $g_μ-2$, is one of the most promising observables to identify hints for physics beyond the Standard Model. QCD contributions are currently responsible for the largest fraction of the overall theoretical uncertainty in $g_μ-2$. The possibility to determine these hadronic contributions from first principles through lattice QCD calculations has triggered a number of recent studies. Recent proposals to improve the accuracy of lattice determinations are reported. We present an update of our studies of the leading-order hadronic contribution to $g_μ-2$ with improved Wilson fermions.

preprint2014arXiv

Study of the Couplings of QED and QCD from the Adler Function

The contribution from hadronic vacuum polarisation effects is responsible for a large fraction of the theoretical uncertainty in the running of the QED coupling. The current level of uncertainty has become a limitation for electroweak precision tests. We use lattice QCD simulations with two flavours of O$(a)$ improved Wilson fermions to determine the Adler function in a broad range of the momentum transfer $Q^2$. The running of the QED coupling, including valence contributions from $u$, $d$, $s$ and $c$ quarks, is compared to phenomenological results at intermediate $Q^2$ values. In the large $Q^2$ regime, the lattice determination of the Adler function is fitted to perturbation theory in order to examine the feasibility of a determination of the strong coupling constant.

preprint2014arXiv

The leading hadronic contribution to (g-2) of the muon: The chiral behavior using the mixed representation method

We extend our analysis of the leading hadronic contribution to the anomalous magnetic moment of the muon using the mixed representation method to study its chiral behavior. We present results derived from local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor ensembles made available to us through the CLS effort with pion masses as low as 190 MeV. The data is analyzed also using the more standard four-momentum method. Both methods are systematically compared as the calculations approach the physical point.

preprint2013arXiv

A new representation of the Adler function for lattice QCD

We address several aspects of lattice QCD calculations of the hadronic vacuum polarization and the associated Adler function. We implement a representation derived previously which allows one to access these phenomenologically important functions for a continuous set of virtualities, irrespective of the flavor structure of the current. Secondly we present a theoretical analysis of the finite-size effects on our particular representation of the Adler function, based on the operator product expansion at large momenta and on the spectral representation of the Euclidean correlator at small momenta. Finally, an analysis of the flavor structure of the electromagnetic current correlator is performed, where a recent theoretical estimate of the Wick-disconnected diagram contributions is rederived independently and confirmed.

preprint2013arXiv

Adler function and hadronic vacuum polarization from lattice vector correlators

We study a representation of the hadronic vacuum polarization based on the time-momentum representation of the vector correlator. This representation suggests a way to compute the hadronic vacuum polarization and the associated Adler function for any value of virtuality, irrespective of the flavor structure of the current. We present results on both of these phenomenologically important functions, derived from local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor ensembles made available to us through the CLS effort.

preprint2013arXiv

Antiscreening of the Ampere force in QED and QCD plasmas

The static forces between electric charges and currents are modified at the loop level by the presence of a plasma. While electric charges are screened, currents are not. The effective coupling constant at long distances is enhanced in both cases as compared to the vacuum, and by different amounts, a clear sign that Lorentz symmetry is broken. We investigate these effects quantitatively, first in a QED plasma and secondly using non-perturbative simulations of QCD with two light degenerate flavors of quarks.

preprint2013arXiv

Lattice QCD and the two-photon decay of the neutral pion

Two-photon decays probe the structure of mesons and represent an important contribution to hadronic light-by-light scattering. For the neutral pion, the decay amplitude tests the effects of the chiral anomaly; for a heavy quarkonium state, it measures the magnitude of its wavefunction at the origin. We rederive the expression of the decay amplitude in terms of a Euclidean correlation function starting from the theory defined on the torus. The derivation shows that for timelike photons the approach to the infinite-volume decay amplitude is exponential in the periodic box size.

preprint2013arXiv

QCD thermodynamics with O(a) improved Wilson fermions at Nf=2

We present an update of our study of the phase diagram of two-flavour QCD at zero baryon density with dynamical $O(a)$ improved Wilson quarks. All simulations are done on lattices with a temporal extent of $N_t=16$ and spatial extent $L=32,48$ and 64, ensuring that discretisation effects are small and finite size effects can be controlled. In the approach to the chiral limit we currently have three scans with pion masses between 540 and 200 MeV. In this proceedings article the focus is on the new scan at $m_π=200$ MeV and the measurement of screening masses. We also present first results concerning a test of scaling in the approach to the chiral limit and the chiral extrapolation of the difference of screening masses in scalar and pseudoscalar channels, which provides a measure for the strength of the anomalous breaking of the $U_A(1)$ symmetry.

preprint2013arXiv

Renormalization of the momentum density on the lattice using shifted boundary conditions

In order to extract transport quantities from energy-momentum-tensor (EMT) correlators in Lattice QCD there is a strong need for a non-perturbative renormalization of these operators. This is due to the fact that the lattice regularization explicitly breaks translational invariance, invalidating the non-renormalization-theorem. Here we present a non-perturbative calculation of the renormalization constant of the off-diagonal components of the EMT in SU(3) pure gauge theory using lattices with shifted boundary conditions. This allows us to induce a non-zero momentum in the system controlled by the shift parameter and to determine the normalization of the momentum density operator.

preprint2013arXiv

Theory Uncertainty in Extracting the Proton's Weak Charge: White Paper

We review the state-of-the-art and address open questions relative to the theory uncertainty of the $γ-Z$ box contribution to the $A_{PV}$ measurement within the QWEAK experiment at Jefferson Lab. This white paper summarizes the contributions by participants and discussion sessions on this topic within the MITP Workshop on Precision Electroweak Physics held in Mainz, Germany, September 23 - October 11, 2013 \url{http://indico.cern.ch/conferenceDisplay.py?confId=248072

preprint2013arXiv

Thermal field theories and shifted boundary conditions

The analytic continuation to an imaginary velocity of the canonical partition function of a thermal system expressed in a moving frame has a natural implementation in the Euclidean path-integral formulation in terms of shifted boundary conditions. The Poincare' invariance underlying a relativistic theory implies a dependence of the free-energy on the compact length L_0 and the shift xi only through the combination beta=L_0(1+xi^2)^(1/2). This in turn implies that the energy and the momentum distributions of the thermal theory are related, a fact which is encoded in a set of Ward identities among the correlators of the energy-momentum tensor. The latter have interesting applications in lattice field theory: they offer novel ways to compute thermodynamic potentials, and a set of identities to renormalize non-perturbatively the energy-momentum tensor. At fixed bare parameters the shifted boundary conditions also provide a simple method to vary the temperature in much smaller steps than with the standard procedure.

preprint2013arXiv

Two-flavour lattice QCD correlation functions in the deconfinement transition region

We report on a lattice QCD calculation with two dynamical flavors of the isovector vector correlator in the high-temperature phase. We analyze the correlator in terms of the associated spectral function by performing a fit for the difference of the thermal and vacuum spectral functions, using also an exact sum rule that constrains this difference. Additonally we carry out a direct fit for the thermal spectral function, and obtain good agreement between the two analyses for frequencies below the two-pion threshold. Under the assumption that the spectral function is smooth in that region, we give an estimate of the electrical conductivity.

preprint2013arXiv

Vector correlator and scale determination in lattice QCD

We implement a proposal made in [arXiv:1107.4388] to determine the lattice spacing by matching the lattice vector correlator at a reference distance scale with the same correlator obtained by a dispersion relation based on the $R$-ratio determined experimentally. We work with the isovector current, requiring a separation of the isovector hadronic final states on the phenomenological side. We also discuss the finite-size effect on the correlator, which must be controlled in order for the method to be applicable.

preprint2012arXiv

Excited state systematics in extracting nucleon electromagnetic form factors

We present updated preliminary results for the nucleon electromagnetic form factors for non-perturbatively $\mathcal{O}(a)$ improved Wilson fermions in $N_f=2$ QCD measured on the CLS ensembles. The use of the summed operator insertion method allows us to suppress the influence of excited states in our measurements. A study of the effect that excited state contaminations have on the $Q^2$ dependence of the extracted nucleon form factors may then be made through comparisons of the summation method to standard plateau fits, as well as to excited state fits.

preprint2012arXiv

Implications of Poincare symmetry for thermal field theories in finite-volume

The analytic continuation to an imaginary velocity $iξ$ of the canonical partition function of a thermal system expressed in a moving frame has a natural implementation in the Euclidean path-integral formulation in terms of shifted boundary conditions. Writing the Boltzmann factor as $\exp[-L_0(H-iξ.P)]$, the Poincare invariance underlying a relativistic theory implies a dependence of the free-energy on $L_0$ and the shift $ξ$ only through the combination $β= L_0 \sqrt{1+ξ^2}$. This in turn implies a set of Ward identities, some of which were previously derived by us, among the correlators of the energy-momentum tensor. In the infinite-volume limit they lead to relations among the cumulants of the total energy distribution and those of the momentum, i.e. they connect the energy and the momentum distributions in the canonical ensemble. In finite volume the Poincare symmetry translates into exact relations among partition functions and correlation functions defined with different sets of (generalized) periodic boundary conditions. They have interesting applications in lattice field theory. In particular, they offer Ward identities to renormalize non-perturbatively the energy-momentum tensor and novel ways to compute thermodynamic potentials. At fixed bare parameters they also provide a simple method to vary the temperature in much smaller steps than with the standard procedure.

preprint2012arXiv

Lattice QCD and the Timelike Pion Form Factor

We present a formula that allows one to calculate the pion form factor in the timelike region 2mpi <= sqrt{s} <= 4mpi in lattice QCD. The form factor quantifies the contribution of two-pion states to the vacuum polarization. It must be known very accurately in order to reduce the theoretical uncertainty on the anomalous magnetic moment of the muon. At the same time, the formula constitutes a rare example where, in a restricted kinematic regime, the spectral function of a conserved current can be determined from Euclidean observables without an explicit analytic continuation.

preprint2012arXiv

QCD thermodynamics with two flavours of Wilson fermions on large lattices

We explore the phase diagram of two flavour QCD at vanishing chemical potential using dynamical O(a)-improved Wilson quarks. In the approach to the chiral limit we use lattices with a temporal extent of N_t=16 and spatial extent L=32,48 and 64 to enable the extrapolation to the thermodynamic limit with small discretisation effects. In addition to an update on the scans at constant κ, reported earlier, we present first results from scans along lines of constant physics at a pion mass of 290 MeV. We probe the transition using the Polyakov loop and the chiral condensate, as well as spectroscopic observables such as screening masses.

preprint2012arXiv

Thermal Correlators in the ρ channel of two-flavor QCD

We present a lattice QCD calculation with two dynamical flavors of the isovector vector correlator in the high-temperature phase. We analyze the correlator in terms of the associated spectral function, for which we review the theoretical expectations. In our main analysis, we perform a fit for the difference of the thermal and vacuum spectral functions, and we use an exact sum rule that constrains this difference. We also perform a direct fit for the thermal spectral function, and obtain good agreement between the two analyses for frequencies below the two-pion threshold. Under the assumption that the spectral function is smooth in that region, we give an estimate of the electrical conductivity.

preprint2012arXiv

Vector Correlators in Lattice QCD: methods and applications

We discuss the calculation of the leading hadronic vacuum polarization in lattice QCD. Exploiting the excellent quality of the compiled experimental data for the e^+e^- --> hadrons cross-section, we predict the outcome of large-volume lattice calculations at the physical pion mass, and design computational strategies for the lattice to have an impact on important phenomenological quantities such as the leading hadronic contribution to (g-2)mu and the running of the electromagnetic coupling constant. First, the R(s) ratio can be calculated directly on the lattice in the threshold region, and we provide the formulae to do so with twisted boundary conditions. Second, the current correlator projected onto zero spatial momentum, in a Euclidean time interval where it can be calculated accurately, provides a potentially critical test of the experimental R(s) ratio in the region that is most relevant for (g-2)mu. This observation can also be turned around: the vector correlator at intermediate distances can be used to determine the lattice spacing in fm, and we make a concrete proposal in this direction. Finally, we quantify the finite-size effects on the current correlator coming from low-energy two-pion states and provide a general parametrization of the vacuum polarization on the torus.

preprint2011arXiv

Hadron Structure on the Lattice

A few chosen nucleon properties are described from a lattice QCD perspective: the nucleon sigma term and the scalar strangeness in the nucleon; the vector form factors in the nucleon, including the vector strangeness contribution, as well as parity breaking effects like the anapole and electric dipole moment; and finally the axial and tensor charges of the nucleon. The status of the lattice calculations is presented and their potential impact on phenomenology is discussed.

preprint2011arXiv

The errant life of a heavy quark in the quark-gluon plasma

In the high-temperature phase of QCD, the heavy quark momentum diffusion constant determines, via a fluctuation-dissipation relation, how fast a heavy quark kinetically equilibrates. This transport coefficient can be extracted from thermal correlators via a Kubo formula. We present a lattice calculation of the relevant Euclidean correlators in the gluon plasma, based on a recent formulation of the problem in heavy-quark effective field theory (HQET). We find a $\approx20%$ enhancement of the Euclidean correlator at maximal time separation as the temperature is lowered from $6T_c$ to $2T_c$, pointing to stronger interactions at lower temperatures. At the same time, the correlator becomes flatter from $6T_c$ down to $2T_c$, indicating a relative shift of the spectral weight to lower frequencies. A recent next-to-leading order perturbative calculation of the correlator agrees with the time dependence of the lattice data at the few-percent level. We estimate how much additional contribution from the $ω\lesssim T$ region of the perturbative spectral function would be required to bring it in agreement with the lattice data at $3.1T_c$.

preprint2011arXiv

Thermal momentum distribution from path integrals with shifted boundary conditions

For a thermal field theory formulated in the grand canonical ensemble, the distribution of the total momentum is an observable characterizing the thermal state. We show that its cumulants are related to thermodynamic potentials. In a relativistic system for instance, the thermal variance of the total momentum is a direct measure of the enthalpy. We relate the generating function of the cumulants to the ratio of (a) a partition function expressed as a Matsubara path integral with shifted boundary conditions in the compact direction, and (b) the ordinary partition function. In this form the generating function is well suited for Monte-Carlo evaluation, and the cumulants can be extracted straightforwardly. We test the method in the SU(3) Yang-Mills theory and obtain the entropy density at three different temperatures.

preprint2011arXiv

Thermodynamic potentials from shifted boundary conditions: the scalar-field theory case

In a thermal field theory, the cumulants of the momentum distribution can be extracted from the dependence of the Euclidean path integral on a shift in the fields built into the temporal boundary condition. When combined with the Ward identities associated with the invariance of the theory under the Poincare' group, thermodynamic potentials such as the entropy or the pressure can be directly inferred from the response of the system to the shift. Crucially the argument holds, up to harmless finite-size and discretization effects, even if translational and rotational invariance are broken to a discrete subgroup of finite shifts and rotations such as in a lattice box. The formulas are thus applicable at finite lattice spacing and volume provided the derivatives are replaced by their discrete counterpart, and no additive or multiplicative ultraviolet-divergent renormalizations are needed to take the continuum limit. In this paper we present a complete derivation of the relevant formulas in the scalar field theory, where several technical complications are avoided with respect to gauge theories. As a by-product we obtain a recursion relation among the cumulants of the momentum distribution, and formulae for finite-volume corrections to several well-known thermodynamic identities.

preprint2011arXiv

Transport Properties of the Quark-Gluon Plasma -- A Lattice QCD Perspective

Transport properties of a thermal medium determine how its conserved charge densities (for instance the electric charge, energy or momentum) evolve as a function of time and eventually relax back to their equilibrium values. Here the transport properties of the quark-gluon plasma are reviewed from a theoretical perspective. The latter play a key role in the description of heavy-ion collisions, and are an important ingredient in constraining particle production processes in the early universe. We place particular emphasis on lattice QCD calculations of conserved current correlators. These Euclidean correlators are related by an integral transform to spectral functions, whose small-frequency form determines the transport properties via Kubo formulae. The universal hydrodynamic predictions for the small-frequency pole structure of spectral functions are summarized. The viability of a quasiparticle description implies the presence of additional characteristic features in the spectral functions. These features are in stark contrast with the functional form that is found in strongly coupled plasmas via the gauge/gravity duality. A central goal is therefore to determine which of these dynamical regimes the quark-gluon plasma is qualitatively closer to as a function of temperature. We review the analysis of lattice correlators in relation to transport properties, and tentatively estimate what computational effort is required to make decisive progress in this field.

preprint2011arXiv

Twisted-mass reweighting for O(a) improved Wilson fermions

We test the reweighting of the quark determinant of O(a) improved Wilson fermions in the domain-decomposed hybrid Monte-Carlo algorithm. Specifically, we implement a reweighting in a twisted-mass parameter proposed by Palombi and Lüscher in $N_{\rm f}=2$ QCD. We find that at equal acceptance rate, the algorithm is significantly more stable on a $32\times64^3$ lattice upon switching on the reweighting parameter. At the same time, the reweighting factor does not fluctuate strongly and hence is under control. At equal statistics, the uncertainty on the pion correlator is comparable to the case of the standard, unreweighted algorithm.

preprint2010arXiv

Lattice Gauge Theory Sum Rule for the Shear Channel

An exact expression is derived for the $(ω,p)=0$ thermal correlator of shear stress in SU($N_c$) lattice gauge theory. I remove a logarithmic divergence by taking a suitable linear combination of the shear correlator and the correlator of the energy density. The operator product expansion shows that the same linear combination has a finite limit when $ω\to\infty$. It follows that the vacuum-subtracted shear spectral function vanishes at large frequencies at least as fast as $α_s^2(ω)$ and obeys a sum rule. The trace anomaly makes a potential contribution to the spectral sum rule which remains to be fully calculated, but which I estimate to be numerically small for $T\gtrsim 3T_c$. By contrast with the bulk channel, the shear channel spectral density is then overall enhanced as compared to the spectral density in vacuo.

preprint2010arXiv

Nucleon structure from mixed action calculations using 2+1 flavors of asqtad sea and domain wall valence fermions

We present high statistics results for the structure of the nucleon from a mixed-action calculation using 2+1 flavors of asqtad sea and domain wall valence fermions. We perform extrapolations of our data based on different chiral effective field theory schemes and compare our results with available information from phenomenology. We discuss vector and axial form factors of the nucleon, moments of generalized parton distributions, including moments of forward parton distributions, and implications for the decomposition of the nucleon spin.

preprint2010arXiv

Quantum Adiabatic Algorithms, Small Gaps, and Different Paths

We construct a set of instances of 3SAT which are not solved efficiently using the simplest quantum adiabatic algorithm. These instances are obtained by picking random clauses all consistent with two disparate planted solutions and then penalizing one of them with a single additional clause. We argue that by randomly modifying the beginning Hamiltonian, one obtains (with substantial probability) an adiabatic path that removes this difficulty. This suggests that the quantum adiabatic algorithm should in general be run on each instance with many different random paths leading to the problem Hamiltonian. We do not know whether this trick will help for a random instance of 3SAT (as opposed to an instance from the particular set we consider), especially if the instance has an exponential number of disparate assignments that violate few clauses. We use a continuous imaginary time Quantum Monte Carlo algorithm in a novel way to numerically investigate the ground state as well as the first excited state of our system. Our arguments are supplemented by Quantum Monte Carlo data from simulations with up to 150 spins.

preprint2010arXiv

The Bulk Channel in Thermal Gauge Theories

We investigate the thermal correlator of the trace of the energy-momentum tensor in the SU(3) Yang-Mills theory. Our goal is to constrain the spectral function in that channel, whose low-frequency part determines the bulk viscosity. We focus on the thermal modification of the spectral function, $ρ(ω,T)-ρ(ω,0)$. Using the operator-product expansion we give the high-frequency behavior of this difference in terms of thermodynamic potentials. We take into account the presence of an exact delta function located at the origin, which had been missed in previous analyses. We then combine the bulk sum rule and a Monte-Carlo evaluation of the Euclidean correlator to determine the intervals of frequency where the spectral density is enhanced or depleted by thermal effects. We find evidence that the thermal spectral density is non-zero for frequencies below the scalar glueball mass $m$ and is significantly depleted for $m\lesssimω\lesssim 3m$.

preprint2010arXiv

The transverse structure of the QCD string

The characterization of the transverse structure of the QCD string is discussed. We formulate a conjecture as to how the stress-energy tensor of the underlying gauge theory couples to the string degrees of freedom. A consequence of the conjecture is that the energy density and the longitudinal-stress operators measure the distribution of the transverse position of the string, to leading order in the string fluctuations, whereas the transverse-stress operator does not. We interpret recent numerical measurements of the transverse size of the confining string and show that the difference of the energy and longitudinal-stress operators is the appropriate probe to use when comparing with the next-to-leading order string prediction. Secondly we derive the constraints imposed by open-closed string duality on the transverse structure of the string. We show that a total of three independent `gravitational' form factors characterize the transverse profile of the closed string, and obtain the interpretation of recent effective string theory calculations: the square radius of a closed string of length β, defined from the slope of its gravitational form factor, is given by (d-1)/(2πσ)\log(β/(4r_0)) in d space dimensions. This is to be compared with the well-known result that the width of the open-string at mid-point grows as (d-1)/(2πσ) log(r/r_0). We also obtain predictions for transition form factors among closed-string states.

preprint2009arXiv

A Quantum Monte Carlo Method at Fixed Energy

In this paper we explore new ways to study the zero temperature limit of quantum statistical mechanics using Quantum Monte Carlo simulations. We develop a Quantum Monte Carlo method in which one fixes the ground state energy as a parameter. The Hamiltonians we consider are of the form $H=H_{0}+λV$ with ground state energy E. For fixed $H_{0}$ and V, one can view E as a function of $λ$ whereas we view $λ$ as a function of E. We fix E and define a path integral Quantum Monte Carlo method in which a path makes no reference to the times (discrete or continuous) at which transitions occur between states. For fixed E we can determine $λ(E)$ and other ground state properties of H.

preprint2009arXiv

Cutoff Effects on Energy-Momentum Tensor Correlators in Lattice Gauge Theory

We investigate the discretization errors affecting correlators of the energy-momentum tensor $T_{μν}$ at finite temperature in SU($N_c$) gauge theory with the Wilson action and two different discretizations of $T_{μν}$. We do so by using lattice perturbation theory and non-perturbative Monte-Carlo simulations. These correlators, which are functions of Euclidean time $x_0$ and spatial momentum ${\bf p}$, are the starting point for a lattice study of the transport properties of the gluon plasma. We find that the correlator of the energy $\int d^3x T_{00}$ has much larger discretization errors than the correlator of momentum $\int d^3x T_{0k}$. Secondly, the shear and diagonal stress correlators ($T_{12}$ and $T_{kk}$) require $\Nt\geq 8$ for the $Tx_0={1/2}$ point to be in the scaling region and the cutoff effect to be less than 10%. We then show that their discretization errors on an anisotropic lattice with $\as/\at=2$ are comparable to those on the isotropic lattice with the same temporal lattice spacing. Finally, we also study finite ${\bf p}$ correlators.

preprint2009arXiv

High-Precision Thermodynamics and Hagedorn Density of States

We compute the entropy density of the confined phase of QCD without quarks on the lattice to very high accuracy. The results are compared to the entropy density of free glueballs, where we include all the known glueball states below the two-particle threshold. We find that an excellent, parameter-free description of the entropy density between 0.7Tc and Tc is obtained by extending the spectrum with the exponential spectrum of the closed bosonic string.

preprint2008arXiv

Energy-momentum tensor correlators and viscosity

Collective flow has been observed in heavy ion collisions, with a large anisotropic component, and ideal hydrodynamic calculations had significant successful in describing the distribution of produced particles at the RHIC experiments. In order to account for this near ideal fluid behavior, the shear and bulk viscosity of the quark gluon plasma (QGP) must be computed from first principles in a regime where the QGP is not weakly coupled. In this talk I describe recent progress in computing energy-momentum tensor correlators on the lattice from which the viscosities can be extracted via Kubo formulas. I also show how to cumulate information from several channels, including at non-vanishing spatial momentum, in order to best constrain the viscosities. These methods should soon yield predictions at the higher temperatures that will be explored at the LHC experiments.

preprint2007arXiv

A calculation of the shear viscosity in SU(3) gluodynamics

We perform a lattice Monte-Carlo calculation of the two-point functions of the energy-momentum tensor at finite temperature in the SU(3) gauge theory. Unprecedented precision is obtained thanks to a multi-level algorithm. The lattice operators are renormalized non-perturbatively and the classical discretization errors affecting the correlators are corrected for. A robust upper bound for the shear viscosity to entropy density ratio is derived, eta/s < 1.0, and our best estimate is eta/s = 0.134(33) at T=1.65Tc under the assumption of smoothness of the spectral function in the low-frequency region.

preprint2005arXiv

Glueballs and the Pomeron -- a Lattice Study

We perform lattice calculations of the lightest J=0,2,4,6 glueball masses in the D=3+1 SU(3) gauge theory and extrapolate to the continuum limit. Assuming that these masses lie on linear Regge trajectories we find a leading glueball trajectory alpha(t) = 0.93(24) + 0.28(2)alpha_R' t, where alpha_R' \simeq 0.9 {GeV}^{-2} is the slope of the usual mesonic Regge trajectories. This glueball trajectory has an intercept and slope similar to that of the Pomeron. We contrast this with the situation in D=2+1 where the leading glueball Regge trajectory is found to have too small an intercept to be important for high-energy cross-sections. We attempt to interpret the observed states and trajectories in terms of open and closed string models of glueballs.

preprint2005arXiv

Vortices on the worldsheet of the QCD string

We investigate the properties of the QCD string in the Euclidean SU(N) pure gauge theory when the space-time dimensions transverse to it are periodic. From the point of view of an effective string theory, the string tension $σ$ and the low-energy constants $c_k$ of the theory are arbitrary functions of the sizes of the transverse dimensions L_p. Since the gauge theory is linearly confining in D=2, 3 and 4 dimensions, we propose an effective string action for the flux-tube energy levels at any choice of $L_p$, given $σ(L_p)$ and $c_k(L_p)$. The Luscher term only depends on the number of massless bosonic degrees of freedom and the effective theory can account for its evolution as a function of $L_p$. As the size of one transverse dimension is varied, we predict a Kosterlitz-Thouless transition of the worldsheet field theory at $σ(L_p)L_p^2 \simeq 1/8π$ driven by vortices, after which the periodic component of the worldsheet displacement vector develops a mass gap and the effective central charge drops by one unit. The universal properties of the transition are emphasised.

preprint2003arXiv

Locality and Statistical Error Reduction on Correlation Functions

We propose a multilevel Monte-Carlo scheme, applicable to local actions, which is expected to reduce statistical errors on correlation functions. We give general arguments to show how the efficiency and parameters of the algorithm are determined by the low-energy spectrum. As an application, we measure the euclidean-time correlation of pairs of Wilson loops in SU(3) pure gauge theory with constant relative errors. In this case the ratio of the new method's efficiency to the standard one increases as exp{m_0t/2}, where m_0 is the mass gap and t the time separation.