Search

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment, $(g-2$)/2, is computed in the infinite volume QED framework with lattice QCD. We report $a_\mu^\text{HLbL}=12.47(1.15)(0.99) \times 10^{-10}$ where the first error is statistical and the second systematic. The result is mainly based on the 2+1 flavor M\"obius domain wall fermion ensemble with inverse lattice spacing $a^{-1} = 1.73~\mathrm{GeV}$, lattice size $L=5.5~\mathrm{fm}$, and $m_\pi = 139~\mathrm{MeV}$, generated by the RBC-UKQCD collaborations. The leading systematic error of this result comes from the lattice discretization. This result is consistent with previous determinations., Comment: 29 pages

A study of two-pion scattering for the isospin channels, I=0 and I=2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of L=4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the I=0 and I=2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for I=0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying K→ππ decay, on physical-mass ensembles using the GEVP.

We present a lattice calculation of the $K\toππ$ matrix elements and amplitudes with both the $ΔI = 3/2$ and 1/2 channels and $\varepsilon'$, the measure of direct $CP$ violation. We use periodic boundary conditions (PBC), where the correct kinematics of $K\toππ$ can be achieved via an excited two-pion final state. To overcome the difficulty associated with the extraction of excited states, our previous work~\cite{Bai:2015nea,RBC:2020kdj} successfully employed G-parity boundary conditions, where pions are forced to have non-zero momentum enabling the $I=0$ two-pion ground state to express the on-shell kinematics of the $K\toππ$ decay. Here instead we overcome the problem using the variational method which allows us to resolve the two-pion spectrum and matrix elements up to the relevant energy where the decay amplitude is on-shell. In this paper we report an exploratory calculation of $K\toππ$ decay amplitudes and $\varepsilon'$ using PBC on a coarser lattice size of $24^3\times64$ with inverse lattice spacing $a^{-1}=1.023$~GeV and the physical pion and kaon masses. The results are promising enough to motivate us to continue our measurements on finer lattice ensembles in order to improve the precision in the near future.

We study the chirality of staggered quarks on the Dirac eigenvalue spectrum using deep learning (DL) techniques. The Kluberg-Stern method to construct staggered bilinear operators conserves continuum property such as recursion relations, uniqueness of chirality, and Ward identities, which leads to a unique and characteristic pattern (we call it "leakage pattern (LP)") in the matrix elements of the chirality operator sandwiched between two quark eigenstates of staggered Dirac operator. DL analysis gives $99.4(2)\%$ accuracy on normal gauge configurations and $0.998$ AUC (Area Under ROC Curve) for classifying non-zero mode octets in the Dirac eigenvalue spectrum. It confirms that the leakage pattern is universal on normal gauge configurations. The multi-layer perceptron (MLP) method turns out to be the best DL model for our study on the LP., Comment: 12 pages, 7 figures, Lattice 2021 proceeding

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

We calculate $K\to\pi\pi$ matrix elements using periodic boundary conditions as an independent calculation from our previous study with G-parity boundary conditions. We present our preliminary results for $K\to\pi\pi$ three-point functions and matrix elements on a $24^3, a^{-1} = 1$~GeV, $2+1$-flavor M\"obius DWF ensemble at physical pion and kaon masses generated by the RBC and UKQCD collaborations and discuss the prospect for high-precision computation of $\varepsilon'$ with periodic boundary conditions., Proceedings of the 38th International Symposium on Lattice Field Theory (Lattice 2021), July 26-30, 2021, Zoom/Gather@MIT

