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

Boram Yoon

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

ResearcherUnclaimed source record
hep-lathep-phnucl-thMachine Learningquant-phhep-exhep-thHuman-Computer Interactionmath.NAnucl-exphysics.data-an

Catalog footprint

What is connected

63works
11topics
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

63 published item(s)

preprint2024arXiv

Update on flavor diagonal nucleon charges from clover fermions

We present a summary of the full calculation of the axial, scalar and tensor flavor diagonal charges of the nucleon carried out using Wilson-clover fermions on eight ensembles generated using 2+1+1-flavors of highly improved staggered quarks (HISQ) by the MILC collaboration. We also give results for the $3\times 3$ matrix of renormalization factors between the RI-sMOM and $\overline{\rm MS}$ scheme for the 2+1 flavor theory that include flavor mixing. Preliminary results for $g_{A,S,T}^{u,d,s}$ are presented in the $\overline{\rm MS}$ scheme at scale 2 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

nEDM from the theta-term and chromoEDM operators

In a previous work, we showed that unresolved excited state contaminations provide a major source of systematic uncertainty in the calculation of the nucleon electric dipole moment due to the QCD topological term theta. Here we extend the calculation to the quark chromo-electric dipole moment operator (qcEDM). We also show quantitatively the impact of the mixing of the qcEDM with lower-dimensional operators on the lattice. Finally, we present preliminary results from a unitary clover-on-clover calculation for the QCD topological term.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Nucleon form factors and the pion-nucleon sigma term

This talk summarizes the progress made since Lattice 2021 in understanding and controlling the contributions of towers of multihadron excited states with mass gaps starting lower than of radial excitations, and in increasing our confidence in the extraction of ground state nucleon matrix elements. The most clear evidence for multihadron excited state contributions (ESC) is in axial/pseudoscalar form factors that are required to satisfy the PCAC relation between them. The talk examines the broader question--which and how many of the theoretically allowed positive parity states $N(\textbf p)π(-\textbf p)$, $N(\textbf 0)π(\textbf 0)π(\textbf 0)$, $N(\textbf p)π(\textbf 0)$, $N(\textbf 0)π(\textbf p),\ \ldots$ make significant contributions to a given nucleon matrix element? New data for the axial, electric and magnetic form factors are presented. They continue to show trends observed in Ref[1]. The N${}^2$LO $χ$PT analysis of the ESC to the pion-nucleon sigma term, $σ_{πN}$, has been extended to include the $Δ$ as an explicit degree of freedom [2]. The conclusion reached in Ref [3] that $N π$ and $N ππ$ states each contribute about 10 MeV to $σ_{πN}$, and the consistency between the lattice result with $N π$ state included and the phenomenological estimate is not changed by this improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Update on Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the sum of the connected and disconnected contributions for the up, down and strange quarks, control over fits to remove excited state contamination, and the simultaneous chiral-continuum fit used to extract the charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

AI and Theoretical Particle Physics

Theoretical particle physicists continue to push the envelope in both high performance computing and in managing and analyzing large data sets. For example, the goals of sub-percent accuracy in predictions of quantum chromodynamics (QCD) using large scale simulations of lattice QCD and in finding signals of rare events and new physics in exabytes of data produced by experiments at the high luminosity large hadron collider (LHC) require new tools beyond just developments in hardware. Machine learning and artificial intelligence offer the promise of dramatically reducing the computational cost and time. This chapter reviews selected areas where AI/ML tools could have a major impact, provides an overview of the challenges, and discusses how new ideas such as normalizing flows can speed up the generation of gauge configurations needed in lattice QCD calculations; the growth of ML in surrogate models and pattern matching to reduce the cost of event generators and in the analysis of experimental data; and in the search for viable vacua in the landscape of string theories. While such approaches transform aspects of particle theory into computational problems, and thus black boxes, we argue that physics-aware development of these tools combined with algorithms that ensure that the results are bias free will continue to require a deep understanding of the physics. We see this broader transformation as akin to formulating and extracting observables from simulations of lattice QCD, a numerical integration of the path integral formulation of QCD that nevertheless requires a deep understanding of the underlying quantum field theory, the standard model of particle physics and effective field theory methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Calculation of neutron electric dipole moment due to the QCD topological term, Weinberg three-gluon operator and the quark chromoelectric moment

