Source author record

G. P. Lepage

G. P. Lepage 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

45works
3topics
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

45 published item(s)

preprint2022arXiv

Properties of low-lying charmonia and bottomonia from lattice QCD + QED

The precision of lattice QCD calculations has been steadily improving for some time and is now approaching, or has surpassed, the 1% level for multiple quantities. At this level QED effects, i.e. the fact that quarks carry electric as well as color charge, come into play. In this report we will summarise results from the first lattice QCD+QED computations of the properties of ground-state charmonium and bottomonium mesons by the HPQCD Collaboration.

preprint2020arXiv

Hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment from four-flavor lattice QCD

We calculate the contribution to the muon anomalous magnetic moment hadronic vacuum polarization from {the} connected diagrams of up and down quarks, omitting electromagnetism. We employ QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks and the physical pion mass, and analyze five ensembles with lattice spacings ranging from $a \approx 0.06$ to~0.15~fm. The up- and down-quark masses in our simulations have equal masses $m_l$. We obtain, in this world where all pions have the mass of the $π^0$, $10^{10} a_μ^{ll}({\rm conn.}) = 637.8\,(8.8)$, in agreement with independent lattice-QCD calculations. We then combine this value with published lattice-QCD results for the connected contributions from strange, charm, and bottom quarks, and an estimate of the uncertainty due to the fact that our calculation does not include strong-isospin breaking, electromagnetism, or contributions from quark-disconnected diagrams. Our final result for the total $\mathcal{O}(α^2)$ hadronic vacuum polarization to the muon's anomalous magnetic moment is~$10^{10}a_μ^{\rm HVP,LO} = 699(15)_{u,d}(1)_{s,c,b}$, where the errors are from the light-quark and heavy-quark contributions, respectively. Our result agrees with both {\it ab-initio} lattice-QCD calculations and phenomenological determinations from experimental $e^+e^-$-scattering data. It is $1.3σ$ below the "no new physics" value of the hadronic-vacuum-polarization contribution inferred from combining the BNL E821 measurement of $a_μ$ with theoretical calculations of the other contributions.

preprint2015arXiv

$B \rightarrow π\ell ν$ at zero recoil from lattice QCD with physical $u/d$ quarks

The exclusive semileptonic decay $B \rightarrow π\ell ν$ is a key process for the determination of the Cabibbo-Kobayashi-Maskawa matrix element $V_{ub}$ from the comparison of experimental rates as a function of $q^2$ with theoretically determined form factors. The sensitivity of the form factors to the $u/d$ quark mass has meant significant systematic uncertainties in lattice QCD calculations at unphysically heavy pion masses. Here we give the first lattice QCD calculations of this process for u/d quark masses going down to their physical values, calculating the $f_0$ form factor at zero recoil to 3\%. We are able to resolve a long-standing controversy by showing that the soft-pion theorem result $f_0(q^2_{max}) = f_B/f_π$ does hold as $m_π \rightarrow 0$. We use the Highly Improved Staggered Quark formalism for the light quarks and show that staggered chiral perturbation theory for the $m_π$ dependence is almost identical to continuum chiral perturbation theory for $f_0$, $f_B$ and $f_π$. We also give results for other processes such as $B_s \rightarrow K \ell ν$.

preprint2015arXiv

An estimate of the hadronic vacuum polarization disconnected contribution to the anomalous magnetic moment of the muon from lattice QCD

The quark-line disconnected diagram is a potentially important ingredient in lattice QCD calculations of the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. It is also a notoriously difficult one to evaluate. Here, for the first time, we give an estimate of this contribution based on lattice QCD results that have a statistically significant signal, albeit at one value of the lattice spacing and an unphysically heavy value of the $u/d$ quark mass. We use HPQCD's method of determining the anomalous magnetic moment by reconstructing the Adler function from time-moments of the current-current correlator at zero spatial momentum. Our results lead to a total (including $u$, $d$ and $s$ quarks) quark-line disconnected contribution to $a_μ$ of $-0.15\%$ of the $u/d$ hadronic vacuum polarization contribution with an uncertainty which is 1\% of that contribution.

preprint2015arXiv

B-meson decay constants: a more complete picture from full lattice QCD

