Improved Wilson Quarks Research Articles

Nucleon form factors in Nf=2+1 lattice QCD at the physical point: Finite lattice spacing effect on the root-mean-square radii

We present results for the nucleon form factors: electric (GE), magnetic (GM), axial (FA), induced pseudoscalar (FP), and pseudoscalar (GP) form factors, using the second PACS10 ensemble that is one of three sets of 2+1 flavor lattice QCD configurations at physical quark masses in large spatial volumes [exceeding (10 fm)3]. The second PACS10 gauge configurations are generated by the PACS Collaboration with the six stout-smeared O(a) improved Wilson quark action and Iwasaki gauge action at the second gauge coupling β=2.00 corresponding to the lattice spacing of a=0.063 fm. We determine the isovector electric, magnetic and axial radii, and magnetic moment from the corresponding form factors, as well as the axial-vector coupling gA. Combining our previous results for the coarser lattice spacing [E. Shintani , ; ], the finite lattice spacing effects on the isovector radii, magnetic moment, and axial-vector coupling are investigated using the difference between the two results. It was found that the effect on gA is kept smaller than the statistical error of 2% while the effect on the isovector radii was observed as a possible discretization error of about 10%, regardless of the channel. We also report the partially conserved axial-vector current relation using a set of nucleon three-point correlation functions in order to verify the effect by O(a) improvement of the axial-vector current. Published by the American Physical Society 2024

Hadronic vacuum polarization in the muon g − 2: the short-distance contribution from lattice QCD

We present results for the short-distance window observable of the hadronic vacuum polarization contribution to the muon g – 2, computed via the time-momentum representation (TMR) in lattice QCD. A key novelty of our calculation is the reduction of discretization effects by a suitable subtraction applied to the TMR kernel function, which cancels the leading \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${x}_{0}^{4}$$\\end{document}-behaviour at short distances. To compensate for the subtraction, one must substitute a term that can be reliably computed in perturbative QCD. We apply this strategy to our data for the vector current collected on ensembles generated with 2 + 1 flavours of O(a)-improved Wilson quarks at six values of the lattice spacing and pion masses in the range 130 – 420 MeV. Our estimate at the physical point contains a full error budget and reads \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\left({a}_{\\mu }^{{\ ext{hvp}}}\\right)}^{{\ ext{SD}}}$$\\end{document} = 68.85(14)stat (42)syst·10−10, which corresponds to a relative precision of 0.7%. We discuss the implications of our result for the observed tensions between lattice and data-driven evaluations of the hadronic vacuum polarization.

Aspects of chiral symmetry in QCD at T=128 MeV

We investigate several aspects of chiral symmetry in QCD at a temperature of $T=128\text{ }\text{ }\mathrm{MeV}$. The study is based on a $24\ifmmode\times\else\texttimes\fi{}9{6}^{3}$ lattice-QCD ensemble with $\mathrm{O}(a)$-improved Wilson quarks and physical up, down and strange quark masses. The pion quasiparticle turns out to be significantly lighter than the zero-temperature pion mass, even though the corresponding static correlation length is shorter. We perform a quantitative comparison of our findings to predictions of chiral perturbation theory. Among several order parameters for chiral symmetry restoration, we compute the difference of the vector- and axial-vector time-dependent correlators and find it to be reduced by a factor $\ensuremath{\sim}2/3$ as compared to its vacuum counterpart.

Nucleon isovector couplings in Nf=2+1 lattice QCD at the physical point

