Source author record

U. M. Heller

U. M. Heller 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

97works
5topics
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

97 published item(s)

preprint2020arXiv

FLAG Review 2019

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to π$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_π$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ as well as those for $D$- and $B$-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $α_s$. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.

preprint2016arXiv

$B\to Kl^+l^-$ decay form factors from three-flavor lattice QCD

We compute the form factors for the $B \to Kl^+l^-$ semileptonic decay process in lattice QCD using gauge-field ensembles with 2+1 flavors of sea quark, generated by the MILC Collaboration. The ensembles span lattice spacings from 0.12 to 0.045 fm and have multiple sea-quark masses to help control the chiral extrapolation. The asqtad improved staggered action is used for the light valence and sea quarks, and the clover action with the Fermilab interpretation is used for the heavy $b$ quark. We present results for the form factors $f_+(q^2)$, $f_0(q^2)$, and $f_T(q^2)$, where $q^2$ is the momentum transfer, together with a comprehensive examination of systematic errors. Lattice QCD determines the form factors for a limited range of $q^2$, and we use the model-independent $z$ expansion to cover the whole kinematically allowed range. We present our final form-factor results as coefficients of the $z$ expansion and the correlations between them, where the errors on the coefficients include statistical and all systematic uncertainties. We use this complete description of the form factors to test QCD predictions of the form factors at high and low $q^2$. We also compare a Standard-Model calculation of the branching ratio for $B \to Kl^+l^-$ with experimental data.

preprint2016arXiv

$B^0_{(s)}$-mixing matrix elements from lattice QCD for the Standard Model and beyond

We calculate---for the first time in three-flavor lattice QCD---the hadronic matrix elements of all five local operators that contribute to neutral $B^0$- and $B_s$-meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral $B$-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio $ξ= 1.206(18)(6)$, where the second error stems from the omission of charm sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from $B$ mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light valence quarks, and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral $B$-meson oscillation frequencies to determine the CKM matrix elements $|V_{td}| = 8.00(34)(8) \times 10^{-3}$, $|V_{ts}| = 39.0(1.2)(0.4) \times 10^{-3}$, and $|V_{td}/V_{ts}| = 0.2052(31)(10)$, which differ from CKM-unitarity expectations by about 2$σ$. These results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.

preprint2016arXiv

Decay constants $f_B$ and $f_{B_s}$ and quark masses $m_b$ and $m_c$ from HISQ simulations

We present a progress report on our calculation of the decay constants $f_B$ and $f_{B_s}$ from lattice-QCD simulations with highly-improved staggered quarks. Simulations are carried out with several heavy valence-quark masses on $(2+1+1)$-flavor ensembles that include charm sea quarks. We include data at six lattice spacings and several light sea-quark masses, including an approximately physical-mass ensemble at all but the smallest lattice spacing, 0.03 fm. This range of parameters provides excellent control of the continuum extrapolation to zero lattice spacing and of heavy-quark discretization errors. Finally, using the heavy-quark effective theory expansion we present a method of extracting from the same correlation functions the charm- and bottom-quark masses as well as some low-energy constants appearing in the heavy-quark expansion.

preprint2016arXiv

Electromagnetic effects on the light pseudoscalar mesons and determination of $m_u/m_d$

The MILC Collaboration has completed production running of electromagnetic effects on light mesons using asqtad improved staggered quarks. In these calculations, we use quenched photons in the noncompact formalism. We study four lattice spacings from $\approx\!0.12\:$fm to $\approx\!0.045\:$fm. To study finite-volume effects, we used six spatial lattice sizes $L/a=12$, 16, 20, 28, 40, and 48, at $a\!\approx\!0.12\:$fm. We update our preliminary values for the correction to Dashen's theorem ($ε$) and the quark-mass ratio $m_u/m_d$.

preprint2016arXiv

Gradient flow and scale setting on MILC HISQ ensembles

We report on a scale determination with gradient-flow techniques on the $N_f=2+1+1$ highly improved staggered quark ensembles generated by the MILC Collaboration. The ensembles include four lattice spacings, ranging from approximately 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ and their tree-level improvements, $\sqrt{t_{0,{\rm imp}}}$ and $w_{0,{\rm imp}}$, are computed on each ensemble using Symanzik flow and the cloverleaf definition of the energy density $E$. Using a combination of continuum chiral-perturbation theory and a Taylor-series ansatz for the lattice-spacing and strong-coupling dependence, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We determine the scales $\sqrt{t_0} = 0.1416({}_{-5}^{+8})$ fm and $w_0 = 0.1714({}_{-12}^{+15})$ fm, where the errors are sums, in quadrature, of statistical and all systematic errors. The precision of $w_0$ and $\sqrt{t_0}$ is comparable to or more precise than the best previous estimates, respectively. We then find the continuum mass dependence of $\sqrt{t_0}$ and $w_0$, which will be useful for estimating the scales of new ensembles. We also estimate the integrated autocorrelation length of $\langle E(t) \rangle$. For long flow times, the autocorrelation length of $\langle E \rangle$ appears to be comparable to that of the topological charge.

preprint2016arXiv

Kaon semileptonic decays with $N_f=2+1+1$ HISQ fermions and physical light-quark masses

We discuss the reduction of errors in the calculation of the form factor $f_+^{K π}(0)$ with HISQ fermions on the $N_f=2+1+1$ MILC configurations from increased statistics on some key ensembles, new data on ensembles with lattice spacings down to 0.042 fm and the study of finite-volume effects within staggered ChPT. We also study the implications for the unitarity of the CKM matrix in the first row and for current tensions with leptonic determinations of $\vert V_{us}\vert$.

preprint2016arXiv

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

preprint2015arXiv

$|V_{ub}|$ from $B\toπ\ellν$ decays and (2+1)-flavor lattice QCD

We present a lattice-QCD calculation of the $B\toπ\ellν$ semileptonic form factors and a new determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad 2+1-flavor lattice configurations at four lattice spacings and light-quark masses down to 1/20 of the physical strange-quark mass. We extrapolate the lattice form factors to the continuum using staggered chiral perturbation theory in the hard-pion and SU(2) limits. We employ a model-independent $z$ parameterization to extrapolate our lattice form factors from large-recoil momentum to the full kinematic range. We introduce a new functional method to propagate information from the chiral-continuum extrapolation to the $z$ expansion. We present our results together with a complete systematic error budget, including a covariance matrix to enable the combination of our form factors with other lattice-QCD and experimental results. To obtain $|V_{ub}|$, we simultaneously fit the experimental data for the $B\toπ\ellν$ differential decay rate obtained by the BaBar and Belle collaborations together with our lattice form-factor results. We find $|V_{ub}|=(3.72\pm 0.16)\times 10^{-3}$ where the error is from the combined fit to lattice plus experiments and includes all sources of uncertainty. Our form-factor results bring the QCD error on $|V_{ub}|$ to the same level as the experimental error. We also provide results for the $B\toπ\ellν$ vector and scalar form factors obtained from the combined lattice and experiment fit, which are more precisely-determined than from our lattice-QCD calculation alone. These results can be used in other phenomenological applications and to test other approaches to QCD.

preprint2015arXiv

$B\toπ\ell\ell$ form factors for new-physics searches from lattice QCD

The rare decay $B\toπ\ell^+\ell^-$ arises from $b\to d$ flavor-changing neutral currents and could be sensitive to physics beyond the Standard Model. Here, we present the first $ab$-$initio$ QCD calculation of the $B\toπ$ tensor form factor $f_T$. Together with the vector and scalar form factors $f_+$ and $f_0$ from our companion work [J. A. Bailey $et~al.$, Phys. Rev. D 92, 014024 (2015)], these parameterize the hadronic contribution to $B\toπ$ semileptonic decays in any extension of the Standard Model. We obtain the total branching ratio ${\text{BR}}(B^+\toπ^+μ^+μ^-)=20.4(2.1)\times10^{-9}$ in the Standard Model, which is the most precise theoretical determination to date, and agrees with the recent measurement from the LHCb experiment [R. Aaij $et~al.$, JHEP 1212, 125 (2012)]. Note added: after this paper was submitted for publication, LHCb announced a new measurement of the differential decay rate for this process [T. Tekampe, talk at DPF 2015], which we now compare to the shape and normalization of the Standard-Model prediction.

