Sea Quark Masses Research Articles

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q2∼4.3 GeV2 and the nucleon up to Q2∼1 GeV2. The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼250 to ∼540 MeV, including the unitary point ∼340 MeV. We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

Lattice calculation of pion form factors with overlap fermions

We present a precise calculation of the pion form factor using overlap fermions on seven ensembles of 2+1-flavor domain-wall configurations with pion masses varying from 139 to 340 MeV. Taking advantage of the fast Fourier transform and other techniques to access many combinations of source and sink momenta, we find the pion mean square charge radius to be $\langle {r_\pi^2} \rangle= 0.430(5)(13)\ {\rm{fm^2}}$, which agrees well with the experimental result, and includes the systematic uncertainties from chiral extrapolation, lattice spacing and finite-volume dependence. We also find that $\langle {r_\pi^2} \rangle$ depends on both the valence and sea quark masses strongly and predict the pion form factor up to $Q^2 = 1.0 \ {\rm{GeV^2}}$ which agrees with experiments very well.

Ratio of flavour non-singlet and singlet scalar density renormalisation parameters in N_mathrm {f}=3 QCD with Wilson quarks

We determine non-perturbatively the normalisation factor r_{mathrm{m}}equiv Z_{mathrm{S}}/Z_{mathrm{S}}^{0}, where Z_{mathrm{S}} and Z_{mathrm{S}}^{0} are the renormalisation parameters of the flavour non-singlet and singlet scalar densities, respectively. This quantity is required in the computation of quark masses with Wilson fermions and for instance the renormalisation of nucleon matrix elements of scalar densities. Our calculation involves simulations of finite-volume lattice QCD with the tree-level Symanzik-improved gauge action, N_{mathrm{f}}= 3 mass-degenerate {mathrm{O}}(a) improved Wilson fermions and Schrödinger functional boundary conditions. The slope of the current quark mass, as a function of the subtracted Wilson quark mass is extracted both in a unitary setup (where nearly chiral valence and sea quark masses are degenerate) and in a non-unitary setup (where all valence flavours are chiral and the sea quark masses are small). These slopes are then combined with Z equiv Z_{mathrm{P}}/(Z_{mathrm{S}}Z_{mathrm{A}}) in order to obtain r_{mathrm{m}}. A novel chiral Ward identity is employed for the calculation of the normalisation factor Z. Our results cover the range of gauge couplings corresponding to lattice spacings below 0.1,fm, for which N_{mathrm{f}}= 2+1 QCD simulations in large volumes with the same lattice action are typically performed.

Correlated Dirac Eigenvalues and Axial Anomaly in Chiral Symmetric QCD.

We introduce novel relations between the derivatives [∂^{n}ρ(λ,m_{l})/∂m_{l}^{n}] of the Dirac eigenvalue spectrum [ρ(λ,m_{l})] with respect to the light sea quark mass (m_{l}) and the (n+1)-point correlations among the eigenvalues (λ) of the massless Dirac operator. Using these relations we present lattice QCD results for ∂^{n}ρ(λ,m_{l})/∂m_{l}^{n} (n=1, 2, 3) for m_{l} corresponding to pion masses m_{π}=160-55 MeV and at a temperature of about 1.6 times the chiral phase transition temperature. Calculations were carried out using (2+1) flavors of highly improved staggered quarks with the physical value of strange quark mass, three lattice spacings a=0.12, 0.08, 0.06fm, and lattices having aspect ratios 4-9. We find that ρ(λ→0,m_{l}) develops a peaked structure. This peaked structure arises due to non-Poisson correlations within the infrared part of the Dirac eigenvalue spectrum, becomes sharper as a→0, and its amplitude is proportional to m_{l}^{2}. We demonstrate that this ρ(λ→0,m_{l}) is responsible for the manifestations of axial anomaly in two-point correlation functions of light scalar and pseudoscalar mesons. After continuum and chiral extrapolations we find that axial anomaly remains manifested in two-point correlation functions of scalar and pseudoscalar mesons in the chiral limit.

Splittings of low-lying charmonium masses at the physical point