We present results for the scalar and tensor isovector-couplings ($g_S$ and $g_T$) of the nucleon measured at the physical point ($M_{\pi}=135$ MeV) with a single lattice spacing of $0.085\ \mathrm{fm}$ in 2+1 flavor QCD. Our calculations are carried out with two ensembles of gauge configurations generated by the PACS Collaboration with nonperturbatively ${\cal O}(a)$ improved Wilson quark action and Iwasaki gauge action on $(10.9\ {\rm fm})^4$ and $(5.5\ {\rm fm})^4$ lattices, where the finite-size effect on the nucleon mass was not shown at the level of the statistical precision less than 0.5%. We compute the nucleon three-point correlation functions in the vector, axial, scalar, and tensor channels. We confirm that our previous result of the nucleon axial coupling on the large spatial volume of $(10.9\ {\rm fm})^4$ has no finite-size effect at the level of the statistical precision of 1.9%. For the renormalization, we first renormalize $g_S$ and $g_T$ nonperturbatively using the RI/SMOM$_{(\gamma_\mu)}$ scheme, a variant of Rome-Southampton RI/MOM scheme with reduced systematic errors, as the intermediate scheme. We evaluate our final results at the renormalization scale of 2 GeV in the $\overline{\rm MS}$ scheme through matching procedure between the RI/SMOM$_{(\gamma_\mu)}$ and $\overline{\rm MS}$ schemes with the help of perturbation theory, and then obtain $g_S=0.927(71)_{\rm stat}(22)_{\rm syst}$ and $g_T=1.036(6)_{\rm stat}(20)_{\rm syst}$.

Kℓ3 form factors at the physical point: Toward the continuum limit

We present updated results for the form factors of the kaon semileptonic $(K_{\ell 3})$ decay process calculated with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on large volumes of more than (10 fm)$^4$. In addition to our previous calculation at the lattice spacing $a = 0.085$ fm, we perform a calculation at the second lattice spacing of $0.063$ fm. Using the results for the form factors extracted from 3-point functions with the local and also conserved vector currents at the two lattice spacings, continuum extrapolation and interpolation of the momentum transfer are carried out simultaneously to obtain the value of the form factor $f_+(0)$ at the zero momentum transfer in the continuum limit. After investigation of stability of $f_+(0)$ against several fit forms and different data, we obtain $f_+(0) = 0.9615(10)(^{+47}_{\ -3})(5)$, where the first, second, and third errors are statistical, systematic errors from choice of the fit forms and isospin breaking effect, respectively. Combining our value of $f_+(0)$ and experimental input of the $K_{\ell 3}$ decay, one of the Cabibbo-Kobayashi-Maskawa matrix elements $|V_{us}|$ is determined as $|V_{us}| = 0.2252(^{\ +5}_{-12})$, whose error contains the experimental one as well as that in the lattice calculation. This value is reasonably consistent with the ones determined from recent lattice QCD results of $f_+(0)$ and also the one determined through the kaon leptonic decay process. We observe some tension between our value and $|V_{us}|$ evaluated from the unitarity of the CKM matrix with $|V_{ud}|$, while it depends on the size of the error of $|V_{ud}|$. It is also found that $|V_{us}|$ determined with our phase space integrals through six $K_{\ell 3}$ decay processes is consistent with the above one using $f_+(0)$.

Nucleon scalar and tensor couplings from lattice QCD at the physical points

We present results for the axial, scalar and tensor isovector-couplings (gA, gS and gT) of the nucleon obtained from 2+1 flavor QCD with the physical light quark masses (Mπ = 135 MeV). Our calculations are performed at a single lattice spacing of 0.085 fm, but with two large volumes of (10.9 fm)4 and (5.5 fm)4. The configurations are generated by the PACS Collaboration with nonperturbatively 𝒪(a) improved Wilson quark action and Iwasaki gauge action. The result of gA is a good indicator for determination of gS and gT with respect to accuracy and precision. Our result of gA well reproduces the experimental value within a statistical error of less than 2%. As for gS and gT, we compute the renormalization constants at the scale of 2 GeV in the ̅M̅S scheme through the RI/SMOM(γμ) intermediate scheme, and then obtain gS = 0.927(83)stat(22)syst and gT = 1.036(6)stat(20)syst.

Calculation of the derivative of nucleon form factors in Nf=2+1 lattice QCD at Mπ=138 MeV on a (5.5 fm)3 volume