We present results for the neutron electric dipole moment due to the dimension 4 and dimension 6 gluonic CP violation, and the isovector quark chromoelectric dipole moment using clover valence quarks on HISQ dynamical ensembles generated by the MILC Collaboration. For the gluonic operators, we use the gradient flow scheme to obtain divergence-free continuum results. For the chromoelectric dipole moment operator, we use the unflowed local operator but discuss how the quadratically divergent mixing with the pseudoscalar operator can be controlled nonperturbatively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effects of Virtual Room Size and Objects on Relative Translation Gain Thresholds in Redirected Walking

This paper investigates how the size of virtual space and objects within it affect the threshold range of relative translation gains, a Redirected Walking (RDW) technique that scales the user's movement in virtual space in different ratios for the width and depth. While previous studies assert that a virtual room's size affects relative translation gain thresholds on account of the virtual horizon's location, additional research is needed to explore this assumption through a structured approach to visual perception in Virtual Reality (VR). We estimate the relative translation gain thresholds in six spatial conditions configured by three room sizes and the presence of virtual objects (3 X 2), which were set according to differing Angles of Declination (AoDs) between eye-gaze and the forward-gaze. Results show that both size and virtual objects significantly affect the threshold range, it being greater in the large-sized condition and furnished condition. This indicates that the effect of relative translation gains can be further increased by constructing a perceived virtual movable space that is even larger than the adjusted virtual movable space and placing objects in it. Our study can be applied to adjust virtual spaces in synchronizing heterogeneous spaces without coordinate distortion where real and virtual objects can be leveraged to create realistic mutual spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Excited states and precision results for nucleon charges and form factors

The exponentially falling signal-to-noise ratio in all nucleon correlation functions, and the presence of towers of multihadron excited states with relatively small mass gaps makes extraction of matrix elements of various operators within the ground state nucleon challenging. Theoretically, the allowed positive parity states with the smallest mass gaps are the $N(\bm p)π(-\bm p)$, $N(\bm 0)π(\bm 0)π(\bm 0)$, $N(\bm p)π(\bm 0)$, $N(\bm 0)π(\bm p),\ \ldots$, states. A priori, the contribution of these states arises at one loop in chiral perturbation theory ($χ$PT), however, in many cases the contributions are enhanced. In this talk, I will review four such cases: the correlation functions from which the axial form factors, electric and magnetic form factors, the $Θ$-term contribution to neutron electric dipole moment (nEDM), and the pion-nucleon sigma term are extracted. Including appropriate multihadron states in the analysis can lead to significantly different results compared to standard analyses with the mass gaps taken from fits to 2-point functions. The $χ$PT case for $N π$ states is the most clear in the axial/pseudoscalar form factors which need to satisfy the PCAC relation between them. Our analyses, supported by $χ$PT, suggests similarly large effects in the calculations of the $Θ$-term and the pion-nucleon sigma term that have significant phenomenological implications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the connected and disconnected contributions, calculation of the renormalization constants and mixing in the RI-sMOM scheme, and control over the simultaneous chiral-continuum-finite-volume fit used to extract the final charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Precision Nucleon Charges and Form Factors Using 2+1-flavor Lattice QCD

We present high statistics results for the isovector nucleon charges and form factors using seven ensembles of 2+1-flavor Wilson-clover fermions. The axial and pseudoscalar form factors obtained on each ensemble satisfy the PCAC relation once the lowest energy $Nπ$ excited state is included in the spectral decomposition of the correlation functions used for extracting the ground state matrix elements. Similarly, we find evidence that the $Nππ$ excited state contributes to the correlation functions with the vector current, consistent with the vector meson dominance model. The resulting form factors are consistent with the Kelly parameterization of the experimental electric and magnetic data. Our final estimates for the isovector charges are $g_{A}^{u-d} = 1.31(06)(05)_{sys}$, $g_{S}^{u-d} = 1.06(10)(06)_{sys}$, and $g_{T}^{u-d} = 0.95(05)(02)_{sys}$, where the first error is the overall analysis uncertainty and the second is an additional combined systematic uncertainty. The form factors yield: (i) the axial charge radius squared, ${\langle r_A^2 \rangle}^{u-d}=0.428(53)(30)_{sys}\ {\rm fm}^2$, (ii) the induced pseudoscalar charge, $g_P^\ast=7.9(7)(9)_{sys}$, (iii) the pion-nucleon coupling $g_{π{\rm NN}} = 12.4(1.2)$, (iv) the electric charge radius squared, ${\langle r_E^2 \rangle}^{u-d} = 0.85(12)(19)_{sys} \ {\rm fm}^2$, (v) the magnetic charge radius squared, ${\langle r_M^2 \rangle}^{u-d} = 0.71(19)(23)_{\rm sys} \ {\rm fm}^2$, and (vi) the magnetic moment $μ^{u-d} = 4.15(22)(10)_{\rm sys}$. All our results are consistent with phenomenological/experimental values but with larger errors. Lastly, we present a Padé parameterization of the axial, electric and magnetic form factors over the range $0.04< Q^2 <1$ GeV${}^2$ for phenomenological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The pion-nucleon sigma term from Lattice QCD

