Source author record

Hartmut Wittig

Hartmut Wittig 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

63works
4topics
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

63 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

Leading isospin breaking effects in the HVP contribution to $a_μ$ and to the running of $α$

The anomalous magnetic moment of the muon $a_μ$ and the running of the electromagnetic coupling $α$ play a fundamental role in beyond Standard Model (SM) physics searches. Non-perturbative hadronic contributions to both quantities, which are related to the hadronic vacuum polarization (HVP) function consisting of two electromagnetic currents, are a main source of uncertainty in the SM prediction. We compute the HVP function in lattice QCD+QED applying the time-momentum representation method. We expand the relevant correlation functions around the isosymmetric limit. In particular, we focus on leading isospin breaking effects taking quark-connected contributions into account, which we evaluate on isosymmetric $N_{\mathrm{f}}=2+1$ QCD gauge ensembles generated by the CLS initiative with non-perturbatively $O(a)$-improved Wilson fermions.

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.

preprint2019arXiv

Rho resonance, timelike pion form factor, and implications for lattice studies of the hadronic vacuum polarisation

We study isospin-1 P-wave $ππ$ scattering in lattice QCD with two flavours of O($a$) improved Wilson fermions. For pion masses ranging from $m_π=265$ MeV to $m_π=437$ MeV, we determine the energy spectrum in the centre-of-mass frame and in three moving frames. We obtain the scattering phase shifts using Lüscher's finite-volume quantisation condition. Fitting the dependence of the phase shifts on the scattering momentum to a Breit-Wigner form allows us to determine the resonance parameters $m_ρ$ and $g_{ρππ}$. By combining the scattering phase shifts with the decay matrix element of the vector current, we calculate the timelike pion form factor, $F_π$, and compare the results to the Gounaris-Sakurai representation of the form factor in terms of the resonance parameters. In addition, we fit our data for the form factor to the functional form suggested by the Omnès representation, which allows for the extraction of the charge radius of the pion. As a further application, we discuss the long-distance behaviour of the vector correlator, which is dominated by the two-pion channel. We reconstruct the long-distance part in two ways: one based on the finite-volume energies and matrix elements and the other based on $F_π$. It is shown that this part can be accurately constrained using the reconstructions, which has important consequences for lattice calculations of the hadronic vacuum polarisation contribution to the muon anomalous magnetic moment.

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

Nucleon matrix elements from lattice QCD with all-mode-averaging and a domain-decomposed solver: an exploratory study

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.

preprint2016arXiv

The hadronic vacuum polarization function with O(a)-improved Wilson fermions - an update

We present an update of our lattice QCD study of the vacuum polarization function using O$(a)$-improved $N_ {\rm f} =2$ Wilson fermions with increased statistics and a large set of momenta. The resulting points are highly correlated and thus require a correlated fitting procedure. We employ an extended frequentist method to estimate the systematic uncertainties due to the momentum dependence and to the continuum and chiral extrapolations. We present preliminary results of the leading order hadronic contribution of the anomalous magnetic moment of the muon $\left(a_μ^{\mathrm{HLO}}\right)$ at the physical point for $u,d,s$ and $c$ valence quarks.

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

Search for the H-Dibaryon in two flavor Lattice QCD