We present high-precision results from lattice QCD for the mass splittings of the low-lying charmonium states. For the valence charm quark, the calculation uses Wilson-clover quarks in the Fermilab interpretation. The gauge-field ensembles are generated in the presence of up, down, and strange sea quarks, based on the improved staggered (asqtad) action, and gluon fields, based on the one-loop, tadpole-improved gauge action. We use five lattice spacings and two values of the light sea quark mass to extrapolate the results to the physical point. An enlarged set of interpolating operators is used for a variational analysis to improve the determination of the energies of the ground states in each channel. We present and implement a continuum extrapolation within the Fermilab interpretation, based on power-counting arguments, and thoroughly discuss all sources of systematic uncertainty. We compare our results for various mass splittings with their experimental values, namely, the 1S hyperfine splitting, the 1P-1S splitting and the P-wave spin-orbit and tensor splittings. Given the uncertainty related to the width of the resonances, we find excellent agreement.

Up-, down-, strange-, charm-, and bottom-quark masses from four-flavor lattice QCD

We calculate the up-, down-, strange-, charm-, and bottom-quark masses using the MILC highly improved staggered-quark ensembles with four flavors of dynamical quarks. We use ensembles at six lattice spacings ranging from $a\ensuremath{\approx}0.15$ to 0.03 fm and with both physical and unphysical values of the two light and the strange sea-quark masses. We use a new method based on heavy-quark effective theory (HQET) to extract quark masses from heavy-light pseudoscalar meson masses. Combining our analysis with our separate determination of ratios of light-quark masses we present masses of the up, down, strange, charm, and bottom quarks. Our results for the $\overline{\mathrm{MS}}$-renormalized masses are ${m}_{u}(2\text{ }\text{ }\mathrm{GeV})=\phantom{\rule{0ex}{0ex}}2.130(41)\text{ }\text{ }\mathrm{MeV}$, ${m}_{d}(2\text{ }\text{ }\mathrm{GeV})=4.675(56)\text{ }\text{ }\mathrm{MeV}$, ${m}_{s}(2\text{ }\text{ }\mathrm{GeV})=92.47(69)\text{ }\text{ }\mathrm{MeV}$, ${m}_{c}(3\text{ }\text{ }\mathrm{GeV})=\phantom{\rule{0ex}{0ex}}983.7(5.6)\text{ }\text{ }\mathrm{MeV}$, and ${m}_{c}({m}_{c})=1273(10)\text{ }\text{ }\mathrm{MeV}$, with four active flavors; and ${m}_{b}({m}_{b})=4195(14)\text{ }\text{ }\mathrm{MeV}$ with five active flavors. We also obtain ratios of quark masses ${m}_{c}/{m}_{s}=11.783(25)$, ${m}_{b}/{m}_{s}=53.94(12)$, and ${m}_{b}/{m}_{c}=4.578(8)$. The result for ${m}_{c}$ matches the precision of the most precise calculation to date, and the other masses and all quoted ratios are the most precise to date. Moreover, these results are the first with a perturbative accuracy of ${\ensuremath{\alpha}}_{s}^{4}$. As byproducts of our method, we obtain the matrix elements of HQET operators with dimension 4 and 5: ${\overline{\mathrm{\ensuremath{\Lambda}}}}_{\mathrm{MRS}}=555(31)\text{ }\text{ }\mathrm{MeV}$ in the minimal renormalon-subtracted (MRS) scheme, ${\ensuremath{\mu}}_{\ensuremath{\pi}}^{2}=0.05(22)\text{ }\text{ }{\mathrm{GeV}}^{2}$, and ${\ensuremath{\mu}}_{G}^{2}({m}_{b})=0.38(2)\text{ }\text{ }{\mathrm{GeV}}^{2}$. The MRS scheme [Phys. Rev. D 97, 034503 (2018)] is the key new aspect of our method.

Determination of quark masses from nf=4 lattice QCD and the RI-SMOM intermediate scheme