We summarize recent evidence, both from lattice QCD and chiral perturbation theory, that suggests that larger-than-expected excited-state contamination could be the reason for the tension between phenomenological determinations and previous direct lattice-QCD calculations of the pion--nucleon sigma term $σ_{πN}$. In addition, we extend the $χ$PT analysis by calculating the corrections due to including the $Δ(1232)$ resonance as an explicit degree of freedom. This correction is found to be small, thereby corroborating the excited-state effects found in the $Δ$-less calculation and the result for $σ_{πN}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A regression algorithm for accelerated lattice QCD that exploits sparse inference on the D-Wave quantum annealer

We propose a regression algorithm that utilizes a learned dictionary optimized for sparse inference on a D-Wave quantum annealer. In this regression algorithm, we concatenate the independent and dependent variables as a combined vector, and encode the high-order correlations between them into a dictionary optimized for sparse reconstruction. On a test dataset, the dependent variable is initialized to its average value and then a sparse reconstruction of the combined vector is obtained in which the dependent variable is typically shifted closer to its true value, as in a standard inpainting or denoising task. Here, a quantum annealer, which can presumably exploit a fully entangled initial state to better explore the complex energy landscape, is used to solve the highly non-convex sparse coding optimization problem. The regression algorithm is demonstrated for a lattice quantum chromodynamics simulation data using a D-Wave 2000Q quantum annealer and good prediction performance is achieved. The regression test is performed using six different values for the number of fully connected logical qubits, between 20 and 64. The scaling results indicate that a larger number of qubits gives better prediction accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Lossy compression of statistical data using quantum annealer

We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3.5 times better compression performance than the algorithms based on a neural-network autoencoder and principal component analysis. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Moments of nucleon isovector structure functions in $2+1+1$-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u-d}$, helicity moment, $\langle x \rangle_{Δu-Δd}$, and the transversity moment, $\langle x \rangle_{δu-δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC collaboration using $2+1+1$-flavors of dynamical highly improved staggered quarks (HISQ). The correlation functions are calculated using the Wilson-Clover action and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI${}^\prime$-MOM scheme. The data have been collected at lattice spacings $a \approx 0.15,\ 0.12,\ 0.09,$ and 0.06 fm and $M_π\approx 310,\ 220$ and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{MS}$ scheme at 2 GeV, are $\langle x \rangle_{u-d} = 0.173(14)(07)$, $\langle x \rangle_{Δu-Δd} = 0.213(15)(22)$ and $\langle x \rangle_{δu-δd} = 0.208(19)(24)$, where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nucleon isovector momentum fraction, helicity and transversity moment using Lattice QCD

We present our recent high precision calculations (Phys. Rev. D102 (2020) no.5, 054512 and JHEP 04 (2021) 044, JHEP 21 (2020) 004) of the first moment of nucleon isovector polarized, unpolarized and transversity distributions, i.e., momentum fraction, helicity and transversity moment, respectively. We use the standard method for the calculation of these moments (via matrix elements of twist two operators), and carry out a detailed analysis of the sources of systematic uncertainty, in particular of excited state contributions. Our calculations have been performed using two different lattice setups (Clover-on-HISQ and Clover-on-Clover), each with several ensembles. They give consistent results that are in agreement with global fit analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Heavy Flavor Program for the Future Electron-Ion Collider

The proposed high-energy and high-luminosity Electron-Ion Collider (EIC) will provide one of the cleanest environments to precisely determine the nuclear parton distribution functions (nPDFs) in a wide $x$-$Q^{2}$ range. Heavy flavor production at the EIC provides access to nPDFs in the poorly constrained high Bjorken-$x$ region, allows us to study the quark and gluon fragmentation processes, and constrains parton energy loss in cold nuclear matter. Scientists at the Los Alamos National Laboratory are developing a new physics program to study heavy flavor production, flavor tagged jets, and heavy flavor hadron-jet correlations in the nucleon/nucleus going direction at the future EIC. The proposed measurements will provide a unique way to explore the flavor dependent fragmentation functions and energy loss in a heavy nucleus. They will constrain the initial-state effects that are critical for the interpretation of previous and ongoing heavy ion measurements at the Relativistic Heavy Ion Collider and the Large Hadron Collider. We show an initial conceptual design of the proposed Forward Silicon Tracking (FST) detector at the EIC, which is essential to carry out the heavy flavor measurements. We further present initial feasibility studies/simulations of heavy flavor hadron reconstruction using the proposed FST.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Machine-Learning Prediction for Quasi-PDF Matrix Elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-$x$ dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken-$x$) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small-$x$ reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_z a)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements (in finer lattice units), increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning (ML) algorithms can make correlator predictions to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different fit quality and systematic errors. The predictions from smaller displacement $z$ to larger ones work better than those for momentum $p$ due to the higher correlation among the data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutron Electric Dipole Moments with Clover Fermions