preprint2015arXiv

Decay constants $f_B$ and $f_{B_s}$ from HISQ simulations

We give a progress report on a project aimed at a high-precision calculation of the decay constants $f_B$ and $f_{B_s}$ from simulations with HISQ heavy and light valence and sea quarks. Calculations are carried out with several heavy valence-quark masses on ensembles with 2+1+1 flavors of HISQ sea quarks at five lattice spacings and several light sea-quark mass ratios $m_{ud}/m_s$, including approximately physical sea-quark masses. This range of parameters provides excellent control of the continuum limit and of heavy-quark discretization errors. We present a preliminary error budget with projected uncertainties of 2.2~MeV and 1.5~MeV for $f_B$ and $f_{B_s}$, respectively.

preprint2015arXiv

Electromagnetic effects on the light hadron spectrum

For some time, the MILC Collaboration has been studying electromagnetic effects on light mesons. These calculations use fully dynamical QCD, but only quenched photons, which suffices to NLO in XPT. That is, the sea quarks are electrically neutral, while the valence quarks carry charge. For the photons we use the non-compact formalism. We have new results with lattice spacing as small as 0.045 fm and a large range of volumes. We consider how well chiral perturbation theory describes these results and the implications for light quark masses

preprint2015arXiv

The $B \to D \ell ν$ form factors at nonzero recoil and $|V_{cb}|$ from $2+1$-flavor lattice QCD

We present the first unquenched lattice-QCD calculation of the hadronic form factors for the exclusive decay $\overline{B} \rightarrow D \ell \overlineν$ at nonzero recoil. We carry out numerical simulations on fourteen ensembles of gauge-field configurations generated with 2+1 flavors of asqtad-improved staggered sea quarks. The ensembles encompass a wide range of lattice spacings (approximately 0.045 to 0.12 fm) and ratios of light (up and down) to strange sea-quark masses ranging from 0.05 to 0.4. For the $b$ and $c$ valence quarks we use improved Wilson fermions with the Fermilab interpretation, while for the light valence quarks we use asqtad-improved staggered fermions. We extrapolate our results to the physical point using rooted staggered heavy-light meson chiral perturbation theory. We then parameterize the form factors and extend them to the full kinematic range using model-independent functions based on analyticity and unitarity. We present our final results for $f_+(q^2)$ and $f_0(q^2)$, including statistical and systematic errors, as coefficients of a series in the variable $z$ and the covariance matrix between these coefficients. We then fit the lattice form-factor data jointly with the experimentally measured differential decay rate from BaBar to determine the CKM matrix element, $|V_{cb}|=(39.6 \pm 1.7_{\rm QCD+exp} \pm 0.2_{\rm QED})\times 10^{-3}$. As a byproduct of the joint fit we obtain the form factors with improved precision at large recoil. Finally, we use them to update our calculation of the ratio $R(D)$ in the Standard Model, which yields $R(D) = 0.299(11)$.

preprint2014arXiv

$B\toπ\ellν$ semileptonic form factors from unquenched lattice QCD and determination of $|V_{ub}|$

We compute the $B\toπ\ellν$ semileptonic form factors and update the determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad ensembles with $N_f=2+1$ sea quarks at four different lattice spacings in the range $a \approx 0.045$~fm to $0.12$~fm. The lattice form factors are extrapolated to the continuum limit using SU(2) staggered chiral perturbation theory in the hard pion limit, followed by an extrapolation in $q^2$ to the full kinematic range using a functional $z$-parameterization. The extrapolation is combined with the experimental measurements of the partial branching fraction to extract $|V_{ub}|$. Our preliminary result is $|V_{ub}|=(3.72\pm 0.14)\times 10^{-3}$, where the error reflects both the lattice and experimental uncertainties, which are now on par with each other.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from four-flavor lattice QCD with physical light quarks

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD using the experimentally determined value of $f_{π^+}$ for normalization. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors---up, down, strange, and charm---and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacings ranging from $a\approx 0.06$ fm to $0.15$ fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately two to four times smaller than those of the most precise previous lattice calculations. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from HISQ simulations

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Determination of $|V_{us}|$ from a lattice-QCD calculation of the $K\toπ\ellν$ semileptonic form factor with physical quark masses

We calculate the kaon semileptonic form factor $f_+(0)$ from lattice QCD, working, for the first time, at the physical light-quark masses. We use gauge configurations generated by the MILC collaboration with $N_f=2+1+1$ flavors of sea quarks, which incorporate the effects of dynamical charm quarks as well as those of up, down, and strange. We employ data at three lattice spacings to extrapolate to the continuum limit. Our result, $f_+(0) = 0.9704(32)$, where the error is the total statistical plus systematic uncertainty added in quadrature, is the most precise determination to date. Combining our result with the latest experimental measurements of $K$ semileptonic decays, one obtains the Cabibbo-Kobayashi-Maskawa matrix element $|V_{us}|=0.22290(74)(52)$, where the first error is from $f_+(0)$ and the second one is from experiment. In the first-row test of Cabibbo-Kobayashi-Maskawa unitarity, the error stemming from $|V_{us}|$ is now comparable to that from $|V_{ud}|$.

preprint2014arXiv

Gradient Flow Analysis on MILC HISQ Ensembles

We report on a preliminary scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ are computed using Symanzik flow and the cloverleaf definition of $\langle E \rangle$ on each ensemble. Then both scales and the meson masses $aM_π$ and $aM_K$ are adjusted for mistunings in the charm mass. Using a combination of continuum chiral perturbation theory and a Taylor series ansatz in the lattice spacing, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. Our preliminary results are $\sqrt{t_0} = 0.1422(7)$fm and $w_0 = 0.1732(10)$fm. We also find the continuum mass-dependence of $w_0$.

preprint2014arXiv

K semileptonic form factor with HISQ fermions at the physical point

We present results for the form factor $f_+^{K π}(0)$, needed to extract the CKM matrix element $|V_{us}|$ from experimental data on semileptonic $K$ decays, on the HISQ $N_f=2+1+1$ MILC configurations. The HISQ action is also used for the valence sector. The data set used for our final result includes three different values of the lattice spacing and data at the physical light quark masses. We discuss the error budget and how this calculation improves on our previous determination of $f_+^{K π}(0)$ on the asqtad $N_f=2+1$ MILC configurations.

preprint2014arXiv

The $D_s$, $D^+$, $B_s$ and $B$ decay constants from $2+1$ flavor lattice QCD

We present a study of the $D$ and $B$ leptonic decay constants on the MILC $N_f=2+1$ asqtad gauge ensembles using asqtad-improved staggered light quarks and clover heavy quarks in the Fermilab interpretation. Our previous analysis \cite{Bazavov:2011aa} computed the decay constants at lattice spacings $a \approx 0.14, 0.11$ and $0.083$ fm. We have extended the simulations to finer $a \approx 0.058$ and $0.043$ fm lattice spacings, and have also increased statistics; this allows us to address many important sources of uncertainty. Technical advances include a two-step two-point fit procedure, better tuning of the heavy quark masses and a better determination of the axial-vector current matching. The present analysis remains blinded, so here we focus on the improvements and their predicted impact on the error budget compared to the prior analysis.

preprint2014arXiv

The equation of state in (2+1)-flavor QCD

We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent $N_τ=6,~8,~10$, and $12$. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range $(130-400)$ MeV. We compare our results with previous calculations, and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, $145~ {\rm MeV} \leq T \leq 163$ MeV, defined by the chiral transition, is $ε_c=(0.18-0.5)~{\rm GeV}/{\rm fm}^3$, $i.e.$, $(1.2-3.1)\ ε_{\rm nuclear}$. At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters $r_0$ from the static quark potential and $w_0$ from the gradient flow.

preprint2014arXiv

Update of $|V_{cb}|$ from the $\bar{B}\to D^*\ell\barν$ form factor at zero recoil with three-flavor lattice QCD

