Source author record

C. McNeile

C. McNeile 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

28works
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

28 published item(s)

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.

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.

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

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

High-precision scale setting in lattice QCD

Scale setting is of central importance in lattice QCD. It is required to predict dimensional quantities in physical units. Moreover, it determines the relative lattice spacings of computations performed at different values of the bare coupling, and this is needed for extrapolating results into the continuum. Thus, we calculate a new quantity, $w_0$, for setting the scale in lattice QCD, which is based on the Wilson flow like the scale $t_0$ (M. Luscher, JHEP 1008 (2010) 071). It is cheap and straightforward to implement and compute. In particular, it does not involve the delicate fitting of correlation functions at asymptotic times. It typically can be determined on the few per-mil level. We compute its continuum extrapolated value in 2+1-flavor QCD for physical and non-physical pion and kaon masses, to allow for mass-independent scale setting even away from the physical mass point. We demonstrate its robustness by computing it with two very different actions (one of them with staggered, the other with Wilson fermions) and by showing that the results agree for physical quark masses in the continuum limit.

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).

preprint2012arXiv

Towards the glueball spectrum from unquenched lattice QCD

We use a variational technique to study heavy glueballs on gauge configurations generated with 2+1 flavours of ASQTAD improved staggered fermions. The variational technique includes glueball scattering states. The measurements were made using 2150 configurations at 0.092 fm with a pion mass of 360 MeV. We report masses for 10 glueball states. We discuss the prospects for unquenched lattice QCD calculations of the oddballs.

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

Precision computation of the kaon bag parameter

Indirect CP violation in K \rightarrow ππ decays plays a central role in constraining the flavor structure of the Standard Model (SM) and in the search for new physics. For many years the leading uncertainty in the SM prediction of this phenomenon was the one associated with the nonperturbative strong interaction dynamics in this process. Here we present a fully controlled lattice QCD calculation of these effects, which are described by the neutral kaon mixing parameter B_K . We use a two step HEX smeared clover-improved Wilson action, with four lattice spacings from a\approx0.054 fm to a\approx0.093 fm and pion masses at and even below the physical value. Nonperturbative renormalization is performed in the RI-MOM scheme, where we find that operator mixing induced by chiral symmetry breaking is very small. Using fully nonperturbative continuum running, we obtain our main result B_K^{RI}(3.5GeV)=0.531(6)_{stat}(2)_{sys}. A perturbative 2-loop conversion yields B_K^{MSbar-NDR}(2GeV)=0.564(6)_{stat}(3)_{sys}(6)_{PT}, which is in good agreement with current results from fits to experimental data.

preprint2011arXiv

Spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states for $N_f=2+ 1$

We study the spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states. We use two sets of $N_f = 2 + 1$ gauge configurations, one set from QCDSF using the SLiNC action, and the other configurations from the Budapest-Marseille-Wuppertal collaboration, using the HEX smeared clover action. To extract information about the excited states, we choose a suitable basis of operators to implement the variational method.

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

Precise determination of the lattice spacing in full lattice QCD

We compare three different methods to determine the lattice spacing in lattice QCD and give results from calculations on the MILC ensembles of configurations that include the effect of $u$, $d$ and $s$ sea quarks. It is useful, for ensemble to ensemble comparison, to express the results as giving a physical value for $r_1$, a parameter from the heavy quark potential. Combining the three methods gives a value for $r_1$ in the continuum limit of 0.3133(23)(3) fm. Using the MILC values for $r_0/r_1$, this corresponds to a value for the $r_0$ parameter of 0.4661(38) fm. We also discuss how to use the $η_s$ for determining the lattice spacing and tuning the $s$-quark mass accurately, by giving values for $m_{η_s}$ (0.6858(40) GeV) and $f_{η_s}$ (0.1815(10) GeV).

preprint2009arXiv

The omega-rho meson mass splitting and mixing from lattice QCD

We compare flavour singlet and non-singlet vector mesons from first principles using lattice QCD. With N_f=2 flavours of light quark, this addresses the omega-rho mass difference. Using maximally twisted-mass lattice QCD, we are able for the first time to determine this mass difference precisely and we compare with experiment. We also discuss omega-rho mixing effects arising within QCD through the u-d quark mass difference.

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.

preprint2007arXiv

Dynamical Twisted Mass Fermions with Light Quarks

We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

preprint2006arXiv

Decay width of light quark hybrid meson from the lattice

Lattice QCD with $N_f=2$ flavours of sea quark is used to explore the spectrum and decay of a $J^{PC}=1^{-+}$ spin-exotic hybrid meson. We test lattice determination of S-wave decay amplitudes at threshold using $b_1 \to πω$ where agreement with data is found. We find a hybrid meson state at 2.2(2) GeV with a partial width to $πb_1$ of 400(120) MeV and to $πf_1$ of 90(60) MeV.

preprint2006arXiv

Properties of light scalar mesons from lattice QCD

Lattice QCD with $N_f=2$ flavours of sea quark is used to explore the spectrum and decay of scalar mesons. We are able to determine the $b_1$ - $a_0$ mass difference and this leads to the conclusion that the lightest non-singlet scalar meson ($a_0$) has a mass of 1.01(4) GeV. We determine from the lattice the coupling strength to KK and $πη$. We compute the leptonic decay constant of the lightest non-singlet scalar meson. We discuss the impact of these lattice results on the interpretation of the $a_0(980)$ state. We also discuss $K^*_0$ states.

preprint2002arXiv

Hybrid meson decay from the lattice

We discuss the allowed decays of a hybrid meson in the heavy quark limit. We deduce that an important decay will be into a heavy quark non-hybrid state and a light quark meson, in other words, the de-excitation of an excited gluonic string by emission of a light quark-antiquark pair. We discuss the study of hadronic decays from the lattice in the heavy quark limit and apply this approach to explore the transitions from a spin-exotic hybrid to $χ_b η$ and $χ_b S$ where $S$ is a scalar meson. We obtain a signal for the transition emitting a scalar meson and we discuss the phenomenological implications.

preprint2001arXiv

Is the up-quark massless?

We report on determinations of the low-energy constants alpha5 and alpha8 in the effective chiral Lagrangian at O(p^4), using lattice simulations with N_f=2 flavours of dynamical quarks. Precise knowledge of these constants is required to test the hypothesis whether or not the up-quark is massless. Our results are obtained by studying the quark mass dependence of suitably defined ratios of pseudoscalar meson masses and matrix elements. Although comparisons with an earlier study in the quenched approximation reveal small qualitative differences in the quark mass behaviour, numerical estimates for alpha5 and alpha8 show only a weak dependence on the number of dynamical quark flavours. Our results disfavour the possibility of a massless up-quark, provided that the quark mass dependence in the physical three-flavour case is not fundamentally different from the two-flavour case studied here.

preprint1992arXiv

Quenched Hadrons using Wilson and O(a)-Improved Fermion Actions at beta=6.2

We present the first study of the light hadron spectrum and decay constants for quenched QCD using an O(a)-improved nearest-neighbour Wilson fermion action at β=6.2. We compare the results with those obtained using the standard Wilson fermion action, on the same set of 18 gauge field configurations of a 24^3 times 48 lattice. For pseudoscalar meson masses in the range 330-800 MeV, we find no significant difference between the results for the two actions. The scales obtained from the string tension and mesonic sector are consistent, but differ from that derived from baryon masses. The ratio of the pseudoscalar decay constant to the vector meson mass is roughly independent of quark mass as observed experimentally, and in approximate agreement with the measured value.