We present preliminary results for the contributions to the neutron EDM arising from the QCD $θ$-term, the Weinberg three-gluon and the quark chromo-EDM operators from our ongoing lattice calculations using clover valence quarks on the MILC HISQ lattices. We use the gradient-flow technique to smooth the lattices and renormalize the gluonic operators, and use the Schwinger source method to incorporate the quark chromo-EDM interactions in the quark propagator. For the QCD $θ$-term and the Weinberg three-gluon operator, we report results in the gradient-flow scheme from 8 ensembles at four lattice spacings and three pion masses, including 2 physical pion mass ensembles. For the quark chromo-EDM, unrenormalized results are presented at two lattice spacings, $a=0.12$ and $0.09$ fm, and two pion masses, $M_π= 310$ MeV and $220$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon Axial Form Factors from Clover Fermion on 2+1+1-flavor HISQ Lattice

The nucleon axial form factors -- axial $G_A$, induced pseudoscalar $\widetilde{G}_P$ and pseudoscalar $G_P$ -- have displayed large systematics in lattice QCD calculations. The major symptoms were the violation of the partially conserved axial current (PCAC) relation between the three form factors, and the underestimation of the induced pseudoscalar coupling $g_P^\ast$ and the axial charge radius $r_A$ compared to phenomenological estimates. The small $g_P^\ast$ was a consequence of the failure of the pion-pole dominance (PPD) hypothesis, especially at low $M_π^2$. The small charge radius $r_A$ and the underestimate of $g_A$ were related. The dominant systematic responsible is the lack of inclusion of low-energy ($N π$) states that are not manifest in the multiexponential fit to the nucleon two-point correlator. We show that this low-energy state can be determined from the three-point correlator $\langle N A_4 N \rangle $ with the insertion of the temporal component of the axial current $A_4$ within the nucleon state, ie, the strategy labeled $S_{A4}$ [1]. Including this low-energy state in fits to control excited-state contamination (ESC) gives results for $g_A$, $r_A$, and $g_P^\ast$ that are consistent with experimental/phenomenological values. However, the systematic uncertainties, especially in data at small $Q^2$, are now much larger.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon charges and form factors using clover and HISQ ensembles