We extend the picture of $B$-meson decay constants obtained in lattice QCD beyond those of the $B$, $B_s$ and $B_c$ to give the first full lattice QCD results for the $B^*$, $B^*_s$ and $B^*_c$. We use improved NonRelativistic QCD for the valence $b$ quark and the Highly Improved Staggered Quark (HISQ) action for the lighter quarks on gluon field configurations that include the effect of $u/d$, $s$ and $c$ quarks in the sea with $u/d$ quark masses going down to physical values. For the ratio of vector to pseudoscalar decay constants, we find $f_{B^*}/f_B$ = 0.941(26), $f_{B^*_s}/f_{B_s}$ = 0.953(23) (both $2σ$ less than 1.0) and $f_{B^*_c}/f_{B_c}$ = 0.988(27). Taking correlated uncertainties into account we see clear indications that the ratio increases as the mass of the lighter quark increases. We compare our results to those using the HISQ formalism for all quarks and find good agreement both on decay constant values when the heaviest quark is a $b$ and on the dependence on the mass of the heaviest quark in the region of the $b$. Finally, we give an overview plot of decay constants for gold-plated mesons, the most complete picture of these hadronic parameters to date.

preprint2015arXiv

The size of the pion from full lattice QCD with physical $u$, $d$, $s$ and $c$ quarks

We present the first calculation of the electromagnetic form factor of the $π$ meson at physical light quark masses. We use configurations generated by the MILC collaboration including the effect of $u$, $d$, $s$ and $c$ sea quarks with the Highly Improved Staggered Quark formalism. We work at three values of the lattice spacing on large volumes and with $u$/$d$ quark masses going down to the physical value. We study scalar and vector form factors for a range in space-like $q^2$ from 0.0 to -0.1 $\mathrm{GeV}^2$ and from their shape we extract mean square radii. Our vector form factor agrees well with experiment and we find $\langle r^2 \rangle_V = 0.403(18)(6) \,\mathrm{fm}^2$. For the scalar form factor we include quark-line disconnected contributions which have a significant impact on the radius. We give the first results for SU(3) flavour-singlet and octet scalar mean square radii, obtaining: $\langle r^2 \rangle_S^{\mathrm{singlet}} = 0.506(38)(53) \mathrm{fm}^2$ and $\langle r^2 \rangle_S^{\mathrm{octet}} = 0.431(38)(46) \mathrm{fm}^2$. We discuss the comparison with expectations from chiral perturbation theory.

preprint2014arXiv

B-meson mixing from full lattice QCD with physical u, d, s and c quarks

We present the first lattice QCD calculation of the $B_s$ and $B_d$ mixing parameters with physical light quark masses. We use MILC gluon field configurations that include $u$, $d$, $s$ and $c$ sea quarks at 3 values of the lattice spacing and with 3 values of the $u/d$ quark mass going down to the physical value. We use improved NRQCD for the valence $b$ quarks. Preliminary results show significant improvements over earlier values.

preprint2014arXiv

High-precision quark masses and QCD coupling from $n_f=4$ lattice QCD

We present a new lattice QCD analysis of heavy-quark pseudoscalar-pseudoscalar correlators, using gluon configurations from the MILC collaboration that include vacuum polarization from $u$, $d$, $s$ and $c$~quarks($n_f=4$). We extract new values for the QCD coupling and for the $c$ quark's $\overline{\mathrm{MS}}$ mass: $α_{\overline{\mathrm{MS}}}(M_Z,n_f=5) = 0.11822(74)$ and $m_c(3\mathrm{GeV}, n_f=4) = 0.9851(63)$GeV. These agree well with our earlier simulations using $n_f=3$ sea quarks, vindicating the perturbative treatment of $c$ quarks in that analysis. We also obtain a new nonperturbative result for the ratio of $c$~and $s$~quark masses: $m_c/m_s=11.652(65)$. This ratio implies $m_s(2\,\mathrm{GeV}, n_f=3)=93.6(8)$MeV when it is combined with our new~$c$~mass. Combining $m_c/m_s$ with our earlier $m_b/m_c$ gives $m_b/m_s=52.55(55)$, which is several standard deviations (but only 4%) away from the Georgi-Jarlskop prediction from certain GUTs. Finally we obtain an $n_f=4$ estimate for $m_b/m_c=4.528(54)$ which agrees well with our earlier $n_f=3$ result. The new ratio implies~$m_b(m_b,n_f=5)=4.162(48)$GeV.

