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

Weonjong Lee

Weonjong Lee 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-thDistributed, Parallel, and Cluster Computingphysics.comp-phPerformancephysics.data-an

Catalog footprint

What is connected

83works
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

83 published item(s)

preprint2022arXiv

2021 Update on $\varepsilon_K$ with lattice QCD inputs

We present recent updates for $\varepsilon_K$ determined directly from the standard model (SM) with lattice QCD inputs such as $\hat{B}_K$, $|V_{cb}|$, $|V_{us}|$, $ξ_0$, $ξ_2$, $ξ_\text{LD}$, $f_K$, and $m_c$. We find that the standard model with exclusive $|V_{cb}|$ and other lattice QCD inputs describes only 66\% of the experimental value of $|\varepsilon_K|$ and does not explain its remaining 34\%, which leads to a strong tension in $|\varepsilon_K|$ at the $4.5σ\sim 3.7σ$ level between the SM theory and experiment. We also find that this tension disappears when we use the inclusive value of $|V_{cb}|$ obtained using the heavy quark expansion based on the QCD sum rule approach.

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
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
preprint2016arXiv

Complete Monopole Dominance of the Yang-Mills Confining Potential

We continue our investigation of quark confinement using a particular variant of the Cho-Duan-Ge gauge independent Abelian decomposition. The decomposition splits the gauge field into a restricted Abelian part and a coloured part in a way that preserves gauge covariance. The restricted part of the gauge field can be divided into a Maxwell term and a topological term. Previously, we showed that by a particular choice of this decomposition we could fully describe the confining potential using only the restricted gauge field. We proposed that various topological objects (a form of magnetic monopole) could arise in the restricted field which would drive confinement. Our mechanism does not explicitly refer to a dual Meissner effect, nor does it use centre vortices. We did not need to gauge fix or introduce any new dynamical fields. We show that if we do gauge fix as well as performing the Abelian decomposition then it is possible to ensure that the topological part of the restricted field fully accounts for the confining potential. Our relationship is exact: there is no approximation or model involved. This isolates the objects responsible for confinement from non-confining contributions to the gauge field, allowing a direct search for our proposed topological objects. Using numerical studies in SU(2), we confirm that our proposed monopoles are present in the field, and the winding number associated with these monopoles is a key factor driving quark confinement. In SU(2), our monopoles are described by two parameters. We show that it is possible to re-parametrise the Yang Mills action and the functional integration measure in terms of these variables (plus the necessary additional parameters). We can thus treat the monopoles as dynamical variables in the functional integral. This might be the first step in a future analytical computation to complement our numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heavy-heavy current improvement for calculation of $\bar{B}\rightarrow D^{(*)}\ell \barν$ semi-leptonic form factors using the Oktay-Kronfeld action

Lattice calculations of the form factors for $\bar{B}\to D^{(*)}\ell\barν$ decays can be used to extract the CKM matrix element $|V_{cb}|$. The Oktay-Kronfeld action is a highly improved version of the Fermilab action, which systematically reduces heavy quark discretization effects through $\mathcal{O}(λ^3)$ in HQET power counting, for heavy-light meson quantities. To calculate $\bar{B}\rightarrow D^{(*)}\ell \barν$ semi-leptonic form factors using Oktay-Kronfeld heavy quarks, we need to improve the heavy quark currents to the same level. We report our progress in calculating the improvement coefficients for currents composed of bottom and charm quarks. Our results presented in this paper are preliminary.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Status report on $\varepsilon_K$ with lattice QCD inputs

We report the current status of $\varepsilon_K$, the indirect CP violation parameter in the neutral kaon system, evaluated using the lattice QCD inputs. We use lattice QCD to fix $\hat{B}_K$, $ξ_0$, $ξ_2$, $|V_{us}|$, $m_c(m_c)$, and $|V_{cb}|$. Since Lattice 2015, FLAG updated $\hat{B}_K$, exclusive $V_{cb}$ has been updated with new lattice data in the $\bar{B}\to D\ellν$ decay channel, and RBC-UKQCD has updated $ξ_0$ and $ξ_2$. Our preliminary results show that the standard model evaluation of $\varepsilon_K$ with exclusive $|V_{cb}|$ (lattice QCD inputs) has $3.2σ$ tension with the experimental value, while that of $\varepsilon_K$ with inclusive $|V_{cb}|$ (heavy quark expansion) shows no tension.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning the hopping parameter in the Oktay-Kronfeld action for charm and bottom quarks on a MILC HISQ ensemble

The first step in the calculation of semi-leptonic form factors in the decay of heavy mesons is the tuning of the hopping parameter $κ$ for the charm and bottom quark masses. Results for the Oktay-Kronfeld (OK) action are presented for one $N_f=2+1+1$ HISQ ensemble generated by the MILC collaboration at $a\approx 0.12\,\mathrm{fm}$ and $M_π\approx 310$ MeV. Estimates of hyperfine splitting of heavy-light and heavy-heavy mesons are presented and the inconsistency parameter is evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Update on Heavy-Meson Spectrum Tests of the Oktay--Kronfeld Action