We present high statistics ($\mathcal{O}(2\times 10^5)$ measurements) preliminary results on (i) the isovector charges, $g^{u-d}_{A,S,T}$, and form factors, $G^{u-d}_E(Q^2)$, $G^{u-d}_M(Q^2)$, $G^{u-d}_A(Q^2)$, $\widetilde G^{u-d}_P(Q^2)$, $G^{u-d}_P(Q^2)$, on six 2+1-flavor Wilson-clover ensembles generated by the JLab/W&M/LANL/MIT collaboration with lattice parameters given in Table 1. Examples of the impact of using different estimates of the excited state spectra are given for the clover-on-clover data, and as discussed in [1], the biggest difference on including the lower energy (close to $Nπ$ and $Nππ$) states is in the axial channel. (ii) Flavor diagonal axial, tensor and scalar charges, $g^{u,d,s}_{A,S,T}$, are calculated with the clover-on-HISQ formulation using nine 2+1+1-flavor HISQ ensembles generated by the MILC collaboration [2] with lattice parameters given in Table 2. Once finished, the calculations of $g^{u,d,s}_{A,T}$ will update the results given in Refs.[3,4]. The estimates for $g^{u,d,s}_{S}$ and $σ_{Nπ}$ are new. Overall, a large part of the focus is on understanding the excited state contamination (ESC), and the results discussed provide a partial status report on developing defensible analyses strategies that include contributions of possible low-lying excited states to individual nucleon matrix elements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recent results of nucleon structure & matrix element calculations

A review of recent lattice calculations of nucleon structure and matrix elements of operators in nucleons is presented. It primarily covers developments in the calculation of the matrix elements of the scalar, tensor, pseudo-scalar, axial-vector and vector operators in the ground state of neutrons and protons in the isospin symmetric limit. Some preliminary calculations of the electric dipole moment, the gravitational moments and stress-energy distribution, and the magnetic polarizability are briefly described.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Axial Vector Form Factors from Lattice QCD that Satisfy the PCAC Relation

Previous lattice QCD calculations of axial vector and pseudoscalar form factors show significant deviation from the partially conserved axial current (PCAC) relation between them. Since the original correlation functions satisfy PCAC, the observed deviations from the operator identity cast doubt on whether all the systematics in the extraction of form factors from the correlation functions are under control. We identify the problematic systematic as a missed excited state, whose energy as a function of the momentum transfer squared, $Q^2$, is determined from the analysis of the 3-point functions themselves. Its mass is much smaller than those of the excited states previously considered and including it impacts the extraction of all the ground state matrix elements. The form factors extracted using these mass/energy gaps satisfy PCAC and other consistency conditions, and validate the pion-pole dominance hypothesis. We also show that the extraction of the axial charge $g_A$ is very sensitive to the value of the mass gaps of the excited states used and current lattice data do not provide an unambiguous determination of these, unlike the $Q^2 \neq 0$ case. To highlight the differences and improvement between the conventional versus the new analysis strategy, we present a comparison of results obtained on a physical pion mass ensemble at $a\approx 0.0871\,\mathrm{fm}$. With the new strategy, we find $g_A = 1.30(6)$. A very significant improvement over previous lattice results is found for the axial charge radius $r_A = 0.74(6)\,\mathrm{fm}$, extracted using the $z$-expansion to parameterize the $Q^2$ behavior of $G_A(Q^2)$, and $g_P^\ast = 8.06(44)$ obtained using the pion pole-dominance ansatz to fit the $Q^2$ behavior of the induced pseudoscalar form factor $\widetilde{G}_P(Q^2)$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nucleon Electromagnetic Form Factors in the Continuum Limit from 2+1+1-flavor Lattice QCD

Results are presented for the nucleon isovector electromagnetic form factors using 11 ensembles generated by the MILC collaboration using the 2+1+1-flavors HISQ action. They span 4 lattice spacings $a \sim$ 0.06, 0.09, 0.12 and 0.15~fm and 3 values of $M_π\sim 135, 225$ and 315 MeV. High-statistics estimates are used to perform a simultaneous extrapolation in the lattice spacing, lattice volume and light-quark masses. The $Q^2$ dependence over the range 0.05-1.4 ${\rm GeV}^2$ is investigated using both the $z$-expansion and the dipole form. Final $z$-expansion estimates for the isovector r.m.s. radius are $r_E = 0.769(27)(30)$ fm $r_M = 0.671(48)(76)$ fm and $μ^{p-n} = 3.939(86)(138)$ Bohr magneton. The first error is the combined uncertainty from the leading-order analysis, and the second is an estimate of the additional uncertainty due to using the leading order chiral-continuum-finite-volume fits. The dipole estimates, $r_E = 0.765(11)(8)$ fm, $r_M = 0.704(21)(29)$ fm and $μ^{p-n} = 3.975(84)(125)$, are consistent with those from the $z$-expansion but with smaller errors. Our analysis highlights three points. First, all data from the eleven ensembles and existing lattice data on, or close to, physical mass ensembles from other collaborations collapses more clearly onto a single curve when plotted versus $Q^2/M_N^2$ as compared to $Q^2$ with the scale set by quantities other than $M_N$. The difference between these two analyses is indicative of discretization errors, some of which presumably cancel when the data are plotted versus $Q^2/M_N^2$. Second, the size of the remaining deviation of this common curve from the Kelly curve is small and can be accounted for by statistical and possible systematic uncertainties. Third, to improve lattice estimates, high statistics data for $Q^2 < 0.1$ ${\rm GeV}^2$ are needed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Algorithms for Disconnected Diagrams in Lattice QCD