preprint2014arXiv

Radial and orbital excitation energies of charmonium

The charmonium system has several excited states below the energy threshold for decay into $D$ and $\bar{D}$ mesons, which can in principle be studied accurately in lattice QCD. Studies that include many states in the spectrum have typically only been done at one value of the lattice spacing and with relatively heavy light quarks in the sea. Here we give preliminary results for radial and orbital excitation energies for charmonium from a calculation on 2+1+1 MILC configurations at multiple lattice spacings and including physical values for $u/d$ quark masses. We use the HISQ formulation for $c$ to obtain small discretisation errors and smeared operators to improve excited state overlap.

preprint2014arXiv

Strange and charm quark contributions to the anomalous magnetic moment of the muon

We describe a new technique to determine the contribution to the anomalous magnetic moment of the muon coming from the hadronic vacuum polarization using lattice QCD. Our method reconstructs the Adler function, using Padé approximants, from its derivatives at $q^2=0$ 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 on large-volume gluon field configurations that include the effect of up and down (at physical masses), strange and charm quarks in the sea at multiple values of the lattice spacing and multiple volumes and show that 1% accuracy is achievable. For the charm quark contributions we use our previously determined moments with up, down and strange quarks in the sea on very fine lattices. We find the (connected) contribution to the anomalous moment from the strange quark vacuum polarization to be $a_μ^s = 53.41(59) \times 10^{-10}$, and from charm to be $a_μ^c = 14.42(39)\times 10^{-10}$. These are in good agreement with flavour-separated results from non-lattice methods, given caveats about the comparison. The extension of our method to the light quark contribution and to that from the quark-line disconnected diagram is straightforward.

preprint2014arXiv

The $Υ$ and $Υ^{\prime}$ Leptonic Widths, $a_μ^b$ and $m_b$ from full lattice QCD

We determine the decay rate to leptons of the ground-state $Υ$ meson and its first radial excitation in lattice QCD for the first time. We use radiatively-improved NRQCD for the $b$ quarks and include $u$, $d$, $s$ and $c$ quarks in the sea with $u/d$ masses down to their physical values. We find $Γ(Υ\rightarrow e^+e^-)$ = 1.19(11) keV and $Γ(Υ^{\prime} \rightarrow e^+e^-)$ = 0.69(9) keV, both in good agreement with experiment. The decay constants we obtain are included in a summary plot of meson decay constants from lattice QCD given in the Conclusions. We also test time-moments of the vector current-current correlator against values determined from the $b$ quark contribution to $σ(e^+e^- \rightarrow \mathrm{hadrons})$ and calculate the $b$-quark piece of the hadronic vacuum polarisation contribution to the anomalous magnetic moment of the muon, $a_μ^b = 0.271(37) \times 10^{-10}$. Finally we determine the $b$-quark mass, obtaining in the $\overline{MS}$ scheme, $\overline{m}_b(\overline{m}_b, n_f=5)$ = 4.196(23) GeV, the most accurate result from lattice QCD to date.

preprint2013arXiv

B, Bs, K and pi weak matrix elements with physical light quarks

Calculations of pseudoscalar decay constants of B, Bs, K and pi mesons with physical light quarks are presented. We use HISQ ensembles that include u,d,s and c sea quarks at three lattice spacings. HISQ is used for the valence light quarks and a radiatively improved NRQCD action for the heavy quarks. The key results are f_{B^+}=0.184(4)$ GeV, f_{B_s}=0.224(4) GeV, f_{B_s}/f_{B^+}=1.217(8), f_{K^+}/f_{pi^+}=1.1916(21), f_{K^+}=155.37(34) MeV, giving a significant improvement over previous results that required chiral extrapolation. We also calculate the Wilson flow scale w_0, finding w_0=0.1715(9) fm.

preprint2013arXiv

Bottomonium and B results from full lattice QCD