We present updated results of a numerical improvement test with heavy-meson spectrum for the Oktay--Kronfeld (OK) action. The OK action is an extension of the Fermilab improvement program for massive Wilson fermions including all dimension-six and some dimension-seven bilinear terms. Improvement terms are truncated by HQET power counting at $\mathrm{O}(Λ^3/m_Q^3)$ for heavy-light systems, and by NRQCD power counting at $\mathrm{O}(v^6)$ for quarkonium. They suffice for tree-level matching to QCD to the given order in the power-counting schemes. To assess the improvement, we generate new data with the OK and Fermilab action that covers both charm and bottom quark mass regions on a MILC coarse $(a \approx 0.12~\text{fm})$ $2+1$ flavor, asqtad-staggered ensemble. We update the analyses of the inconsistency quantity and the hyperfine splittings for the rest and kinetic masses. With one exception, the results clearly show that the OK action significantly reduces heavy-quark discretization effects in the meson spectrum. The exception is the hyperfine splitting of the heavy-light system near the $B_s$ meson mass, where statistics are too low to draw a firm conclusion, despite promising results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Current status of $\varepsilon_K$ with lattice QCD inputs

We present the Standard Model evaluation of the indirect CP violation parameter $\varepsilon_K$ using inputs determined from lattice QCD together with experiment: $|V_{us}|$, $|V_{cb}|$, $ξ_0$, and $\hat{B}_K$. We use the Wolfenstein parametrization ($|V_{cb}|$, $λ$, $\barρ$, $\barη$) for the CKM matrix elements. For the central value, we take the angle-only fit of the UTfit collaboration, and use $|V_{us}|$ from the $K_{\ell 3}$ and $K_{μ2}$ decays as an independent input to fix $λ$. For the error estimate, we use results of the global unitarity triangle fits from the CKMfitter and UTfit collaborations. We find that the Standard Model (SM) prediction of $\varepsilon_K$ with exclusive $V_{cb}$ (lattice QCD results) is lower than the experimental value by $3.6(2)σ$. However, with inclusive $V_{cb}$ (results of the heavy quark expansion), the tension between the SM prediction of $\varepsilon_K$ and its experimental value disappears.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Determination of $\varepsilon_K$ using lattice QCD inputs

We present results for the indirect CP violation parameter $\varepsilon_K$ determined directly from the standard model using lattice QCD to fix the inputs $\hat{B}_K$, $ξ_0$, $|V_{us}|$, and $|V_{cb}|$. We use the FLAG and SWME results for $\hat{B}_K$. We use the RBC-UKQCD result for $ξ_0$ determined using the experimental value of $\varepsilon'/\varepsilon$ and the lattice result of $\mathrm{Im}\,A_2$. To set the Wolfenstein parameter $λ$, we use $|V_{us}|$, which is determined from $K_{\ell3}$ and $K_{\mu2}$ decays combined with lattice evaluations of the $K \to π\ell ν$ vector form factor and $f_K$. To set the Wolfenstein parameter $A$, we use the FNAL/MILC results for $|V_{cb}|$, which are determined from the exclusive decay $\bar{B} \to D^* \ell \barν$ and the axial form factor at zero recoil. We also use the inclusive $|V_{cb}|$ obtained using the heavy quark expansion based on QCD sum rules and the OPE. We compare the results with those for exclusive $|V_{cb}|$. We find that the standard model prediction of $\varepsilon_K$ with exclusive $|V_{cb}|$ (lattice QCD results) is lower than the experimental value by 3.4$σ$. However, we observe no tension in $\varepsilon_K$ determined from inclusive $|V_{cb}|$.

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

Non-perturbative Renormalization of Bilinear Operators on Fine Lattice

We present results of the wave function renormalization factor $Z_q$ and mass renormalization factor $Z_m$ obtained using non-perturbative renormalization (NPR) method in the RI-MOM scheme with HYP improved staggered quarks. We use fine ensembles of MILC asqtad lattices ($N_f = 2+1$) with $28^3 \times 96$ geometry, $a \approx 0.09$\,fm, and $am_\ell/am_s = 0.0062/0.031 $. We also study on scalability of $Z_q$ and $Z_m$ by comparing the results on the coarse and fine ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonperturbative Renormalization in the RI-SMOM Scheme and Gribov Uncertainty in the RI-MOM Scheme for Staggered Bilinears

We present results of renormalization factors for bilinear operators obtained using the nonperturbative renormalization method (NPR) in the RI-SMOM schemes. The operators are constructed using HYP staggered quarks on the MILC asqtad lattice ($N_f=2+1$). We compare results in the RI-SMOM schemes with those in the RI-MOM scheme for the $V\otimes S$ and $S\otimes S$ operators. Since we use Landau gauge fixing, we study the effect of Gribov ambiguity on the wave function renormalization $Z_q$ in the RI-MOM scheme. We find that the Gribov uncertainty is negligibly small for $Z_q$ in the RI-MOM scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Performance of GTX Titan X GPUs and Code Optimization

Recently Nvidia has released a new GPU model: GTX Titan X (TX) in a linage of the Maxwell architecture. We use our conjugate gradient code and non-perturbative renormalization code to measure the performance of TX. The results are compared with those of GTX Titan Black (TB) in a lineage of the Kepler architecture. We observe a significant gain in the single and double precision calculations much greater than the theoretical expectation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Standard Model evaluation of $\varepsilon_K$ using lattice QCD inputs for $\hat{B}_K$ and $V_{cb}$