Computing disconnected diagrams in Lattice QCD (operator insertion in a quark loop) entails the computationally demanding problem of taking the trace of the all to all quark propagator. We first outline the basic algorithm used to compute a quark loop as well as improvements to this method. Then, we motivate and introduce an algorithm based on the synergy between hierarchical probing and singular value deflation. We present results for the chiral condensate using a 2+1-flavor clover ensemble and compare estimates of the nucleon charges with the basic algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Estimation of matrix trace using machine learning

We present a new trace estimator of the matrix whose explicit form is not given but its matrix multiplication to a vector is available. The form of the estimator is similar to the Hutchison stochastic trace estimator, but instead of the random noise vectors in Hutchison estimator, we use small number of probing vectors determined by machine learning. Evaluation of the quality of estimates and bias correction are discussed. An unbiased estimator is proposed for the calculation of the expectation value of a function of traces. In the numerical experiments with random matrices, it is shown that the precision of trace estimates with $\mathcal{O}(10)$ probing vectors determined by the machine learning is similar to that with $\mathcal{O}(10000)$ random noise vectors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Iso-vector and Iso-scalar Tensor Charges of the Nucleon from Lattice QCD

We present results for the iso-vector and flavor diagonal tensor charges $g^{u-d}_T$, $g^{u}_T$, $g^{d}_T$, and $g^{s}_T$ needed to probe novel tensor interactions at the TeV scale in neutron and nuclear $β$-decays and the contribution of the quark electric dipole moment (EDM) to the neutron EDM. The lattice QCD calculations were done using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. These ensembles span three lattice spacings $a \approx 0.06, 0.09$ and $0.12 $ fm and three quark masses corresponding to the pion masses $M_π\approx 130, 220$ and $310 $ MeV. Using estimates from these ensembles, we quantify all systematic uncertainties and perform a simultaneous extrapolation in the lattice spacing, volume and light quark masses for the connected contributions. The final estimates of the connected nucleon (proton) tensor charge for the iso-vector combination is $g_T^{u-d} = 1.020(76) $ in the $\bar{\text{MS}}$ scheme at $2$ GeV. The additional disconnected quark loop contributions needed for the flavor-diagonal matrix elements are calculated using a stochastic estimator employing the truncated solver method with the all-mode-averaging technique. We find that the size of the disconnected contribution is smaller than the statistical error in the connected contribution. This allows us to bound the disconnected contribution and include it as an additional uncertainty in the flavor-diagonal charges. After a continuum extrapolation, we find $g_T^{u} = 0.774(66) $, $g_T^{d} = -0.233(28) $ and $g_T^{u+d} = 0.541(67) $. The strangeness tensor charge, that can make a significant contribution to the neutron EDM due to the large ratio $m_s/m_{u,d}$, is $g_T^{s}=0.008(9)$ in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD

We present Lattice QCD results on the neutron tensor charges including, for the first time, a simultaneous extrapolation in the lattice spacing, volume, and light quark masses to the physical point in the continuum limit. We find that the "disconnected" contribution is smaller than the statistical error in the "connected" contribution. Our estimates in the $\bar{\text{MS}}$ scheme at $2$ GeV, including all systematics, are $g_T^{d-u}=1.020(76)$, $g_T^d = 0.774(66)$, $g_T^u = - 0.233(28)$, and $g_T^s = 0.008(9)$. The flavor diagonal charges determine the size of the neutron electric dipole moment (EDM) induced by quark EDMs that are generated in many new scenarios of CP-violation beyond the Standard Model (BSM). We use our results to derive model-independent bounds on the EDMs of light quarks and update the EDM phenomenology in split Supersymmetry with gaugino mass unification, finding a stringent upper bound of $d_n < 4 \times 10^{-28} \, e$ cm for the neutron EDM in this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment from quark Chromoelectric Dipole Moment

