Source author record

Steven Gottlieb

Steven Gottlieb 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

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

100 published item(s)

preprint2023arXiv

Two-link Staggered Quark Smearing in QUDA

Gauge covariant smearing based on the 3D lattice Laplacian can be used to create extended operators that have better overlap with hadronic ground states. For staggered quarks, we make use of two-link parallel transport to preserve taste properties. We have implemented the procedure in QUDA. We present the performance of this code on the NVIDIA A100 GPUs in Indiana University's Big Red 200 supercomputer and on the AMD MI250X GPUs in Oak Ridge Leadership Computer Facility's (OLCF's) Crusher and discuss its scalability. We also study the performance improvement from using NVSHMEM on OLCF's Summit. Reusing precomputed two-link products for all sources and sinks, it reduces the total smearing time for a baryon correlator measurement by a factor of 100-120 as compared with the original MILC code and reduces the overall time by 60-70%.

preprint2022arXiv

A lattice QCD perspective on weak decays of b and c quarks Snowmass 2022 White Paper

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics.

preprint2022arXiv

Lattice QCD and the Computational Frontier

The search for new physics requires a joint experimental and theoretical effort. Lattice QCD is already an essential tool for obtaining precise model-free theoretical predictions of the hadronic processes underlying many key experimental searches, such as those involving heavy flavor physics, the anomalous magnetic moment of the muon, nucleon-neutrino scattering, and rare, second-order electroweak processes. As experimental measurements become more precise over the next decade, lattice QCD will play an increasing role in providing the needed matching theoretical precision. Achieving the needed precision requires simulations with lattices with substantially increased resolution. As we push to finer lattice spacing we encounter an array of new challenges. They include algorithmic and software-engineering challenges, challenges in computer technology and design, and challenges in maintaining the necessary human resources. In this white paper we describe those challenges and discuss ways they are being dealt with. Overcoming them is key to supporting the community effort required to deliver the needed theoretical support for experiments in the coming decade.

preprint2020arXiv

Computing Nucleon Charges with Highly Improved Staggered Quarks

This work continues our program of lattice-QCD baryon physics using staggered fermions for both the sea and valence quarks. We present a proof-of-concept study that demonstrates, for the first time, how to calculate baryon matrix elements using staggered quarks for the valence sector. We show how to relate the representations of the continuum staggered flavor-taste group $\text{SU}(8)_{FT}$ to those of the discrete lattice symmetry group. The resulting calculations yield the normalization factors relating staggered baryon matrix elements to their physical counterparts. We verify this methodology by calculating the isovector vector and axial-vector charges $g_V$ and $g_A$. We use a single ensemble from the MILC Collaboration with 2+1+1 flavors of sea quark, lattice spacing $a\approx 0.12$ fm, and a pion mass $M_π\approx305$ MeV. On this ensemble, we find results consistent with expectations from current conservation and neutron beta decay. Thus, this work demonstrates how highly-improved staggered quarks can be used for precision calculations of baryon properties, and, in particular, the isovector nucleon charges.

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.

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.

preprint2020arXiv

Lattice QCD Impact on Determination of CKM Matrix: Status and Prospects

Lattice QCD is an important tool for theoretical input for flavor physics. There have been four reviews by the Flavour Lattice Averaging Group (FLAG). This talk will review the current status of the magnitude of eight of the nine CKM matrix elements, borrowing heavily from the most recent FLAG review (co-authored by the speaker). Future prospects for improving the determination of the CKM matrix will be discussed.

preprint2019arXiv

$B$-meson semileptonic form factors on (2+1+1)-flavor HISQ ensembles

We report updates to an ongoing lattice-QCD calculation of the form factors for the semileptonic decays $B \to π\ell ν$, $B_s \to K \ell ν$, $B \to π\ell^+ \ell^-$, and $B \to K \ell^+ \ell^-$. The tree-level decays $B_{(s)} \to π(K) \ell ν$ enable precise determinations of the CKM matrix element $|V_{ub}|$, while the flavor-changing neutral-current interactions $B \to π(K) \ell^+ \ell^-$ are sensitive to contributions from new physics. This work uses MILC's (2+1+1)-flavor HISQ ensembles at approximate lattice spacings between $0.057$ and $0.15$ fm, with physical sea-quark masses on four out of the seven ensembles. The valence sector is comprised of a clover $b$ quark (in the Fermilab interpretation) and HISQ light and $s$ quarks. We present preliminary results for the form factors $f_0$, $f_+$, and $f_T$, including studies of systematic errors.

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

Lattice field theory applications in high energy physics

Lattice gauge theory was formulated by Kenneth Wilson in 1974. In the ensuing decades, improvements in actions, algorithms, and computers have enabled tremendous progress in QCD, to the point where lattice calculations can yield sub-percent level precision for some quantities. Beyond QCD, lattice methods are being used to explore possible beyond the standard model (BSM) theories of dynamical symmetry breaking and supersymmetry. We survey progress in extracting information about the parameters of the standard model by confronting lattice calculations with experimental results and searching for evidence of BSM effects.

preprint2016arXiv

MILC staggered conjugate gradient performance on Intel KNL

We review our work done to optimize the staggered conjugate gradient (CG) algorithm in the MILC code for use with the Intel Knights Landing (KNL) architecture. KNL is the second gener- ation Intel Xeon Phi processor. It is capable of massive thread parallelism, data parallelism, and high on-board memory bandwidth and is being adopted in supercomputing centers for scientific research. The CG solver consumes the majority of time in production running, so we have spent most of our effort on it. We compare performance of an MPI+OpenMP baseline version of the MILC code with a version incorporating the QPhiX staggered CG solver, for both one-node and multi-node runs.

preprint2016arXiv

Phenomenology of semileptonic B-meson decays with form factors from lattice QCD

We study the exclusive semileptonic $B$-meson decays $B\to K(π)\ell^+\ell^-$, $B\to K(π)ν\barν$, and $B\toπτν$, computing observables in the Standard model using the recent lattice-QCD results for the underlying form factors from the Fermilab Lattice and MILC Collaborations. These processes provide theoretically clean windows into physics beyond the Standard Model because the hadronic uncertainties are now under good control for suitably binned observables. For example, the resulting partially integrated branching fractions for $B\toπμ^+μ^-$ and $B\to Kμ^+μ^-$ outside the charmonium resonance region are 1-2$σ$ higher than the LHCb Collaboration's recent measurements, where the theoretical and experimental errors are commensurate. The combined tension is 1.7$σ$. Combining the Standard-Model rates with LHCb's measurements yields values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{td}|=7.45{(69)}\times10^{-3}$, $|V_{ts}|=35.7(1.5)\times10^{-3}$, and $|V_{td}/V_{ts}|=0.201{(20)}$, which are compatible with the values obtained from neutral $B_{(s)}$-meson oscillations and have competitive uncertainties. Alternatively, taking the CKM matrix elements from unitarity, we constrain new-physics contributions at the electroweak scale. The constraints on the Wilson coefficients ${\rm Re}(C_9)$ and ${\rm Re}(C_{10})$ from $B\toπμ^+μ^-$ and $B\to Kμ^+μ^-$ are competitive with those from $B\to K^* μ^+μ^-$, and display a 2.0$σ$ tension with the Standard Model. Our predictions for $B\to K(π)ν\barν$ and $B\toπτν$ are close to the current experimental limits.

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

Dynamical QCD+QED simulation with staggered quarks

Electromagnetic effects play an important role in many phenomena such as isospin-symmetry breaking in the hadron spectrum and the hadronic contributions to g-2. We have generalized the MILC QCD code to include the electromagnetic field. In this work, we focus on simulations including charged sea quarks using the RHMC algorithm. We show details of the dynamical QCD+QED simulation algorithm with compact QED. We analyze the code performance and results for hadron-spectrum observables.

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

Staggered Dslash Performance on Intel Xeon Phi Architecture

The conjugate gradient (CG) algorithm is among the most essential and time consuming parts of lattice calculations with staggered quarks. We test the performance of CG and dslash, the key step in the CG algorithm, on the Intel Xeon Phi, also known as the Many Integrated Core (MIC) architecture. We try different parallelization strategies using MPI, OpenMP, and the vector processing units (VPUs).

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

$B\toπ\ellν$ and $B\toπ\ell^+\ell^-$ semileptonic form factors from unquenched lattice QCD

We update the lattice calculation of the $B\toπ$ semileptonic form factors, which have important applications to the CKM matrix element $|V_{ub}|$ and the $B\toπ\ell^+\ell^-$ rare decay. We use MILC asqtad ensembles with $N_f=2+1$ sea quarks and over a range of lattice spacings $a \approx 0.045$--$0.12$ fm. We perform a combined chiral and continuum extrapolation of our lattice data using SU(2) staggered chiral perturbation theory in the hard pion limit. To extend the results for the form factors to the full kinematic range, we take a functional approach to parameterize the form factors using the Bourrely-Caprini-Lellouch formalism in a model-independent way. Our analysis is still blinded with an unknown off-set factor which will be disclosed when we present the final results.

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

Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model

We calculate the quadratic fluctuations of net baryon number, electric charge and strangeness as well as correlations among these conserved charges in (2+1)-flavor lattice QCD at zero chemical potential. Results are obtained using calculations with tree level improved gauge and the highly improved staggered quark (HISQ) actions with almost physical light and strange quark masses at three different values of the lattice cut-off. Our choice of parameters corresponds to a value of 160 MeV for the lightest pseudo scalar Goldstone mass and a physical value of the kaon mass. The three diagonal charge susceptibilities and the correlations among conserved charges have been extrapolated to the continuum limit in the temperature interval 150 MeV <T < 250 MeV. We compare our results with the hadron resonance gas (HRG) model calculations and find agreement with HRG model results only for temperatures T<= 150 MeV. We observe significant deviations in the temperature range 160 MeV < T < 170 MeV and qualitative differences in the behavior of the three conserved charge sectors. At T < 160 MeV quadratic net baryon number fluctuations in QCD agree with HRG model calculations while, the net electric charge fluctuations in QCD are about 10% smaller and net strangeness fluctuations are about 20% larger. These findings are relevant to the discussion of freeze-out conditions in relativistic heavy ion collisions.

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

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.

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

Form factors for $B$ to $Kll$ semileptonic decay from three-flavor lattice QCD

We study the $B \to Kl^+l^-$ semileptonic decay process in three-flavor lattice QCD. We analyze several ensembles generated by the MILC collaboration at different lattice spacings and sea-quark masses. We use the asqtad improved staggered action for the light quarks and the clover action with the Fermilab interpretation for the heavy $b$ quark. We present preliminary results for the vector current induced form factors for a range of kaon energies. Our analysis includes chiral and continuum extrapolations based on SU(2) staggered χPT.

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

Results for light pseudoscalar mesons

We present the current status of the MILC collaboration's calculations of the properties of the light pseudoscalar meson sector. We use asqtad staggered ensembles with 2+1 dynamical flavors down to $a \approx 0.045$ fm and light quark mass down to 0.05 $m_s$. Here we describe fits to the data using chiral forms from SU(3) chiral perturbation theory, including all staggered taste violations at NLO and the continuum NNLO chiral logarithms. We emphasize issues of convergence of the chiral expansion.

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

Cell processor implementation of a MILC lattice QCD application

We present results of the implementation of one MILC lattice QCD application-simulation with dynamical clover fermions using the hybrid-molecular dynamics R algorithm-on the Cell Broadband Engine processor. Fifty-four individual computational kernels responsible for 98.8% of the overall execution time were ported to the Cell's Synergistic Processing Elements (SPEs). The remaining application framework, including MPI-based distributed code execution, was left to the Cell's PowerPC processor. We observe that we only infrequently achieve more than 10 GFLOPS with any of the kernels, which is just over 4% of the Cell's peak performance. At the same time, many of the kernels are sustaining a bandwidth close to 20 GB/s, which is 78% of the Cell's peak. This indicates that the application performance is limited by the bandwidth between the main memory and the SPEs. In spite of this limitation, speedups of 8.7x (for 8x8x16x16 lattice) and 9.6x (for 16x16x16x16 lattice) were achieved when comparing a 3.2 GHz Cell processor to a single core of a 2.33 GHz Intel Xeon processor. When comparing the code scaled up to execute on a dual-Cell blade and a quad-core dual-chip Intel Xeon blade, the speedups are 1.5x (8x8x16x16 lattice) and 4.1x (16x16x16x16 lattice).

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

MILC results for light pseudoscalars

We present the latest preliminary results of the MILC collaboration's analysis of the light pseudoscalar meson sector. The analysis includes data from new ensembles with smaller lattice spacings, smaller light quark masses and lighter-than-physical strange quark masses. Both SU(2) and SU(3) chiral fits, including NNLO chiral logarithms, are shown. We give results for decay constants, quark masses, Gasser-Leutwyler low energy constants, and condensates in the two- and three-flavor chiral limits.

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

Quarkonium mass splittings in three-flavor lattice QCD

We report on calculations of the charmonium and bottomonium spectrum in lattice QCD. We use ensembles of gauge fields with three flavors of sea quarks, simulated with the asqtad improved action for staggered fermions. For the heavy quarks we employ the Fermilab interpretation of the clover action for Wilson fermions. These calculations provide a test of lattice QCD, including the theory of discretization errors for heavy quarks. We provide, therefore, a careful discussion of the results in light of the heavy-quark effective Lagrangian. By and large, we find that the computed results are in agreement with experiment, once parametric and discretization errors are taken into account.

preprint2009arXiv

Quarkonium mass splittings with Fermilab heavy quarks and 2+1 flavors of improved staggered sea quarks

We present results from an ongoing lattice study of the lowest lying charmonium and bottomonium level splittings using the Fermilab heavy quark formalism. Our objective is to test the performance of this action on MILC-collaboration ensembles of (2+1) flavors of light improved staggered (asqtad) quarks. Measurements are done on 16 ensembles with degenerate up and down quarks of various masses, thus permitting a chiral extrapolation, and over lattice spacings ranging from 0.09 fm to 0.18 fm, thus permitting study of lattice-spacing dependence. We examine combinations of the mass splittings that are sensitive to components of the effective quarkonium potential.

preprint2009arXiv

Results from the MILC collaboration's SU(3) chiral perturbation theory analysis

We present the status of the MILC collaboration's analysis of the light pseudoscalar meson sector with SU(3) chiral fits. The analysis includes data from new ensembles with smaller lattice spacing, smaller light quark masses and lighter than physical strange quark masses. Our fits include the NNLO chiral logarithms. We present results for decay constants, quark masses, Gasser-Leutwyler low energy constants, and condensates in the two- and three-flavor chiral limits.

preprint2009arXiv

SU(2) chiral fits to light pseudoscalar masses and decay constants

We present the results of fits to recent asqtad data in the light pseudoscalar sector using SU(2) partially-quenched staggered chiral perturbation theory. Superfine (a~0.06 fm) and ultrafine (a~0.045 fm) ensembles are used, where light sea quark masses and taste splittings are small compared to the strange quark mass. Our fits include continuum NNLO chiral logarithms and analytic terms. We give preliminary results for the pion decay constant, SU(2) low-energy constants and the chiral condensate in the two-flavor chiral limit.

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.

preprint2008arXiv

Heavy baryon mass spectrum from lattice QCD with 2+1 dynamical sea quark flavors

We investigate the mass differences of charm and bottom heavy baryons, using MILC lattice gauge configurations with 2+1 dynamical sea quark flavors. We extend our previous work to three lattice spacings: fine (a~0.09), coarse (a~0.12), and medium-coarse (a~0.15) ensembles. For extrapolations and interpolations, we apply simultaneous quadratic fits and simple linear fits with the full QCD data points for the valence quark masses as well as the sea quark masses.

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

Onium Masses with Three Flavors of Dynamical Quarks

We have greatly extended an earlier calculation of the charmonium spectrum on three flavor dynamical quark ensembles by using more recent ensembles generated by the MILC collaboration. The heavy quarks are treated using the Fermilab formulation. The charmonium state masses are in reasonable agreement with the observed spectrum; however, some of the spin splittings may still be too small.

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

First determination of the strange and light quark masses from full lattice QCD

We compute the strange quark mass $m_s$ and the average of the $u$ and $d$ quark masses $\hat m$ using full lattice QCD with three dynamical quarks combined with experimental values for the pion and kaon masses. The simulations have degenerate $u$ and $d$ quarks with masses $m_u=m_d\equiv \hat m$ as low as $m_s/8$, and two different values of the lattice spacing. The bare lattice quark masses obtained are converted to the $\msbar$ scheme using perturbation theory at $O(alpha_s)$. Our results are: $m_s^\msbar$(2 GeV) = 76(0)(3)(7)(0) MeV, $\hat m^\msbar$(2 GeV) = 2.8(0)(1)(3)(0) MeV and $m_s/\hat m$ = 27.4(1)(4)(0)(1), where the errors are from statistics, simulation, perturbation theory, and electromagnetic effects, respectively.

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.

preprint2004arXiv

The scaling dimension of low lying Dirac eigenmodes and of the topological charge density

As a quantitative measure of localization, the inverse participation ratio of low lying Dirac eigenmodes and topological charge density is calculated on quenched lattices over a wide range of lattice spacings and volumes. Since different topological objects (instantons, vortices, monopoles, and artifacts) have different co-dimension, scaling analysis provides information on the amount of each present and their correlation with the localization of low lying eigenmodes.

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.

preprint2001arXiv

Benchmarking and tuning the MILC code on clusters and supercomputers

Recently, we have benchmarked and tuned the MILC code on a number of architectures including Intel Itanium and Pentium IV (PIV), dual-CPU Athlon, and the latest Compaq Alpha nodes. Results will be presented for many of these, and we shall discuss some simple code changes that can result in a very dramatic speedup of the KS conjugate gradient on processors with more advanced memory systems such as PIV, IBM SP and Alpha.

preprint1999arXiv

Scale Determination Using the Static Potential with Two Dynamical Quark Flavors

We study the static potential using gauge configurations that include the effects of two flavors of dynamical Kogut-Susskind quarks. The configurations, generated by the MILC collaboration, and used to study the spectrum and heavy-light decay constants, cover a range $5.3 \le 6/g^2 \le 5.6$. There are at least four quark masses for each coupling studied. Determination of $r_0$ from the potential can be used to set a scale. This alternative scale is useful to study systematic errors on the spectrum and decay constants.

preprint1997arXiv

QCD Thermodynamics with an Improved Lattice Action

We have investigated QCD with two flavors of degenerate fermions using a Symanzik-improved lattice action for both the gauge and fermion actions. Our study focuses on the deconfinement transition on an $N_t=4$ lattice. Having located the thermal transition, we performed zero temperature simulations nearby in order to compute hadronic masses and the static quark potential. We find that the present action reduces lattice artifacts present in thermodynamics with the standard Wilson (gauge and fermion) actions. However, it does not bring studies with Wilson-type quarks to the same level as those using the Kogut--Susskind formulation.

preprint1995arXiv

Resonance Scattering Phase Shifts on a Non-Rest Frame Lattice

Many low energy hadrons, such as the rho, can be observed as resonances in scattering experiments. A proposal by Lüscher enables one to determine infinite volume elastic scattering phases from the two-particle energy spectrum measured from finite periodic lattices. In this work, we generalize the formalism to the case where the total momentum of the particles is non-zero; i.e. the lattice frame is not the center-of-mass frame of the scattering particles. There are several advantages to this procedure including making a wider variety of center of mass energies accessible with a fixed lattice volume, and making the avoided level crossing in a P-wave decay occur with a smaller volume. The formalism is tested with a simple lattice model of two fields with different masses and a 3-point coupling in 3+1 dimensions. We find remarkable agreement between the rest-frame and non-rest-frame scattering.

preprint1995arXiv

Status and prospect for determining $f_B$, $f_{B_s}$, $f_{B_s} / f_B$ on the lattice

Preliminary results from the MILC collaboration for $f_B$, $f_{B_s}$, $f_D$, $f_{D_s}$ and their ratios are presented. We compute in the quenched approximation at $β=6.3$, 6.0 and 5.7 with Wilson light quarks and static and Wilson heavy quarks. We attempt to quantify all systematic errors other than quenching, and have a first indication of the size of quenching errors.

preprint1992arXiv

Finite Size Effects on the QCD Spectrum Revisited

We have continued our study of finite size effects in the QCD spectrum on lattices ranging in size from $8^3\times 24$ to $16^3\times24$. We have increased our statistics for quark mass $am_q=0.025$ for the smallest lattice size. In addition, we have studied quark mass 0.0125 for lattice sizes $12^3\times24$ and $16^3\times24$. These lattice sizes correspond to a box 1.8--3.6 fm on a side when the rho mass at zero quark mass is used to set the scale. We discuss the nucleon to rho mass ratio at a smaller value of $m_π/m_ρ$ than previously studied with two dynamical flavors.