Source author record

Sungwoo Park

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

30 published item(s)

preprint2026arXiv

The Spectrum and Scale Setting on 2+1-flavor NME Lattices

This paper describes the thirteen ensembles, named NME, generated with 2+1-flavor Wilson-clover fermions by the JLab/W\&M/LANL/MIT/Marseille collaborations, and presents an analysis of the meson and baryon spectrum, decay constants $f_π$ and $f_K$, flow scales $t_0$ and $w_0$, and time histories of the $Θ$ and Weinberg operators under gradient flow. Using these quantities, the physical point values of the two flow scales, ${t_0^{\rm Phy}}$ and ${w_0^{\rm Phy}}$, and the ratio $\mathop{f_K / f_π}^{\rm Phy}$ are determined. The masses of the octet and decuplet baryons are analyzed using both the next-to-leading order (NLO) and the next-next-to-leading order (NNLO) ansatz from heavy baryon chiral perturbation theory (HB$χ$PT). The NNLO fit to the octet baryons, $M_N$, $M_Σ$, $M_Λ$ and $M_Ξ$, is preferred while the corresponding fits to the decuplet Omega mass, $M_Ω$, are not distinguished. We also present a study of the autocorrelations in the data and show that there is no evidence, even at large flow time, of the freezing of the topological charge or the Weinberg three-gluon operator.

preprint2024arXiv

Update on flavor diagonal nucleon charges from clover fermions

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

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.

preprint2023arXiv

Update on Flavor Diagonal Nucleon Charges

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

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

ETF Portfolio Construction via Neural Network trained on Financial Statement Data

Recently, the application of advanced machine learning methods for asset management has become one of the most intriguing topics. Unfortunately, the application of these methods, such as deep neural networks, is difficult due to the data shortage problem. To address this issue, we propose a novel approach using neural networks to construct a portfolio of exchange traded funds (ETFs) based on the financial statement data of their components. Although a number of ETFs and ETF-managed portfolios have emerged in the past few decades, the ability to apply neural networks to manage ETF portfolios is limited since the number and historical existence of ETFs are relatively smaller and shorter, respectively, than those of individual stocks. Therefore, we use the data of individual stocks to train our neural networks to predict the future performance of individual stocks and use these predictions and the portfolio deposit file (PDF) to construct a portfolio of ETFs. Multiple experiments have been performed, and we have found that our proposed method outperforms the baselines. We believe that our approach can be more beneficial when managing recently listed ETFs, such as thematic ETFs, of which there is relatively limited historical data for training advanced machine learning methods.

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

Flavor Diagonal Nucleon Charges

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

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.

preprint2022arXiv

Reinforcement Learning Portfolio Manager Framework with Monte Carlo Simulation

Asset allocation using reinforcement learning has advantages such as flexibility in goal setting and utilization of various information. However, existing asset allocation methods do not consider the following viewpoints in solving the asset allocation problem. First, State design without considering portfolio management and financial market characteristics. Second, Model Overfitting. Third, Model training design without considering the statistical structure of financial time series data. To solve the problem of the existing asset allocation method using reinforcement learning, we propose a new reinforcement learning asset allocation method. First, the state of the portfolio managed by the model is considered as the state of the reinforcement learning agent. Second, Monte Carlo simulation data are used to increase training data complexity to prevent model overfitting. These data can have different patterns, which can increase the complexity of the data. Third, Monte Carlo simulation data are created considering various statistical structures of financial markets. We define the statistical structure of the financial market as the correlation matrix of the assets constituting the financial market. We show experimentally that our method outperforms the benchmark at several test intervals.

preprint2022arXiv

Shai-am: A Machine Learning Platform for Investment Strategies

The finance industry has adopted machine learning (ML) as a form of quantitative research to support better investment decisions, yet there are several challenges often overlooked in practice. (1) ML code tends to be unstructured and ad hoc, which hinders cooperation with others. (2) Resource requirements and dependencies vary depending on which algorithm is used, so a flexible and scalable system is needed. (3) It is difficult for domain experts in traditional finance to apply their experience and knowledge in ML-based strategies unless they acquire expertise in recent technologies. This paper presents Shai-am, an ML platform integrated with our own Python framework. The platform leverages existing modern open-source technologies, managing containerized pipelines for ML-based strategies with unified interfaces to solve the aforementioned issues. Each strategy implements the interface defined in the core framework. The framework is designed to enhance reusability and readability, facilitating collaborative work in quantitative research. Shai-am aims to be a pure AI asset manager for solving various tasks in financial markets.

preprint2022arXiv

The pion-nucleon sigma term from Lattice QCD

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

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.

preprint2016arXiv

Dynamic Subarrays for Hybrid Precoding in Wideband mmWave MIMO Systems

Hybrid analog/digital precoding architectures can address the trade-off between achievable spectral efficiency and power consumption in large-scale MIMO systems. This makes it a promising candidate for millimeter wave systems, which require deploying large antenna arrays at both the transmitter and receiver to guarantee sufficient received signal power. Most prior work on hybrid precoding focused on narrowband channels and assumed fully-connected hybrid architectures. MmWave systems, though, are expected to be wideband with frequency selectivity. In this paper, a closed-form solution for fully-connected OFDM-based hybrid analog/digital precoding is developed for frequency selective mmWave systems. This solution is then extended to partially-connected but fixed architectures in which each RF chain is connected to a specific subset of the antennas. The derived solutions give insights into how the hybrid subarray structures should be designed. Based on them, a novel technique that dynamically constructs the hybrid subarrays based on the long-term channel characteristics is developed. Simulation results show that the proposed hybrid precoding solutions achieve spectral efficiencies close to that obtained with fully-digital architectures in wideband mmWave channels. Further, the results indicate that the developed dynamic subarray solution outperforms the fixed hybrid subarray structures in various system and channel conditions.

preprint2016arXiv

Optimal User Loading in Massive MIMO Systems with Regularized Zero Forcing Precoding

We consider a downlink multiuser multiple-input multiple output (MIMO) system employing regularized zero-forcing (RZF) precoding. We derive the asymptotic signal-to-leakage-plus-noise ratio (SLNR) as both the number of antennas and the number of users go to infinity at a fixed ratio. Focusing on the symmetric uncorrelated channels, we show that the SLNR is asymptotically equal to signal-to-interference-plus-noise ratio (SINR) which allows us to optimize the user loading for spectral efficiency. The results show that the optimal user loading varies depending on the channel signal-to-noise ratio (SNR) but is equal to one in both the low or high SNR regimes.

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.

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

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.

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.

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

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

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.

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