We determine the charm and strange quark masses in the $\overline{\text{MS}}$ scheme, using $n_f=2+1+1$ lattice QCD calculations with highly improved staggered quarks (HISQ) and the RI-SMOM intermediate scheme to connect the bare lattice quark masses to continuum renormalisation schemes. Our study covers analysis of systematic uncertainties from this method, including nonperturbative artefacts and the impact of the non-zero physical sea quark masses. We find $m_c^{\overline{\text{MS}}}(3 \text{GeV}) = 0.9896(61)$ GeV and $m_s^{\overline{\text{MS}}}(3 \text{GeV}) = 0.08536(85)$ GeV, where the uncertainties are dominated by the tuning of the bare lattice quark masses. These results are consistent with, and of similar accuracy to, those using the current-current correlator approach coupled to high-order continuum QCD perturbation theory, implemented in the same quark formalism and on the same gauge field configurations. This provides a strong test of the consistency of methods for determining the quark masses to high precision from lattice QCD. We also give updated lattice QCD world averages for $c$ and $s$ quark masses.

Effects of nonequilibrated topological charge distributions on pseudoscalar meson masses and decay constants

We study the effects of failure to equilibrate the squared topological charge $Q^2$ on lattice calculations of pseudoscalar masses and decay constants. The analysis is based on chiral perturbation theory calculations of the dependence of these quantities on the QCD vacuum angle $\theta$. For the light-light partially quenched case, we rederive the known chiral perturbation theory results of Aoki and Fukaya, but using the nonperturbatively-valid chiral theory worked out by Golterman, Sharpe and Singleton, and by Sharpe and Shoresh. We then extend these calculations to heavy-light mesons. Results when staggered taste-violations are important are also presented. The derived $Q^2$ dependence is compared to that of simulations using the MILC collaboration's ensembles of lattices with four flavors of HISQ dynamical quarks. We find agreement, albeit with large statistical errors. These results can be used to correct for the leading effects of unequilibrated $Q^2$, or to make estimates of the systematic error coming from the failure to equilibrate $Q^2$. In an appendix, we show that the partially quenched chiral theory may be extended beyond a lower bound on valence masses discovered by Sharpe and Shoresh. Subtleties occurring when a sea-quark mass vanishes are discussed in another appendix.

Numerical experiments using deflation with the HISQ action

We report on numerical experiments using deflation to compute quark propagators for the highly improved staggered quark (HISQ) action. The method is tested on HISQ gauge configurations, generated by the MILC collaboration, with lattice spacings of 0.15 fm, with a range of volumes, and sea quark masses down to the physical quark mass.

Charmed states and flavour symmetry breaking

Extending the SU(3) flavour symmetry breaking expansion from up, down and strange sea quark masses to partially quenched valence quark masses allows an extrapolation to the charm quark mass. This approach leads to a determination of charmed quark hadron masses and decay constants. We describe our recent progress and give preliminary results in particular with regard to the recently discovered doubly charmed baryon (the [see formula in PDF]) by the LHCb Collaboration.

Semileptonic B-meson decays to light pseudoscalar mesons on the HISQ ensembles

We report the status of an ongoing lattice-QCD calculation of form factors for exclusive semileptonic decays of B mesons with both charged currents (B → πlv, Bs → Klv) and neutral currents (B → πl+l-, B → Kl+l-). The results are important for constraining or revealing physics beyond the Standard Model. This work uses MILC’s (2+1 + 1)-flavor ensembles with the HISQ action for the sea and light valence quarks and the clover action in the Fermilab interpretation for the b quark. Simulations are carried out at three lattice spacings down to 0.088 fm, with both physical and unphysical sea-quark masses. We present preliminary results for correlation-function fits.

Pseudoscalar meson electromagnetic form factor at high Q2 from full lattice QCD

