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

Chulwoo Jung

Chulwoo Jung 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-thphysics.comp-phhep-exnucl-thphysics.data-an

Catalog footprint

What is connected

47works
7topics
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

47 published item(s)

preprint2022arXiv

$K\toππ$ decay matrix elements at the physical point with periodic boundary conditions

We calculate $K\toππ$ 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ππ$ three-point functions and matrix elements on a $24^3, a^{-1} = 1$~GeV, $2+1$-flavor Möbius 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deep learning study on the Dirac eigenvalue spectrum of staggered quarks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discovering new physics in rare kaon decays

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lattice QCD and the Computational Frontier

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chiral Ward identities for Dirac eigenmodes with staggered fermions

We study chiral properties of eigenvalue spectrum for staggered quarks. We present a new method to identify would-be zero modes and nonzero modes using their symmetry and chiral properties. Here, we review the traditional method with HYP improved staggered quarks, and extend it to a completely new method which uses the chiral Ward identities and leakage patterns to achieve the goal.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Light quark vacuum polarization at the physical point and contribution to the muon $g-2$

We report on the computation of the connected light quark vacuum polarization with 2+1+1 flavors of HISQ fermions at the physical point and its contribution to the muon anomalous magnetic moment. Three ensembles, generated by the MILC collaboration, are used to take the continuum limit. The finite volume correction to this result is computed in the (Euclidean) time-momentum representation to NNLO in chiral perturbation theory. We find $a_μ^{ll}({\rm HVP})=(659\pm 20\pm 5\pm 5\pm 4)\times 10^{-10}$, where the errors are statistical and estimates of residual uncertainties from taking the continuum limit, scale setting, and truncation of chiral perturbation theory at NNLO. We compare our results with recent ones in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon mass and isovector couplings in 2+1-flavor dynamical domain-wall lattice QCD near physical mass

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

We report the first result for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment with all errors systematically controlled. Several ensembles using 2+1 flavors of physical mass Möbius 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_μ^{\rm HLbL} = 7.87(3.06)_\text{stat}(1.77)_\text{sys}\times 10^{-10}$. Our value is consistent with previous model results and leaves little room for this notoriously difficult hadronic contribution to explain the difference between the Standard Model and the BNL experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The connected and leading disconnected diagrams of the hadronic light-by-light contribution to muon $g - 2$

We report our recent lattice calculation of hadronic light-by-light contribution to muon $g-2$ using our recently developed moment method. The connected diagrams and the leading disconnected diagrams are included. The calculation is performed on a $48^3 \times 96$ lattice with physical pion mass and 5.5 fm box size. We expect sizable finite volume and finite lattice spacing corrections to the results of these calculations which will be estimated in calculations to be carried out over the next 1-2 years.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Covariant approximation averaging

We present a new class of statistical error reduction techniques for Monte-Carlo simulations. Using covariant symmetries, we show that correlation functions can be constructed from inexpensive approximations without introducing any systematic bias in the final result. We introduce a new class of covariant approximation averaging techniques, known as all-mode averaging (AMA), in which the approximation takes account of contributions of all eigenmodes through the inverse of the Dirac operator computed from the conjugate gradient method with a relaxed stopping condition. In this paper we compare the performance and computational cost of our new method with traditional methods using correlation functions and masses of the pion, nucleon, and vector meson in $N_f=2+1$ lattice QCD using domain-wall fermions. This comparison indicates that AMA significantly reduces statistical errors in Monte-Carlo calculations over conventional methods for the same cost.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kaon BSM B-parameters using improved staggered fermions from $N_f=2+1$ unquenched QCD

We present results for the matrix elements of the additional $Δ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 $ΔM_K$ and $ε_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, $μ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimizing the domain wall fermion Dirac operator using the R-Stream source-to-source compiler