We present a direct calculation for the first derivative of the isovector nucleon form factors with respect to the momentum transfer ${q}^{2}$ using the lower moments of the nucleon 3-point function in the coordinate space. Our numerical simulations are performed using the ${N}_{f}=2+1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action near the physical point, corresponding to the pion mass ${M}_{\ensuremath{\pi}}=138\text{ }\text{ }\mathrm{MeV}$, on a $(5.5\text{ }\text{ }\mathrm{fm}{)}^{4}$ lattice at a single lattice spacing of $a=0.085\text{ }\text{ }\mathrm{fm}$. In the momentum derivative approach, we can directly evaluate the mean square radii for the electric, magnetic, and axial-vector form factors, and also the magnetic moment without the ${q}^{2}$ extrapolation to the zero momentum point. These results are compared with the ones determined by the standard method, where the ${q}^{2}$ extrapolations of the corresponding form factors are carried out by fitting models. We find that the new results from the momentum derivative method are obtained with a larger statistical error than the standard method, but with a smaller systematic error associated with the data analysis. Within the total error range of the statistical and systematic errors combined, the two results are in good agreement. On the other hand, two variations of the momentum derivative of the induced pseudoscalar form factor at the zero momentum point show some discrepancy. It seems to be caused by a finite volume effect, since a similar trend is not observed on a large volume, but seen on a small volume in our pilot calculations at a heavier pion mass of ${M}_{\ensuremath{\pi}}=510\text{ }\text{ }\mathrm{MeV}$. Furthermore, we discuss an equivalence between the momentum derivative method and the similar approach with the point splitting vector current.