We have developed two methods for handling $b$ quarks in lattice QCD. One uses NRQCD (now improved to include radiative corrections) and the other uses Highly Improved Staggered Quarks (HISQ), extrapolating to the $b$ quark from lighter masses and using multiple lattice spacings to control discretisation errors. Comparison of results for the two different methods gives confidence in estimates of lattice QCD systematic errors, since they are very different in these two cases. Here we show results for heavyonium hyperfine splittings and vector current-current correlator moments using HISQ quarks, to add to earlier results testing the heavy HISQ method with pseudoscalar mesons. We also show the form factor for $B \rightarrow πl ν$ decay at zero recoil using NRQCD $b$ quarks and $u/d$ quarks with physical masses. This allows us to test the soft pion theorem relation ($f_0(q^2_{max})=f_B/f_π$) accurately and we find good agreement as $M_π \rightarrow 0$. }

preprint2013arXiv

Direct determination of strange and light quark condensates from full lattice QCD

We determine the strange and light quark condensates in full lattice QCD for the first time. This is done by direct calculation of the expectation value of the trace of the quark propagator followed by subtraction of the appropriate perturbative contribution to convert to a value for the condensate in the MS-bar scheme at 2 GeV. We use lattice QCD configurations including u, d, s and c quarks in the sea with u/d quark masses going down to the physical value. We find the ratio of the strange to the light quark condensate to be 1.08(16).

preprint2013arXiv

Precision tests of the J/psi from full lattice QCD: mass, leptonic width and radiative decay rate to eta_c

We show results from calculations in full lattice QCD of the mass, leptonic width and radiative decay rate to eta_c of the J/psi meson. These provide few % tests of QCD. Another (1.5%) test comes from comparison of time-moments of the vector charmonium correlator with results derived from the experimental values of R(e+e- to hadrons) in the charm region.

preprint2013arXiv

Prediction of the $D_s^*$ width from a calculation of its radiative decay in full lattice QCD

We determine the rate for $D_s^* \rightarrow D_s γ$ for the first time from lattice QCD and include the full effect of $u$, $d$ and $s$ sea quarks. The valence quarks are implemented using the Highly Improved Staggered Quark (HISQ) formalism and we normalise the vector current nonperturbatively. We obtain $M(D_s^*)-M(D_s)$ of 148(4) MeV, in good agreement with experiment. The value of the decay constant of the $D_s^*$ we find to be 274(6) MeV, so that $f_{D_s^*}/f_{D_s} = 1.10(2)$. For the radiative decay we find $Γ(D_s^* \rightarrow D_s γ)$ = 0.066(26) keV. Given the experimental branching fraction for this decay we predict a total width for the $D_s^*$ of 0.070(28) keV, making this the longest lived charged vector meson.

preprint2013arXiv

The shape of the D -> K semileptonic form factor from full lattice QCD and V_cs

We present a new study of the form factors for D -> K semileptonic decay from lattice QCD that allows us to compare the shape of the vector form factor to experiment and, for the first time, to extract V_cs using results from all experimental q^2 bins. The valence quarks are implemented with the Highly Improved Staggered Quark action on MILC configurations that include u, d and s sea quarks. The scalar and vector currents are nonperturbatively normalised and, using phased boundary conditions, we are able to cover the full q^2 range accessible to experiment. Our result is V_cs = 0.963(5)_{expt}(14)_{lattice}. We also demonstrate that the form factors are insensitive to whether the spectator quark is u/d or s, which has implications for other decay channels.

preprint2013arXiv

V_cs from D_s to ϕlν semileptonic decay and full lattice QCD

We determine the complete set of axial and vector form factors for the Ds to ϕlν decay from full lattice QCD for the first time. The valence quarks are implemented using the Highly Improved Staggered Quark action and we normalise the appropriate axial and vector currents fully nonperturbatively. The q^2 and angular distributions we obtain for the differential rate agree well with those from the BaBar experiment and, from the total branching fraction, we obtain Vcs = 1.017(63), in good agreement with that from D to Klν semileptonic decay. We also find the mass and decay constant of the ϕ meson in good agreement with experiment, showing that its decay to K{\bar{K}} (which we do not include here) has at most a small effect. We include an Appendix on nonperturbative renormalisation of the complete set of staggered vector and axial vector bilinears needed for this calculation.

preprint2013arXiv

Vus from pi and K decay constants in full lattice QCD with physical u, d, s and c quarks