We present a lattice QCD calculation of the $��I=1/2$, $K\to����$ decay amplitude $A_0$ and $\varepsilon'$, the measure of direct CP-violation in $K\to����$ decay, improving our 2015 calculation of these quantities. Both calculations were performed with physical kinematics on a $32^3\times 64$ lattice with an inverse lattice spacing of $a^{-1}=1.3784(68)$ GeV. However, the current calculation includes nearly four times the statistics and numerous technical improvements allowing us to more reliably isolate the $����$ ground-state and more accurately relate the lattice operators to those defined in the Standard Model. We find ${\rm Re}(A_0)=2.99(0.32)(0.59)\times 10^{-7}$ GeV and ${\rm Im}(A_0)=-6.98(0.62)(1.44)\times 10^{-11}$ GeV, where the errors are statistical and systematic, respectively. The former agrees well with the experimental result ${\rm Re}(A_0)=3.3201(18)\times 10^{-7}$ GeV. These results for $A_0$ can be combined with our earlier lattice calculation of $A_2$ to obtain ${\rm Re}(\varepsilon'/\varepsilon)=21.7(2.6)(6.2)(5.0) \times 10^{-4}$, where the third error represents omitted isospin breaking effects, and Re$(A_0)$/Re$(A_2) = 19.9(2.3)(4.4)$. The first agrees well with the experimental result of ${\rm Re}(\varepsilon'/\varepsilon)=16.6(2.3)\times 10^{-4}$. A comparison of the second with the observed ratio Re$(A_0)/$Re$(A_2) = 22.45(6)$, demonstrates the Standard Model origin of this "$��I = 1/2$ rule" enhancement., Updated to published version. 95 pages, 12 figures

We report preliminary results for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment. Several ensembles using 2+1 flavors of M\"obius domain-wall fermions, generated by the RBC/UKQCD collaborations, are employed to take the continuum and infinite volume limits of finite volume lattice QED+QCD. We find $a_\mu^{\rm HLbL} = (7.41\pm6.33)\times 10^{-10}$, Comment: 8 pages, 19 figures, to appear in the proceedings of Rencontres de Moriond (EW)

We report nucleon mass, isovector vector and axial-vector charges, and tensor and scalar couplings, calculated using two recent 2+1-flavor dynamical domain-wall fermions lattice-QCD ensembles generated jointly by the RIKEN-BNL-Columbia and UKQCD collaborations. These ensembles were generated with Iwasaki $\times$ dislocation-suppressing-determinant-ratio gauge action at inverse lattice spacing of 1.378(7) GeV and pion mass values of 249.4(3) and 172.3(3) MeV. The nucleon mass extrapolates to a value $m_N = 0.950(5)$ GeV at physical point. The isovector vector charge renormalizes to unity in the chiral limit, narrowly constraining excited-state contamination in the calculation. The ratio of the isovector axial-vector to vector charges shows a deficit of about ten percent. The tensor coupling no longer depends on mass and extrapolates to 1.04(5) in $\overline {\rm MS}$ 2-GeV renormalization at physical point, in a good agreement with the value obtained at the lightest mass in our previous calculations and other calculations that followed. The scalar charge, though noisier, does not show mass dependence and is in agreement with other calculations., Comment: Accepted for publication in Physical Review D, 11 pages, 13 figures and 7 tables. arXiv admin note: text overlap with arXiv:1710.06656

We present the first calculation of the kaon semileptonic form factor with sea and valence quark masses tuned to their physical values in the continuum limit of 2+1 flavour domain wall lattice QCD. We analyse a comprehensive set of simulations at the phenomenologically convenient point of zero momentum transfer in large physical volumes and for two different values of the lattice spacing. Our prediction for the form factor is f+(0)=0.9685(34)(14) where the first error is statistical and the second error systematic. This result can be combined with experimental measurements of K->pi decays for a determination of the CKM-matrix element for which we predict |Vus|=0.2233(5)(9) where the first error is from experiment and the second error from the lattice computation., 21 pages, 7 figures, 6 tables

We present results for the matrix elements of the additional $\Delta S=2$ operators that appear in models of physics beyond the Standard Model (BSM), expressed in terms of four BSM $B$-parameters. Combined with experimental results for $\Delta M_K$ and $\epsilon_K$, these constrain the parameters of BSM models. We use improved staggered fermions, with valence HYP-smeared quarks and $N_f=2+1$ flavors of "asqtad" sea quarks. The configurations have been generated by the MILC collaboration. The matching between lattice and continuum four-fermion operators and bilinears is done perturbatively at one-loop order. We use three lattice spacings for the continuum extrapolation: $a\approx 0.09$, $0.06$ and $0.045\;$fm. Valence light-quark masses range down to $\approx m_s^{\rm phys}/13$ while the light sea-quark masses range down to $\approx m_s^{\rm phys}/20$. Compared to our previous published work, we have added four additional lattice ensembles, leading to better controlled extrapolations in the lattice spacing and sea-quark masses. We report final results for two renormalization scales, $\mu=2\;\text{GeV}$ and $3\;\text{GeV}$, and compare them to those obtained by other collaborations. Agreement is found for two of the four BSM $B$-parameters ($B_2$ and $B_3^\text{SUSY}$). The other two ($B_4$ and $B_5$) differ significantly from those obtained using RI-MOM renormalization as an intermediate scheme, but are in agreement with recent preliminary results obtained by the RBC-UKQCD collaboration using RI-SMOM intermediate schemes., Comment: 17 pages, 24 figures, 13 tables, and updated

We present updated results for kaon B-parameters for operators arising in models of new physics. We use HYP-smeared staggered quarks on the $N_f = 2+1$ MILC asqtad lattices. During the last year we have added new ensembles, which has necessitated chiral-continuum fitting with more elaborate fitting functions. We have also corrected an error in a two-loop anomalous dimension used to evolve results between different scales. Our results for the beyond-the-Standard-Model B-parameters have total errors of $5-10$\%. We find that the discrepancy observed last year between our results and those of the RBC/UKQCD and ETM collaborations for some of the B-parameters has been reduced from $4\!-\!5\,\sigma$ to $2\!-\!3\,\sigma$., Comment: 7 pages, 8 figures, Lattice 2014 proceeding

We calculate the kaon mixing B-parameters for operators arising generically in theories of physics beyond the standard model. We use HYP-smeared improved staggered fermions on the $N_f = 2+1$ MILC asqtad lattices. Operator matching is done perturbatively at one-loop order. Chiral extrapolations are done using "golden combinations" in which one-loop chiral logarithms are absent. For the combined sea-quark mass and continuum extrapolation, we use three lattice spacings: $a \approx 0.045, 0.06$ and $0.09 \text{fm}$. Our results have a total error of 5-6%, which is dominated by the systematic error from matching and continuum extrapolation. For two of the BSM $B$-parameters, we agree with results obtained using domain-wall and twisted-mass dynamical fermions, but we disagree by $(4-5)\sigma$ for the other two., Comment: 7 pages, 5 figures, Lattice 2013 Proceeding

We report on a study of the finite-temperature QCD transition region for temperatures between 139 and 196 MeV, with a pion mass of 200 MeV and two space-time volumes: $24^3\times8$ and $32^3\times8$, where the larger volume varies in linear size between 5.6 fm (at T=139 MeV) and 4.0 fm (at T=195 MeV). These results are compared with the results of an earlier calculation using the same action and quark masses but a smaller, $16^3\times8$ volume. The chiral domain wall fermion formulation with a combined Iwasaki and dislocation suppressing determinant ratio gauge action are used. This lattice action accurately reproduces the $\sua$ and $\ua$ symmetries of the continuum. Results are reported for the chiral condensates, connected and disconnected susceptibilities and the Dirac eigenvalue spectrum. We find a pseudo-critical temperature, $T_c$, of approximately 165 MeV consistent with previous results and strong finite volume dependence below $T_c$. Clear evidence is seen for $\ua$ symmetry breaking above $T_c$ which is quantitatively explained by the measured density of near-zero modes in accordance with the dilute instanton gas approximation., Comment: 58 pages, 9 figures and 8 tables

We update our results for $B_K$ obtained using HYP-smeared staggered valence quarks on the MILC asqtad lattices. In the last year, we have added 5 new measurments on the fine ($a\approx 0.09\;$fm) ensembles, and 2 new measurements on the superfine ($a\approx 0.06\;$fm) ensembles. These allow a simultaneous extrapolation in $a^2$ and sea quark masses, reducing the corresponding systematic error significantly. Our updated result is $\hat{B}_K = 0.738 \pm 0.005 (\text{stat}) \pm 0.034 (\text{sys})$., Comment: 7 pages, 4 figures, 2 tables, Lattice 2013 proceeding

We present results for matrix elements of $\Delta S=2$ four-fermion operators arising generically in models of new physics. These are needed to constrain such models using the measured values of $\varepsilon_K$ and $\Delta M_K$. We use lattice QCD with 2+1 flavors of improved staggered fermions on lattices generated by the MILC collaboration. We extrapolate to the continuum from three lattice spacings ranging down to $a\approx 0.045\;$fm. Total errors are $\sim 5-6%$, arising primarily from our use of one-loop matching between lattice and continuum operators. For two of the matrix elements, our results disagree significantly from those obtained using different fermion discretizations., Comment: 5 pages, 3 figures, added references

We report on recent progress in the calculation of $B_K$ using HYP-smeared improved staggered fermions on the MILC asqtad lattices. We have added measurements on fine ($a\sim 0.09 $fm) and superfine ($a\sim 0.06 $fm) ensembles at different values of the light sea quark mass ($a m_\ell$), as well as increased the statistics on some other ensembles. We find that the results on the fine lattices show a significantly stronger $am_\ell$ dependence than those on the superfine and coarse ($a\sim 0.12 $fm) lattices. We discuss different methods for accounting for these new results when doing the $a m_\ell$ and continuum extrapolations., 7 pages, 6 figures

We calculate the B-parameters for operators arising in theories of new physics beyond the standard model (BSM) using HYP-smeared improved staggered fermions on the MILC asqtad lattices with N_f = 2+1 flavors. We use three different lattice spacings (a ~ 0.045, 0.06 and 0.09 fm) at a fixed ratio of light to strange quarks, m_l/m_s = 1/5, to obtain the continuum results. Operator matching is done using perturbative matching at one-loop order, and results are run to 2 or 3 GeV using two-loop running in the MS-bar scheme. For the chiral and continuum extrapolations, we use SU(2) staggered chiral perturbation theory. We present preliminary results with only statistical errors., 7 pages, 4 figures, Contribution to proceedings of 30th International Symposium on Lattice Field Theory (Lattice 2012), June 24-29, 2012; Cairns, Australia

We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_\tau = 8 and spatial extent N_\sigma = 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_\pi = 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations., Comment: 80 pages, 14 figures, 4 appendices

We present a calculation of the kaon $B$-parameter, $B_K$, using lattice QCD. We use improved staggered valence and sea fermions, the latter generated by the MILC collaboration with $N_f=2+1$ light flavors. To control discretization errors, we use four different lattice spacings ranging down to $a\approx 0.045\;$fm. The chiral and continuum extrapolations are done using SU(2) staggered chiral perturbation theory. Our final result is $\hat{B}_K = 0.727 \pm 0.004 (\text{stat}) \pm 0.038 (\text{sys})$, where the dominant systematic error is from our use of truncated (one-loop) matching factors., 4 pages, 4 figures, 2 tables, updated references, updated grant numbers

We report a possible solution to the trouble that the covariance fitting fails when the data is highly correlated and the covariance matrix has small eigenvalues. As an example, we choose the data analysis of highly correlated $B_K$ data on the basis of the SU(2) staggered chiral perturbation theory. Basically, the essence of the problem is that we do not have an accurate fitting function so that we cannot fit the highly correlated and precise data. When some eigenvalues of the covariance matrix are small, even a tiny error of fitting function can produce large chi-square and spoil the fitting procedure. We have applied a number of prescriptions available in the market such as diagonal approximation and cutoff method. In addition, we present a new method, the eigenmode shift method which fine-tunes the fitting function while keeping the covariance matrix untouched., 7 pages, 4 figures

We discuss finite volume errors in our calculations of $B_K$ using improved staggered fermions on the MILC asqtad lattices. Using GPUs, we are now able to extrapolate using next-to-leading order (NLO) staggered SU(2) chiral perturbation theory including the finite volume corrections arising from pion loops. We find that the impact of FV fitting is very small, giving a 0.5% shift in the continuum limit., Comment: 7 pages, 6 figures, Lattice 2011 proceeding

We present recent progress in our calculation of $B_K$ with improved staggered fermions using chiral extrapolations based on SU(3) staggered chiral perturbation theory. We have accumulated significantly higher statistics on the coarse, fine, and ultrafine MILC asqtad lattices. This leads to a reduction in statistical error and an improved continuum extrapolation. Our updated result is $\hat{B}_K = B_K(\text{RGI}) = 0.737 \pm 0.003(\text{stat}) \pm 0.046 (\text{sys})$. This is consistent with the result obtained using chiral extrapolations based on SU(2) staggered chiral perturbation theory, although the total error is somewhat larger with the SU(3) analysis., Comment: 7 pages, 6 figures

We report on recent progress in the calculation of $B_K$ using HYP-smeared staggered fermions on the MILC asqtad lattices. Our main focus is on the continuum extrapolation, which is done using (up to) four different lattice spacings---$a\approx$ 0.12, 0.09, 0.06 and 0.045 fm. Since Lattice 2010, we have reduced the statistical errors on the $a\approx 0.09\;$fm lattices by a factor of $\sim 3$, and roughly doubled the size of the $a\approx0.045\;$fm ensemble. We find that these improvements have a very significant impact on the continuum extrapolation, with the $a\approx 0.12\;$fm data lying outside the range of applicability of simple functional forms. Hence we use only the three smallest lattice spacings to perform the extrapolation, finding $\hat{B}_K = B_K(\text{RGI}) = 0.725 \pm 0.004(\text{stat}) \pm 0.038(\text{sys}) $. This value is consistent with our published value from 2010 (based the three coarsest lattice spacings), but has smaller errors., Comment: 7 pages, 6 figures, references updated

We address a frequently asked question on the covariance fitting of the highly correlated data such as our $B_K$ data based on the SU(2) staggered chiral perturbation theory. Basically, the essence of the problem is that we do not have an accurate fitting function enough to fit extremely precise data. When eigenvalues of the covariance matrix are small, even a tiny error of fitting function yields large chi-square and spoils the fitting procedure. We have applied a number of prescriptions available in the market such as the cut-off method, modified covariance matrix method, and Bayesian method. We also propose a brand new method, the eigenmode shift method which allows a full covariance fitting without modifying the covariance matrix at all. In our case, the eigenmode shift (ES) method and Bayesian method turn out to be the best prescription to the problem. We also provide a pedagogical example of data analysis in which the diagonal approximation and the cut-off method fail in fitting manifestly, but the ES method and the Bayesian approach work well., Comment: 33 pages, 8 figures, updated appendix

We report updated results for $B_K$ using HYP-smeared staggered valence quarks on MILC asqtad lattices based on an analysis using SU(3) staggered chiral perturbation theory. The most important new feature of our data sample is the inclusion of a fourth ("ultrafine") lattice spacing. This improves the control over the continuum extrapolation and errors due our use of one-loop perturbative matching. We present a complete updated error budget, which leads to $B_K(\text{NDR}, \mu = 2 \text{GeV}) = 0.5309 \pm 0.0051 \pm 0.0424$ and $\hat{B}_K = B_K(\text{RGI}) = 0.727 \pm 0.07 \pm 0.058$. The results of the SU(3) analysis are inferior to those based on SU(2) staggered chiral perturbation theory, primarily because of the dependence on the Bayesian priors we use in the SU(3) fits., Comment: 7 pages, 6 figures, Lattice 2010 Proceeding, references updated

We study three sources of error in our calculation of $B_K$ using HYP-smeared staggered fermions on the MILC asqtad lattices. These are (1) dependence on the light sea quark mass; (2) finite volume effects; and (3) the impact of an order of magnitude increase in the number of measurements. Our main results are (1) the dependence on the light sea-quark mass is weaker than expected by naive dimensional analysis, (2) including finite volume effects in SU(2) staggered chiral perturbation theory fits leads to a very small change in $B_K$, of size $\approx 0.1%$, and (3) increasing the statistics on one of the coarse MILC lattices resolves a potential discrepancy with other coarse results., 7 pages, 8 figures, Lattice 2010 Proceeding, references updated

We report updated results for $B_K$ calculated using HYP-smeared staggered fermions on the MILC asqtad 2+1 flavor lattices. We use four different lattice spacings ($a \approx$ 0.12, 0.09, 0.06 and 0.045 fm) to control the continuum extrapolation. We use SU(2) staggered chiral perturbation theory to do the data analysis. We find that $B_K(\text{NDR}, \mu=2 \text{GeV}) = 0.526 \pm 0.007 \pm 0.024$ and $\hat{B}_K = B_K(\text{RGI}) = 0.720 \pm 0.010 \pm 0.033$. Here the first error is statistical and the second systematic. The dominant source of error is that due to our use of a truncated (one-loop) matching factor., Comment: 7 pages, 7 figures, Lattice 2010 proceeding, references updated

We present preliminary results for $B_K$ from the MILC ultrafine lattices, based on a partial ensemble of 305 configurations. We use HYP-smeared improved staggered valence quarks. The analysis is done using fitting forms based on both SU(2) and SU(3) staggered chiral perturbation thery. For the SU(2) analysis, we find that the result using the NLO fit function is consistent with that from a partial NNLO fit. For the SU(3) analysis, where we have to use partially constrained fits due to the number of fit parameters, we find that our two preferred fits ("N-BB1" and "N-BB2") are also consistent, both with each other and with the results of the SU(2) fits. These results are used in companion proceedings to improve the control over the continuum extrapolation., Comment: 7 pages, 9 figures, Lattice 2010 Proceeding, references updated

We present results for kaon mixing parameter $B_K$ calculated using HYP-smeared improved staggered fermions on the MILC asqtad lattices. We use three lattice spacings ($a\approx 0.12$, $0.09$ and $0.06\;$fm), ten different valence quark masses ($m\approx m_s/10-m_s$), and several light sea-quark masses in order to control the continuum and chiral extrapolations. We derive the next-to-leading order staggered chiral perturbation theory (SChPT) results necessary to fit our data, and use these results to do extrapolations based both on SU(2) and SU(3) SChPT. The SU(2) fitting is particularly straightforward because parameters related to taste-breaking and matching errors appear only at next-to-next-to-leading order. We match to the continuum renormalization scheme (NDR) using one-loop perturbation theory. Our final result is from the SU(2) analysis, with the SU(3) result providing a (less accurate) cross check. We find $B_K(\text{NDR}, \mu = 2 \text{GeV}) = 0.529 \pm 0.009 \pm 0.032$ and $\hat{B}_K =B_K(\text{RGI})= 0.724 \pm 0.012 \pm 0.043$, where the first error is statistical and the second systematic. The error is dominated by the truncation error in the matching factor. Our results are consistent with those obtained using valence domain-wall fermions on lattices generated with asqtad or domain-wall sea quarks., Comment: 37 pages, 31 figures, most updated version

We present the results for $B_K$ calculated using HYP-smeared staggered fermions using lattices generated by the MILC collaboration using the asqtad staggered action. We have done the calculation on 8 ensembles of these lattices, including three different lattice spacings ($a=0.12, 0.09, 0.06$ fm). We fit the data to forms based on those predicted by SU(3) mixed-action partially-quenched staggered chiral perturbation theory. Our preliminary result is $ B_K(\text{NDR},\mu=2\text{GeV}) = 0.528 \pm 0.011 \pm 0.048 $, where the first error is statistical and the second systematic. The error turns out to be larger than that from an analysis using SU(2) chiral perturbation theory., Comment: 7 pages, 2 figures, Lattice 2009 proceeding

We present results for $B_K$ calculated using HYP-smeared improved staggered fermions on the MILC asqtad lattices. In this report, the data is analyzed using the results of SU(2) staggered chiral perturbation theory (SChPT). We outline the derivation of the NLO SU(2) SChPT result, explain our fitting procedure, and outline how we estimate systematic errors. We also show the light sea-quark mass and lattice spacing dependence for both SU(2) and SU(3)-based analyses. Our preliminary result from the SU(2) analysis is $B_K(\text{NDR}, \mu = 2 \text{GeV}) = 0.512 \pm 0.014 \pm 0.034$ and $\hat{B}_K = B_K(\text{RGI})= 0.701 \pm 0.019 \pm 0.047$. This is somewhat more accurate than our result from the SU(3) analysis. It is consistent with results obtained using valence domain-wall fermions, Comment: 7 pages, 6 figures, Lattice 2009 proceeding

We study the finite-volume effects in our calculation of $B_K$ using HYP-smeared improved staggered valence fermions. We calculate the predictions of both SU(3) and SU(2) staggered chiral perturbation theory at one-loop order. We compare these to the results of a direct calculation, using MILC coarse lattices with two different volumes: $20^3$ and $28^3$. From the direct calculation, we find that the finite volume effect is $\approx 2%$ for the SU(3) analysis and $\approx 0.9%$ for the SU(2) analysis. We also show how the statistical error depends on the number of measurements made per configuration, and make a first study of autocorrelations., Comment: 7 pages, 10 figures, Lattice 2009 proceeding

We give an update on our calculation of $B_K$ using HYP-smeared valence staggered quarks. We have results for $B_K$ at tree-level on several coarse MILC lattices ($a\approx 0.12 $fm) and one of the fine lattices ($a\approx 0.09 $fm), using 10 light valence quarks ranging down to $m_s^{\rm phys}/10$. We have generalized staggered chiral perturbation theory to our mixed action setup, and outline the results. We explain our present fitting strategy, and give some preliminary results., 7 pages, 4 figures, Lattice 2008 proceeding

We study the pion spectrum (and in particular taste-symmetry breaking within it) using HYP-smeared valence staggered fermions on the coarse and fine MILC lattices (which have asqtad staggered sea quarks). We focus on the dependence on lattice spacing and sea-quark mass. We also update our results on source dependence. Our main conclusion is that on the MILC fine lattices the appropriate power-counting for SU(3) staggered chiral perturbation theory may have discretization errors entering at next-to-leading order rather than at leading-order., 7 pages, 8 figures, Lattice 2008 proceeding

We study the impact of hypercubic (HYP) smearing on the size of taste breaking for staggered fermions, comparing to unimproved and to asqtad-improved staggered fermions. As in previous studies, we find a substantial reduction in taste-breaking compared to unimproved staggered fermions (by a factor of 4-7 on lattices with spacing $a\approx 0.1 $fm). In addition, we observe that discretization effects of next-to-leading order in the chiral expansion (${\cal O}(a^2 p^2)$) are markedly reduced by HYP smearing. Compared to asqtad valence fermions, we find that taste-breaking in the pion spectrum is reduced by a factor of 2.5-3, down to a level comparable to the expected size of generic ${\cal O}(a^2)$ effects. Our results suggest that, once one reaches a lattice spacing of $a\approx 0.09 $fm, taste-breaking will be small enough after HYP smearing that one can use a modified power counting in which ${\cal O}(a^2) \ll {\cal O}(p^2)$, simplify fitting to phenomenologically interesting quantities., 14 pages, 13 figures, references updated, minor changes

We describe details of 1/a ~ 2.2Gev, L ~ 3 fm dynamical domain wall fermion simulations which will allow us to do a more systematic continuum extrapolation in combination with existing simu- lations. Details of the simulations such as algorithm choices and machine performance, as well as results of basic measurements are presented. These configurations are presently being generated on the QCDOC machine at Edinburgh and the DOE QCDOC machine at Brookhaven as part of a joint project with LHPC., 7 pages, presented at the 25th International Symposium on Lattice Field Theory, 30 Jul - 4 Aug 2007, Regensburg, Germany

We extend our previous study of taste-symmetry breaking using HYP-smeared staggered fermions in two ways. First, we improve the statistics of a comparison of unimproved and HYP-smeared staggered fermions on quenched lattices with $a\approx 0.1 $fm. This allows us to obtain a signal for all pion tastes, rather than just a subset, and thus to make a complete comparison. In addition, it allows us to differentiate between wall and local sources. Second, we compare HYP-smeared valence quarks to asqtad valence quarks on 2+1 flavor unquenched MILC lattices. We find that taste breaking is substantially reduced by HYP-smearing, bringing the size of this discretization effect (which is the dominant such effect with staggered fermions) down to the size expected generically for any fermion type., 7 pages, 7 figures

We present results for the dependence of the residual mass of domain wall fermions (DWF) on the size of the fifth dimension and its relation to the density and localization properties of low-lying eigenvectors of the corresponding hermitian Wilson Dirac operator relevant to simulations of 2+1 flavor domain wall QCD. Using the DBW2 and Iwasaki gauge actions, we generate ensembles of configurations with a $16^3\times 32$ space-time volume and an extent of 8 in the fifth dimension for the sea quarks. We demonstrate the existence of a regime where the degree of locality, the size of chiral symmetry breaking and the rate of topology change can be acceptable for inverse lattice spacings $a^{-1} \ge 1.6$ GeV., Comment: 59 Pages, 23 figures, 1 MPG linked

The decays and mixing of $K$ mesons are remarkably sensitive to the weak interactions of quarks and leptons at high energies. They provide important tests of the standard model at both first and second order in the Fermi constant $G_F$ and offer a window into possible new phenomena at energies as high as 1,000 TeV. These possibilities become even more compelling as the growing capabilities of lattice QCD make high-precision standard model predictions possible. Here we discuss and attempt to forecast some of these capabilities., Comment: submitted to the Rare Processes and Precision, Theory and Computational Frontiers for the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

The QCDSP machine at Columbia University has grown to 2,048 nodes achieving a peak speed of 100 Gigaflops. Software for quenched and Hybrid Monte Carlo (HMC) evolution schemes has been developed for staggered fermions, with support for Wilson and clover fermions under development. We provide an overview of the runtime environment, the current status of the QCDSP construction program and preliminary results not presented elsewhere in these proceedings., Comment: 3 pages, 1 figure, LaTeX with espcrc2.sty, talk presented at LATTICE97(Poster) by G. R. Fleming