We compute the zero-recoil form factor for the semileptonic decay $\bar{B}^0\to D^{*+}\ell^-\barν$ (and modes related by isospin and charge conjugation) using lattice QCD with three flavors of sea quarks. We use an improved staggered action for the light valence and sea quarks (the MILC \asqtad\ configurations), and the Fermilab action for the heavy quarks. Our calculations incorporate higher statistics, finer lattice spacings, and lighter quark masses than our 2008 work. As a byproduct of tuning the new data set, we obtain the $D_s$ and $B_s$ hyperfine splittings with few-MeV accuracy. For the zero-recoil form factor, we obtain $\mathcal{F}(1)=0.906(4)(12)$, where the first error is statistical and the second is the sum in quadrature of all systematic errors. With the latest HFAG average of experimental results and a cautious treatment of QED effects, we find $|V_{cb}| = (39.04 \pm 0.49_\text{expt} \pm 0.53_\text{QCD} \pm 0.19_\text{QED})\times10^{-3}$. The QCD error is now commensurate with the experimental error.

preprint2013arXiv

Charmed and strange pseudoscalar meson decay constants from HISQ simulations

We update our determinations of $f_{D^+}$, $f_{D_s}$, $f_K$, and quark mass ratios from simulations with four flavors of HISQ dynamical quarks. The availability of ensembles with light quarks near their physical mass means that we can extract physical results with only small corrections for valence- and sea-quark mass mistunings instead of a chiral extrapolation. The adjusted valence-quark masses and lattice spacings may be determined from an ensemble-by-ensemble analysis, and the results for the quark mass ratios then extrapolated to the continuum limit. Our central values of the charmed meson decay constants, however, come from an alternative analysis, which uses staggered chiral perturbation theory for the heavy-light mesons, and allows us to incorporate data at unphysical quark masses where statistical errors are often smaller. A jackknife analysis propagated through all of these steps takes account of the correlations among all the quantities used in the analysis. Systematic errors from the finite spatial size and EM effects are estimated by varying the parameters in the analysis, and systematic errors from the assumptions in the continuum extrapolation are estimated from the spread of values from different extrapolations.

preprint2013arXiv

Electromagnetic contributions to pseudoscalar masses

We report on the calculation by the MILC Collaboration of the electromagnetic effects on kaon and pion masses. These masses are computed in QCD with dynamical (asqtad staggered) quarks plus quenched photons at three lattice spacings varying from 0.12 to 0.06 fm. The masses are fit to staggered chiral perturbation theory with NLO electromagnetic terms, as well as analytic terms at higher order. We extrapolate the results to physical light-quark masses and to the continuum limit. At the current stage of the analysis, most, but not all, of the systematic errors have been estimated. The main goal is the comparison of kaon electromagnetic splittings to those of the pion, i.e., an evaluation of the corrections to "Dashen's theorem." This in turn will allow us to significantly reduce the systematic errors in our determination of m_u/m_d.

preprint2013arXiv

Heavy-meson semileptonic decays for the Standard Model and beyond

We calculate the form factors for the semileptonic decays $B_s\to K\ellν$ and $B\to K\ell\ell$ with lattice QCD. We work at several lattice spacings and a range of light quark masses, using the MILC 2+1-flavor asqtad ensembles. We use the Fermilab method for the $b$ quark. We obtain chiral-continuum extrapolations for $E_K$ up to $\sim1.2$ GeV and then extend to the entire kinematic range with the model-independent $z$ expansion.

preprint2013arXiv

Kaon semileptonic vector form factor and determination of |V_{us}| using staggered fermions

Using staggered fermions and twisted boundary conditions, we calculate the K meson semileptonic decay vector form factor at zero momentum transfer. The HISQ formulation is used for the valence quarks, while the sea quarks are simulated with the asqtad action (MILC N_f=2+1 configurations). For the chiral and continuum extrapolation we use two-loop continuum CHPT, supplemented by partially quenched staggered CHPT at one loop. Our result is f_+^{Kπ}(0) = 0.9667+-0.0023+-0.0033, where the first error is statistical and the second is the sum in quadrature of the systematic uncertainties. This result is the first N_f=2+1 calculation with two lattice spacings and a controlled continuum extrapolation. It is also the most precise result to date for the vector form factor and, although the central value is larger than previous unquenched lattice calculations, it is compatible with them within errors. Combining our value for f_+^{Kπ}(0) with the latest experimental measurements of K semileptonic decays, we obtain |V_{us}| = 0.2238+-0.0009+-0.0005, where the first error is from f_+^{Kπ}(0) and the second one is experimental. As a byproduct of our calculation, we obtain the combination of low-energy constants [C_{12}^r+C_{34}^r-(L_5^r)^2](M_ρ) = (3.62+-1.00)x10^{-6}.

preprint2013arXiv

Lattice QCD ensembles with four flavors of highly improved staggered quarks

We present results from our simulations of quantum chromodynamics (QCD) with four flavors of quarks: u, d, s, and c. These simulations are performed with a one-loop Symanzik improved gauge action, and the highly improved staggered quark (HISQ) action. We are generating gauge configurations with four values of the lattice spacing ranging from 0.06 fm to 0.15 fm, and three values of the light quark mass, including the value for which the Goldstone pion mass is equal to the physical pion mass. We discuss simulation algorithms, scale setting, taste symmetry breaking, and the autocorrelations of various quantities. We also present results for the topological susceptibility which demonstrate the improvement of the HISQ configurations relative to those generated earlier with the asqtad improved staggered action.

preprint2013arXiv

Leptonic decay-constant ratio f_{K^+}/f_{pi^+} from lattice QCD with physical light quarks

A calculation of the ratio of leptonic decay constants f_{K^+}/f_{π^+} makes possible a precise determination of the ratio of CKM matrix elements |V_{us}|/|V_{ud}| in the Standard Model, and places a stringent constraint on the scale of new physics that would lead to deviations from unitarity in the first row of the CKM matrix. We compute f_{K^+}/f_{π^+} numerically in unquenched lattice QCD using gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, strange, and charm. We analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with simulated pion masses down to the physical value 135 MeV. We obtain f_{K^+}/f_{π^+} = 1.1947(26)(37), where the errors are statistical and total systematic, respectively. This is our first physics result from our N_f = 2+1+1 ensembles, and the first calculation of f_{K^+}/f_{π^+} from lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD determination of f_{K^+}/f_{π^+}, with an error comparable to the current world average. When combined with experimental measurements of the leptonic branching fractions, it leads to a precise determination of |V_{us}|/|V_{ud}| = 0.2309(9)(4) where the errors are theoretical and experimental, respectively.

preprint2013arXiv

Symanzik flow on HISQ ensembles

We report on a scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The lattice scale $w_0/a$, originally proposed by the BMW collaboration, is computed using Symanzik flow at four lattice spacings ranging from 0.15 to 0.06 fm. With a Taylor series ansatz, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We give a preliminary determination of the scale $w_0$ in physical units, along with associated systematic errors, and compare with results from other groups. We also present a first estimate of autocorrelation lengths as a function of flowtime for these ensembles.

preprint2013arXiv

Update on the 2+1+1 flavor QCD equation of state with HISQ

We present recent results on the QCD equation of state with 2+1+1 flavors of highly improved staggered quarks (HISQ). We focus on three sets of ensembles with temporal extent 6, 8 and 10, that reach up to temperatures of 967, 725 and 580 MeV, respectively. The strange and charm quark masses are tuned to the physical values and the light quarks mass is set to one fifth of the strange. This corresponds to a Goldstone pion of about 300 MeV.

preprint2012arXiv

B_s->D_s/B->D Semileptonic Form-Factor Ratios and Their Application to BR(B^0_s->μ^+μ^-)