We determine the decay constants of the pi and K mesons on gluon field configurations from the MILC collaboration including u, d, s and c quarks. We use three values of the lattice spacing and u/d quark masses going down to the physical value. We use the w_0 parameter to fix the relative lattice spacing and f_pi to fix the overall scale. This allows us to obtain a value for f{K^+}/f{pi^+} = 1.1916(21). Comparing to the ratio of experimental leptonic decay rates gives |Vus| = 0.22564(28){Br(K^+)}(20){EM}(40){latt}(5){Vud} and the test of unitarity of the first row of the Cabibbo-Kobayashi-Maskawa matrix: |Vud|^2+|Vus|^2+|Vub|^2 - 1 = 0.00009(51).

preprint2012arXiv

Direct determination of the strange and light quark condensates from full lattice QCD

We determine the strange quark condensate from lattice QCD for the first time and compare its value to that of the light quark and chiral condensates. The results come from a direct calculation of the expectation value of the trace of the quark propagator followed by subtraction of the appropriate perturbative contribution, derived here, to convert the non-normal-ordered $m\barψψ$ to the $\bar{MS}$ scheme at a fixed scale. This is then a well-defined physical `nonperturbative' condensate that can be used in the Operator Product Expansion of current-current correlators. The perturbative subtraction is calculated through $\mathcal{O}(α_s)$ and estimates of higher order terms are included through fitting results at multiple lattice spacing values. The gluon field configurations used are `second generation' ensembles from the MILC collaboration that include 2+1+1 flavors of sea quarks implemented with the Highly Improved Staggered Quark action and including $u/d$ sea quarks down to physical masses. Our results are : $<\bar{s}{s}>^{\bar{MS}}(2 \mathrm{GeV})= -(290(15) \mathrm{MeV})^3$, $<\bar{l}{l}>^{\bar{MS}}(2\, \mathrm{GeV})= -(283(2) \mathrm{MeV})^3$, where $l$ is a light quark with mass equal to the average of the $u$ and $d$ quarks. The strange to light quark condensate ratio is 1.08(16). The light quark condensate is significantly larger than the chiral condensate in line with expectations from chiral analyses. We discuss the implications of these results for other calculations.

preprint2012arXiv

Heavy meson masses and decay constants from relativistic heavy quarks in full lattice QCD

We determine masses and decay constants of heavy-heavy and heavy-charm pseudoscalar mesons as a function of heavy quark mass using a fully relativistic formalism known as Highly Improved Staggered Quarks for the heavy quark. We are able to cover the region from the charm quark mass to the bottom quark mass using MILC ensembles with lattice spacing values from 0.15 fm down to 0.044 fm. We obtain f_{B_c} = 0.427(6) GeV; m_{B_c} = 6.285(10) GeV and f_{η_b} = 0.667(6) GeV. Our value for f_{η_b} is within a few percent of f_Υ confirming that spin effects are surprisingly small for heavyonium decay constants. Our value for f_{B_c} is significantly lower than potential model values being used to estimate production rates at the LHC. We discuss the changing physical heavy-quark mass dependence of decay constants from heavy-heavy through heavy-charm to heavy-strange mesons. A comparison between the three different systems confirms that the B_c system behaves in some ways more like a heavy-light system than a heavy-heavy one. Finally we summarise current results on decay constants of gold-plated mesons.

preprint2012arXiv

Precision tests of the J/ψ from full lattice QCD: mass, leptonic width and radiative decay rate to η_c

We calculate the J/ψ mass, leptonic width and radiative decay rate to γη_c from lattice QCD including u, d and s quarks in the sea for the first time. We use the Highly Improved Staggered Quark formalism and nonperturbatively normalised vector currents for the leptonic and radiative decay rates. Our results are: M_{J/ψ} -M_{η_c} = 116.5(3.2) MeV; Γ(J/ψ to e^+e^-) = 5.48(16) keV; Γ(J/ψ to γη_c) = 2.49(19) keV. The first two are in good agreement with experiment, with Γ(J/ψ to e^+e^-) providing a test of a decay matrix element in QCD, independent of CKM uncertainties, to 2%. At the same time results for the time moments of the correlation function can be compared to values from the charm contribution to Γ(e^+e^- to hadrons), giving a 1.5% test of QCD. Our results show that an improved experimental error would enable a similarly strong test from Γ(J/ψ to γη_c).