We give an accurate determination of the vector (electromagnetic) form factor, $F(Q^2)$, for a light pseudoscalar meson up to squared momentum transfer $Q^2$ values of 6 $\mathrm{GeV}^2$ for the first time from full lattice QCD, including $u$, $d$, $s$ and $c$ quarks in the sea at multiple values of the lattice spacing. Our results show good control of lattice discretisation and sea quark mass effects. We study a pseudoscalar meson made of valence $s$ quarks but the qualitative picture obtained applies also to the $\pi$ meson, relevant to upcoming experiments at Jefferson Lab. We find that $Q^2F(Q^2)$ becomes flat in the region between $Q^2$ of 2 $\mathrm{GeV}^2$ and 6 $\mathrm{GeV}^2$, with a value well above that of the asymptotic perturbative QCD expectation, but well below that of the vector-meson dominance pole form appropriate to low $Q^2$ values. Our calculations show that we can reach higher $Q^2$ values in future to shed further light on where the perturbative QCD result emerges.

Quark chiral condensate from the overlap quark propagator**Supported by National Natural Science Foundation of China (11575197, 11575196, 11335001, 11405178), joint funds of NSFC (U1632104, U1232109), YC and ZL acknowledge the support of NSFC and DFG (CRC110)

From the overlap lattice quark propagator calculated in the Landau gauge, we determine the quark chiral condensate by fitting operator product expansion formulas to the lattice data. The quark propagators are computed on domain wall fermion configurations generated by the RBC-UKQCD Collaborations with Nf = 2+1 flavors. Three ensembles with different light sea quark masses are used at one lattice spacing 1/a = 1.75(4) GeV. We obtain in the SU(2) chiral limit.

Strange Quark Magnetic Moment of the Nucleon at the Physical Point.

We report a lattice QCD calculation of the strange quark contribution to the nucleon's magnetic moment and charge radius. This analysis presents the first direct determination of strange electromagnetic form factors including at the physical pion mass. We perform a model-independent extraction of the strange magnetic moment and the strange charge radius from the electromagnetic form factors in the momentum transfer range of 0.051 GeV^{2}≲Q^{2}≲1.31 GeV^{2}. The finite lattice spacing and finite volume corrections are included in a global fit with 24 valence quark masses on four lattices with different lattice spacings, different volumes, and four sea quark masses including one at the physical pion mass. We obtain the strange magnetic moment G_{M}^{s}(0)=-0.064(14)(09)μ_{N}. The four-sigma precision in statistics is achieved partly due to low-mode averaging of the quark loop and low-mode substitution to improve the statistics of the nucleon propagator. We also obtain the strange charge radius ⟨r_{s}^{2}⟩_{E}=-0.0043(16)(14) fm^{2}.

Stochastic calculation of the Dirac spectrum on the lattice and a determination of chiral condensate in 2 + 1-flavor QCD

We compute the chiral condensate in 2+1-flavor QCD through the spectrum of low-lying eigenmodes of Dirac operator. The number of eigenvalues of the Dirac operator is evaluated using a stochastic method with an eigenvalue filtering technique on the background gauge configurations generated by lattice QCD simulations including the effects of dynamical up, down and strange quarks described by the Mobius domain-wall fermion formulation. The low-lying spectrum is related to the chiral condensate, which is one of the leading order low-energy constants in chiral effective theory, as dictated by the Banks-Casher relation. The spectrum shape and its dependence on the sea quark masses calculated in numerical simulations are consistent with the expectation from one-loop chiral perturbation theory. After taking the chiral limit as well as the continuum limit using the data at three lattice spacings ranging 0.080-0.045 fm, we obtain $\Sigma^{1/3}$(2 GeV) = 270.0(4.9) MeV, with the error combining those from statistical and from various sources of systematic errors. Finite volume effect is confirmed to be under control by a direct comparison of the results from two different volumes at the lightest available sea quarks corresponding to 230 MeV pions.

Diquark mass differences from unquenched lattice QCD**Supported by National Science Foundation of China (11575197, 10835002, 11405178, 11335001), joint funds of NSFC (U1232109), MG and ZL are partially supported by the Youth Innovation Promotion Association of CAS (2015013, 2011013), YC and ZL acknowledge support of NSFC and DFG (CRC110)