We calculate form-factor ratios between the semileptonic decays \bar{B}->D^+\ell^-\barν and \bar{B}_s->D_s^+\ell^-\barν with lattice QCD. These ratios are a key theoretical input in a new strategy to determine the fragmentation fractions of the neutral B decays, which are needed for measurements of BR(B^0_s-> μ^+μ^-). We use the MILC ensembles of gauge configurations with 2+1 flavors of sea quarks at two lattice spacings of approximately 0.12 fm and 0.09 fm. We use the model-independent z-parametrization to extrapolate our simulation results at small recoil toward maximum recoil. Our results for the form-factor ratios are $f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_K) =1.046(44)_{stat.}(15)_{syst.}$ and $f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_π)=1.054(47)_{stat.}(17)_{syst.}$. In contrast to a QCD sum-rule calculation, no significant departure from U-spin (d<->s) symmetry is observed.

preprint2012arXiv

Critical analysis of topological charge determination in the background of center vortices in SU(2) lattice gauge theory

We analyze topological charge contributions from classical SU(2) center vortices with shapes of planes and spheres using different topological charge definitions, namely the center vortex picture of topological charge, a discrete version of FF in the plaquette and hypercube definitions and the lattice index theorem. For the latter the zeromodes of the Dirac operator in the fundamental and adjoint representations using both the overlap and asqtad staggered fermion formulations are investigated. We find several problems for the individual definitions and discuss the discrepancies between the different topological charge definitions. Our results show that the interpretation of topological charge in the background of center vortices is rather subtle.

preprint2012arXiv

Finite-Volume Scaling of the Wilson-Dirac Operator Spectrum

The microscopic spectral density of the Hermitian Wilson-Dirac operator is computed numerically in quenched lattice QCD. We demonstrate that the results given for fixed index of the Wilson-Dirac operator can be matched by the predictions from Wilson chiral perturbation theory. We test successfully the finite volume and the mass scaling predicted by Wilson chiral perturbation theory at fixed lattice spacing.

preprint2012arXiv

Kaon semileptonic decay form factors with HISQ valence quarks

We report on the status of our kaon semileptonic form factor calculations using the highly-improved staggered quark (HISQ) formulation to simulate the valence fermions. We present results for the form factor f_+^{K π}(0) on the asqtad N_f=2+1 MILC configurations, discuss the chiral-continuum extrapolation, and give a preliminary estimate of the total error. We also present a more preliminary set of results for the same form factor but with the sea quarks also simulated with the HISQ action; these results include data at the physical light quark masses. The improvements that we expect to achieve with the use of the HISQ configurations and simulations at the physical quark masses are briefly discussed.

preprint2012arXiv

Neutral B-meson mixing from three-flavor lattice QCD: Determination of the SU(3)-breaking ratio ξ

We study SU(3)-breaking effects in the neutral B_d-\bar B_d and B_s-\bar B_s systems with unquenched N_f=2+1 lattice QCD. We calculate the relevant matrix elements on the MILC collaboration's gauge configurations with asqtad-improved staggered sea quarks. For the valence light-quarks (u, d, and s) we use the asqtad action, while for b quarks we use the Fermilab action. We obtain ξ=f_{B_s}\sqrt{B_{B_s}}/f_{B_d}\sqrt{B_{B_d}}=1.268+-0.063. We also present results for the ratio of bag parameters B_{B_s}/B_{B_d} and the ratio of CKM matrix elements |V_{td}|/|V_{ts}|. Although we focus on the calculation of ξ, the strategy and techniques described here will be employed in future extended studies of the B mixing parameters ΔM_{d,s} and ΔΓ_{d,s} in the Standard Model and beyond.

preprint2012arXiv

New Ways to Determine Low-Energy Constants with Wilson Fermions

We show how the leading physical and Wilson low-energy constants associated with Wilson fermions in lattice gauge theory can be determined individually by using spectral information of the Wilson Dirac operator with fixed index at finite volume. The methods are demonstrated in simulations with leading-order improved Wilson fermions. In addition to the expected suppression of the leading term in Wilson chiral perturbation theory we observe a substantial reduction also of the higher-order Wilson low-energy constants.

preprint2012arXiv

Refining new-physics searches in B -> D tau nu decay with lattice QCD

The semileptonic decay channel B -> D tau nu is sensitive to the presence of a scalar current, such as that mediated by a charged-Higgs boson. Recently the BaBar experiment reported the first observation of the exclusive semileptonic decay B -> D tau nu, finding an approximately 2-sigma disagreement with the Standard-Model prediction for the ratio R(D)=BR(B->D tau nu)/BR(B->D l nu), where l=e,mu. We compute this ratio of branching fractions using hadronic form factors computed in unquenched lattice QCD and obtain R(D) = 0.316(12)(7), where the errors are statistical and total systematic, respectively. This result is the first Standard-Model calculation of R(D) from ab initio full QCD. Its error is smaller than that of previous estimates, primarily due to the reduced uncertainty in the scalar form factor f_0(q^2). Our determination of R(D) is approximately 1-sigma higher than previous estimates and, thus, reduces the tension with experiment. We also compute R(D) in models with electrically charged scalar exchange, such as the type II two-Higgs doublet model. Once again, our result is consistent with, but approximately 1-sigma higher than, previous estimates for phenomenologically relevant values of the scalar coupling in the type II model. As a byproduct of our calculation, we also present the Standard-Model prediction for the longitudinal polarization ratio P_L (D)= 0.325(4)(3).

preprint2012arXiv

Status of the MILC calculation of electromagnetic contributions to pseudoscalar masses

We calculate pseudoscalar masses on gauge configurations containing the effects of 2+1 flavors of dynamical asqtad quarks and quenched electromagnetism. The lattice spacings vary from 0.12 to 0.06 fm. The masses are fit with staggered chiral perturbation theory including NLO electromagnetic terms. We attempt to extract the fit parameters for the electromagnetic contributions, while taking into account the finite volume effects, and extrapolate them to the physical limit.

preprint2012arXiv

The chiral and deconfinement aspects of the QCD transition

We present results on the chiral and deconfinement properties of the QCD transition at finite temperature. Calculations are performed with 2+1 flavors of quarks using the p4, asqtad and HISQ/tree actions. Lattices with temporal extent N_tau=6, 8 and 12 are used to understand and control discretization errors and to reliably extrapolate estimates obtained at finite lattice spacings to the continuum limit. The chiral transition temperature is defined in terms of the phase transition in a theory with two massless flavors and analyzed using O(N) scaling fits to the chiral condensate and susceptibility. We find consistent estimates from the HISQ/tree and asqtad actions and our main result is T_c=154 +/- 9 MeV.

preprint2012arXiv

U_A(1) breaking at finite temperature from the Dirac spectrum with the dynamical HISQ action

We investigate $U_A(1)$ breaking above $T_c$ in terms of the Dirac spectrum on configurations with (2+1)-flavors, using the HISQ action. The strange quark mass is at its physical value. We use several light quark masses corresponding to the Goldstone pion masses in the range of about 115 -- 230 MeV on lattices of size 32$^3 \times$8 and 48$^3 \times$8. We calculate the 100 lowest-lying Dirac eigenvalues at temperatures below and above $T_c$. We investigate the volume dependence of the Dirac eigenvalue density to determine whether there is a gap around zero, which can appear if $U_A(1)$ symmetry is restored in the chiral symmetric phase. We also investigate the quark mass dependence of the Dirac eigenvalue density at zero and check whether there is a linear behavior that would signal the $U_A(1)$ breaking above $T_c$.

preprint2011arXiv

B- and D-meson decay constants from three-flavor lattice QCD

We calculate the leptonic decay constants of B_{(s)} and D_{(s)} mesons in lattice QCD using staggered light quarks and Fermilab bottom and charm quarks. We compute the heavy-light meson correlation functions on the MILC asqtad-improved staggered gauge configurations which include the effects of three light dynamical sea quarks. We simulate with several values of the light valence- and sea-quark masses (down to ~m_s/10) and at three lattice spacings (a ~ 0.15, 0.12, and 0.09 fm) and extrapolate to the physical up and down quark masses and the continuum using expressions derived in heavy-light meson staggered chiral perturbation theory. We renormalize the heavy-light axial current using a mostly nonperturbative method such that only a small correction to unity must be computed in lattice perturbation theory and higher-order terms are expected to be small. We obtain f_{B^+} = 196.9(8.9) MeV, f_{B_s} = 242.0(9.5) MeV, f_{D^+} = 218.9(11.3) MeV, f_{D_s} = 260.1(10.8) MeV, and the SU(3) flavor-breaking ratios f_{B_s}/f_{B} = 1.229(26) and f_{D_s}/f_{D} = 1.188(25), where the numbers in parentheses are the total statistical and systematic uncertainties added in quadrature.