Erratum: Nf=2+1 QCD thermodynamics with gradient flow using two-loop matching coefficients [Phys. Rev. D 102 , 014510 (2020)

We study thermodynamic properties of Nf=2+1 QCD on the lattice adopting O(a)-improved Wilson quark action and Iwasaki gauge action. To cope with the problems due to explicit violation of the Poincare and chiral symmetries, we apply the Small Flow-time eXpansion (SFtX) method based on the gradient flow, which is a general method to correctly calculate any renormalized observables on the lattice. In this method, the matching coefficients in front of operators in the small flow-time expansion are calculated by perturbation theory. In a previous study using one-loop matching coefficients, we found that the SFtX method works well for the equation of state, chiral condensates and susceptibilities. In this paper, we study the effect of two-loop matching coefficients by Harlander et al. We also test the influence of the renormalization scale in the SFtX method. We find that, by adopting the mu_0 renormalization scale of Harlander et al. instead of the conventional mu_d=1/sqrt{8t} scale, the linear behavior at large t is improved so that we can perform the t -> 0 extrapolation of the SFtX method more confidently. In the calculation of the two-loop matching coefficients by Harlander et al., the equation of motion for quark fields was used. For the entropy density in which the equation of motion has no effects, we find that the results using the two-loop coefficients agree well with those using one-loop coefficients. On the other hand, for the trace anomaly which is affected by the equation of motion, we find discrepancies between the one- and two-loop results at high temperatures. By comparing the results of one-loop coefficients with and without using the equation of motion, the main origin of the discrepancies is suggested to be attributed to O((aT)^2)=O(1/N_t^2) discretization errors in the equation of motion at N_t =< 10.

Nf=2+1 QCD thermodynamics with gradient flow using two-loop matching coefficients

We study thermodynamic properties of Nf=2+1 QCD on the lattice adopting O(a)-improved Wilson quark action and Iwasaki gauge action. To cope with the problems due to explicit violation of the Poincare and chiral symmetries, we apply the Small Flow-time eXpansion (SFtX) method based on the gradient flow, which is a general method to correctly calculate any renormalized observables on the lattice. In this method, the matching coefficients in front of operators in the small flow-time expansion are calculated by perturbation theory. In a previous study using one-loop matching coefficients, we found that the SFtX method works well for the equation of state, chiral condensates and susceptibilities. In this paper, we study the effect of two-loop matching coefficients by Harlander et al. We also test the influence of the renormalization scale in the SFtX method. We find that, by adopting the mu_0 renormalization scale of Harlander et al. instead of the conventional mu_d=1/sqrt{8t} scale, the linear behavior at large t is improved so that we can perform the t -> 0 extrapolation of the SFtX method more confidently. In the calculation of the two-loop matching coefficients by Harlander et al., the equation of motion for quark fields was used. For the entropy density in which the equation of motion has no effects, we find that the results using the two-loop coefficients agree well with those using one-loop coefficients. On the other hand, for the trace anomaly which is affected by the equation of motion, we find discrepancies between the one- and two-loop results at high temperatures. By comparing the results of one-loop coefficients with and without using the equation of motion, the main origin of the discrepancies is suggested to be attributed to O((aT)^2)=O(1/N_t^2) discretization errors in the equation of motion at N_t =< 10.

Kl3 form factors at the physical point on a (10.9 fm)3 volume

We present the calculation of the $K_{l3}$ form factors with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on a large volume of (10.9 fm)$^3$ at one lattice spacing of $a = 0.085$ fm. We extract the form factors from 3-point functions with three different time separations between the source and sink operators to confirm suppression of excited state contributions. The form factors are calculated in very close to the zero momentum transfer, $q^2 = 0$, thanks to the large volume, so that stable interpolations to $q^2 = 0$ are carried out. Using our form factors, we obtain the form factor at $q^2 = 0$, $f_+(0) = 0.9603(16)(^{+14}_{\ -4})(44)(19)(1)$, where the first, second, and fifth errors are statistical, systematic errors from fit functions and the isospin breaking effect, respectively. The third and fourth errors denote the finite lattice spacing effects estimated from the renormalization factor and contribution beyond the leading order SU(3) chiral perturbation theory (ChPT). The result of $f_+(0)$ yields the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, $|V_{us}| = 0.2255(13)(4)$, where the first error comes from our calculation and the second from the experiment. This value is consistent with the ones determined from the unitarity of the CKM matrix and the $K_{l2}$ decay within one standard deviation, while it is slightly larger than recent lattice calculations by at most 1.5 $\sigma$. Furthermore, we evaluate the shape of the form factors and the phase space integral from our results. We confirm that those results are consistent with the experiment, and also $|V_{us}|$ determined with our phase space integral agrees with the one in the above.

Lattice calculation of the hadronic leading order contribution to the muon g − 2

The persistent discrepancy of about 3.5 standard deviations between the experimental measurement and the Standard Model prediction for the muon anomalous magnetic moment, aµ, is one of the most promising hints for the possible existence of new physics. Here we report on our lattice QCD calculation of the hadronic vacuum polarisation contribution $ a_\mu ^{{\rm{hvp}}} $, based on gauge ensembles with Nf = 2 + 1 flavours of O(a) improved Wilson quarks. We address the conceptual and numerical challenges that one encounters along the way to a sub-percent determination of the hadronic vacuum polarisation contribution. The current status of lattice calculations of $ a_\mu ^{{\rm{hvp}}} $ is presented by performing a detailed comparison with the results from other groups.

Charmonium resonances with JPC=1−− and 3−− from D¯D scattering on the lattice

We present a lattice QCD study of charmonium resonances and bound states with $J^{PC}=1^{--}$ and $3^{--}$ near the open-charm threshold, taking into account their strong transitions to $\bar DD$. Vector charmonia are the most abundant in the experimentally established charmonium spectrum, while recently LHCb reported also the first discovery of a charmonium with likely spin three. The $\bar DD$ scattering amplitudes for partial waves $l=1$ and $l=3$ are extracted on the lattice by means of the L\"uscher formalism, using multiple volumes and inertial frames. Parameterizations of the scattering amplitudes provide masses and widths of the resonances, as well as the masses of bound states. CLS ensembles with 2+1 dynamical flavors of non-perturbatively $O(a)$ improved Wilson quarks are employed with $m_\pi\simeq 280$ MeV, a single lattice spacing of $a\simeq0.086$ fm and two lattice spatial extents of $L=24$ and $32$. Two values of the charm quark mass are considered to examine the influence of the position of the $\bar{D}D$ threshold on the hadron masses. For the lighter charm quark mass we find the vector resonance $\psi(3770)$ with mass $m=3780(7)$ MeV and coupling $g=16.0(^{+2.1}_{-0.2})$ (related to the width), both consistent with their experimental values. The vector $\psi(2S)$ appears as a bound state with $m=3666(10)$ MeV. The charmonium resonance with $J^{PC}=3^{--}$ is found at $m=3831(^{+10}_{-16})$ MeV, consistent with the $X(3842)$ recently discovered by LHCb. At our heavier charm-quark mass the $\psi(2S)$ as well as the $\psi(3770)$ are bound states and the $X(3842)$ remains a resonance. We stress that all quoted uncertainties are only statistical, while lattice spacing effects and the approach to the physical point still need to be explored. This study of conventional charmonia sets the stage for more challenging future studies of unconventional charmonium-like states.

How much do charm sea quarks affect the charmonium spectrum?

The properties of charmonium states are or will be intensively studied by the B-factories Belle II and BESIII, the LHCb and PANDA experiments and at a future Super-c-tau Factory. Precise lattice calculations provide valuable input and several results have been obtained by simulating up, down and strange quarks in the sea. We investigate the impact of a charm quark in the sea on the charmonium spectrum, the renormalization group invariant charm–quark mass M_{mathrm{c}} and the scalar charm–quark content of charmonium. The latter is obtained by the direct computation of the mass-derivatives of the charmonium masses. We do this investigation in a model, QCD with two degenerate charm quarks. The absence of light quarks allows us to reach very small lattice spacings down to 0.023~hbox {fm}. By comparing to pure gauge theory we find that charm quarks in the sea affect the hyperfine splitting at a level around 2%. The most significant effects are 5% in M_c and 3% in the value of the charm quark content of the eta _c meson. Given that we simulate two charm quarks these estimates are upper bounds for the contribution of a single charm quark. We show that lattice spacings <0.06~hbox {fm} are needed for safe continuum extrapolations of the charmonium spectrum with O(a) improved Wilson quarks. A useful relation for the projection to the desired parity of operators in two-point functions computed with twisted mass fermions is proven.

Nucleon form factors on a large volume lattice near the physical point in 2+1 flavor QCD

We present results for the isovector nucleon form factors measured on a $96^4$ lattice at almost the physical pion mass with a lattice spacing of 0.085 fm in 2+1 flavor QCD. The configurations are generated with the stout-smeared $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $\beta$=1.82. The pion mass at the simulation point is about 146 MeV. A large spatial volume of $(8.1~{\rm fm})^3$ allows us to investigate the form factors in the small momentum transfer region. We determine the isovector electric radius and magnetic moment from nucleon electric ($G_E$) and magnetic ($G_M$) form factors as well as the axial-vector coupling $g_A$. We also report on the results of the axial-vector ($F_A$), induced pseudoscalar ($F_P$) and pseudoscalar ($G_P$) form factors in order to verify the axial Ward-Takahashi identity in terms of the nucleon matrix elements, which may be called the generalized Goldberger-Treiman relation.

Non-perturbative quark mass renormalisation and running in N_{scriptstyle mathrm{f}}=3 QCD

We determine from first principles the quark mass anomalous dimension in N_{scriptstyle mathrm{f}}=3 QCD between the electroweak and hadronic scales. This allows for a fully non-perturbative connection of the perturbative and non-perturbative regimes of the Standard Model in the hadronic sector. The computation is carried out to high accuracy, employing massless text{ O }(a)-improved Wilson quarks and finite-size scaling techniques. We also provide the matching factors required in the renormalisation of light quark masses from lattice computations with text{ O }(a)-improved Wilson fermions and a tree-level Symanzik improved gauge action. The total uncertainty due to renormalisation and running in the determination of light quark masses in the SM is thus reduced to about 1%.

Leptonic decay constants for D-mesons from 3-flavour CLS ensembles

e report on the status of an ongoing effort by the RQCD and ALPHA Collaborations, aimed at determining leptonic decay constants of charmed mesons. Our analysis is based on large-volume ensembles generated within the CLS effort, employing Nf = 2 + 1 non-perturbatively O(a) improved Wilson quarks, tree-level Symanzik-improved gauge action and open boundary conditions. The ensembles cover lattice spac-ings from a ≈ 0.09 fm to a ≈ 0.05 fm, with pion masses varied from 420 to 200 MeV. To extrapolate to the physical masses, we follow both the (2ml + ms) = const. and the ms = const. lines in parameter space.

Baryon interactions from lattice QCD with physical quark masses – Nuclear forces and ΞΞ forces –

We present the latest lattice QCD results for baryon interactions obtained at nearly physical quark masses. Nf = 2 + 1 nonperturbatively O(a)-improved Wilson quark action with stout smearing and Iwasaki gauge action are employed on the lattice of (96a)4 ≃(8.1fm)4 with a-1 ≃2.3 GeV, where mπ ≃146 MeV and mK ≃525 MeV. In this report, we study the two-nucleon systems and two-Ξ systems in 1S0 channel and 3S1-3D1 coupled channel, and extract central and tensor interactions by the HAL QCD method. We also present the results for the NΩ interaction in 5S2 channel which is relevant to the NΩ pair-momentum correlation in heavy-ion collision experiments.

Nucleon structure from 2+1 flavor lattice QCD near the physical point

We present an update on our results of nucleon form factors measured on a large-volume lattice (8:1fm)4 at almost the physical point in 2+1 flavor QCD. The configurations are generated with the stout-smeared O(a) improved Wilson quark action and Iwasaki gauge action at β = 1:82, which corresponds to the lattice spacing of 0.085 fm. The pion mass at the simulation point is about 145 MeV. We determine the isovector electric radius and magnetic moment from nucleon electric (GE) and magnetic (GM) form factors. We also report on preliminary results of the axial-vector (FA), induced pseudo-scalar (FP) and pseudo-scalar (GP) form factors in order to verify the axial Ward-Takahashi identity in terms of the nucleon matrix elements, which may be called as the generalized Goldberger-Treiman relation.

Equation of state in (2+1)-flavor QCD at physical point with improved Wilson fermion action using gradient flow

We study the energy-momentum tensor and the equation of state as well as the chiral condensate in (2+1)-flavor QCD at the physical point applying the method of Makino and Suzuki based on the gradient flow. We adopt a nonperturbatively O(a)- improved Wilson quark action and the renormalization group-improved Iwasaki gauge action. At Lattice 2016, we have presented our preliminary results of our study in (2+1)- flavor QCD at a heavy u; d quark mass point. We now extend the study to the physical point and perform finite-temperature simulations in the range T ≃ 155.544 MeV (Nt = 4-14 including odd Nt’s) at a ≃ 0:09 fm. We show our final results of the heavy QCD study and present some preliminary results obtained at the physical point so far.

Hadroquarkonium from lattice QCD

The hadro-quarkonium picture [S. Dubinskiy and M.B. Voloshin, Phys. Lett. B 666, 344 (2008)] provides one possible interpretation for the pentaquark candidates with hidden charm, recently reported by the LHCb Collaboration, as well as for some of the charmonium-like "X,Y,Z" states. In this picture, a heavy quarkonium core resides within a light hadron giving rise to four- or five-quark/antiquark bound states. We test this scenario in the heavy quark limit by investigating the modification of the potential between a static quark-antiquark pair induced by the presence of a hadron. Our lattice QCD simulations are performed on a Coordinated Lattice Simulations (CLS) ensemble with $N_f = 2+1$ flavours of non-perturbatively improved Wilson quarks at a pion mass of about 223 MeV and a lattice spacing of about 0.0854 fm. We study the static potential in the presence of a variety of light mesons as well as of octet and decuplet baryons. In all these cases, the resulting configurations are favoured energetically. The associated binding energies between the quarkonium in the heavy quark limit and the light hadron are found to be smaller than a few MeV, similar in strength to deuterium binding. It needs to be seen if the small attraction survives in the infinite volume limit and supports bound states or resonances.