We calculate diquark correlation functions in the Landau gauge on the lattice using overlap valence quarks and 2+1-flavor domain wall fermion configurations. Quark masses are extracted from the scalar part of quark propagators in the Landau gauge. The scalar diquark quark mass difference and axial vector scalar diquark mass difference are obtained for diquarks composed of two light quarks and of a strange and a light quark. The light sea quark mass dependence of the results is examined. Two lattice spacings are used to check the discretization effects. The coarse and fine lattices are of sizes 243 × 64 and 323 × 64 with inverse spacings 1/a = 1.75(4) GeV and 2.33(5) GeV, respectively.

The strange and charm quark contributions to the anomalous magnetic moment of the muon from lattice QCD

We describe a new technique (published in [1]) to determine the contribution to the anomalous magnetic moment of the muon coming from the hadronic vacuum polarisation using lattice QCD. Our method uses Padé approximants to reconstruct the Adler function from its derivatives at q2=0. These are obtained simply and accurately from time-moments of the vector current-current correlator at zero spatial momentum. We test the method using strange quark correlators calculated on MILC Collaboration's nf=2+1+1 HISQ ensembles at multiple values of the lattice spacing, multiple volumes and multiple light sea quark masses (including physical pion mass configurations). We find the (connected) contribution to the anomalous moment from the strange quark vacuum polarisation to be aμs=53.41(59)×10−10, and the contribution from charm quarks to be aμc=14.42(39)×10−10−1% accuracy is achieved for the strange quark contribution. The extension of our method to the light quark contribution and to that from the quark-line disconnected diagram is straightforward.

Heavy flavour precision physics from Nf = 2 + 1 + 1 lattice simulations

We present precision lattice calculations of the pseudoscalar decay constants of the charmed sector as well as determinations of the bottom quark mass and its ratio to the charm quark mass. We employ Nf=2+1+1 dynamical quark gauge configurations generated by the European Twisted Mass Collaboration, using data at three values of the lattice spacing and pion masses as low as 210 MeV. Strange and charm sea quark masses are close to their physical values.

Lattice calculation of the leading strange quark-connected contribution to the muon g − 2

We present results for the leading hadronic contribution to the muon anomalous magnetic moment due to strange quark-connected vacuum polarisation effects. Simulations were performed using RBC--UKQCD's $N_f=2+1$ domain wall fermion ensembles with physical light sea quark masses at two lattice spacings. We consider a large number of analysis scenarios in order to obtain solid estimates for residual systematic effects. Our final result in the continuum limit is $a_\mu^{(2)\,{\rm had},\,s}=53.1(9)\left(^{+1}_{-3}\right)\times10^{-10}$.

Kaon BSMB-parameters using improved staggered fermions fromNf=2+1unquenched QCD

We present results for the matrix elements of the additional $\Delta S=2$ operators that appear in models of physics beyond the Standard Model (BSM), expressed in terms of four BSM $B$-parameters. Combined with experimental results for $\Delta M_K$ and $\epsilon_K$, these constrain the parameters of BSM models. We use improved staggered fermions, with valence HYP-smeared quarks and $N_f=2+1$ flavors of "asqtad" sea quarks. The configurations have been generated by the MILC collaboration. The matching between lattice and continuum four-fermion operators and bilinears is done perturbatively at one-loop order. We use three lattice spacings for the continuum extrapolation: $a\approx 0.09$, $0.06$ and $0.045\;$fm. Valence light-quark masses range down to $\approx m_s^{\rm phys}/13$ while the light sea-quark masses range down to $\approx m_s^{\rm phys}/20$. Compared to our previous published work, we have added four additional lattice ensembles, leading to better controlled extrapolations in the lattice spacing and sea-quark masses. We report final results for two renormalization scales, $\mu=2\;\text{GeV}$ and $3\;\text{GeV}$, and compare them to those obtained by other collaborations. Agreement is found for two of the four BSM $B$-parameters ($B_2$ and $B_3^\text{SUSY}$). The other two ($B_4$ and $B_5$) differ significantly from those obtained using RI-MOM renormalization as an intermediate scheme, but are in agreement with recent preliminary results obtained by the RBC-UKQCD collaboration using RI-SMOM intermediate schemes.