We report the Standard Model evaluation of the indirect CP violation parameter $\varepsilon_K$ using inputs determined from lattice QCD: the kaon bag parameter $\hat{B}_K$, $ξ_0$, $|V_{us}|$ from the $K_{\ell 3}$ and $K_{μ2}$ decays, and $|V_{cb}|$ from the axial current form factor for the exclusive decay $\bar{B} \to D^* \ell \barν$ at zero-recoil. The theoretical expression for $\varepsilon_K$ is thoroughly reviewed to give an estimate of the size of the neglected corrections, including long distance effects. The Wolfenstein parametrization $(|V_{cb}|, λ, \barρ, \barη)$ is adopted for CKM matrix elements which enter through the short distance contribution of the box diagrams. For the central value, we take the Unitarity Triangle apex $(\barρ, \barη)$ from the angle-only fit of the UTfit collaboration and use $V_{us}$ as an independent input to fix $λ$. We find that the Standard Model prediction of $\varepsilon_K$ with exclusive $V_{cb}$ (lattice QCD results) is lower than the experimental value by $3.4σ$. However, with inclusive $V_{cb}$ (results of the heavy quark expansion), there is no gap between the Standard Model prediction of $\varepsilon_K$ and its experimental value. For the calculation of $\varepsilon_K$, we perform the renormalization group running to obtain $η_{cc}$ at next-to-next-to-leading-order; we find $η_{cc}^\mathrm{NNLO}=1.72(27)$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The conserved axial current in the presence of multiple chiral symmetries

In response to a recent work by Mandula, we investigate whether there are any ambiguities in the expression for the pion mass resulting from multiple chiral symmetries. If the conserved current for Ginsparg Wilson chiral symmetries is calculated in the usual way, different expressions of the chiral symmetry lead to different currents. This implies an ambiguity in the definition of the pion and pion decay constant for all Ginsparg-Wilson expressions of the Dirac operator, including the overlap operator on the lattice (although all these currents would have the same continuum limit). We use a renormalisation group mapping procedure to consider local chiral symmetry transformations for a continuum Ginsparg-Wilson "Dirac-operator." We find that this naturally leads to an expression for the conserved current which is independent of which of the Ginsparg-Wilson symmetries is chosen. We recover the standard expressions for the massive Dirac operator, propagator, and chiral condensate. Our main conclusion is that, when the currents are properly constructed and consistently applied, no observable depends on which Mandula symmetry is used; at least in these continuum Ginsparg-Wilson theories. We will consider whether the same argument applies to lattice theories in a subsequent paper.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The eigSUMR inverter for overlap fermion

We discuss the usage and applicability of deflation methods for the overlap lattice Dirac operator, focussing on calculating the eigenvalues using a method similar to the eigCG algorithm used for other Dirac operators. The overlap operator, which contains several theoretical advantages over other formulations of lattice Quantum Chromodynamics, is more computationally expensive because it requires the computation of the matrix sign function. The principle change made compared to deflation methods for other formulations of lattice QCD is that it is necessary for best performance to tune or relax the accuracy of the matrix sign function as the computation proceeds. We adapt the eigCG algorithm for two inversion algorithms for overlap fermions, GMRESR(relCG) and GMRESR(relSUMR). Before deflation, the rate of convergence of these routines in terms of iterations is similar, but, since the Shifted Unitary Minimal Residual (SUMR) algorithm only requires one call to the matrix sign function compared to the two calls required for Conjugate Gradient (CG), SUMR is usually preferred for single inversions of the Dirac operator. We construct bounds for the required accuracy of the matrix sign function during the eigenvalue calculation. For the SUMR algorithm, we use a Galerkin projection to perform the deflation; while for the CG algorithm, we are able to use a considerably superior spectral pre-conditioner. The superior performance of the spectral preconditioner, and its need for less accurate eigenvalues, almost erodes SUMR's advantage over CG as an inversion algorithm. We see factor of three gains for the inversion algorithm from the deflation on our small test lattices. There is, however, a significant cost in the eigenvalue calculation because we cannot relax the accuracy of the matrix sign function as aggressively when calculating the eigenvalues as we do while performing the inversions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Static Quark Potential from the Gauge Independent Abelian Decomposition

We investigate the relationship between colour confinement and the gauge independent Cho-Duan-Ge Abelian decomposition. The decomposition is defined in terms of a colour field $n$; the principle novelty of our study is that we have defined this field in terms of the eigenvectors of the Wilson Loop. This establishes an equivalence between the path ordered integral of the non-Abelian gauge fields with an integral over an Abelian restricted gauge field which is tractable both theoretically and numerically in lattice QCD. We circumvent path ordering without needing an additional path integral. By using Stokes' theorem, we can compute the Wilson Loop in terms of a surface integral over a restricted field strength, and show that the restricted field strength may be dominated by certain structures, which occur when one of the quantities parametrising the colour field $n$ winds itself around a non-analyticity in the colour field. If they exist, these structures will lead to a area law scaling for the Wilson Loop and provide a mechanism for quark confinement. We search for these structures in quenched lattice QCD. We perform the Abelian decomposition, and compare the electric and magnetic fields with the patterns expected theoretically. We find that the restricted field strength is dominated by objects which may be peaks a single lattice spacing in size or extended string-like lines of electromagnetic flux. The objects are not isolated monopoles, as they generate electric fields in addition to magnetic fields, and the fields are not spherically symmetric, but may be either caused by a monopole/anti-monopole condensate, some other types of topological objects or a combination of these. Removing these peaks removes the area law scaling of the string tension, suggesting that they are responsible for confinement.

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