preprint2011arXiv

Eigenvalue distribution of the Dirac operator at finite temperature with (2+1)-flavor dynamical quarks using the HISQ action

We report on the behavior of the eigenvalue distribution of the Dirac operator in (2+1)-flavor QCD at finite temperature, using the HISQ action. We calculate the eigenvalue density at several values of the temperature close to the pseudocritical temperature. For this study we use gauge field configurations generated on lattices of size $32^3 \times 8$ with two light quark masses corresponding to pion masses of about 160 and 115 MeV. We find that the eigenvalue density below $T_c$ receives large contributions from near-zero modes which become smaller as the temperature increases or the light quark mass decreases. Moreover we find no clear evidence for a gap in the eigenvalue density up to 1.1$T_c$. We also analyze the eigenvalue density near $T_c$ where it appears to show a power-law behavior consistent with what is expected in the critical region near the second order chiral symmetry restoring phase transition in the massless limit.

preprint2011arXiv

Properties of light pseudoscalars from lattice QCD with HISQ ensembles

We fit lattice-QCD data for light-pseudoscalar masses and decay constants, from HISQ configurations generated by MILC, to SU(3) staggered chiral perturbation theory. At present such fits have rather high values of chi^2/d.o.f., possibly due to the lack of ensembles with lighter-than-physical sea strange-quark masses. We propose solutions to this problem for future work. We also perform simple linear interpolations near the physical point on two ensembles with different lattice spacings, and obtain the preliminary result (f_K / f_pi)^phys = 1.1872(41) in the continuum limit.

preprint2011arXiv

Tuning Fermilab Heavy Quarks in 2+1 Flavor Lattice QCD with Application to Hyperfine Splittings

We report the non-perturbative tuning of parameters--- kappa_c, kappa_b, and kappa_crit ---that determine the heavy-quark mass in the Fermilab action. This requires the computation of the masses of Ds^(*) and Bs^(*) mesons comprised of a Fermilab heavy quark and a staggered light quark. Additionally, we report the hyperfine splittings for Ds and Bs mesons as a cross-check of our simulation and analysis methods. We find a splitting of 145 +/- 15 MeV for the Ds system and 40 +/- 9 MeV for the Bs system. These are in good agreement with the Particle Data Group average values of 143.9 +/- 0.4 MeV and 46.1 +/- 1.5 MeV, respectively. The calculations are carried out with the MILC 2+1 flavor gauge configurations at three lattice spacings $a$ approximately 0.15, 0.12, and 0.09 fm.

preprint2010arXiv

B -> D* l nu at zero recoil: an update

We present an update of our calculation of the form factor for B -> D* l nu at zero recoil, with higher statistics and finer lattices. As before, we use the Fermilab action for b and c quarks, the asqtad staggered action for light valence quarks, and the MILC ensembles for gluons and light quarks (Lüscher-Weisz married to 2+1 rooted staggered sea quarks). In this update, we have reduced the total uncertainty on F(1) from 2.6% to 1.7%. At Lattice2010 we presented a still-blinded result, but this writeup includes the unblinded result from the September 2010 CKM workshop.

preprint2010arXiv

QCD thermodynamics with nonzero chemical potential at $N_t=6$ and effects from heavy quarks

We extend our work on QCD thermodynamics with 2+1 quark flavors at nonzero chemical potential to finer lattices with $N_t=6$. We study the equation of state and other thermodynamic quantities, such as quark number densities and susceptibilities, and compare them with our previous results at $N_t=4$. We also calculate the effects of the addition of the charm and bottom quarks on the equation of state at zero and nonzero chemical potential. These effects are important for cosmological studies of the early Universe.

preprint2010arXiv

Scaling studies of QCD with the dynamical HISQ action

We study the lattice spacing dependence, or scaling, of physical quantities using the highly improved staggered quark (HISQ) action introduced by the HPQCD/UKQCD collaboration, comparing our results to similar simulations with the asqtad fermion action. Results are based on calculations with lattice spacings approximately 0.15, 0.12 and 0.09 fm, using four flavors of dynamical HISQ quarks. The strange and charm quark masses are near their physical values, and the light-quark mass is set to 0.2 times the strange-quark mass. We look at the lattice spacing dependence of hadron masses, pseudoscalar meson decay constants, and the topological susceptibility. In addition to the commonly used determination of the lattice spacing through the static quark potential, we examine a determination proposed by the HPQCD collaboration that uses the decay constant of a fictitious "unmixed s bar s" pseudoscalar meson. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses.

preprint2010arXiv

Semileptonic decays of K and D mesons in 2+1 flavor QCD

The experimentally measured rates of the semileptonic decays K -> pi l nu and D -> K(pi) l nu can be combined with lattice calculations of the associated form factors to precisely extract the CKM matrix elements |V_{us}| and |V_{cs(d)}|. We report on the status of form factor calculations with Fermilab charm quarks and staggered light quarks on the 2+1 flavor asqtad staggered MILC ensembles. Analysis of data for the D -> pi l nu form factor provides a nontrivial test of our methods via comparison with CLEO data. We discuss the use of HISQ valence quarks to calculate the K -> pi l nu form factor f_+^{K pi}(0) and describe tests of our method.

preprint2010arXiv

Simulations with dynamical HISQ quarks

We report on the status of a program of generating and using configurations with four flavors of dynamical quarks, using the HISQ action. We study the lattice spacing dependence of physical quantities in these simulations, using runs at several lattice spacings, but with the light quark mass held fixed at two tenths of the strange quark mass. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses. We also discuss methods for setting the scale, or assigning a lattice spacing to ensembles run at unphysical parameters.

preprint2010arXiv

Staggered chiral perturbation theory in the two-flavor case and SU(2) analysis of the MILC data

In the light pseudoscalar sector, we study rooted staggered chiral perturbation theory in the two-flavor case. The pion mass and decay constant are calculated through NLO for a partially-quenched theory. In the limit where the strange quark mass is large compared to the light quark masses and the taste splittings, we show that the SU(2) staggered chiral theory emerges from the SU(3) staggered chiral theory, as expected. Explicit relations between SU(2) and SU(3) low energy constants and taste-violating parameters are given. A brief summary of updated SU(2) chiral fits to the MILC lattice data is presented.

preprint2010arXiv

Topological susceptibility with the asqtad action

Chiral perturbation theory predicts that in quantum chromodynamics (QCD), light dynamical quarks suppress the gauge-field topological susceptibility of the vacuum. The degree of suppression depends on quark multiplicity and masses. It provides a strong consistency test for fermion formulations in lattice QCD. Such tests are especially important for staggered fermion formulations that lack a full chiral symmetry and use the "fourth-root" procedure to achieve the desired number of sea quarks. Over the past few years we have measured the topological susceptibility on a large database of 18 gauge field ensembles, generated in the presence of 2+1 flavors of dynamical asqtad quarks with up and down quark masses ranging from 0.05 to 1 in units of the strange quark mass and lattice spacings ranging from 0.045 fm to 0.12 fm. Our study also includes three quenched ensembles with lattice spacings ranging from 0.06 to 0.12 fm. We construct the topological susceptibility from the integrated point-to-point correlator of the discretized topological charge density F-Fdual. To reduce its variance, we model the asymptotic tail of the correlator. The continuum extrapolation of our results for the topological susceptibility agrees nicely at small quark mass with the predictions of lowest-order SU(3) chiral perturbation theory, thus lending support to the validity of the fourth-root procedure.

preprint2009arXiv

Full nonperturbative QCD simulations with 2+1 flavors of improved staggered quarks