The application of the Dirac operator on a spinor field, the Dslash operation, is the most computation-intensive part of the lattice QCD simulations. It is often the key kernel to optimize to achieve maximum performance on various platforms. Here we report on a project to optimize the domain wall fermion Dirac operator in Columbia Physics System (CPS) using the R-Stream source-to-source compiler. Our initial target platform is the Intel PC clusters. We discuss the optimization strategies involved before and after the automatic code generation with R-Stream and present some preliminary benchmark results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The kaon semileptonic form factor in Nf=2+1 domain wall lattice QCD with physical light quark masses

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Calculation of BSM Kaon B-parameters using Staggered Quarks

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\,σ$ to $2\!-\!3\,σ$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved determination of $B_K$ with staggered quarks

We present results for the kaon mixing parameter $B_K$ obtained using improved staggered fermions on a much enlarged set of MILC asqtad lattices. Compared to our previous publication, which was based largely on a single ensemble at each of the three lattice spacings $a\approx 0.09\;$fm, $0.06\;$fm and $0.045\;$fm, we have added seven new fine and four new superfine ensembles, with a range of values of the light and strange sea-quark masses. We have also increased the number of measurements on one of the original ensembles. This allows us to do controlled extrapolations in the light and strange sea-quark masses, which we do simultaneously with the continuum extrapolation. This reduces the extrapolation error and improves the reliability of our error estimates. Our final result is $\hat{B}_K = 0.7379 \pm 0.0047 (\text{stat}) \pm 0.0365 (\text{sys})$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD chiral transition, $\ua$ symmetry and the Dirac spectrum using domain wall fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD phase transition with physical-mass, chiral quarks

We report on the first lattice calculation of the QCD phase transition using chiral fermions at physical values of the quark masses. This calculation uses 2+1 quark flavors, spatial volumes between (4 fm$)^3$ and (11 fm$)^3$ and temperatures between 139 and 196 MeV . Each temperature was calculated using a single lattice spacing corresponding to a temporal Euclidean extent of $N_t=8$. The disconnected chiral susceptibility, $χ_{\rm disc}$ shows a pronounced peak whose position and height depend sensitively on the quark mass. We find no metastability in the region of the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm. Each result is strong evidence that the QCD ``phase transition'' is not first order but a continuous cross-over for $m_π=135$ MeV. The peak location determines a pseudo-critical temperature $T_c = 155(1)(8)$ MeV. Chiral $SU(2)_L\times SU(2)_R$ symmetry is fully restored above 164 MeV, but anomalous $U(1)_A$ symmetry breaking is non-zero above $T_c$ and vanishes as $T$ is increased to 196 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Beyond the Standard Model B-parameters with improved staggered fermions in $N_f=2+1$ QCD

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)σ$ for the other two.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutral kaon mixing from new physics: matrix elements in $N_f=2+1$ QCD

We present results for matrix elements of $Δ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 $Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Update on $B_K$ and $\varepsilon_K$ with staggered quarks

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})$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Beyond the Standard Model corrections to $K^0-\bar{K}^0$ mixing

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Covariance fitting of highly correlated data in lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Full QED+QCD Low-Energy Constants through Reweighting

The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor dynamical domain-wall quark ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kaon $B$-parameter from improved staggered fermions in $N_f=2+1$ QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Light Quark Mass Reweighting

We present a systematic study of the effectiveness of light quark mass reweighting. This method allows a single lattice QCD ensemble, generated with a specific value of the dynamical light quark mass, to be used to determine results for other, nearby light dynamical quark masses. We study two gauge field ensembles generated with 2+1 flavors of dynamical domain wall fermions with light quark masses m_l=0.02 (m_π=620 MeV) and m_l=0.01 (m_π=420 MeV). We reweight each ensemble to determine results which could be computed directly from the other and check the consistency of the reweighted results with the direct results. The large difference between the 0.02 and 0.01 light quark masses suggests that this is an aggressive application of reweighting as can be seen from fluctuations in the magnitude of the reweighting factor by four orders of magnitude. Never-the-less, a comparison of the reweighed topological charge, average plaquette, residual mass, pion mass, pion decay constant, and scalar correlator between these two ensembles shows agreement well described by the statistical errors. The issues of the effective number of configurations and finite sample size bias are discussed. An examination of the topological charge distribution implies that it is more favorable to reweight from heavier mass to lighter quark mass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Recent progress in calculation of $B_K$ using staggered fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The chiral transition and U(1)_A symmetry restoration from lattice QCD using Domain Wall Fermions

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_τ= 8 and spatial extent N_σ= 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_π= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

