BZPEER

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.

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.

preprint2022arXiv

Discovering new physics in rare kaon decays

The decays and mixing of $K$ mesons are remarkably sensitive to the weak interactions of quarks and leptons at high energies. They provide important tests of the standard model at both first and second order in the Fermi constant $G_F$ and offer a window into possible new phenomena at energies as high as 1,000 TeV. These possibilities become even more compelling as the growing capabilities of lattice QCD make high-precision standard model predictions possible. Here we discuss and attempt to forecast some of these capabilities.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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

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.

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

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\,σ$ .

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.

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.

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.

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.

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

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

Yong-Chull Jang

Research overview

Is this your work?

41 publications

Nucleon form factors and the pion-nucleon sigma term

2021 Update on εK\varepsilon_KεK​ with lattice QCD inputs

Discovering new physics in rare kaon decays

Excited states and precision results for nucleon charges and form factors

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

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

Nucleon charges and form factors using clover and HISQ ensembles

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

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

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

Status report on εK\varepsilon_KεK​ with lattice QCD inputs

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

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

Current status of εK\varepsilon_KεK​ with lattice QCD inputs

Determination of εK\varepsilon_KεK​ using lattice QCD inputs

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

Standard Model evaluation of εK\varepsilon_KεK​ using lattice QCD inputs for B^K\hat{B}_KB^K​ and VcbV_{cb}Vcb​

Calculation of BSM Kaon B-parameters using Staggered Quarks

Code Optimization on Kepler GPUs and Xeon Phi

Current status of εK\varepsilon_KεK​ calculated with lattice QCD inputs

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

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

Improved determination of BKB_KBK​ with staggered quarks

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

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

Optimization of the Oktay-Kronfeld Action Conjugate Gradient Inverter

Update on BKB_KBK​ and εK\varepsilon_KεK​ with staggered quarks

Beyond the Standard Model corrections to K0−Kˉ0K^0-\bar{K}^0K0−Kˉ0 mixing

Covariance fitting of highly correlated data in lattice QCD

Current Status of Indirect CP Violation in Neutral Kaon System

Kaon BBB-parameter from improved staggered fermions in Nf=2+1N_f=2+1Nf​=2+1 QCD

Recent progress in calculation of BKB_KBK​ using staggered fermions

BKB_KBK​ from improved staggered fermions using SU(3) chiral perturbation theory

Continuum extrapolation of BKB_KBK​ with staggered fermions

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

Trouble shooting for covariance fitting in highly correlated data

BKB_KBK​ using HYP-smeared staggered fermions in Nf=2+1N_f=2+1Nf​=2+1 unquenched QCD

BKB_KBK​ with improved staggered fermions: analysis using SU(2) staggered chiral perturbation theory

BKB_KBK​ with improved staggered fermions: analysis using SU(3) staggered chiral perturbation theory

First results for BKB_KBK​ on the ultrafine (a=0.045a=0.045a=0.045 fm) ensemble

Systematic Uncertainties in BKB_KBK​ with Improved Staggered Fermions

See the researcher in context

Building this map preview

2021 Update on $\varepsilon_K$ with lattice QCD inputs

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

Status report on $\varepsilon_K$ with lattice QCD inputs

Current status of $\varepsilon_K$ with lattice QCD inputs

Determination of $\varepsilon_K$ using lattice QCD inputs

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

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

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

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

Improved determination of $B_K$ with staggered quarks

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

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

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

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

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

Recent progress in calculation of $B_K$ using staggered fermions

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

Continuum extrapolation of $B_K$ with staggered fermions

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

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

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

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

Systematic Uncertainties in $B_K$ with Improved Staggered Fermions