Dramatic progress has been made over the last decade in the numerical study of quantum chromodynamics (QCD) through the use of improved formulations of QCD on the lattice (improved actions), the development of new algorithms and the rapid increase in computing power available to lattice gauge theorists. In this article we describe simulations of full QCD using the improved staggered quark formalism, ``asqtad'' fermions. These simulations were carried out with two degenerate flavors of light quarks (up and down) and with one heavier flavor, the strange quark. Several light quark masses, down to about 3 times the physical light quark mass, and six lattice spacings have been used. These enable controlled continuum and chiral extrapolations of many low energy QCD observables. We review the improved staggered formalism, emphasizing both advantages and drawbacks. In particular, we review the procedure for removing unwanted staggered species in the continuum limit. We then describe the asqtad lattice ensembles created by the MILC Collaboration. All MILC lattice ensembles are publicly available, and they have been used extensively by a number of lattice gauge theory groups. We review physics results obtained with them, and discuss the impact of these results on phenomenology. Topics include the heavy quark potential, spectrum of light hadrons, quark masses, decay constant of light and heavy-light pseudoscalar mesons, semileptonic form factors, nucleon structure, scattering lengths and more. We conclude with a brief look at highly promising future prospects.

preprint2009arXiv

HISQ action in dynamical simulations

We report on recent progress in employing the Highly Improved Staggered Quark (HISQ) action introduced by the HPQCD/UKQCD collaboration in simulations with dynamical fermions. The HISQ action is an order $a^2$ Symanzik-improved action with further suppressed taste symmetry violations. The improvement in taste symmetry is achieved by introducing Fat7 smearing of the original gauge links and reunitarization (projection to an element of U(3) or SU(3)) followed by Asq-type smearing. Major challenges for calculating the fermion force are related to the reunitarization step. We present a preliminary study of the HISQ action on two 2+1+1 flavor ensembles with the lattice spacing roughly equivalent to the MILC asqtad a=0.125 and 0.09 fm ensembles.

preprint2009arXiv

Progress on charm semileptonic form factors from 2+1 flavor lattice QCD

Lattice calculations of the form factors for the charm semileptonic decays D to K l nu and D to pi l nu provide inputs to direct determinations of the CKM matrix elements |V(cs)| and |V(cd)| and can be designed to validate calculations of the form factors for the bottom semileptonic decays B to pi l nu and B to K l l-bar. We are using Fermilab charm (bottom) quarks and asqtad staggered light quarks on the 2+1 flavor asqtad MILC ensembles to calculate the charm (bottom) form factors. We outline improvements to the previous calculation of the charm form factors and detail our progress. We expect our current round of data production to allow us to reduce the theoretical uncertainties in |V(cs)| and |V(cd)| from 10.5% and 11%, respectively, to about 7%.

preprint2009arXiv

Progress on four flavor QCD with the HISQ action

We describe recent progress on generation of gauge configurations using the Highly Improved Staggered Quark (HISQ) action that was designed by the HPQCD/UKQCD collaboration. The HISQ action requires two levels of smearing with a reunitarization of the links before the second smearing. We describe how we deal with the occurrence of occasional large forces arising from the reunitarization step. The MILC collaboration is currently generating ensembles with approximate lattice spacings of 0.15, 0.12, 0.09, and 0.06 fm, with the strange and charm quark masses close to their physical values and the mass of the light quarks m_l set to 0.2 m_s. We present recent results for pion taste splittings, light hadron masses, the static potential, the eta_c dispersion relation and the topological susceptibility.

preprint2009arXiv

The Ds and D+ Leptonic Decay Constants from Lattice QCD

We present the leptonic decay constants fDs and fD+ computed on the MILC collaboration's 2+1 flavor asqtad gauge ensembles. We use clover heavy quarks with the Fermilab interpretation and improved staggered light quarks. The simultaneous chiral and continuum extrapolation, which determines both decay constants, includes partially-quenched lattice results at lattice spacings a ~ 0:09, 0:12 and 0:15 fm. We have made several recent improvements in our analysis: a) we include terms in the fit describing leading order heavy-quark discretization effects, b) we have adopted a more precise input r1 value consistent with our other D and B meson studies, c) we have retuned the input bare charm masses based upon the new r1. Our preliminary results are fDs = 260 +/-10 MeV and fD+ = 217 +/-10 MeV.

preprint2008arXiv

Electromagnetic splittings of hadrons from improved staggered quarks in full QCD

We present our initial study of the electromagnetic splittings of charged and neutral mesons, and the violation of Dashen's theorem. Hadron masses are calculated on MILC N_f=2+1 QCD ensembles at lattice spacing \approx 0.15fm, together with quenched non-compact U(1) configurations. The O(a^2) tadpole improved staggered quark (asqtad) action is used both for the sea quarks and for six different valence quark masses. Chiral extrapolations are performed using partially quenched chiral perturbation theory incorporating electromagnetic corrections.

preprint2007arXiv

Baryon masses with improved staggered quarks

The MILC collaboration's simulations with improved staggered quarks are being extended with runs at a lattice spacing of 0.06 fm with quark masses down to one tenth the strange quark mass. We give a brief introduction to these new simulations and the determination of the lattice spacing. Then we combine these new runs with older results to study the masses of the nucleon and the Omega minus in the continuum and chiral limits.

preprint2007arXiv

Status of the MILC light pseudoscalar meson project

We discuss the current status of our calculation of the physics of pi and K mesons using three dynamical flavors of improved staggered quarks. This year, we have a new ensemble with a lattice spacing of 0.06 fm and a light sea mass of 0.2 m_s, as well as significant increases in statistics at several coarser lattice spacings and/or heavier sea masses. Results for decay constants, quark masses, low energy constants, condensates, and V_{us} are presented.

preprint2007arXiv

The 2+1 flavor topological susceptibility from the asqtad action at 0.06 fm

We report new data for the topological susceptibility computed on 2+1 flavor dynamical configurations with lattice spacing 0.06 fm, generated with the asqtad action. The topological susceptibility is computed by HYP smearing and compared with rooted staggered chiral perturbation theory as the pion mass goes to zero. At 0.06 fm, the raw data is already quite close to the continuum extrapolated values obtained from coarser lattices. These results provide a further test of the asqtad action with rooted staggered flavors.

preprint2006arXiv

The QCD equation of state with asqtad staggered fermions

We report on our result for the equation of state (EOS) with a Symanzik improved gauge action and the asqtad improved staggered fermion action at $N_t=4$ and 6. In our dynamical simulations with 2+1 flavors we use the inexact R algorithm and here we estimate the finite step-size systematic error on the EOS. Finally we discuss the non-zero chemical potential extension of the EOS and give some preliminary results.

preprint2006arXiv

Update on the physics of light pseudoscalar mesons

We present an update of the MILC investigation of the properties of light pseudoscalar mesons using three flavors of improved staggered quarks. Results are presented for the $π$ and $K$ leptonic decay constants, the CKM matrix element $V_{us}$, the up, down and strange quark masses, and the coefficients of the $O(p^4)$ chiral lagrangian. We have new data for lattice spacing $a \approx 0.15$ fm with several values of the light quark mass down to one-tenth the strange quark mass, higher statistics for $a \approx 0.09$ fm with the light quark mass equal to one-tenth the strange quark mass, and initial results for our smallest lattice spacing, $a \approx 0.06$ fm with light quark mass two-fifths of the strange quark mass.

preprint2005arXiv

Group-theoretical construction of extended baryon operators in lattice QCD

The design and implementation of large sets of spatially-extended, gauge-invariant operators for use in determining the spectrum of baryons in lattice QCD computations are described. Group-theoretical projections onto the irreducible representations of the symmetry group of a cubic spatial lattice are used in all isospin channels. The operators are constructed to maximize overlaps with the low-lying states of interest, while minimizing the number of sources needed in computing the required quark propagators. Issues related to the identification of the spin quantum numbers of the states in the continuum limit are addressed.

preprint2005arXiv

The Equation of State for QCD with 2+1 Flavors of Quarks