Code Optimization on Kepler GPUs and Xeon Phi

Kepler GTX Titan Black and Kepler Tesla K40 are still the best GPUs for high performance computing, although Maxwell GPUs such as GTX 980 are available in the market. Hence, we measure the performance of our lattice QCD codes using the Kepler GPUs. We also upgrade our code to use the latest CPS (Columbia Physics System) library along with the most recent QUDA (QCD CUDA) library for lattice QCD. These new libraries improve the performance of our conjugate gradient (CG) inverter so that it runs twice faster than before. We also investigate the performance of Xeon Phi 7120P coprocessor. It has similar computing power with the Kepler GPUs in principle. However, its performance for our CG code is significantly inferior to that of the GTX Titan Black GPUs at present.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Confinement, the Abelian Decomposition, and the Contribution of Topology to the Static Quark Potential

In the past few years, we have presented a new way of considering quark confinement. Through a careful choice of a Cho-Duan-Ge Abelian Decomposition, we can construct the QCD Wilson Loop in terms of an Abelian restricted field. The relationship between the QCD and restricted string tensions is exact; and we do not need to gauge fix, apply any path ordering of gauge links, or additional path integrals. This hints at why mesons are colour neutral. Furthermore, the Abelian restricted field contains two parts: a Maxwell term, and a topological term. The topological term can describe magnetic monopoles and other topological objects, which can be studied both numerically and theoretically. By examining the topological part of the restricted field strength we have found evidence suggesting that these objects, which will contribute to confinement if present, are indeed there. Previous studies have used simplifications, breaking the exact relationship between the restricted and QCD string tensions, but it was found that the topological term dominated the restricted string tension. Here we remove those simplifications, and show that the Abelian restricted field does indeed fully explain confinement. However, our results for how much of the restricted string tension arises from the topological objects show strong dependence on the lattice spacing and level of smearing, so we are not yet able to draw a definitive conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Current status of $\varepsilon_K$ calculated with lattice QCD inputs

We present results for $\varepsilon_K$, the indirect CP violation parameter, calculated in the Standard Model using inputs from lattice QCD: the kaon bag parameter $\hat{B}_K$, and the CKM matrix element $V_{cb}$ from the axial current form factor for the exclusive decay $\bar{B}\to D^*\ell\barν$ at zero-recoil. In addition, we take the coordinates of the unitarity triangle apex $(\barρ,\barη)$ from the angle-only fit of the UTfit Collaboration and use $V_{us}$ to fix $λ$. In order to estimate the systematic error, we also use Wolfenstein parameters from the CKMfitter and UTfit. We find a $3.3(2)σ$ difference between $\varepsilon_K$ and experiment with exclusive $V_{cb}$. We report details of this preliminary result.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heavy-Meson Spectrum Tests of the Oktay--Kronfeld Action

We present heavy-meson spectrum results obtained using the Oktay--Kronfeld (OK) action on MILC asqtad lattices. The OK action was designed to improve the heavy-quark action of the Fermilab formulation, such that heavy-quark discretization errors are reduced. The OK action includes dimension-6 and -7 operators necessary for tree-level matching to QCD through order $\mathrm{O}(Λ^3/m_Q^3)$ for heavy-light mesons and $\mathrm{O}(v^6)$ for quarkonium, or, equivalently, through $\mathrm{O}(a^2)$ with some $\mathrm{O}(a^3)$ terms with Symanzik power counting. To assess the improvement, we extend previous numerical tests with heavy-meson masses by analyzing data generated on a finer ($a \approx 0.12\;$fm) lattice with the correct tadpole factors for the $c_5$ term in the action. We update the analyses of the inconsistency parameter and the hyperfine splittings for the rest and kinetic masses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved currents for $\bar{B}\to D^{(*)}\ell\barν$ form factors from Oktay-Kronfeld heavy quarks

The CKM matrix element $|V_{cb}|$ can be extracted by combining experimentally determined branching fractions for $\bar{B}\to D^{(*)}\ell\barν$ decays with form factors from the lattice. While successful, the precision of this approach has been limited by heavy-quark discretization effects. An improved version of the Fermilab action, the Oktay-Kronfeld action, can be used to reduce heavy-quark discretization effects in calculations performed at the physical bottom and charm quark masses. Treating charm and bottom quarks as massive, we are carrying out improvement of the flavor-changing currents through third order in the momentum (HQET) expansion.

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

Non-perturbative Renormalization of Four-Fermion Operators Relevant to $B_K$ with Staggered Quarks

We present preliminary results of matching factors of the four-fermion operators relevant to $B_K$, which are obtained using the non-perturbative renormalization (NPR) method in the RI-MOM scheme with HYP-smeared improved staggered fermions. We use the MILC asqtad coarse ($a \cong 0.12\,$fm) ensembles with $20^3 \times 64$ geometry and $am_{\ell}/am_s = 0.01/0.05$. We compare NPR results with those of one-loop perturbative matching.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The static quark potential from the gauge invariant Abelian decomposition