preprint2011arXiv

Fast Fits for Lattice QCD Correlators

We illustrate a technique for fitting lattice QCD correlators to sums of exponentials that is significantly faster than traditional fitting methods --- 10--40 times faster for the realistic examples we present. Our examples are drawn from a recent analysis of the Upsilon spectrum, and another recent analysis of the D -> pi semileptonic form factor. For single correlators, we show how to simplify traditional effective-mass analyses.

preprint2011arXiv

High-Precision $f_{B_s}$ and HQET from Relativistic Lattice QCD

We present a new determination of the $B_s$ leptonic decay constant from lattice QCD simulations that use gluon configurations from MILC and a highly improved discretization of the relativistic quark action for both valence quarks. Our result, $f_{B_s} = 0.225(4)$\,GeV, is almost three times more accurate than previous determinations. We analyze the dependence of the decay constant on the heavy quark's mass and obtain the first empirical evidence for the leading $1/\sqrt{m_h}$ dependence predicted by Heavy Quark Effective Theory (HQET). As a check, we use our analysis technique to calculate the $m_{B_s}-m_{η_b}/2$ mass difference. Our result agrees with experiment to within errors of $11\,\mathrm{MeV}$ (better than 2%). We discuss how to extend our analysis to other quantities in $B_s$ and $B$ physics, making 2%-precision possible for the first time.

preprint2011arXiv

The Upsilon spectrum and the determination of the lattice spacing from lattice QCD including charm quarks in the sea