We report results for the interaction measure, pressure and energy density for nonzero temperature QCD with 2+1 flavors of improved staggered quarks. In our simulations we use a Symanzik improved gauge action and the Asqtad $O(a^2)$ improved staggered quark action for lattices with temporal extent $N_t=4$ and 6. The heavy quark mass $m_s$ is fixed at approximately the physical strange quark mass and the two degenerate light quarks have masses $m_{ud} =0.1m_s$ or $0.2m_s$. The calculation of the thermodynamic observables employs the integral method where energy density and pressure are obtained by integration over the interaction measure.

preprint2005arXiv

Update on pi and K Physics

We present an update of the MILC studies of the physics of light pseudoscalars using improved staggered fermions. New runs at lighter quark mass, as well as increased statistics for older sets, are enabling us to improve the results for decay constants in full QCD. In addition, we have analyzed quenched runs at two different lattice spacings. This makes possible a test of the applicability of staggered chiral perturbation theory in a different context.

preprint2004arXiv

Analysis of $N^{\star}$ spectra using matrices of correlation functions based on irreducible baryon operators

We present results for ground and excited-state nucleon masses in quenched lattice QCD using anisotropic lattices. Group theoretical constructions of local and nonlocal straight-link irreducible operators are used to obtain suitable sources and sinks. Matrices of correlation functions are diagonalized to determine the eigenvectors. Both chi-square fitting and Bayesian inference with an entropic prior are used to extract masses from the correlation functions in a given channel. We observe clear separation of the excited state masses from the ground state mass. States of spin >= 5/2 have been isolated by use of G_2 operators.

preprint2004arXiv

Baryonic sources using irreducible representations of the double-covered octahedral group

Irreducible representations (IRs) of the double-covered octahedral group are used to construct lattice source and sink operators for three-quark baryons. The goal is to achieve a good coupling to higher spin states as well as ground states. Complete sets of local and nonlocal straight-link operators are explicitly shown for isospin 1/2 and 3/2 baryons. The orthogonality relations of the IR operators are confirmed in a quenched lattice simulation.

preprint2004arXiv

Leptonic decay constants f_Ds and f_D in three flavor lattice QCD

We determine the leptonic decay constants in three flavor unquenched lattice QCD. We use O(a^2)-improved staggered light quarks and O(a)-improved charm quarks in the Fermilab heavy quark formalism. Our preliminary results, based upon an analysis at a single lattice spacing, are f_Ds = 263(+5-9)(+/-24) MeV and f_D = 225(+11-13)(+/-21) MeV. In each case, the first reported error is statistical while the is the combined systematic uncertainty.

preprint2003arXiv

Baryon operators and spectroscopy in lattice QCD

The construction of the operators and correlators required to determine the excited baryon spectrum is presented, with the aim of exploring the spatial and spin structure of the states while minimizing the number of propagator inversions. The method used to construct operators that transform irreducibly under the symmetries of the lattice is detailed, and the properties of example operators are studied using domain-wall fermion valence propagators computed on MILC asqtad dynamical lattices.

preprint2003arXiv

Baryonic Operators for Lattice Simulations

The construction of baryonic operators for determining the N* excitation spectrum is discussed. The operators are designed with one eye towards maximizing overlaps with the low-lying states of interest, and the other eye towards minimizing the number of sources needed in computing the required quark propagators. Issues related to spin identification are outlined. Although we focus on tri-quark baryon operators, the construction method is applicable to both mesons and penta-quark operators.

preprint2003arXiv

Spectroscopy using the Anisotropic Clover Action

The calculation of the light-hadron spectrum in the quenched approximation to QCD using an anisotropic clover fermion action is presented. The tuning of the parameters of the action is discussed, using the pion and rho dispersion relation. The adoption of an anisotropic lattice provides clear advantages in the determination of the baryonic resonances, and in particular that of the so-called Roper resonance, the lightest radial excitation of the nucleon.

preprint2003arXiv

The Phase Diagram of High Temperature QCD with Three Flavors of Improved Staggered Quarks

We report on progress in our study of high temperature QCD with three flavors of improved staggered quarks. Simulations are being carried out with three degenerate quarks with masses less than or equal to the strange quark mass, $m_s$, and with degenerate up and down quarks with masses in the range $0.1 m_s \leq m_{u,d}\leq 0.6 m_s$, and the strange quark mass fixed near its physical value. For the quark masses studied to date we find rapid crossovers, which sharpen as the quark mass is reduced, rather than bona fide phase transitions.

preprint2002arXiv

Gluon Propagator on Coarse Lattices in Laplacian Gauges

The Laplacian gauge is a nonperturbative gauge fixing that reduces to Landau gauge in the asymptotic limit. Like Landau gauge, it respects Lorentz invariance, but it is free of Gribov copies; the gauge fixing is unambiguous. In this paper we study the infrared behavior of the lattice gluon propagator in Laplacian gauge by using a variety of lattices with spacings from $a = 0.125$ to 0.35 fm, to explore finite volume and discretization effects. Three different implementations of the Laplacian gauge are defined and compared. The Laplacian gauge propagator has already been claimed to be insensitive to finite volume effects and this is tested on lattices with large volumes.

preprint2002arXiv

Lattice quark propagator with staggered quarks in Landau and Laplacian gauges

We report on the lattice quark propagator using standard and improved Staggered quark actions, with the standard, Wilson gauge action. The standard Kogut-Susskind action has errors of \oa{2} while the ``Asqtad'' action has \oa{4}, \oag{2}{2} errors. The quark propagator is interesting for studying the phenomenon of dynamical chiral symmetry breaking and as a test-bed for improvement. Gauge dependent quantities from lattice simulations may be affected by Gribov copies. We explore this by studying the quark propagator in both Landau and Laplacian gauges. Landau and Laplacian gauges are found to produce very similar results for the quark propagator.

preprint2002arXiv

Modelling the quark propagator

The quark propagator is at the core of lattice hadron spectrum calculations as well as studies in other nonperturbative schemes. We investigate the quark propagator with an improved staggered action (Asqtad) and an improved gluon action, which provides good quality data down to small quark masses. This is used to construct ansätze suitable for model hadron calculations as well as adding to our intuitive understanding of QCD.

preprint2002arXiv

Moments of Nucleon Light Cone Quark Distributions Calculated in Full Lattice QCD

Moments of the quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD. Calculations of proton matrix elements of operators corresponding to these moments through the operator product expansion have been performed on $16^3 \times 32$ lattices for Wilson fermions at $β= 5.6$ using configurations from the SESAM collaboration and at $β= 5.5$ using configurations from SCRI. One-loop perturbative renormalization corrections are included. At quark masses accessible in present calculations, there is no statistically significant difference between quenched and full QCD results, indicating that the contributions of quark-antiquark excitations from the Dirac Sea are small. Close agreement between calculations with cooled configurations containing essentially only instantons and the full gluon configurations indicates that quark zero modes associated with instantons play a dominant role. Naive linear extrapolation of the full QCD calculation to the physical pion mass yields results inconsistent with experiment. Extrapolation to the chiral limit including the physics of the pion cloud can resolve this discrepancy and the requirements for a definitive chiral extrapolation are described.

preprint2001arXiv

Lattice Quark Propagator in Landau and Laplacian Gauges

We present results for the lattice quark propagator in both Landau and Laplacian gauges using standard and improved staggered quark actions. The standard Kogut-Susskind action has errors of \oa{2} while the improved ``Asqtad'' action has \oa{4}, \oag{2}{2} errors. This improvement is seen in the quark propagator. We demonstrate the application of tree-level corrections to these actions and see that Landau and Laplacian gauges produce very similar results.

preprint2001arXiv

Quark propagator from an improved staggered action in Laplacian and Landau gauges

Studies of gauge dependent quantities are afflicted with Gribov copies, but Laplacian gauge fixing provides one possible solution to this problem. We present results for the lattice quark propagator in both Landau and Laplacian gauges using standard and improved staggered quark actions. The standard Kogut-Susskind action has errors of \oa{2} while the improved ``Asqtad'' action has \oa{4}, \oag{2}{2} errors and this improvement is seen in the quark propagator. We demonstrate the application of tree-level corrections to these actions and see that Landau and Laplacian gauges produce very similar results. In addition, we test an ansatz for the quark mass function, with promising results. In the chiral limit, the infrared quark mass, $M(q^2 = 0)$ is found to be $260\pm 20$ MeV.