We present preliminary results from a lattice QCD calculation of the H-dibaryon using two flavors of $\mathcal{O}(a)$ improved Wilson fermions. We employ local six-quark interpolating operators at the source with a combination of local six-quark and two-baryon operators at the sink with the appropriate quantum numbers of the H-dibaryon and its coupling to the two-baryon channels. We find that the two-baryon operators provide an improved overlap onto the ground state in comparison to the local six-quark operators. We also apply Lüscher's finite volume formalism to obtain information on the nature of the infinite-volume interaction of two particles. Further, the momentum projection to three moving frames enables the isolation of the pole in the infinite-volume scattering amplitude. Preliminary results at pion masses of 450 MeV and 1 GeV clearly indicate the presence of states below the $ΛΛ$ threshold while a finite-volume analysis fails to conclusively show the existence of an infinite-volume bound state.

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 scalar radius of the pion from Lattice QCD in the continuum limit

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^π_{\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.

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

Search for a bound H-dibaryon using local six-quark interpolating operators

We present early results from a lattice QCD study seeking a bound $H$-dibaryon using $N_f=2$ flavors of $O(a)$ improved Wilson fermions and a quenched strange quark. We compute a matrix of two-point functions using operators consisting of the two independent local products of six positive-parity-projected quarks with the appropriate quantum numbers, which belong to the singlet and 27-plet irreducible representations of flavor SU(3). To expand this basis, we also independently vary the quark-field smearing, and apply a new scheme to reduce the noise caused by smearing. We then find the ground-state mass by solving the generalized eigenvalue problem. We show results from ensembles with pion masses 451 MeV and 1 GeV, and compare with other lattice calculations.

preprint2014arXiv

Study of the anomalous magnetic moment of the muon computed from the Adler function

We compute the Adler function on the lattice from vacuum polarization data with twisted boundary conditions using numerical derivatives. The study is based on CLS ensembles with two flavours of $O(a)$ improved Wilson fermions. We extrapolate the lattice data for the Adler function to the continuum limit and to the physical pion mass and analyze its dependence on the momentum transfer. We discuss the application of this method to the extraction of the $u,d$ contribution to $a_μ^{\mathrm{HLO}}$.

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 disconnected contribution to the anomalous magnetic moment of the muon

The hadronic vacuum polarization can be determined from the vector correlator in a mixed time-momentum representation. We explicitly calculate the disconnected contribution to the vector correlator, both in the $N_f = 2$ theory and with an additional quenched strange quark, using non-perturbatively $O(a)$-improved Wilson fermions. All-to-all propagators are computed using stochastic sources and a generalized hopping parameter expansion. Combining the result with the dominant connected contribution, we are able to estimate an upper bound for the systematic error that arises from neglecting the disconnected contribution in the determination of $(g-2)_μ$.

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.

preprint2014arXiv

The scalar pion form factor in two-flavor lattice QCD

We calculate the scalar form factor of the pion using two dynamical flavors of non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions, including both the connected and the disconnected contribution to the relevant correlation functions. We employ the calculation of all-to-all propagators using stochastic sources and a generalized hopping parameter expansion. From the form factor data at vanishing momentum transfer, $Q^2=0$, and two non-vanishing $Q^2$ we obtain an estimate for the scalar radius $\left<r^2\right>^π_{_{\rm S}}$ of the pion at one value of the lattice spacing and for five different pion masses. Using Chiral Perturbation Theory at next-to-leading order, we find $\left<r^2\right>^π_{_{\rm S}}=0.635\pm0.016$ fm$^2$ at the physical pion mass (statistical error only). This is in good agreement with the phenomenological estimate from $ππ$-scattering. The inclusion of the disconnected contribution is essential for achieving this level of agreement.

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

Computing the Adler function from the vacuum polarization function

We use a lattice determination of the hadronic vacuum polarization tensor to study the associated Ward identities and compute the Adler function. The vacuum polarization tensor is computed from a combination of point-split and local vector currents, using two flavours of O($a$)-improved Wilson fermions. Partially twisted boundary conditions are employed to obtain a fine momentum resolution. The modifications of the Ward identities by lattice artifacts and by the use of twisted boundary conditions are monitored. We determine the Adler function from the derivative of the vacuum polarization function over a large region of momentum transfer $q^2$. As a first account of systematic effects, a continuum limit scaling analysis is performed in the large $q^2$ regime.

preprint2013arXiv

Numerical studies of Minimally Doubled Fermions

We have performed the first numerical study of minimally doubled fermions of the Karsten-Wilczek class in the quenched approximation. This requires fixing the counterterms, which arise due to hypercubic symmetry breaking induced by the Karsten-Wilczek term. Non-perturbative renormalisation criteria are formulated after a detailed study of the parameter dependence of mesonic observables. Minimisation of the mass anisotropy of the pseudoscalar ground state fixes non-perturbative renormalisation conditions for the counterterm coefficients. These anisotropies are mapped out by probing different euclidean components of the transfer matrix through calculations of the pseudoscalar ground state mass in different directions. The chiral behaviour of the pseudoscalar ground state is studied with the tuned Karsten-Wilczek action for multiple lattice spacings. Light pseudoscalar masses ($ M_{PS} \lesssim 250\,MeV $) were achieved in the quenched approximation without encountering exceptional configurations. The presence of quenched chiral logarithms is studied under the tentative assumption of Goldstone Boson-like behaviour.

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

The electromagnetic form factor of the pion in two-flavour lattice QCD

We present the current status of our lattice calculation of the electromagnetic form factor of the pion with two flavours of non-perturbatively O(a)-improved Wilson fermions. Using twisted boundary conditions and stochastic sources we obtain accurate results with a fine momentum resolution near $q^2=0$. This enables the computation of the charge radius without model dependence. The ensembles cover various lattice spacings and pion masses, ranging down to 250 MeV. This allows to compare the data to continuum chiral perturbation theory to NNLO including corrections of finite lattice spacing to perform a simultaneous chiral and continuum extrapolation. An estimate for the systematic error resulting from the extrapolation can be obtained by looking at the spread of results obtained from other functional forms such as polynomials.

preprint2013arXiv

The pion vector form factor from lattice QCD and NNLO chiral perturbation theory

We present a comprehensive study of the electromagnetic form factor, the decay constant and the mass of the pion computed in lattice QCD with two degenerate O(a)-improved Wilson quarks at three different lattice spacings in the range 0.05-0.08fm and pion masses between 280 and 630MeV at m_pi L >~ 4. Using partially twisted boundary conditions and stochastic estimators, we obtain a dense set of precise data points for the form factor at very small momentum transfers, allowing for a model-independent extraction of the charge radius. Chiral Perturbation Theory (ChPT) augmented by terms which model lattice artefacts is then compared to the data. At next-to-leading order the effective theory fails to produce a consistent description of the full set of pion observables but describes the data well when only the decay constant and mass are considered. By contrast, using the next-to-next-to-leading order expressions to perform global fits result in a consistent description of all data. We obtain <r^2_pi>=0.481(33)(13)fm^2 as our final result for the charge radius at the physical point. Our calculation also yields estimates for the pion decay constant in the chiral limit, F_pi/F=1.080(16)(6), the quark condensate, Sigma^{1/3}_MSbar(2GeV)=261(13)(1)MeV and several low-energy constants of SU(2) ChPT.

preprint2013arXiv

The Shape of Covariantly Smeared Sources in Lattice QCD

Covariantly smeared sources are commonly used in lattice QCD to enhance the projection onto the ground state. Here we investigate the dependence of their shape on the gauge field background and find that the presence of localized concentrations of magnetic field can lead to strong distortions which reduce the smearing radii achievable by iterative smearing prescriptions. In particular, as $a\to 0$, iterative procedures like Jacobi smearing require increasingly large iteration counts in order to reach physically-sized smearing radii $r_{sm}\sim$ 0.5 fm, and the resulting sources are strongly distorted. To bypass this issue, we propose a covariant smearing procedure (``free-form smearing'') that allows us to create arbitrarily shaped sources, including in particular Gaussians of arbitrary radius.

preprint2013arXiv

Two-Baryon Correlation Functions in 2-flavour QCD

We present first results for two-baryon correlation functions, computed using $N_f=2$ flavours of O($a$) improved Wilson quarks, with the aim of explaining potential dibaryon bound states, specifically the H-dibaryon. In particular, we use a GEVP to isolate the groundstate using two-baryon (hyperon-hyperon) correlation functions $\big(\langle C_{XY}(t)C_{XY}(0) \rangle$, where $XY=ΛΛ, ΣΣ, NΞ, \cdots\big)$, each of which has an overlap with the H-dibaryon. We employ a `blocking' algorithm to handle the large number of contractions, which may easily be extended to N-baryon correlation functions. We also comment on its application to the analysis of single baryon masses ($n$, $Λ$, $Ξ$, $\cdots$). This study is performed on an isotropic lattice with $m_π= 460$ MeV, $m_πL = 4.7$ and $a = 0.063$ fm.

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.

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

Lattice calculations of the leading hadronic contribution to (g-2)_mu

We report on our ongoing project to calculate the leading hadronic contribution to the anomalous magnetic moment of the muon a_mu^HLO using two dynamical flavours of non-perturbatively O(a) improved Wilson fermions. In this study, we changed the vacuum polarisation tensor to a combination of local and point-split currents which significantly reduces the numerical effort. Partially twisted boundary conditions allow us to improve the momentum resolution of the vacuum polarisation tensor and therefore the determination of the leading hadronic contribution to (g-2)_mu. We also extended the range of ensembles to include a pion mass below 200 MeV which allows us to check the non-trivial chiral behaviour of a_mu^HLO.

preprint2012arXiv

Low-energy particle physics and chiral extrapolations

In this review I discuss the role of chiral extrapolations for the determination of several phenomenologically relevant quantities, including light quark masses, meson decay constants and the axial charge of the nucleon. In particular, I investigate whether chiral extrapolations are sufficiently controlled in order to rightfully claim the accuracy which is quoted in recent compilations of these quantities. While this is the case for the masses of the light quarks and the ratio fK/fpi of decay constants, small inconsistencies in the chiral and continuum behaviour of individual decay constants fK and fpi, as well as the hadronic radii r0, r1 remain and must be clarified. In the case of the nucleon axial charge, gA, the chiral behaviour is still poorly understood due to the presence of other systematic effects.

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

Towards a precise lattice determination of the leading hadronic contribution to (g-2)_mu

We report on our computation of the leading hadronic contribution to the anomalous magnetic moment of the muon using two dynamical flavours of non-perturbatively O(a) improved Wilson fermions. The strange quark is introduced in the quenched approximation. Partially twisted boundary conditions are applied to improve the momentum resolution in the relevant integral. Our results, obtained at three different values of the lattice spacing, allow for a preliminary study of discretization effects. We explore a wide range of lattice volumes, namely 2 fm < L < 3 fm, with pion masses from 600 to 280 MeV and discuss different chiral extrapolations to the physical point. We observe a non-trivial dependence of a_mu(HLO) on m_pi especially for small pion masses. The final result, a_mu(HLO)=618(64)*10^(-10), is obtained by considering only the quark connected contribution to the vacuum polarization. We present a detailed analysis of systematic errors and discuss how they can be reduced in future simulations.

preprint2011arXiv

Calculation of the pion electromagnetic form factor from lattice QCD

We present a lattice calculation of the vector form factor of the pion for two flavours of non-perturbatively O(a) improved Wilson fermions. For the measurements we utilise the CLS ensembles which include various lattice spacings and pion masses down to about 250 MeV. To obtain a fine momentum resolution near zero momentum transfer (q^2) partially twisted boundary conditions are employed using several twist angles. Due to the fine resolution around q^2=0 we are able to determine the slope of the form factor and, in turn, extract the charge radius of the pion without any model dependence. The results for the form factor and the charge radius are then compared to chiral perturbation theory and phenomenological models which are used to extrapolate the results to the physical point.

preprint2011arXiv

Lattice Determination of the Anomalous Magnetic Moment of the Muon

We compute the leading hadronic contribution to the anomalous magnetic moment of the muon a_mu^HLO using two dynamical flavours of non-perturbatively O(a) improved Wilson fermions. By applying partially twisted boundary conditions we are able to improve the momentum resolution of the vacuum polarisation, an important ingredient for the determination of the leading hadronic contribution. We check systematic uncertainties by studying several ensembles, which allows us to discuss finite size effects and lattice artefacts. The chiral behavior of a_mu^HLO turns out to be non-trivial, especially for small pion masses.

preprint2011arXiv

Review of lattice results concerning low energy particle physics

We review lattice results relevant for pion and kaon physics with the aim of making them easily accessible to the particle physics community. Specifically, we review the determination of the light-quark masses, the form factor f_+(0), relevant for the semileptonic K -> pi transition at zero momentum transfer as well as the ratio f_K/f_pi of decay constants and discuss the consequences for the elements V_{us} and V_{ud} of the CKM matrix. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)_LxSU(2)_R and SU(3)_LxSU(3)_R Chiral Perturbation Theory and review the determination of the B_K parameter of neutral kaon mixing. We introduce quality criteria and use these when forming averages. Although subjective and imperfect, these criteria may help the reader to judge different aspects of current lattice computations. Our main results are summarized in section 1.2, but we stress the importance of the detailed discussion that underlies these results and constitutes the bulk of the present review.

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

Implementation of the Neuberger-Dirac operator on GPUs

Recent developments have shown that a lot can be gained for QCD simulations from GPU hardware. This can be exploited especially in the case of Ginsparg-Wilson fermions when the com putational costs are particularly high. In this work, we use the Neuberger-Dirac operator as our realisation of Ginsparg-Wilson fermions, which greatly facilitate lattice investigations of decays like $K \to ππ$. We report on the ongoing study of our GPU implementation of the Neuberger-Dirac operator including the exact treatment of the low lying eigenmodes of the Wilson-Dirac operator. Our benchmarks show that we achieve speed-up factors of around 23 and 16 in single and double precision, respectively.

preprint2010arXiv

Minimally doubled fermions and their renormalization

Minimally doubled fermions have been proposed as a strictly local discretization of the QCD quark action, which also preserves chiral symmetry at finite cut-off. We study the renormalization and mixing properties of two particular realizations of minimally doubled fermions in lattice perturbation theory at one loop, and we construct conserved axial currents which have a simple form involving only nearest-neighbours sites. We also introduce a notation which allows a unified description of the renormalization properties of both actions.

preprint2010arXiv

Renormalization of minimally doubled fermions

We investigate the renormalization properties of minimally doubled fermions, at one loop in perturbation theory. Our study is based on the two particular realizations of Borici-Creutz and Karsten-Wilczek. A common feature of both formulations is the breaking of hyper-cubic symmetry, which requires that the lattice actions are supplemented by suitable counterterms. We show that three counterterms are required in each case and determine their coefficients to one loop in perturbation theory. For both actions we compute the vacuum polarization of the gluon. It is shown that no power divergences appear and that all contributions which arise from the breaking of Lorentz symmetry are cancelled by the counterterms. We also derive the conserved vector and axial-vector currents for Karsten-Wilczek fermions. Like in the case of the previously studied Borici-Creutz action, one obtains simple expressions, involving only nearest-neighbour sites. We suggest methods how to fix the coefficients of the counterterms non-perturbatively and discuss the implications of our findings for practical simulations.

preprint2010arXiv

Systematic errors in extracting nucleon properties from lattice QCD

Form factors of the nucleon have been extracted from experiment with high precision. However, lattice calculations have failed so far to reproduce the observed dependence of form factors on the momentum transfer. We have embarked on a program to thoroughly investigate systematic effects in lattice calculation of the required three-point correlation functions. Here we focus on the possible contamination from higher excited states and present a method which is designed to suppress them. Its effectiveness is tested for several baryonic matrix elements, different lattice sizes and pion masses.

preprint2010arXiv

Towards the N_f=2 deconfinement transition temperature with O(a) improved Wilson fermions

A lot of effort in lattice simulations over the last years has been devoted to studies of the QCD deconfinement transition. Most state-of-the-art simulations use rooted staggered fermions, while Wilson fermions are affected by large systematic uncertainties, such as coarse lattices or heavy sea quarks. Here we report on an ongoing study of the transition, using two degenerate flavours of nonperturbatively $\Ord(a)$ improved Wilson fermions. We start with $N_{t}=12$ and 16 lattices and pion masses of 600 to 450 MeV, aiming at chiral and continuum limits with light quarks.

preprint2010arXiv

Towards the N_f=2 deconfinement transition temperature with O(a) improved Wilson fermions: An update

We give an update on our current project to determine the transition temperature and the order of the deconfinement transition in the chiral limit of two flavour QCD. We use nonperturbatively O(a) improved Wilson fermions of the Sheikholeslami-Wohlert type, employing the efficient deflation accelerated DDHMC algorithm. We start at lattices with N_t>=12 and pion masses below 600 MeV, aiming at chiral and continuum limits with light quarks.

preprint2009arXiv

Minimally doubled fermions at one loop

Minimally doubled fermions have been proposed as a cost-effective realization of chiral symmetry at non-zero lattice spacing. Using lattice perturbation theory at one loop, we study their renormalization properties. Specifically, we investigate the consequences of the breaking of hyper-cubic symmetry, which is a typical feature of this class of fermionic discretizations. Our results for the quark self-energy indicate that the four-momentum undergoes a renormalization which contains a linearly divergent piece. We also compute renormalization factors for quark bilinears, construct the conserved vector and axial-vector currents and verify that at one loop the renormalization factors of the latter are equal to one.

preprint2009arXiv

Minimally doubled fermions at one-loop level

Single fermionic degrees of freedom together with standard chiral symmetry at finite lattice spacing, correct continuum limit and local interactions only are precluded by the Nielsen-Ninomiya no-go theorem. The class of minimally doubled fermion actions exhibits exactly two chiral modes. Recent interest in these actions has been sparked by the investigation of fermionic actions defined on "hyperdiamond" lattices. Due to the necessity of breaking hypercubic symmetry explicitly, radiative corrections generate operator mixings with relevant and marginal operators that should vanish in continuum QCD. These cannot be avoided and must be taken into account in particular by a peculiar wave-function renormalisation and additive momentum renormalisation. Renormalisation properties at one-loop level of the self-energy, local bilinears and conserved vector and axial-vector currents are presented for Borici-Creutz and Karsten-Wilczek actions. Distinct differences and similarities between both actions are elucidated.

preprint2003arXiv

Status of lattice calculations of B-meson decays and mixing

The present status of lattice calculations of f_B, B_B and SU(3) flavour breaking ratios such as f_Bs/f_B is reviewed. Particular attention is devoted to systematic uncertainties, such as those arising from the lack of simulation data for dynamical quarks with realistic masses, and the related difficulties associated with chiral extrapolations. Global averages for decay constants and mixing parameters are presented (Table 1), and the procedures to obtain them are discussed in detail.

preprint2002arXiv

Panel discussion on chiral extrapolation of physical observables

This is an approximate reconstruction of the panel discussion on chiral extrapolation of physical observables. The session consisted of brief presentations from panelists, followed by responses from the panel, and concluded with questions and comments from the floor with answers from panelists. In the following, the panelists have summarized their statements, and the ensuing discussion has been approximately reconstructed from notes.

preprint2001arXiv

Non-perturbative renormalization of the quark condensate in Ginsparg-Wilson regularizations

We present a method to compute non-perturbatively the renormalization constant of the scalar density for Ginsparg-Wilson fermions. It relies on chiral symmetry and is based on a matching of renormalization group invariant masses at fixed pseudoscalar meson mass, making use of results previously obtained by the ALPHA Collaboration for O(a)-improved Wilson fermions. Our approach is quite general and enables the renormalization of scalar and pseudoscalar densities in lattice regularizations that preserve chiral symmetry and of fermion masses in any regularization. As an application we compute the non-perturbative factor which relates the renormalization group invariant quark condensate to its bare counterpart, obtained with overlap fermions at beta=5.85 in the quenched approximation.

preprint2001arXiv

Scalar condensate and light quark masses from overlap fermions

We have studied pseudoscalar correlation functions computed using the overlap operator. Within the accuracy of our calculation we find that the quark mass dependence agrees with the prediction of lowest-order Chiral Perturbation Theory (ChPT) for quark masses in the range of m ~ m_s/2-2m_s. We present the results of an analysis which assumes lowest-order ChPT to be valid to extract the low-energy constants Sigma and f_P, as well as the strange quark mass. Non-perturbative renormalization is implemented via a matching procedure with data obtained using Wilson fermions in the Schroedinger functional set-up. We find that the scalar condensate computed here agrees with the one obtained previously through a finite-size scaling analysis.

preprint2000arXiv

Effective Chiral Lagrangians and Lattice QCD

We propose a general method to obtain accurate estimates for some of the "low-energy constants" in the one-loop effective chiral Lagrangian by means of simulating lattice QCD. In particular, the method is sensitive to those constants whose values are required to test the hypothesis of a massless up-quark. Initial tests performed in the quenched approximation confirm that good statistical precision can be achieved. As a byproduct we obtain an accurate estimate for the ratio of pseudoscalar decay constants, F_K/F_pi, in the quenched approximation, which lies 10% below the experimental result. The quantities that serve to extract the low-energy constants also allow a test of the scaling behaviour of different discretizations of QCD and a search for the effects of dynamical quarks.

preprint1999arXiv

Hadron masses and matrix elements from the QCD Schr"odinger functional

We explain how masses and matrix elements can be computed in lattice QCD using Schr"odinger functional boundary conditions. Numerical results in the quenched approximation demonstrate that good precision can be achieved. For a statistical sample of the same size, our hadron masses have a precision similar to what is achieved with standard methods, but for the computation of matrix elements such as the pseudoscalar decay constant the Schr"odinger functional technique turns out to be much more efficient than the known alternatives.

preprint1996arXiv

Non-perturbative determination of the axial current normalization constant in O(a) improved lattice QCD

A finite-size technique is employed to compute the normalization constant $Z_A$ of the isovector axial current in lattice QCD. The calculation is carried out in the quenched approximation for values of the bare gauge coupling $g_0$ ranging from 0 to 1. In the lattice action and the lattice expression for the axial current we include the counterterms required for O(a) improvement, with non-perturbatively determined coefficients. With little additional work the normalization constant $Z_V$ of the improved isospin current is also obtained.