BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

L. Levkova

L. Levkova appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phhep-exnucl-th

Catalog footprint

What is connected

81works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

81 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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

0 citations0 reviewsNo code linkedOpen access
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)$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}|$.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quark-gluon plasma in an external magnetic field

Using numerical simulations of lattice QCD we calculate the effect of an external magnetic field on the equation of state of the quark-gluon plasma. The results are obtained using a Taylor expansion of the pressure with respect to the magnetic field for the first time. The coefficients of the expansion are computed to second order in the magnetic field. Our setup for the external magnetic field avoids complications arising from toroidal boundary conditions, making a Taylor series expansion straightforward. This study is exploratory and is meant to serve as a proof of principle.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charmonium mass splittings at the physical point

We present results from an ongoing study of mass splittings of the lowest lying states in the charmonium system. We use clover valence charm quarks in the Fermilab interpretation, an improved staggered (asqtad) action for sea quarks, and the one-loop, tadpole-improved gauge action for gluons. This study includes five lattice spacings, 0.15, 0.12, 0.09, 0.06, and 0.045 fm, with two sets of degenerate up- and down-quark masses for most spacings. We use an enlarged set of interpolation operators and a variational analysis that permits study of various low-lying excited states. The masses of the sea quarks and charm valence quark are adjusted to their physical values. This large set of gauge configurations allows us to extrapolate results to the continuum physical point and test the methodology.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charm annihilation effects on the hyperfine splitting in charmonium

In calculations of the hyperfine splitting in charmonium, the contributions of the disconnected diagrams are considered small and are typically ignored. We aim to estimate nonperturbatively the size of the resulting correction, which may eventually be needed in high precision calculations of the charmonium spectrum. We study this problem in the quenched and unquenched QCD cases. On dynamical ensembles the disconnected charmonium propagators contain light modes which complicate the extraction of the signal at large distances. In the fully quenched case, where there are no such light modes, the interpretation of the signal is simplified. We present results from lattices with $a\approx 0.09$ fm and $a\approx 0.06$ fm.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effects of the charm quark on the QCD equation of state

We study the effects of the addition of the charm quark on the QCD equation of state at zero and nonzero chemical potential on lattices with $N_t=6$. Our ensembles are quenched with respect to charm and the charm quark is a valence staggered quark. Along the trajectory of constant physics the ratio $m_s/m_c$ is kept constant after tuning the charm quark mass at a lattice spacing of about 0.09 fm. We find that the charm quark has a significant contribution to the equation of state at zero chemical potential already at temperatures between about $1.2T_c$ and $2T_c$. The additional contribution at nonzero chemical potential vanishes within the current statistical uncertainty.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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%.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Contributions of charm anihilation to the hyperfine splitting in charmonium

In calculations of the hyperfine splitting in charmonium, the contributions of the disconnected diagrams is considered small and is typically ignored. We aim to estimate nonperturbatively the size of the resulting error, which could potentially affect the high precision calculations of the charmonium spectrum. Following our work on the effects of the disconnected diagrams in unquenched QCD presented at Lattice 2007, we study the same problem in the quenched case. On dynamical ensembles the disconnected charmonium propagators contain light modes which complicate the extraction of the signal at large distances. In the fully quenched case, where there are no such light modes, the interpretation of the signal is simplified. We present results from lattices with $a\approx 0.09$ fm and $a\approx 0.063$ fm.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Effects of the disconnected flavor singlet corrections on the hyperfine splitting in charmonium

Experimentally the charmonium hyperfine splitting is $M_{J/ψ}-M_{η_c}=117$ MeV and current lattice results are generally below this value. The difference could be due to the effects of the disconnected flavor singlet diagrams which have not been included in these calculations. Previous attempts to determine the disconnected flavor singlet corrections have led just to rough estimates in the range of $\pm 20$ MeV. We present preliminary results for these corrections calculated on fine ($a\approx 0.09$ fm) Asqtad 2+1 flavor lattices provided by the MILC Collaboration.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Domain-Wall Fermions at Strong Coupling

The DWF formulation becomes increasingly problematic at gauge couplings for which $a^{-1}<2$ GeV, where the roughness of the gauge field leads to increased explicit chiral symmetry breaking ($\mres$). This problem becomes especially severe for sufficiently strong coupling where the underlying 4-dimensional Wilson theory is in the Aoki phase. We review our attempts to find a suitable modification of the gauge and/or the fermion action which would allow the DWF method to work reliably at stronger coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Hardware and software status of QCDOC

QCDOC is a massively parallel supercomputer whose processing nodes are based on an application-specific integrated circuit (ASIC). This ASIC was custom-designed so that crucial lattice QCD kernels achieve an overall sustained performance of 50% on machines with several 10,000 nodes. This strong scalability, together with low power consumption and a price/performance ratio of $1 per sustained MFlops, enable QCDOC to attack the most demanding lattice QCD problems. The first ASICs became available in June of 2003, and the testing performed so far has shown all systems functioning according to specification. We review the hardware and software status of QCDOC and present performance figures obtained in real hardware as well as in simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Improving Dynamical Domain-Wall Fermion Simulations

We report on studies of the chiral properties of dynamical domain wall fermions combined with the DBW2 gauge action for different gauge couplings and fermion masses. For quenched theories, the DBW2 action gives a residual chiral symmetry breaking much smaller than what was found with more traditional choices for the gauge action. Our goal is to investigate the possibilities which this and further improvements provide for the study of QCD thermodynamics and other simulations at stronger couplings.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Staggered Fermion Thermodynamics using Anisotropic Lattices

Numerical simulations of full QCD on anisotropic lattices provide a convenient way to study QCD thermodynamics with fixed physics scales and reduced lattice spacing errors. We report results from calculations with 2-flavors of dynamical fermions where all bare parameters and hence the physics scales are kept constant while the temperature is changed in small steps by varying only the number of the time slices. The results from a series of zero-temperature scale setting simulations are used to determine the Karsch coefficients and the equation of state at finite temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Anisotropic Lattices and Dynamical Fermions

We report results from full QCD calculations with two flavors of dynamical staggered fermions on anisotropic lattices. The physical anisotropy as determined from spatial and temporal masses, their corresponding dispersion relations, and spatial and temporal Wilson loops is studied as a function of the bare gauge anisotropy and the bare velocity of light appearing in the Dirac operator. The anisotropy dependence of staggered fermion flavor symmetry breaking is also examined. These results will then be applied to the study of 2-flavor QCD thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

QCDOC: A 10-teraflops scale computer for lattice QCD

The architecture of a new class of computers, optimized for lattice QCD calculations, is described. An individual node is based on a single integrated circuit containing a PowerPC 32-bit integer processor with a 1 Gflops 64-bit IEEE floating point unit, 4 Mbyte of memory, 8 Gbit/sec nearest-neighbor communications and additional control and diagnostic circuitry. The machine's name, QCDOC, derives from ``QCD On a Chip''.

0 citations0 reviewsNo code linkedOpen access