$B_K$ from improved staggered fermions using SU(3) chiral perturbation theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuum extrapolation of $B_K$ with staggered fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite volume effects in $B_K$ with improved staggered fermions

We extend our recent unquenched ($N_f=2+1$ flavor) calculation of $B_K$ using improved staggered fermions by including in the fits the finite volume shift predicted by one-loop staggered chiral perturbation theory. The net result is to lower the result in the continuum limit by 0.6%. This shift is slightly smaller than our previous estimate of finite volume effects based on a direct comparison between different volumes. To include the finite volume effects in a reasonable time, we found it necessary to calculate them using Graphics Processing Units.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Trouble shooting for covariance fitting in highly correlated data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

$B_K$ using HYP-smeared staggered fermions in $N_f=2+1$ unquenched QCD

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}, μ= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

$B_K$ with improved staggered fermions: analysis using SU(2) staggered chiral perturbation theory

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}, μ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

$B_K$ with improved staggered fermions: analysis using SU(3) staggered chiral perturbation theory

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}, μ= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First results for $B_K$ on the ultrafine ($a=0.045$ fm) ensemble

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Systematic Uncertainties in $B_K$ with Improved Staggered Fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Determination of $B_K$ using improved staggered fermions (I): SU(3) chiral perturbation theory fit

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},μ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Determination of $B_K$ using improved staggered fermions (II) SU(2) chiral perturbation theory fit

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}, μ= 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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Determination of $B_K$ using improved staggered fermions (III) Finite volume effects

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Calculating $B_K$ using HYP staggered fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Scaling behavior and sea quark dependence of pion spectrum with HYP-smeared staggered fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Thermodynamics using p4-improved staggered fermion action on QCDOC

We present an exploratory study of the thermodynamics of $N_f=3$ QCD with an improved staggered fermions using the QCDOC supercomputer. We use a p4 action with MILC-style smeared links (Fat 7). Some details of the implementation of the p4 action on QCDOC are discussed and performance benchmarks are given. We show preliminary results for the quark mass dependence of the pseudo-critical temperature $T_c$ from several lattice volumes . We also make a comparison between p4fat7 and the old p4 action.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Relating U(1) monopole configurations to SU(2) saddle-point configurations

We have studied field configurations of the 3-dimensional Georgi-Glashow model which interpolate between the $U(1)$ and the $SU(2)$ limits. In the intermediate region, these configurations contain 't-Hooft--Polyakov monopoles. We use cooling and extremization to find these configurations and investigate their evolution as we adiabatically move towards the $U(1)$ and the $SU(2)$ limits. We also evolve an $SU(2)$ saddle point solution towards the $U(1)$ limit to see the relation between the unstable solutions in the $SU(2)$ theory and the stable ones in the $U(1)$ theory.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Analysis of saddle-point configurations in 3-dimensional SU(2) gauge theory

We discuss the properties of a class of saddle point solutions in SU(2) in three dimensions (SU$(2)_3$), exhibiting localized peaks in the action. These configurations are generated by deterministic cooling and extremizing algorithms from analytic configurations. They share some characteristics with cooled and extremized Monte Carlo generated lattices. We have investigated physical behavior such as the string tension by averaging over this class of saddle point configurations. We have also measured the eigenvalues for harmonic fluctuations around these configurations.

0 citations0 reviewsNo code linkedOpen access