Researcher profile

C. Bernard

C. Bernard contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
59works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

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

K semileptonic form factor with HISQ fermions at the physical point

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

preprint2014arXiv

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

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

preprint2014arXiv

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

Matrix Elements for $D$- and $B$-Mixing from 2+1 Flavor Lattice QCD

We present the status of our calculation of hadronic matrix elements for $D$- and $B$-meson mixing. We use a large set of the MILC collaboration's $N_f=2+1$ asqtad ensembles, which includes lattice spacings in the range $a\approx0.12$-0.045 fm, and up/down to strange quark mass ratios as low as 0.05. The asqtad action is also employed for the light valence quarks. For the heavy quarks we use the Sheikholeslami-Wohlert action with the Fermilab interpretation. Our calculation covers the complete set of five local operators needed to describe $B$-meson mixing in the Standard Model and Beyond. In the charm sector, our calculation of local mixing matrix elements may be used to constrain new physics models. We present final correlator fit results on the full data set for the $B$-meson mixing project and preliminary fit results for the $D$-meson mixing project.

preprint2013arXiv

Neutral B Mixing in Staggered Chiral Perturbation Theory

I calculate, at one loop in staggered chiral perturbation theory, the matrix elements of the complete set of five local operators that may contribute to B mixing both in the Standard Model and in beyond-the-Standard-Model theories. Lattice computations of these matrix elements by the Fermilab Lattice/MILC collaborations (and earlier by the HPQCD collaboration) convert a light staggered quark into a naive quark, and construct the relevant 4-quark operators as local products of two local bilinears, each involving the naive light quark and the heavy quark. This particular representation of the operators turns out to be important in the chiral calculation, and it results in the presence of "wrong-spin" operators, whose contributions however vanish in the continuum limit. If the matrix elements of all five operators are computed on the lattice, then no additional low energy constants are required to describe wrong-spin chiral effects.

preprint2013arXiv

Semileptonic Kaon Decay in Staggered Chiral Perturbation Theory

The determination of $\vert V_{us}\vert$ from kaon semileptonic decays requires the value of the form factor $f_+(q^2=0)$, which can be calculated precisely on the lattice. We provide the one-loop partially quenched staggered chiral perturbation theory expressions that may be employed to analyze staggered simulations of $f_+(q^2)$ with three light flavors. We consider both the case of a mixed action, where the valence and sea sectors have different staggered actions, and the standard case where these actions are the same. The momentum transfer $q^2$ of the form factor is allowed to have an arbitrary value. We give results for the generic situation where the $u$, $d$, and $s$ quark masses are all different, $N_f=1+1+1$, and for the isospin limit, $N_f=2+1$. The expression we obtain for $f_+(q^2)$ is independent of the mass of the (valence) spectator quark. In the limit of vanishing lattice spacing, our results reduce to the one-loop continuum partially quenched expression for $f_+(q^2)$, which has not previously been reported in the literature for the $N_f=1+1+1$ case. Our expressions have already been used in staggered lattice analyses of $f_+(0)$, and should prove useful in future calculations as well.

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

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$ mixing from 2+1 flavor lattice-QCD: the Standard Model and beyond

We report on the status of our lattice-QCD calculation of the hadronic contribution to $B_d^0$ and $B^0_s$ mixing, with 2+1 flavors of dynamical sea quarks. Preliminary results for hadronic mixing matrix elements are given for a basis of five four-quark, dimension-six, $ΔB=2$ mixing operators that spans the space of all possible hadronic mixing contributions in the Standard Model and beyond. At the intermediate stage of analysis reported on in this work, our errors are competitive with published Standard Model matrix element results. For beyond the Standard Model matrix elements, this is the first unquenched calculation and the first new lattice-QCD calculation in ten years.

preprint2012arXiv

Neutral B mixing from 2+1 flavor lattice QCD

We present an update of the Fermilab-MILC Collaboration&#39;s calculation of hadronic matrix elements for B^0-\bar{B^0} mixing. This work is a more extended analysis than our recent publication of the SU(3)-breaking ratio xi [arXiv:1205.7013]. We use the asqtad staggered action for light valence quarks in combination with the Fermilab interpretation of the Sheikoleslami-Wohlert action for heavy quarks. The calculations use MILC&#39;s 2+1 flavor asqtad ensembles. Ensembles include four lattice spacings from approximately 0.125 fm to 0.045 fm and up/down to strange quark mass ratios as low as 0.05. Our calculation covers the complete set of five operators needed to describe B mixing in the Standard Model and beyond. In addition to an update including a fuller set of analyzed data, we comment on the form of the staggered ChPT extrapolation function.

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&#39;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

Staggered Chiral Perturbation Theory for Neutral B Mixing

I describe a calculation of B meson mixing at one-loop in staggered chiral perturbation theory, for the complete set of Standard Model and beyond-the-Standard Model operators. The particular lattice representation of the continuum operators used by the Fermilab Lattice/MILC collaborations (and earlier by the HPQCD collaboration) turns out to be important, and results in the presence of &#34;wrong-spin&#34; operators, whose contributions however vanish in the continuum limit. The relation between staggered and naive fermions also plays a key role.

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.

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

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

Electromagnetic Corrections in Staggered Chiral Perturbation Theory

To reduce errors in light-quark mass determinations, it is now necessary to consider electromagnetic contributions to light-meson masses. Calculations using staggered quarks and quenched photons are currently underway. Suitably-extended chiral perturbation theory is necessary to extrapolate the lattice data to the physical limit. Here we give (preliminary) results for light-meson masses using staggered chiral perturbation theory including electromagnetism, and discuss the extent to which quenched-photon simulations can improve quark-mass calculations.

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 &#34;unmixed s bar s&#34; 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

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.

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&#39;&#39; fermions. These simulations were carried out with two degenerate flavors of light quarks (up and down) and with one heavier flavor, the strange quark. Several light quark masses, down to about 3 times the physical light quark mass, and six lattice spacings have been used. These enable controlled continuum and chiral extrapolations of many low energy QCD observables. We review the improved staggered formalism, emphasizing both advantages and drawbacks. In particular, we review the procedure for removing unwanted staggered species in the continuum limit. We then describe the asqtad lattice ensembles created by the MILC Collaboration. All MILC lattice ensembles are publicly available, and they have been used extensively by a number of lattice gauge theory groups. We review physics results obtained with them, and discuss the impact of these results on phenomenology. Topics include the heavy quark potential, spectrum of light hadrons, quark masses, decay constant of light and heavy-light pseudoscalar mesons, semileptonic form factors, nucleon structure, scattering lengths and more. We conclude with a brief look at highly promising future prospects.

preprint2009arXiv

HISQ action in dynamical simulations

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

preprint2009arXiv

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

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

preprint2009arXiv

Progress on four flavor QCD with the HISQ action

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

preprint2009arXiv

Results from the MILC collaboration&#39;s SU(3) chiral perturbation theory analysis

We present the status of the MILC collaboration&#39;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&#39;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&#39;s theorem. Hadron masses are calculated on MILC N_f=2+1 QCD ensembles at lattice spacing \approx 0.15fm, together with quenched non-compact U(1) configurations. The O(a^2) tadpole improved staggered quark (asqtad) action is used both for the sea quarks and for six different valence quark masses. Chiral extrapolations are performed using partially quenched chiral perturbation theory incorporating electromagnetic corrections.

preprint2007arXiv

Baryon masses with improved staggered quarks

The MILC collaboration&#39;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

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.

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.