The connection between a regularization-independent symmetric momentum substraction (RI-$\tilde{\rm S}$MOM) and the $\overline{\rm MS}$ scheme for the quark chromo EDM operators is discussed. A method for evaluating the neutron EDM from quark chromoEDM is described. A preliminary study of the signal in the matrix element using clover quarks on a highly improved staggered quark (HISQ) ensemble is shown.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nucleon Charges, Form-factors and Neutron EDM

We present an update of our analysis of statistical and systematic errors in the calculation of iso-vector scalar, axial and tensor charges of the nucleon. The calculations are done using $N_f=2+1+1$ flavor HISQ ensembles generated by the MILC Collaboration at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV); and clover fermions for calculating the correlation functions, i.e., we use a clover-on-HISQ lattice formulation. The all-mode-averaging method allows us to increase the statistics by a factor of eight for the same computational cost leading to a better understanding of and control over excited state contamination. Our current results, after extrapolation to the continuum limit and physical pion mass are $g_A^{u-d} = 1.21(3)$, $g_T^{u-d} = 1.005(59)$ and $g_S^{u-d} = 0.95(12) $. Further checks of control over all systematic errors, especially in $g_A^{u-d}$, are still being performed. Using results for the flavor-diagonal charges, $g_T^{u} $, $g_T^{d} $ and $g_T^{s}$, we analyze contributions of the quark electric dipole moment to the neutron EDM and the consequences for split SUSY model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quark Chromoelectric Dipole Moment Contribution to the Neutron Electric Dipole Moment

The quark chromo-electric dipole moment operator and the pseudo-scalar fermion bilinear with which it mixes under renormalization can both be included in a calculation of the electromagnetic form factor of the nucleon using the Schwinger source method. A preliminary calculation of these operators using clover quarks on HISQ lattices generated by MILC collaboration will be presented showing the quality of the signal in the correlators necessary for calculating the neutron electric dipole moment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disconnected Quark Loop Contributions to Nucleon Structure

We calculate the disconnected contribution to isoscalar nucleon charges for scalar, axial and tensor channels of light and strange quarks. The calculation has been done with the Clover valence quarks on the MILC $N_f=2+1+1$ HISQ lattices whose light quark masses corresponding to the pion masses of 305 MeV and 217 MeV at $a \approx 0.12$ fm and 312 MeV at $a \approx 0.09$ fm. All-mode-averaging technique is used for the evaluation two-point functions. Disconnected quark loops are estimated by using the truncated solver method with Gaussian random noise sources. Contamination from the excited states is removed by fitting the results of various source-sink separations and operator insertions to the formula including up to the first excited state, simultaneously.

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

Precision calculations of nucleon charges $g_A$, $g_S$, $g_T$

We present a detailed analysis of statistical and systematic errors in the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These analyses are being done on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. Using ensembles at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV) we find that the estimates of the tensor charge are stable and it can be extracted with $5\%$ precision with O(10,000) measurements. We also find that higher statistics are needed to resolve the various uncertainties in the calculation of $g_A$ and improve the signal in $g_S$, which with present data has large errors. A brief status report on the mixing and renormalization of novel operators contributing to nEDM is also given.

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

Nucleon Charges and Electromagnetic Form Factors from 2+1+1-Flavor Lattice QCD

We present lattice-QCD results on the nucleon isovector axial, scalar and tensor charges, the isovector electromagnetic Dirac and Pauli form factors, and the connected parts of the isoscalar charges. The calculations have been done using two ensembles of HISQ lattices generated by the MILC Collaboration with 2+1+1 dynamical flavors at a lattice spacing of 0.12 fm and with light-quark masses corresponding to pions with masses 310 and 220 MeV. We perform a systematic study including excited-state degrees of freedom and examine the dependence of the extracted nucleon matrix elements on source-sink separation. This study demonstrates with high-statistics data that including excited-state contributions and generating data at multiple separations is necessary to remove contamination that would otherwise lead to systematic error. We also determine the renormalization constants of the associated quark bilinear operators in the RI-sMOM scheme and make comparisons of our renormalized results with previous dynamical-lattice calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Toolkit for staggered $ΔS=2$ matrix elements