We investigate the relationship between colour confinement and topological structures derived from the gauge invariant Abelian (Cho-Duan-Ge) decomposition. This Abelian decomposition is made imposing an isometry on a colour field $n$ which selects the Abelian direction; the principle novelty of our study is that we have defined this field in terms of the eigenvectors of the Wilson Loop. This allows us to establish an equivalence between the path ordered integral of the non-Abelian gauge fields with an integral over an Abelian restricted gauge field which is tractable both theoretically and numerically in lattice QCD. By using Stokes' theorem, we can relate the Wilson Loop in terms of a surface integral over a restricted field strength, and show that the restricted field strength may be dominated by topological structures, which occur when one of the parameters parametrising the colour field $n$ winds itself around a non-analyticity in the colour field. If they exist, these objects will lead to an area law scaling for the Wilson Loop and provide a mechanism for quark confinement. We search for these structures in quenched lattice QCD. We perform the Abelian decomposition, and find that the restricted field strength is dominated by peaks on the lattice. Wilson Loops containing these peaks show a stronger area-Law and thus provide the dominant contribution to the string tension.

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

Confinement From The Gauge Invariant Abelian Decomposition

A common approach while considering confinement is to study the dominance of an Abelian subgroup of the SU(3) gauge Links. A good way to find the Abelian component of the field is through the Cho-Guan-De gauge invariant Abelian Decomposition, which uses a carefully chosen direction vector $n$ to split the gauge field into an Abelian restricted field and a remnant coloured field. The restricted field can be further subdivided into topological and non-topological terms. We show that there is a choice of $n$ which allows us to exactly represent the Wilson Loop of full QCD as a function of only the restricted Abelian field without requiring any path ordering or additional path integrals. We present numerical evidence showing that the topological part of the restricted field dominates the string tension. We also show that $n$ contains certain topological objects, which, if they exist, will be at least partially responsible for confinement. These leave distinctive patterns in the restricted field strength, and we search for these structures in quenched lattice QCD.

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

Optimization of the Oktay-Kronfeld Action Conjugate Gradient Inverter

Improving the Fermilab action to third order in heavy quark effective theory yields the Oktay-Kronfeld action, a promising candidate for precise calculations of the spectra of heavy quark systems and weak matrix elements relevant to searches for new physics. We have optimized the bi-stabilized conjugate gradient inverter in the SciDAC QOPQDP library and are developing a GPU code. The action is rewritten and the needed gauge-link combinations are precalculated. In tests with a MILC coarse lattice, this procedure accelerates the inverter by a factor of four. The remaining floating-point operations are mostly simple matrix multiplications between gauge links and fermion vectors, which we accelerate by more than an order of magnitude by using CUDA. Further gains could be achieved by using QUDA.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Performance of Kepler GTX Titan GPUs and Xeon Phi System

NVIDIA's new architecture, Kepler improves GPU's performance significantly with the new streaming multiprocessor SMX. Along with the performance, NVIDIA has also introduced many new technologies such as direct parallelism, hyper-Q and GPU Direct with RDMA. Apart from other usual GPUs, NVIDIA also released another Kepler 'GeForce' GPU named GTX Titan. GeForce GTX Titan is not only good for gaming but also good for high performance computing with CUDA. Nevertheless, it is remarkably cheaper than Kepler Tesla GPUs. We investigate the performance of GTX Titan and find out how to optimize a CUDA code appropriately for it. Meanwhile, Intel has launched its new many integrated core (MIC) system, Xeon Phi. A Xeon Phi coprocessor could provide similar performance with NVIDIA Kepler GPUs theoretically but, in reality, it turns out that its performance is significantly inferior to GTX Titan.

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

Chiral extrapolation of matrix elements of BSM kaon operators

Models of new physics induce K_0-K_0bar mixing through operators having Dirac structures other than the "left-left" form of the Standard Model. To carry out the chiral-continuum extrapolation of results from numerical simulations, one needs to know the quark mass and lattice spacing dependence of the corresponding B-parameters in the partially quenched theory at least at next-to-leading order. For simulations using staggered fermions (such as that we are doing with HYP-smeared valence fermions on the MILC asqtad lattices) one must determine this dependence using staggered chiral perturbation theory (SChPT). We have calculated the required dependence in both SU(3) and SU(2) SChPT, working at next-to-leading order, and we give here an overview of the methodology and results. The SU(3) SChPT result turns out to be much simpler than that for the Standard Model B_K operator, due to the absence of chiral suppression for the new operators. The SU(2) SChPT result turns out to be closely related to that for B_K: the chiral logarithms are identical, up to an operator-dependent sign. Our results are also useful for fermions with chiral symmetry as they provide, in the continuum limit, the partially quenched generalization of existing continuum results.

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

Current Status of Indirect CP Violation in Neutral Kaon System

In the standard model (SM), the CP violation is introduced through a single phase in the CKM matrix. The neutral kaon system is one of the most precise channels to test how the SM theory describes the experiment data such as $ε_K$ accurately. The indirect CP violation is parametrized into $ε_{K}$, which can be calculated directly using lattice QCD. In this calculation, the largest uncertainty comes from two sources: one is $\hat{B}_K$ and the other is $V_{cb}$. We use the lattice results of $\hat{B}_K$ and exclusive $V_{cb}$ to calculate the theoretical estimate of $ε_K$, which turns out to be $3.1σ$ away from its experimental value. Here, the error is evaluated using the standard error propagation method.

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