We give results for the Upsilon spectrum from lattice QCD using an improved version of the NRQCD action for b quarks which includes radiative corrections to kinetic terms at O(v^4) in the velocity expansion. We also include for the first time the effect of up, down, strange and charm quarks in the sea using 'second generation' gluon field configurations from the MILC collaboration. Using the Υ 2S - 1S splitting to determine the lattice spacing, we are able to obtain the 1P - 1S splitting to 1.4% and the 3S - 1S splitting to 2.4%. Our improved result for M (Upsilon) - M (eta_b) is 70(9) MeV and we predict M(Upsilon) - M(eta_b') = 35(3) MeV. We also calculate pi, K and eta_s correlators using the Highly Improved Staggered Quark action and perform a chiral and continuum extrapolation to give values for M(eta_s) (0.6893(12) GeV) and f_{eta_s} (0.1819(5) GeV) that allow us to tune the strange quark mass as well as providing an independent and consistent determination of the lattice spacing. Combining the NRQCD and HISQ analyses gives mb/ms = 54.7(2.5) and a value for the heavy quark potential parameter of r_1 = 0.3209(26) fm.

preprint2010arXiv

High-Precision c and b Masses, and QCD Coupling from Current-Current Correlators in Lattice and Continuum QCD

We extend our earlier lattice-QCD analysis of heavy-quark correlators to smaller lattice spacings and larger masses to obtain new values for the c mass and QCD coupling, and, for the first time, values for the b mass: m_c(3GeV,n_f=4)=0.986(6)GeV, alpha_msb(M_Z,n_f=5)=0.1183(7), and m_b(10GeV,n_f=5)=3.617(25)GeV. These are among the most accurate determinations by any method. We check our results using a nonperturbative determination of the mass ratio m_b(mu,n_f)/m_c(mu,n_f); the two methods agree to within our 1% errors and taken together imply m_b/m_c=4.51(4). We also update our previous analysis of alpha_msb from Wilson loops to account for revised values for r_1 and r_1/a, finding a new value alpha_\msb(M_Z,n_f=5)=0.1184(6); and we update our recent values for light-quark masses from the ratio m_c/m_s. Finally, in the Appendix, we derive a procedure for simplifying and accelerating complicated least-squares fits.

preprint2010arXiv

Precise charm to strange mass ratio and light quark masses from full lattice QCD

By using a single formalism to handle charm, strange and light valence quarks in full lattice QCD for the first time, we are able to determine ratios of quark masses to 1%. For $m_c/m_s$ we obtain 11.85(16), an order of magnitude more precise than the current PDG average. Combined with 1% determinations of the charm quark mass now possible this gives $\bar{m}_s(2{\rm GeV}) =$ 92.4(1.5) MeV. The MILC result for $m_s/m_l = 27.2(3) $ yields $\bar{m}_l(2{\rm GeV})$ = 3.40(7) MeV for the average of $u$ and $d$ quark masses.

preprint2010arXiv

Update: Precision D_s decay constant from full lattice QCD using very fine lattices

We update our previous determination of both the decay constant and the mass of the $D_s$ meson using the Highly Improved Staggered Quark formalism. We include additional results at two finer values of the lattice spacing along with improved determinations of the lattice spacing and improved tuning of the charm and strange quark masses. We obtain $m_{D_s}$ = 1.9691(32) GeV, in good agreement with experiment, and $f_{D_s}$ = 0.2480(25) GeV. Our result for $f_{D_s}$ is 1.6$σ$ lower than the most recent experimental average determined from the $D_s$ leptonic decay rate and using $V_{cs}$ from CKM unitarity. Combining our $f_{D_s}$ with the experimental rate we obtain a direct determination of $V_{cs} = 1.010(22)$, or alternatively $0.990 {+0.013 \atop -0.016}$ using a probability distribution for statistical errors for this quantity which vanishes above 1. We also include an accurate prediction of the decay constant of the $η_c$, $f_{η_c}$ = 0.3947(24) GeV, as a calibration point for other lattice calculations.

preprint2009arXiv

A Prediction of the B*_c mass in full lattice QCD

By using the Highly Improved Staggered Quark formalism to handle charm, strange and light valence quarks in full lattice QCD, and NRQCD to handle bottom valence quarks we are able to determine accurately ratios of the B meson vector-pseudoscalar mass splittings, in particular, (m(B*_c)-m(B_c))/(m(B*_s)-m(B_s)). We find this ratio to be 1.15(15), showing the `light' quark mass dependence of this splitting to be very small. Hence we predict m(B_c*) = 6.330(7)(2)(6) GeV where the first two errors are from the lattice calculation and the third from existing experiment. This is the most accurate prediction of a gold-plated hadron mass from lattice QCD to date.

preprint2009arXiv

Recent results on B mixing and decay constants from HPQCD

We review recent results for B_d and B_s mixing parameters using MILC N_f=2+1 lattices, NRQCD b-quarks and AsqTad light quarks. Latest numbers for decay constants f_B and f_{B_s} are also presented. Combining our lattice results with experimental determinations of the mass differences ΔM_d and ΔM_s leads to an important ratio of elements of the CKM matrix, |V_{td}|/|V_{ts}| = 0.214(1)(5) and an updated Standard Model number for the branching fraction Br(B_s -> mu^+ mu^-) = 3.19(19) x 10^{-9}. Preliminary new results for f_{B_s} based on other actions are also described.

preprint2008arXiv

Precision charm physics, m_c and alpha_s from lattice QCD

We present an update of results from the HPQCD collaboration on charm physics using the Highly Improved Staggered Quark action. This includes a precise determination of m_c using moments of current-current correlators combined with high-order continuum QCD perturbation theory. We also include an update on the determination of alpha_s from lattice QCD, preliminary results on the determination of m_b and a summary plot of the status of the gold-plated meson spectrum. There is an appendix on tackling systematic errors in fitting using the Bayesian approach.

preprint2002arXiv

Semileptonic B Decays from an NRQCD/D234 Action

Semileptonic, B --> pi, l nubar, decays are studied on quenched anisotropic lattices using tree-level tadpole improved Symanzik glue, NRQCD heavy quark and D234 light quark actions. Constrained fitting methods are applied to extract groundstate contributions to two-point and three-point correlators. We agree with previous lattice determinations of the form factors. The major source of systematic error here, as in previous work, comes from the chiral extrapolation to the physical pion mass. Future calculations must work at lighter quark masses to resolve this.

preprint2002arXiv

The Upsilon Spectrum from Lattice QCD with 2+1 Flavors of Dynamical Quarks

We describe the bottomonium spectrum obtained on the MILC configurations which incorporate 2+1 flavors of dynamical quarks. We compare to quenched and 2 flavor results also on MILC configurations. We show that the lattice spacing determination using different quantities shows clear signs of convergence with 2+1 flavors and give results for the leptonic width and hyperfine splitting, in the form of the ratio of the 1st excited state of the Upsilon to that of the ground state.

preprint2000arXiv

Scaling and Further Tests of Heavy Meson Decay Constant Determinations from Nonrelativistic QCD

We present results for the B_s meson decay constant f_{B_s} from simulations at three lattice spacings in the range a^{-1}=1.1 to 2.6 GeV using NRQCD heavy quarks and clover light quarks in the quenched approximation. We study scaling of this quantity and check the consistency between mesons decaying from rest and from a state with nonzero spatial momentum. The cancellation of power law contributions that arise in the NRQCD formulation of heavy-light currents is discussed. On the coarsest lattice the D_s meson decay constant f_{D_s} is calculated. Our best values for the decay constants are given by f_{B_s} = 187(4)(4)(11)(2)(7)(6) MeV and f_{D_s} = 223(6)(31)(38)(23)(9)(^{+3}_{-1}) MeV.

preprint1999arXiv

Perturbative two- and three-loop coefficients from large beta Monte Carlo

Perturbative coefficients for Wilson loops and the static quark self-energy are extracted from Monte Carlo simulations at large beta on finite volumes, where all the lattice momenta are large. The Monte Carlo results are in excellent agreement with perturbation theory through second order. New results for third order coefficients are reported. Twisted boundary conditions are used to eliminate zero modes and to suppress Z_3 tunneling.

preprint1998arXiv

Bottom and Charm quark masses from lattice NRQCD

We present new values for the $\bar{MS}$ masses of $b$ and $c$ quarks based on lattice NRQCD simulations of the $Υ(b\bar{b})$ and $ψ(c\bar{c})$ systems. These include three measurements of the $b$ mass based on quenched simulations with lattice spacings ranging from 0.05fm to 0.15fm, which we find to be largely independent of lattice spacing. In addition, we find a consistent value from an unquenched simulation at 0.08fm.

preprint1997arXiv

Further Precise Determinations of $α_s$ from Lattice QCD

We present a new determination of the strong coupling constant from lattice QCD simulations. We use four different short-distance quantities to obtain the coupling, three different (infrared) meson splittings to tune the simulation parameters, and a wide range of lattice spacings, quark masses, and lattice volumes to test for systematic errors. Our final result consists of ten different determinations of $α^{(3)}_{P}(8.2 GeV)$, which agree well with each other and with our previous results. The most accurate of these, when evolved perturbatively to the $Z^0$ mass, gives $α^{(5)}_{\msbar}(M_Z) = .1174(24)$. We compare our results with those obtained from other recent lattice simulations.

preprint1996arXiv

Towards Highly Improved Quark Actions

We describe two ideas useful in the construction of highly improved quark actions for simulations on coarse lattices: (1) Field transformations to solve the doubler problem without destroying tree-level improvement for on- or off-shell quantities. (2) Anisotropic lattices with smaller temporal than spatial lattice spacing. We present first simulation results for the quenched light hadron and charmonium spectra obtained with a D234 action on anisotropic lattices.

preprint1994arXiv

Precision Upsilon Spectroscopy from Nonrelativistic Lattice QCD

The spectrum of the Upsilon system is investigated using the Nonrelativistic Lattice QCD approach to heavy quarks and ignoring light quark vacuum polarization. We find good agreement with experiment for the Upsilon(1S), Upsilon(2S), Upsilon(3S) and for the center of mass and fine structure of the chi_b states. The lattice calculations predict b-bbar D-states with center of mass at (10.20 +/- 0.07 +/- 0.03)GeV. Fitting procedures aimed at extracting both ground and excited state energies are developed. We calculate a nonperturbative dispersion mass for the Upsilon(1S) and compare with tadpole-improved lattice perturbation theory.

preprint1992arXiv

Improved Nonrelativistic QCD for Heavy Quark Physics

We construct an improved version of nonrelativistic QCD for use in lattice simulations of heavy quark physics, with the goal of reducing systematic errors from all sources to below 10\%. We develop power counting rules to assess the importance of the various operators in the action and compute all leading order corrections required by relativity and finite lattice spacing. We discuss radiative corrections to tree level coupling constants, presenting a procedure that effectively resums the largest such corrections to all orders in perturbation theory. Finally, we comment on the size of nonperturbative contributions to the coupling constants.