A recent numerical lattice calculation of the kaon mixing matrix elements of general $ΔS=2$ four-fermion operators using staggered fermions relied on two auxiliary theoretical calculations. Here we describe the methodology and present the results of these two calculations. The first concerns one-loop matching coefficients between staggered lattice operators and the corresponding continuum operators. Previous calculations with staggered fermions have used a non-standard regularization scheme for the continuum operators, and here we provide the additional matching factors needed to connect to the standard regularization scheme. This is the scheme in which two-loop anomalous dimensions are known. We also observe that all previous calculations of this operator matching using staggered fermions have overlooked one matching step in the continuum. This extra step turns out to have no impact on three of the five operators (including that relevant for $B_K$), but does affect the other two operators. The second auxiliary calculation concerns the two-loop renormalization group (RG) evolution equations for the $B$-parameters of the $ΔS=2$ operators. For one pair of operators, the standard analytic solution to the two-loop RG equations fails due to a spurious singularity introduced by the approximations made in the calculation. We give a non-singular expression derived using analytic continuation, and check the result using a numerical solution to the RG equations. We also describe the RG evolution for "golden" combinations of $B$-parameters, and give numerical results for RG evolution matrices needed in the companion lattice calculation.

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

Flavored pions and kaons at next-to-leading order in mixed-action staggered chiral perturbation theory

Different versions of improved staggered fermions can be used as valence quarks to reduce discretization effects in lattice QCD calculations while increasing statistics on existing staggered gauge ensembles. Such mixed-action simulations can be used to improve determinations of light quark masses, Gasser-Leutwyler couplings, decay constants, and other parameters relevant to particle phenomenology. We recall the generalization of ordinary, unmixed staggered chiral perturbation theory required to describe data from lattice calculations with a mixed action such as with HYP staggered valence quarks and asqtad sea quarks. We calculate the next-to-leading order loop diagrams contributing to the masses and decay constants of the flavored pseudo-Goldstone bosons of all tastes and here report results for the decay constants and valence-valence masses.

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

Non-perturbative Renormalization of Bilinear Operators with Improved Staggered Quarks

We present renormalization factors for the bilinear operators obtained using the non-perturbative renormalization method (NPR) in the RI-MOM scheme with improved staggered fermions on the MILC asqtad lattices ($N_f = 2+1$). We use the MILC coarse ensembles with $20^3 \times 64$ geometry and $am_{\ell}/am_s = 0.01/0.05$. We obtain the wave function renormalization factor $Z_q$ from the conserved vector current and the mass renormalization factor $Z_m$ from the scalar bilinear operator. We also present preliminary results of renormalization factors for other bilinear operators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Recent progress in staggered chiral perturbation theory

We present a review on recent progress in staggered chiral perturbation theory (SChPT). In the last decade, the scope of the application of SChPT has been extended beyond the level of calibration into the region of prediction with high precision. SChPT becomes an essential tool to do the data analysis reliably for physical observables calculated using improved staggered fermions. Here, we focus on the following examples: pion spectrum, pion decay constants, $\varepsilon_K$, and $π-π$ scattering amplitudes. In each subject, we review the recent progress and future prospects.

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

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

Non-Perturbative Renormalization for Staggered Fermions (Self-energy Analysis)

We present preliminary results of data analysis for the non-perturbative renormalization (NPR) on the self-energy of the quark propagators calculated using HYP improved staggered fermions on the MILC asqtad lattices. We use the momentum source to generate the quark propagators. In principle, using the vector projection operator of $(\bar{\bar{γ_μ\otimes 1}})$ and the scalar projection operator $(\bar{\bar{1 \otimes 1}})$, we should be able to obtain the wave function renormalization factor $Z_q'$ and the mass renormalization factor $Z_q \cdot Z_m$. Using the MILC coarse lattice, we obtain a preliminary but reasonable estimate of $Z_q'$ and $Z_q \cdot Z_m$ from the data analysis on the self-energy.

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

Taste non-Goldstone pion decay constants in staggered chiral perturbation theory

We calculate the next-to-leading order axial current decay constants of taste non-Goldstone pions and kaons in staggered chiral perturbation theory. This is an extension of the taste Goldstone decay constants calculation to that of the non-Goldstone tastes. We present results for the partially quenched case in the SU(3) and SU(2) staggered chiral perturbation theories and discuss the difference between the taste Goldstone and non-Goldstone cases.

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

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