Kaon mixing matrix elements from beyond-the-Standard-Model operators in staggered chiral perturbation theory

Models of new physics induce K-Kbar mixing operators having Dirac structures other than the "left-left" form of the Standard Model. We calculate the functional form of the corresponding B-parameters at next-to-leading order in both SU(3) and SU(2) staggered chiral perturbation theory (SChPT). Numerical results for these matrix elements are being generated using improved staggered fermions; our results can be used to extrapolate these matrix elements to the physical light and strange quark masses. The SU(3) SChPT results turn out to be much simpler than that for the Standard Model B_K operator, due to the absence of chiral suppression in the new operators. The SU(2) SChPT result is of similar simplicity to that for B_K. In fact, in the latter case, the chiral logarithms for two of the new B-parameters are identical to those for B_K, while those for the other two new B-parameters are of opposite sign. In addition to providing results for the 2+1 flavor theory in SU(3) SChPT and the 1+1+1 flavor theory in SU(2) SChPT, we present the corresponding continuum partially quenched results, as these are not available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modifying the molecular dynamics action to increase topological tunnelling rate for dynamical overlap fermions

We describe a new Hybrid Monte Carlo (HMC) algorithm for dynamical overlap fermions, which improves the rate of topological index changes by adding an additional (intensive) term to the action for the molecular dynamics part of the algorithm. The metropolis step still uses the exact action, so that the Monte Carlo algorithm still generates the correct ensemble. By tuning this new term, we hope to be able to balance the acceptance rate of the HMC algorithm and the rate of topological index changes. We also describe how suppressing, but not eliminating, the small eigenvalues of the kernel operator may improve the volume scaling of the cost per trajectory for overlap HMC while still allowing topological index changes. We test this operator on small lattices, comparing our new algorithm with an old overlap HMC algorithm with a slower rate of topological charge changes, and an overlap HMC algorithm which fixes the topology. Our new HMC algorithm more than doubles the rate of topological index changes compared to the previous state of the art, while maintaining the same metropolis acceptance rate. We investigate the effect of topological index changes on the local topological charge density, measured using an improved field theoretic operator after heavy smearing. We find that the creation and annihilation of large lumps of topological charge is increased with the new algorithm.

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

Performance of SSE and AVX Instruction Sets

SSE (streaming SIMD extensions) and AVX (advanced vector extensions) are SIMD (single instruction multiple data streams) instruction sets supported by recent CPUs manufactured in Intel and AMD. This SIMD programming allows parallel processing by multiple cores in a single CPU. Basic arithmetic and data transfer operations such as sum, multiplication and square root can be processed simultaneously. Although popular compilers such as GNU compilers and Intel compilers provide automatic SIMD optimization options, one can obtain better performance by a manual SIMD programming with proper optimization: data packing, data reuse and asynchronous data transfer. In particular, linear algebraic operations of vectors and matrices can be easily optimized by the SIMD programming. Typical calculations in lattice gauge theory are composed of linear algebraic operations of gauge link matrices and fermion vectors, and so can adopt the manual SIMD programming to improve the performance.

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

String tension from gauge invariant Magnetic Monopoles

We investigate the relationship between colour confinement and the monopoles derived from the Cho-Duan-Ge decomposition. These monopoles, unlike Dirac and 't Hooft monopoles, do not require a singular gauge field and are defined for any choice of gauge (and are not just restricted to, for example, the maximum Abelian gauge). The Abelian decomposition is defined in terms of a colour field $n$; the principle novelty of our study is that we have used a unique definition of this field in terms of the eigenvectors of the Wilson Loop. This allows us to investigate the relationship between the gauge invariant monopoles and confinement both analytically and numerically, as well as retaining the maximal possible symmetry within the colour field so that it is able to see all the monopoles in an SU($N_C$) calculation. We describe how the Abelian decomposition is related to the Wilson Loop, so that the string tension may be calculated from the field strength related to the decomposed (or restricted) Abelian field. We discuss the structures in the colour field which may cause an area law in the Wilson Loop, which turn out to be magnetic monopoles. If these monopoles are present, they will lead to an area law scaling of the Wilson Loop and thus be at least partially responsible for confinement. We search for these monopoles in quenched lattice QCD. We show that the string tension is dominated by peaks in the restricted field strength, at least some of which are located close to structures in the colour field consistent with with theoretical expectations for the monopoles. We show that the string tension extracted from the monopole contribution to the restricted field is close to that of the entire original field; again suggesting that confinement can at least partially be explained in terms of these monopoles.

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

A further extension to the group of Ginsparg-Wilson (overlap) chiral symmetries