preprint2000arXiv

Domain Wall Fermions with Exact Chiral Symmetry

We show how the standard domain wall action can be simply modified to allow arbitrarily exact chiral symmetry at finite fifth dimensional extent. We note that the method can be used for both quenched and dynamical calculations. We test the method using smooth and thermalized gauge field configurations. We also make comparisons of the performance (cost) of the domain wall operator for spectroscopy compared to other methods such as the overlap-Dirac operator and find both methods are comparable in cost.

preprint1999arXiv

Improved Staggered Fermion Actions for QCD Thermodynamics

We analyze the cut-off dependence of the fermion contribution to the finite temperature free energy density in ${\cal O}(g^2)$ lattice perturbation theory for several improved staggered fermion actions. Cut-off effects are drastically reduced for the Naik action and an action with improved rotational symmetry of the quark propagator. We show that improvement of rotational symmetry at ${\cal O}(g^2)$ further reduces cut-off effects in thermodynamic observables. We also show that the introduction of fat-links does not have a significant influence on cut-off distortions at ${\cal O}(g^2)$.

preprint1999arXiv

Looking for Effects of Topology in the Dirac Spectrum of Staggered Fermions

We classify SU(3) gauge field configurations in different topological sectors by the smearing technique. In each sector we compute the distribution of low lying eigenvalues of the staggered Dirac operator. In all sectors we find perfect agreement with the predictions for the sector of topological charge zero. The smallest Dirac operator eigenvalues of staggered fermions at presently realistic lattice couplings are thus insensitive to gauge field topology. On the smeared configurations, $4ν$ eigenvalues go to zero in agreement with the index theorem.

preprint1999arXiv

Sea quark effects in B Spectroscopy and Decay Constants

We present comprehensive results for the spectrum and decay constants of hadrons containing a single b quark. The heavy quark is simulated using an $O(1/M)$ NRQCD action and the light quark using the $O(a)$ tadpole-improved clover action on gauge configurations containing two degenerate flavours of sea quarks at $β^{n_f=2}=5.6$ provided by the HEMCGC collaboration. We present detailed results for the lower lying $S$ and $P$ wave $B$ meson states and the $Λ_b$ baryon. We find broad agreement with experiment. In addition, we present results for the pseudoscalar and, for the first time, the vector decay constants fully consistent to $O(α/M): f_B = 186(5)(stat)(19)(pert)(9)(disc)(13)(NRQCD)(+50)(a^{-1})MeV, f_B^* = 181(6) (stat)(18)(pert)(9)(disc)(13)(NRQCD)(+55)(a^{-1})MeV$ and $f_{B_s}/f_B = 1.14(2)(stat)(-2)(κ_s)$. We present an investigation of sea quark effects in the $B$ spectrum and decay constants. We compare our results with those from similar quenched simulations at $β^{n_f=0}=6.0$. For the spectrum, the quenched results reproduce the experimental spectrum well and there is no significant difference between the quenched and $n_f=2$ results. For the decay constants, our results suggest that sea quark effects may be large. We find that $f_B$ increases by approximately 25% between $n_f=0$ and $n_f=2$.

preprint1996arXiv

$B$ Decay Constants from NRQCD with Dynamical Fermions

We present a lattice investigation of the heavy-light meson decay constants using Wilson light quarks and NRQCD heavy quarks, partially including the effects of dynamical sea quarks. We calculate the pseudoscalar and vector decay constants over a wide range in heavy quark mass and are able to perform a detailed analysis of heavy quark symmetry. We find consistency between the extrapolation of the NRQCD results and the static case, as expected. We find the slope of the decay constants with $1/M$ is significantly larger than naive expectations and the results of previous lattice calculations. For the first time we extract the non-perturbative coefficients of the slope arising from the $O(1/M)$ heavy quark interactions separately and show the kinetic energy of the heavy quark is dominant and responsible for the large slope. In addition, we find that significant systematic errors remain in the decay constant extracted around the $B$ meson mass due to truncating the NRQCD series at $O(1/M)$. We estimate the higher order contributions to $f_B$ are approximately $20\%$; roughly the same size as the systematic errors introduced by using the Wilson action for light quarks.

preprint1996arXiv

$B$ Spectroscopy from NRQCD with Dynamical Fermions

We present a lattice investigation, partially including the effects of dynamical quarks, of the heavy-light mesons using NRQCD for the heavy quark and the Wilson action for the light quark. We performed an extensive calculation of the spectrum employing a multi-state, multi-exponential fitting analysis which enabled us to extract the $2S{-}1S$ splitting as well as the $^1P_1{-}\bar{S}$ and hyperfine splittings. The Wilson action introduces a large systematic error into the calculation, and within this uncertainty we obtain agreement with experiment. We performed a comprehensive calculation of the heavy-light meson mass, investigating three methods and their range of validity. The agreement found between these methods confirms that Lorentz invariance can be restored at this order in NRQCD by a constant shift to all masses. We calculated spectroscopic quantities over a wide range in heavy quark mass and were able to perform a detailed investigation of heavy quark symmetry around the $b$ quark mass. In particular, we extracted the nonperturbative coefficients of terms in the heavy quark expansion of the meson binding energy. We demonstrate the importance of using tadpole-improved operators in such a calculation.

preprint1996arXiv

Thermodynamics of lattice QCD with two light quark flavours on a 16^3 x 8 lattice II

We have extended our earlier simulations of the high temperature behaviour of lattice QCD with two light flavours of staggered quarks on a $16^3 \times 8$ lattice to lower quark mass (m_q=0.00625). The transition from hadronic matter to a quark-gluon plasma is observed at $6/g^2=5.49(2)$ corresponding to a temperature of $T_c \approx 140$MeV. We present measurements of observables which probe the nature of the quark-gluon plasma and serve to distinguish it from hadronic matter. Although the transition is quite abrupt, we have seen no indications that it is first order.

preprint1995arXiv

From Spectroscopy to the Strong Coupling Constant with Heavy Wilson Quarks

In this work we present lattice calculations of the masses of P-wave mesons using Monte Carlo simulations. Our valence fermions are defined by the Wilson action. Our gauge fields are generated with both dynamical staggered fermions at a lattice coupling $β\equiv 6/g^2=5.6$ for sea quark masses of $am_q=0.010$ and 0.025, and in the quenched approximation at $β=6.0$. We present results for charm and charmonium spectroscopy and use them to compute the strong coupling constant $α_s$. We compare our results to those of other recent lattice calculations and experiments.

preprint1994arXiv

The heavy quark potential in QCD with 2 flavors of dynamical quarks

We compute the heavy quark potential on configurations generated by the HEMCGC collaboration with dynamical staggered fermions at $6/g^2 = 5.6$ and with dynamical Wilson fermions at $6/g^2 = 5.3$. The computations are done on $16^3 \times 32$ lattices, corresponding to physical sizes of about 1.6 and 2.3 fm, respectively. Up to the distances probed no sign of string breaking is detectable. We also compute the recently proposed scale $r_0$ defined by $r_0^2 F(r_0) = 1.65$.

preprint1993arXiv

Computation of the Spatial String Tension in High Temperature SU(2) Gauge Theory

A detailed investigation of the temperature dependence of the spatial string tension $σ_s$ in $SU(2)$ gauge theory is presented. A sustained performance of 3~GFLOPS on a 64K Connection Machine CM-2 equivalent has been achieved. Scaling of $σ_s$ between $β=2.5115$ and $β=2.74$, on large lattices, is demonstrated. Below the critical temperature, $T_c$, $σ_s$ remains constant. For temperatures larger than $2T_c$ the temperature dependence can be parametrized by $σ_s(T) = (0.369\pm 0.014)^2 g^4(T)T^2$, where $g(T)$ is a 2-loop running coupling constant with the scale parameter determined as $Λ_T = (0.076\pm 0.013)T_c$.