As shown by Mandula, the Ginsparg-Wilson lattice realisation of chiral symmetry has a possible ambiguity: there is no unique lattice chiral symmetry, but an infinite group of symmetries with non-commuting generators. The physical implications of this abundance of symmetry remain unclear. In recent work, it has been shown how these chiral symmetries for overlap fermions can be derived from a renormalisation group blocking in the continuum, transforming the action from the standard continuum action to an equivalent to the lattice overlap action. There is no unique blocking, and different blockings lead to different chiral symmetries. The group of symmetries found by Mandula immediately follows. In this way, the excess chiral symmetry on the lattice can be explained in terms of different renormalisation schemes. The previous work suffered from one technical challenge: there is no continuum analogue of the lattice chiral eigenvectors at eigenvalue $2/a$. As the construction of the overlap operator required a mapping between lattice and continuum eigenvalues, the lack of a continuum counterpart to the doublers of the zero modes creates an ambiguity in the construction. Although the lattice chiral symmetry can still be defined, this leads to difficulties when considering $\mathcal{CP}$ symmetry on the lattice. In this work, we investigate the possibility of resolving this ambiguity by adding a second fermion field to the original continuum action used as a basis of the renormalisation group blockings. This second fermion field has a mass of the order of the momentum cut-off, to simulate the effects of the fermion doublers. Working through the same renormalisation group procedure to map this action to the lattice overlap action yields additional Ginsparg-Wilson relations satisfied by the overlap operator, and more (non-commuting) lattice chiral symmetries.

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

Gell-Mann-Oakes-Renner relation for multiple chiral symmetries

As a first step towards considering a chiral perturbation theory for overlap fermions, we investigate whether there are any ambiguities in the expression for the pion mass resulting from multiple chiral symmetries. The concern is that, calculating the conserved current for Ginsparg Wilson chiral symmetries in the usual way, different expressions of the chiral symmetries lead to different currents. This implies an ambiguity in the definition of the pion and pion decay constant for all Ginsparg-Wilson expressions of the Dirac operator, including the overlap operator. We use a renormalisation group mapping procedure to consider local chiral symmetry transformations for a continuum Ginsparg-Wilson "Dirac-operator." We find that this naturally leads to an expression for the conserved current that differs from the standard expression by cut-off artefacts, but is independent of which of the Ginsparg-Wilson symmetries is chosen. We recover the standard expressions for the massive Dirac operator, propagator, and chiral condensate. With this in place, we proceed to calculate the pion mass in the mapped theory as a function of the quark mass, and discover a unique expression for $F_π$ and $m_π$, recovering the usual Gell-Mann-Oakes-Renner relation, baring the substitution of the chiral condensate with its modified value. We hypothesise that the argument can be carried directly over to the lattice theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi GPU Performance of Conjugate Gradient Solver with Staggered Fermions in Mixed Precision

GPU has a significantly higher performance in single-precision computing than that of double precision. Hence, it is important to take a maximal advantage of the single precision in the CG inverter, using the mixed precision method. We have implemented mixed precision algorithm to our multi GPU conjugate gradient solver. The single precision calculation use half of the memory that is used by the double precision calculation, which allows twice faster data transfer in memory I/O. In addition, the speed of floating point calculations is 8 times faster in single precision than in double precision. The overall performance of our CUDA code for CG is 145 giga flops per GPU (GTX480), which does not include the infiniband network communication. If we include the infiniband communication, the overall performance is 36 giga flops per GPU (GTX480).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Goldstone pion masses with NLO in Staggered Chiral Perturbation Theory

We present results of the masses of taste non-Goldstone $(F \ne ξ_5)$ pions and kaons calculated up to the next-to-leading order in the SU(3) staggered chiral perturbation theory (SChPT). The results can be used to fit data and to understand taste symmetry breaking effect quantitatively. The final expressions for the non-Goldstone masses contain 20 low energy constants unique to the non-Goldstone sector. We have calculated the several cases such as the full QCD, partially quenched QCD, and quenched QCD in the $N_f=1+1+1$ flavor and $N_f=2+1$ flavor cases in the SU(3) and SU(2) SChPT. In this paper, we present only the SU(3) part.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

One-loop matching of improved four-fermion staggered operators with an improved gluon action

We present results for one-loop matching factors of four-fermion operators composed of HYP-smeared staggered fermions. We generalize previous calculations by using the tree-level improved Symanzik gauge action. These results are needed for our companion numerical calculation of $B_K$ and related matrix elements. We find that the impact on one-loop matching factors of using the improved gluon action is much smaller than that from the use of either HYP smearing or mean-field improvement. The one-loop coefficients for mean-field improved, HYP-smeared operators with the Symanzik gauge action have a maximum magnitude of $O(1)\times α_s$, indicating that perturbation theory is reasonably convergent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Taste non-Goldstone, flavor-charged pseudo-Goldstone boson masses in staggered chiral perturbation theory

We calculate the masses of taste non-Goldstone pions and kaons in staggered chiral perturbation theory through next-to-leading order in the standard power counting. The results can be used to quantitatively understand taste violations in existing lattice data generated with staggered fermions and to extract the $u$, $d$, and $s$ quark masses and Gasser-Leutwyler parameters from the experimentally observed spectrum. The expressions for the non-Goldstone masses contain low-energy couplings unique to the non-Goldstone sector. With two exceptions these enter as coefficients of analytic terms; all the new couplings can be fixed by performing spectrum calculations. We report one-loop results for the quenched case and the fully dynamical and partially quenched 1+1+1 and 2+1 flavor cases in the chiral SU(3) and SU(2) theories.

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

Multi GPU Performance of Conjugate Gradient Algorithm with Staggered Fermions

We report results of the performance test of GPUs obtained using the conjugate gradient (CG) algorithm for staggered fermions on the MILC fine lattice ($28^3 \times 96$). We use GPUs of nVIDIA GTX 295 model for the test. When we turn off the MPI communication and use only a single GPU, the performance is 35 giga flops in double precision, which corresponds to 47% of the peak. When we turn on the MPI communication and use multi-GPUs, the performance is reduced down to 12.3 giga flops. The data transfer through the infiniband network and PCI-E bus I/O is a main bottle neck. We suggest two potential solutions of how to optimize the data transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

One loop matching factors for improved staggered four-fermion operators with improved glue

We present results for matching factors for staggered four-fermion operators constructed using HYP-smeared fat links both in the action and the operators. We use perturbation theory to calculate the matching factors and work to one-loop order. The new feaure of this work is the use of the Symanzik-improved gauge action, as opposed to the Wilson gauge action. Our results are needed for our ongoing calculation of weak matrix elements using HYP-smeared staggered valence quarks and operators on MILC lattices. We give explicit results for matching factors of the operator needed to calculate $B_K$. We compare the impact of the improvement of the gauge action on one-loop coefficients with that of mean-field improvement of the operators.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

One-loop matching factors for staggered bilinear operators with improved gauge actions

We present results for one-loop perturbative matching factors using bilinear operators composed of improved staggered fermions, using unimproved (Wilson) and improved (Symanzik, Iwasaki, and DBW2) gluon actions. We consider two fermions actions---HYP/$\bar{\text{Fat7}}$-smeared and "asqtad". The former is being used in calculations of electroweak matrix elements, while the latter have been used extensively by the MILC collaboration. We observe that using the improved gluon action leads to small reductions in the perturbative corrections, but that these reductions are smaller than those obtained when moving from the tadpole-improved naive staggered action to either HYP-smeared or asqtad action.

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
preprint2009arXiv

One loop matching factors for staggered bilinear operators with improved glue

We present results for matching factors for bilinear operators composed of HYP-smeared staggered fermions and constructed using HYP-smeared fat links. The matching factors are calculated perturbatively at one-loop order. The new feature of our calculation compared to previous work on HYP-smeared staggered fermions is the use of the Symanzik-improved gluon propagator, which allows our results to be applied to our ongoing simulations based on configurations generated by the MILC collaboration. We address the issue of the relative efficiency of various improvement schemes in reducing one-loop corrections to the matching factors.

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

Testing improved staggered fermions with $m_s$ and $B_K$

We study the improvement of staggered fermions using hypercubically smeared (HYP) links. We calculate the strange quark mass and the kaon B-parameter, $B_K$, in quenched QCD on a $16^3 \times 64$ lattice at $β=6.0$. We find $m_s(\bar{\rm MS},2 {\rm GeV})=101.2\pm1.3\pm4 $MeV and $B_K(\bar{\rm MS},2 {\rm GeV}) = 0.578 \pm 0.018\pm 0.042$, where the first error is from statistics and fitting, and the second from using one-loop matching factors. The scale ($1/a=1.95$GeV) is set by $M_ρ$, and $m_s$ is determined using the kaon mass. Comparing to quenched results obtained using unimproved staggered fermions and other discretizations, we argue that the size of discretization errors in $B_K$ is substantially reduced by improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Perturbative Matching of the staggered four-fermion operators for e'/e

Using staggered fermions, we calculate the perturbative corrections to the bilinear and four-fermion operators that are used in the numerical study of weak matrix elements for $ε'/ε$. We present results for one-loop matching coefficients between continuum operators, calculated in the Naive Dimensional Regularization (NDR) scheme, and gauge invariant staggered fermion operators. These results, combined with existing results for penguin diagrams, provide the complete one-loop renormalization of the staggered four-fermion operators.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Scaling behavior of improvement and renormalization constants

This talk summarizes results for all the scale independent renormalization constants for bilinear currents ($Z_A$, $Z_V$, and $Z_S/Z_P$), the improvement constants ($c_A$, $c_V$, and $c_T$), the quark mass dependence of $Z_O$, and the coefficients of the equation of motion operators for $O(a)$ improved lattice QCD. Using data at $β=6.0$, 6.2 and 6.4 we study the scaling behavior of these quantities and quantify residual discretization errors.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Non-perturbative Renormalization Constants using Ward Identities

We extend the application of axial Ward identities to calculate $b_A, b_P$ and $b_T$, coefficients that give the mass dependence of the renormalization constants of the corresponding bilinear operators in the quenched theory. The extension relies on using operators with non-degenerate quark masses. It allows a complete determination of the O(a) improvement coefficients for bilinears in the quenched approximation using Ward Identities alone. Only the scale dependent normalization constants $Z_P^0$ (or $Z_S^0$) and $Z_T$ are undetermined. We present results of a pilot numerical study using hadronic correlators.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

K^0--\bar{K}^0 mixing in full lattice QCD

There are at least two methods to calculate $ B_K $ with staggered fermions: one is the two spin trace formalism and the other is the one spin trace formalism. We have performed numerical simulations on a $ 16^3 \times 40 $ lattice in full QCD with $ β= 5.7 $ and a dynamical quark mass 0.01 in lattice units. We try various sources to select only the pseudo-Goldstone bosons and compare the various results.

0 citations0 reviewsNo code linkedOpen access