Source author record

Huey-Wen Lin

Huey-Wen Lin 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

67works
8topics
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

67 published item(s)

preprint2026arXiv

First Nucleon Gluon PDF from Large Momentum Effective Theory

We report the first nucleon gluon parton distribution function (PDF) using Large-Momentum Effective Theory (LaMET). We focus on the gluon operator which was demonstrated to have the best signal-to-noise in the previous attempt [1] in computing gluon PDFs using LaMET. We compute the corresponding Wilson coefficients needed for the hybrid-renormalized matrix elements and the matching kernel to convert the quasi-PDF to the lightcone one at the one-loop level. We demonstrate that with the proper Wilson coefficients in place, the counterterms for the renormalization are independent of the hadron and mass within statistical error. Using the resulting renormalization, we then compute the nucleon PDF using a HISQ ensemble generated by the MILC collaboration with $N_f=2+1+1$, $a \approx 0.12$ fm, with valence pion masses of 310 and 690 MeV and two gauge link smearing techniques. Despite the physics effects of the heavier than physical pion masses and gauge link smearing, this calculation provides excellent proof of principle and compares reasonably with selected global fit results.

preprint2025arXiv

Three-dimensional imaging of hadrons with hard exclusive reactions: advances in experiment, theory, phenomenology, and lattice QCD

Generalized Parton Distributions (GPDs) have emerged as a powerful framework for exploring the internal structure of hadrons in terms of their partonic constituents. Over the past three decades, the field has witnessed significant theoretical and experimental advancements. The interpretation of GPDs in impact parameter space offers a vivid three-dimensional visualization of hadron structure, correlating longitudinal momentum and transverse spatial distributions, thereby enabling tomographic imaging of hadrons. Furthermore, the link between GPDs and the matrix elements of the QCD energy-momentum tensor provides access to fundamental properties of hadrons, including spin decomposition and internal pressure distributions. Notably, recent analyses of Deeply Virtual Compton Scattering (DVCS) data have enabled the empirical extraction of the quark pressure profile inside the proton. This white paper presents an overview of recent developments in GPD theory and phenomenology, as well as progress in lattice QCD studies. It outlines the prospects for advancing our understanding of hadron structure through the next generation of dedicated experiments, including the extension of the Jefferson Lab 12~GeV program (and its potential 22~GeV upgrade), J-PARC, COMPASS/AMBER, LHC ultra-peripheral collisions, and the future electron-ion colliders EIC and EicC.

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

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

Benchmarks for a Global Extraction of Information from Deeply Virtual Exclusive Scattering

We develop a framework to establish benchmarks for machine learning and deep neural networks analyses of exclusive scattering cross sections (FemtoNet). Within this framework we present an extraction of Compton form factors for deeply virtual Compton scattering from an unpolarized proton target. Critical to this effort is a study of the effects of physics constraint built into machine learning (ML) algorithms. We use the Bethe-Heitler process, which is the QED radiative background to deeply virtual Compton scattering, to test our ML models and, in particular, their ability to generalize information extracted from data. We then use these techniques on the full cross section and compare the results to analytic model calculations. We propose a quantification technique, the random targets method, to begin understanding the separation of aleatoric and epistemic uncertainties as they are manifest in exclusive scattering analyses. We propose a set of both physics driven and machine learning based benchmarks providing a stepping stone towards applying explainable machine learning techniques with controllable uncertainties in a wide range of deeply virtual exclusive processes.

preprint2022arXiv

Bubble Chamber Detectors with Light Nuclear Targets: A Snowmass 2021 White Paper

Neutrino cross sections are a critical ingredient in experiments that depend on neutrino scattering to reconstruct event kinematics and infer neutrino characteristics, like NOvA and T2K. An opportunity exists to reduce the 5-10% broad uncertainty on neutrino cross sections by producing more measurements of neutrino scattering from light nuclear targets at the relevant energies. Bubble chambers with light nuclear targets would be ideal for these measurements but the most recent device designed for use with an accelerator neutrino source is at least fifty years old. A new bubble chamber with light nuclear targets could be designed by observing how the technology has progressed for use in dark matter experiments and producing smaller modular devices that use more efficient cooling systems. A smaller modular device could also be designed for deployment to all functioning neutrino beams, but an investigation of the proper operating characteristics is necessary to adapt newer detectors to the structure of contemporary neutrino beams.

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

Impact of lattice $s(x)-\bar{s}(x)$ data in the CTEQ-TEA global analysis

We study the impact of the lattice data on the determination of the strangeness asymmetry distribution $s_-(x) \equiv s(x) - {\bar s}(x)$ in the general CTEQ-TEA global analysis of parton distribution functions (PDFs) of the proton. Firstly, we find that allowing a non-vanishing $s_-(x)$, at the initial $Q_0=1.3$ GeV scale, in a global PDF analysis leads to the quality of the CT18As fit similar to CT18A. Secondly, including the lattice data in the CT18As\_Lat fit greatly reduces the $s_-$-PDF error band size in the large $x$ region. To further reduce its error would require more precise lattice data, extended to smaller $x$ values.

preprint2022arXiv

Lattice QCD Calculations of Parton Physics

In this document, we summarize the status and challenges of calculating parton physics in lattice QCD for the US Particle Physics Community Planning Exercise (a.k.a. "Snowmass"). While PDF-moments calculations have been very successful and been continuously improved, new methods have been developed to calculate distributions directly in $x$-space. Many recent lattice studies have been focused on calculating isovector PDFs of the pion and nucleon, learning to control systematics associated with excited-state contamination, renormalization and continuum extrapolations, pion-mass and finite-volume effects, etc. Although in some cases, the lattice results are already competitive with experimental data, to reach the level of precision in a wide range of $x$ for unpolarized nucleon PDFs impactful for future collider physics remains a challenge, and may require exascale supercomputing power. The new theoretical methods open the door for calculating other partonic observables which will be the focus of the experimental program in nuclear physics, including generalized parton distributions and transverse-momentum dependent PDFs. A fruitful interplay between experimental data and lattice-QCD calculations will usher in a new era for parton physics and hadron structure.

preprint2022arXiv

Neutrino Scattering Measurements on Hydrogen and Deuterium: A Snowmass White Paper

Neutrino interaction uncertainties are a limiting factor in current and next-generation experiments probing the fundamental physics of neutrinos, a unique window on physics beyond the Standard Model. Neutrino-nucleon scattering amplitudes are an important part of the neutrino interaction program. However, since all modern neutrino detectors are composed primarily of heavy nuclei, knowledge of elementary neutrino-nucleon amplitudes relies heavily on experiments performed in the 1970s and 1980s, whose statistical and systematic precision are insufficient for current needs. In this white paper, we outline the motivation for attempting measurements on hydrogen and deuterium that would improve this knowledge, and we discuss options for making these measurements either with the DUNE near detector or with a dedicated facility.

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

Machine-Learning Prediction for Quasi-PDF Matrix Elements

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

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.

preprint2020arXiv

Pion generalized parton distribution from lattice QCD

We present the first lattice calculation of the valence-quark generalized parton distribution (GPD) of the pion using the large-momentum effective theory (LaMET) approach. We focus on the zero-skewness limit, where the GPD has a probability-density interpretation in the longitudinal Bjorken $x$ and the transverse impact-parameter distributions. Our calculation is done using clover valence fermions on an ensemble of gauge configurations with $2+1+1$ flavors (degenerate up/down, strange and charm) of highly improved staggered quarks (HISQ) with lattice spacing $a \approx 0.12$ fm, box size $L \approx 3$ fm and pion mass $m_π\approx 310$ MeV. The parton distribution function and the form factor are reproduced as special limits of the GPD as expected. Due to the large errors, this exploratory study does not show a clear preference among different model assumptions about the kinematic dependence of the GPD. To discriminate between these assumptions, future studies using higher-statistics data will be crucial.

preprint2019arXiv

Hadrons and Nuclei

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD calculations related to the structure and spectroscopy of hadrons and nuclei. An overview of recent lattice calculations of the structure of the proton and other hadrons is presented along with prospects for future extensions. Progress and prospects of hadronic spectroscopy and the study of resonances in the light, strange and heavy quark sectors is summarized. Finally, recent advances in the study of light nuclei from lattice QCD are addressed, and the scope of future investigations that are currently envisioned is outlined.

preprint2019arXiv

Lattice QCD and Neutrino-Nucleus Scattering

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in neutrino-oscillation physics, which inevitably entails nucleon and nuclear structure. In addition to discussing pertinent lattice-QCD calculations of nucleon and nuclear matrix elements, the interplay with models of nuclei is discussed. This program of lattice- QCD calculations is relevant to current and upcoming neutrino experiments, becoming increasingly important on the timescale of LBNF/DUNE and HyperK.

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.

preprint2019arXiv

Pion and Kaon Structure at the Electron-Ion Collider

Understanding the origin and dynamics of hadron structure and in turn that of atomic nuclei is a central goal of nuclear physics. This challenge entails the questions of how does the roughly 1 GeV mass-scale that characterizes atomic nuclei appear; why does it have the observed value; and, enigmatically, why are the composite Nambu-Goldstone (NG) bosons in quantum chromodynamics (QCD) abnormally light in comparison? In this perspective, we provide an analysis of the mass budget of the pion and proton in QCD; discuss the special role of the kaon, which lies near the boundary between dominance of strong and Higgs mass-generation mechanisms; and explain the need for a coherent effort in QCD phenomenology and continuum calculations, in exa-scale computing as provided by lattice QCD, and in experiments to make progress in understanding the origins of hadron masses and the distribution of that mass within them. We compare the unique capabilities foreseen at the electron-ion collider (EIC) with those at the hadron-electron ring accelerator (HERA), the only previous electron-proton collider; and describe five key experimental measurements, enabled by the EIC and aimed at delivering fundamental insights that will generate concrete answers to the questions of how mass and structure arise in the pion and kaon, the Standard Model's NG modes, whose surprisingly low mass is critical to the evolution of our Universe.

preprint2018arXiv

Lattice-QCD Determination of the Hyperon Axial Couplings in the Continuum Limit

We present the first continuum extrapolation of the hyperon octet axial couplings ($g_{ΣΣ}$ and $g_{ΞΞ}$) from $N_f=2+1+1$ lattice QCD. These couplings are important parameters in the low-energy effective field theory description of the octet baryons and fundamental to the nonleptonic decays of hyperons and to hyperon-hyperon and hyperon-nucleon scattering with application to neutron stars. We use clover lattice fermion action for the valence quarks with sea quarks coming from configurations of $N_f=2+1+1$ highly improved staggered quarks (HISQ) generated by MILC Collaboration. Our work includes the first calculation of $g_{ΣΣ}$ and $g_{ΞΞ}$ directly at the physical pion mass on the lattice, and a full account of systematic uncertainty, including excited-state contamination, finite-volume effects and continuum extrapolation, all addressed for the first time. We find the continuum-limit hyperon coupling constants to be $g_{ΣΣ}=0.4455(55)_\text{stat}(65)_\text{sys}$ and $g_{ΞΞ} =-0.2703(47)_\text{stat}(13)_\text{sys}$, which correspond to low-energy constants of $D = 0.708(10)_\text{stat}(6)_\text{sys}$ and $F = 0.438(7)_\text{stat}(6)_\text{sys}$. The corresponding SU(3) symmetry breaking is 9\% which is about a factor of 2 smaller than the earlier lattice estimate.

preprint2016arXiv

From C to Parton Sea: Bjorken-x Dependence of the PDFs

Studying the structure of nucleons is not only important to understanding the strong interactions of quarks and gluons, but also to improving the precision of new-physics searches. Since a broad class of experiments, including the LHC and dark-matter detection, require Standard-Model backgrounds with parton distribution functions (PDFs) as inputs for disentangling SM contributions from potential new physics. For a long time, lattice calculations of the PDFs (as well as many hadron structures) has been limited to the first few moments. In this talk, we present a first direct calculation of the Bjorken-x dependence of the PDFs using Large-Momentum Effective Theory (LaMET). An exploratory study of the antiquark/sea flavor asymmetry of these distributions will be discussed. This breakthrough opens an exciting new frontier calculating more complicated quantities, such as gluon structure and transverse-momentum dependence, which will complement existing theoretical programs for the upcoming Electron-Ion Collider (EIC) or Large Hadron-Electron Collider (LHeC).

preprint2016arXiv

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

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

preprint2016arXiv

Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD

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

preprint2016arXiv

Nucleon Charges, Form-factors and Neutron EDM

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

preprint2016arXiv

Nucleon Helicity and Transversity Parton Distributions from Lattice QCD

We present the first lattice-QCD calculation of the isovector polarized parton distribution functions (both helicity and transversity) using the large-momentum effective field theory (LaMET) approach for direct Bjorken-$x$ dependence. We first review the detailed steps of the procedure in the unpolarized case, then generalize to the helicity and transversity cases. We also derive a new mass-correction formulation for all three cases. We then compare the effects of each finite-momentum correction using lattice data calculated at $M_π\approx 310$ MeV. Finally, we discuss the implications of these results for the poorly known antiquark structure and predict the sea-flavor asymmetry in the transversely polarized nucleon.

preprint2016arXiv

Some Statistics on Women in Lattice QCD

We present a sampling of analyses concerning the gender ratio of plenary speakers during the years 2000--2016 and make comparisons with other conferences, such as the APS April meeting. We hope this will invite discussion of ideas for how to make our field more accessible to women and minorities. We are preparing for an in-depth survey of the lattice field and welcome any ideas or suggestions. To leave post-conference comments and read about news affecting women in our field, see our Facebook page: https://www.facebook.com/WLQCD

preprint2015arXiv

Flavor Structure of the Nucleon Sea from Lattice QCD

We present the first direct lattice calculation of the isovector sea-quark parton distributions using the formalism developed recently by one of the authors. We use $N_f=2+1+1$ HISQ lattice gauge ensembles (generated by MILC Collaboration) and clover valence fermions with pion mass 310 MeV. We are able to obtain the qualitative features of the nucleon sea flavor structure even at this large pion mass: We observe violation of the Gottfried sum rule, indicating $\overline{d}(x) > \overline{u}(x)$; the helicity distribution obeys $Δ\overline{u}(x) > Δ\overline{d}(x)$, which is consistent with the STAR data at large and small leptonic pseudorapidity.

preprint2015arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Probing TeV scale physics in precision UCN decays

We present the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions to the Standard Model at the TeV scale. Results are presented at one value of the lattice spacing, $a=0.12$ fm, and two values of light quarks corresponding to $M_π=310$ and $220$ MeV. We discuss two sources of systematic errors, contribution of excited states to these matrix elements and the renormalization constants, and the efficacy of methods used to control them.

preprint2013arXiv

Neutral Meson Decays into Two Photons from Lattice QCD

A precision determination of the neutral-pion width would improve determinations of the splitting between the up- and down-quark masses, and matrix elements for the decay of neutral mesons into two photons could play a role in the attempt to probe beyond-the-Standard Model physics in muon $g-2$ experiments. The theoretical error is dominated by hadronic light-by-light diagrams, and since direct measurements are extremely difficult, model calculations factorize it into two-photon diagrams connected by the lightest hadrons. We employ perturbative techniques to express the photon as a superposition of QCD eigenstates accessible in lattice-QCD calculations and found that vector-meson dominance is a poor description of the two-photon decay process when both photons are off shell.

preprint2012arXiv

Charmed-Baryon Spectroscopy from Lattice QCD with N_f=2+1+1 Flavors

We present the results of a calculation of the positive-parity ground-state charmed-baryon spectrum using 2+1+1 flavors of dynamical quarks. The calculation uses a relativistic heavy-quark action for the valence charm quark, clover-Wilson fermions for the valence light and strange quarks, and HISQ sea quarks. The spectrum is calculated with a lightest pion mass around 220 MeV, and three lattice spacings (a \approx 0.12 fm, 0.09 fm, and 0.06 fm) are used to extrapolate to the continuum. The light-quark mass extrapolation is performed using heavy-hadron chiral perturbation theory up to O(m_pi^3) and at next-to-leading order in the heavy-quark mass. For the well-measured charmed baryons, our results show consistency with the experimental values. For the controversial J=1/2 Xi_{cc}, we obtain the isospin-averaged value M_{Xi_{cc}}=3595(39)(20)(6) MeV (the three uncertainties are statistics, fitting-window systematic, and systematics from other lattice artifacts, such as lattice scale setting and pion-mass determination), which shows a 1.7 sigma deviation from the experimental value. We predict the yet-to-be-discovered doubly and triply charmed baryons Xi_{cc}^*, Omega_{cc}, Omega_{cc}^* and Omega_{ccc} to have masses 3648(42)(18)(7) MeV, 3679(40)(17)(5) MeV, 3765(43)(17)(5) MeV and 4761(52)(21)(6) MeV, respectively.

preprint2012arXiv

Lattice Hadron Structure: Applications within and beyond QCD

Study of the hadronic matrix elements can provide not only tests of the QCD sector of the Standard Model (in comparing with existing experiments) but also reliable low-energy hadronic quantities applicable to a wide range of beyond-the-Standard Model scenarios where experiments or theoretical calculations are limited or difficult. On the QCD side, progress has been made in the notoriously difficult problem of addressing gluonic structure inside the nucleon, reaching higher-$Q^2$ region of the form factors, and providing a complete picture of the proton spin. However, even further study and improvement of systematic uncertainties are needed. There are also proposed calculations of higher-order operators in the neutron electric dipole moment Lagrangian, which would be useful when combined with effective theory to probe BSM. Lattice isovector tensor and scalar charges can be combined with upcoming neutron beta-decay measurements of the Fierz interference term and neutrino asymmetry parameter to probe new interactions in the effective theory, revealing the scale of potential new TeV particles. Finally, I revisit the systematic uncertainties in recent calculations of $g_A$ and review prospects for future calculations.

preprint2012arXiv

Probing novel TeV physics through precision calculations of scalar and tensor charges of the nucleon

We present an update on the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions of the Standard Model at the TeV scale. Our calculations are being done using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. We provide preliminary estimates of the dependence of these matrix elements on the light quark masses, lattice spacing, and the time separation between the source and sink of the nucleons. We also find that the renormalization constants calculated using the RI-sMOM scheme are close to unity for the HYP smeared HISQ lattices.

preprint2012arXiv

Probing TeV scale physics via ultra cold neutron decays and calculating non-standard baryon matrix elements

We motivate undertaking precision analyses of neutron decays to look for signatures of new scalar and tensor interactions that can arise in extensions of the Standard Model at the TeV scale. The key ingrediant needed to connect experimental data with theoretical analysis are high-precision calculations of matrix elements of isovector bilinear operators between the decaying neutron and final state proton. We describe the status of our Lattice QCD program of using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. On the theoretical side we use the effective field theory method and provide both model independent and dependent analyses to obtain bounds on possible scalar and tensor interactions, both from low energy experiments and LHC data.

preprint2012arXiv

Theoretical Bounds on New Four-Fermion Interactions and TeV Scale Physics

The standard model weak interactions can be described by four-fermion V-A operators at low energies. New physics at the TeV scale can, however, generate the other Lorentz structures. In this talk, we review the constraints on such interactions from nuclear and hadronic decays, as well as from collider searches. Currently the most stringent bounds come from the analysis of the 0+ to 0+ nuclear and the pi to e nu gamma radiative pion decays. In the near future, the ultracold neutron beta decay experiments and the direct LHC measurements will compete in setting the most stringent bounds, provided, however, that the neutron-to-proton non-perturbative transition matrix elements can be calculated to a level of 10-20% accuracy.

preprint2011arXiv

Beyond the Standard Model with Precision Nucleon Matrix Elements on the Lattice

Precision measurements of nucleons provide constraints on the Standard Model and can discern the signatures predicted for particles beyond the Standard Model (BSM). Knowing the Standard Model inputs to nucleon matrix elements will be necessary to constrain the couplings of dark matter candidates such as the neutralino, to relate the neutron electric dipole moment to the CP-violating theta parameter, or to search for new TeV-scale particles though non-$V-A$ interactions in neutron beta decay. However, these matrix elements derive from the properties of quantum chromodynamics (QCD) at low energies, where perturbative treatments fail. Using lattice gauge theory, we can nonperturbatively calculate the QCD path integral on a supercomputer. In this proceeding, I will discuss a few representative areas in which lattice QCD (LQCD) can contribute to the search for BSM physics, emphasizing suppressed operators in neutron decay, and outline prospects for future development.

preprint2011arXiv

Charmed Baryon Spectroscopy from Lattice QCD with $N_f=2+1+1$ flavors

We present preliminary results of the calculation of the positive-parity ground state charm baryon spectrum using of N_f=2+1+1$ dynamical quarks. The calculation uses a relativistic heavy-quark action for the valence charm quark, clover-Wilson fermions for the valence light quarks and HISQ sea quarks. The spectrum is calculated with a lightest pion mass around 220 MeV and two lattice spacings ($a\sim0.12$ fm and 0.09 fm) are used to extrapolate to continuum limit. Our preliminary results are consistent with the currently measured baryon spectrum, except for the isospin-averaged $J=1/2$ $Ξ_{cc}$ which is approximately $2 σ$ above the SELEX observed value. We predict the yet-to-be-discovered double and triple-charm baryons $Ξ_{cc}^*$, $Ω_{cc}$, $Ω_{cc}^*$, $Ω_{ccc}$ to have masses 3665(42)(29) MeV, 3694(40)(45) MeV, 3739(35)(21) MeV and 4782(24)(28) MeV, respectively.

preprint2011arXiv

Comment on `Controversy concerning the definition of quark and gluon angular momentum' by Elliot Leader (arXiv:1101.5956, PRD 83, 096012 (2011))

It is argued by the author that the canonical form of the quark energy-momentum tensor with a partial derivative instead of the covariant derivative is the correct definition for the quark momentum and angular momentum fraction of the nucleon in covariant quantization. Although it is not manifestly gauge invariant, its matrix elements in the nucleon will be non-vanishing and are gauge invariant. We test this idea in the path-integral quantization by calculating correlation functions on the lattice with a gauge-invariant nucleon interpolation field and replacing the gauge link in the quark lattice momentum operator with unity, which corresponds to the partial derivative in the continuum. We find that the ratios of three-point to two-point functions are zero within errors for both the u and d quarks, contrary to the case without setting the gauge links to unity.

preprint2011arXiv

Lattice QCD for Precision Nucleon Matrix Elements

Precision measurements on nucleons provide constraints on the Standard Model and can also discern the signatures predicted for particles beyond the Standard Model. Knowing the Standard Model inputs to nucleon matrix elements will be necessary to constrain the couplings of dark matter candidates such as the neutralino, to relate the neutron electric dipole moment to the CP-violating theta parameter, or to search for new TeV-scale particles though non-$V-A$ interactions in neutron beta decay. However, these matrix elements derive from the properties of quantum chromodynamics (QCD) at low energies, where the coupling is strong and thus perturbative treatments fail. Using lattice gauge theory, we can nonperturbatively calculate the QCD path integral on a supercomputer. In this proceeding, I will review a few representative areas in which lattice QCD can contribute to understanding the structure inside nucleon and how they can contribute to the search for beyond-the-Standard Model physics, with discussions of the difficulties and prospects for future development.

preprint2011arXiv

Medium Effects in Parton Distributions

Understanding the effects of a background hadronic medium on hadronic observables is important in the context of hadron structure. Many experiments probing nucleon structure make use of nuclear targets and unraveling the modifications that ensue is a complex task. Using lattice QCD, we investigate the ab initio computation of hadron structure in a medium, focusing on the structure of the pion in a Bose-condensed gas of pions.

preprint2011arXiv

Nucleon and Pion Form Factors from $N_f=2+1$ Anisotropic Lattices

We report a recent lattice-QCD calculation of nucleon and pion electromagnetic form factors and nucleon axial form factors, with special emphasis on large $Q^2$. Conventional lattice form-factor calculations can only reach about 2.5 GeV$^2$, but in this work the transfer momentum is pushed as large as 6 GeV$^2$. Here, we demonstrate the results on 2+1-flavor anisotropic clover lattices for the nucleon and pion, comparing with low-$Q^2$ quantities, such as Dirac and Pauli radii, anomalous magnetic moments, $g_A$ and $M_A$. Our approach can be applied to isotropic lattices and lattices with smaller lattice spacing to achieve even larger-$Q^2$ form factors. The form factors are processed to obtain transverse charge and magnetization densities across 2-dimensional impact-parameter space. These measurements could give important theoretical input to experiments, such as those of JLab's 12-GeV program, and provide insight into hadronic structure.

preprint2011arXiv

Probing Novel Scalar and Tensor Interactions from (Ultra)Cold Neutrons to the LHC

Scalar and tensor interactions were once competitors to the now well-established V-A structure of the Standard Model weak interactions. We revisit these interactions and survey constraints from low-energy probes (neutron, nuclear, and pion decays) as well as collider searches. Currently, the most stringent limit on scalar and tensor interactions arise from 0+ -> 0+ nuclear decays and the radiative pion decay pi -> e nu gamma, respectively. For the future, we find that upcoming neutron beta decay and LHC measurements will compete in setting the most stringent bounds. For neutron beta decay, we demonstrate the importance of lattice computations of the neutron-to-proton matrix elements to setting limits on these interactions, and provide the first lattice estimate of the scalar charge and a new average of existing results for the tensor charge. Data taken at the LHC is currently probing these interactions at the 10^-2 level (relative to the standard weak interactions), with the potential to reach the 10^-3 level. We show that, with some theoretical assumptions, the discovery of a charged spin-0 resonance decaying to an electron and missing energy implies a lower limit on the strength of scalar interactions probed at low energy.

preprint2011arXiv

Review of Baryon Spectroscopy in Lattice QCD

The complex patterns of the hadronic spectrum have puzzled physicists since the early discovery of the "particle zoo" in the 1960s. Today, the properties of these myriad particles are understood to be the result of quantum chromodynamics (QCD) with some modification by the electroweak interactions. Despite the discovery of this fundamental theory, the description of the hadronic spectrum has long been dominated by phenomenological models, due to the difficulties of addressing QCD in the strong-coupling regime, where nonperturbative effects are essential. By making numerical calculations in discretized spacetime, lattice gauge theory enables the ab initio study of many low-energy properties of QCD. Significant efforts are underway internationally to use lattice QCD to directly compute properties of ground and excited-state baryons. Detailed knowledge of the hadronic spectrum will provide insight into the character of these states beyond what can be extracted from models. In this review, I will focus on the latest progress in lattice calculations of the $P_{11}(1440)$, the poorly known hyperon spectrum and the energies of highly-excited states of the nucleon, Delta and other light-flavor baryons. In the heavy-flavor sector, I will concentrate on recent lattice-QCD calculations of baryon masses, particularly those that make predictions concerning yet-to-be-discovered baryons, such as $Ω_{cc}$, $Ξ^\prime_b$ or triply-heavy baryons.

preprint2010arXiv

Heavy-Baryon Spectroscopy from Lattice QCD

We use a four-dimensional lattice calculation of the full-QCD (quantum chromodynamics, the non-abliean gauge theory of the strong interactions of quarks and gluons) path integrals needed to determine the masses of the charmed and bottom baryons. In the charm sector, our results are in good agreement with experiment within our systematics, except for the spin-1/2 $Ξ_{cc}$, for which we found the isospin-averaged mass to be $Ξ_{cc}$ to be $3665\pm17\pm14^{+0}_{-78}$ MeV. We predict the mass of the (isospin-averaged) spin-1/2 $Ω_{cc}$ to be $3763\pm19\pm26^{+13}_{-79}$ {MeV}. In the bottom sector, our results are also in agreement with experimental observations and other lattice calculations within our statistical and systematic errors. In particular, we find the mass of the $Ω_b$ to be consistent with the recent CDF measurement. We also predict the mass for the as yet unobserved $Ξ^\prime_b$ to be 5955(27) MeV.

preprint2010arXiv

Lattice Calculations of Nucleon Electromagnetic Form Factors at Large Momentum Transfer

In this work, we report a novel technique in lattice QCD for studying the high momentum-transfer region of nucleon form factors. These calculations could give important theoretical input to experiments, such as those of JLab's 12-GeV program and studies of nucleon deformation. There is an extensive history of form-factor calculations on the lattice, primarily with ground states for both the initial and final state. However, determining form factors at large momentum transfer ($Q^2$) has been difficult due to large statistical and systematic errors in this regime. We study the nucleon form factors using three pion masses with both quenched and 2+1-flavor anisotropic lattice configurations with $Q^2$ as large as $6 GeV^2$. These form factors are further processed to obtain transverse charge and magnetization densities across 2-dimensional impact-parameter space. Our approach can be applied to isotropic lattices and lattices with smaller lattice spacing to calculate even larger-$Q^2$ form factors.

preprint2010arXiv

Nucleon isovector structure functions in (2+1)-flavor QCD with domain wall fermions

We report on numerical lattice QCD calculations of some of the low moments of the nucleon structure functions. The calculations are carried out with gauge configurations generated by the RBC and UKQCD collaborations with (2+1)-flavors of dynamical domain wall fermions and the Iwasaki gauge action ($β= 2.13$). The inverse lattice spacing is $a^{-1} = 1.73$ GeV, and two spatial volumes of ((2.7{\rm fm})^3) and ((1.8 {\rm fm})^3) are used. The up and down quark masses are varied so the pion mass lies between 0.33 and 0.67 GeV while the strange mass is about 12 % heavier than the physical one. The structure function moments we present include fully non-perturbatively renormalized iso-vector quark momentum fraction, (< x >_{u-d}), helicity fraction, (< x >_{Δu - Δd}), and transversity, (< 1 >_{δu - δd}), as well as an unrenormalized twist-3 coefficient, (d_1). The ratio of the momentum to helicity fractions, (< x >_{u-d}/< x >_{Δu - Δd}), does not show dependence on the light quark mass and agrees well with the value obtained from experiment. Their respective absolute values, fully renormalized, show interesting trends toward their respective experimental values at the lightest quark mass. A prediction for the transversity, (0.7 < < 1 >_{δu -δd} < 1.1), in the (\bar{\rm MS}) scheme at 2 GeV is obtained. The twist-3 coefficient, (d_1), though yet to be renormalized, supports the perturbative Wandzura-Wilczek relation.

preprint2010arXiv

Nucleon structure from mixed action calculations using 2+1 flavors of asqtad sea and domain wall valence fermions

We present high statistics results for the structure of the nucleon from a mixed-action calculation using 2+1 flavors of asqtad sea and domain wall valence fermions. We perform extrapolations of our data based on different chiral effective field theory schemes and compare our results with available information from phenomenology. We discuss vector and axial form factors of the nucleon, moments of generalized parton distributions, including moments of forward parton distributions, and implications for the decomposition of the nucleon spin.

preprint2010arXiv

Nucleon, Delta and Omega excited state spectra at three pion mass values

The energies of the excited states of the Nucleon, Delta and Omega are computed in lattice QCD, using two light quarks and one strange quark on anisotropic lattices. The calculations are performed at three values of the pion mass: 392(4), 438(3) and 521(3) MeV. We employ the variational method with a basis of about ten interpolating operators enabling six energies to be distinguished clearly in each irreducible representation of the octahedral group. We compare our calculations of nucleon excited states with the low-lying experimental spectrum. There is reasonable agreement for the pattern of states.

preprint2009arXiv

Excited-Nucleon Spectroscopy with 2+1 Fermion Flavors

We present progress made by the Hadron Spectrum Collaboration (HSC) in determining the tower of excited nucleon states using 2+1-flavor anisotropic clover lattices. The HSC has been investigating interpolating operators projected into irreducible representations of the cubic group in order to better calculate two-point correlators for nucleon spectroscopy; results are published for quenched and 2-flavor anisotropic Wilson lattices. In this work, we present the latest results using a new technique, distillation, which allows us to reach higher statistics than before. Future directions will be outlined at the end.

preprint2009arXiv

High Statistics Analysis using Anisotropic Clover Lattices: (III) Baryon-Baryon Interactions

Low-energy baryon-baryon interactions are calculated in a high-statistics lattice QCD study on a single ensemble of anisotropic clover gauge-field configurations at a pion mass of m_π~ 390 MeV, a spatial volume of L^3 ~ (2.5 fm)^3, and a spatial lattice spacing of b~0.123 fm. Lüscher's method is used to extract nucleon-nucleon, hyperon-nucleon and hyperon-hyperon scattering phase shifts at one momentum from the one- and two-baryon ground-state energies in the lattice volume. The isospin-3/2 NΣinteractions are found to be highly spin-dependent, and the interaction in the ^3S_1 channel is found to be strong. In contrast, the NΛinteractions are found to be spin-independent, within the uncertainties of the calculation, consistent with the absence of one-pion-exchange. The only channel for which a negative energy-shift is found is ΛΛ, indicating that the ΛΛinteraction is attractive, as anticipated from model-dependent discussions regarding the H-dibaryon. The NN scattering lengths are found to be small, clearly indicating the absence of any fine-tuning in the NN-sector at this pion mass. This is consistent with our previous Lattice QCD calculation of NN interactions. The behavior of the signal-to-noise ratio in the baryon-baryon correlation functions, and in the ratio of correlation functions that yields the ground-state energy splitting is explored. In particular, focus is placed on the window of time slices for which the signal-to-noise ratio does not degrade exponentially, as this provides the opportunity to extract quantitative information about multi-baryon systems.

preprint2009arXiv

Nucleon form factors with 2+1 flavor dynamical domain-wall fermions

We report our numerical lattice QCD calculations of the isovector nucleon form factors for the vector and axialvector currents: the vector, induced tensor, axialvector, and induced pseudoscalar form factors. The calculation is carried out with the gauge configurations generated with N_f=2+1 dynamical domain wall fermions and Iwasaki gauge actions at beta = 2.13, corresponding to a cutoff 1/a = 1.73 GeV, and a spatial volume of (2.7 fm)^3. The up and down quark masses are varied so the pion mass lies between 0.33 and 0.67 GeV while the strange quark mass is about 12% heavier than the physical one. We calculate the form factors in the range of momentum transfers, 0.2 < q^2 < 0.75 GeV^2. The vector and induced tensor form factors are well described by the conventional dipole forms and result in significant underestimation of the Dirac and Pauli mean-squared radii and the anomalous magnetic moment compared to the respective experimental values. We show that the axialvector form factor is significantly affected by the finite spatial volume of the lattice. In particular in the axial charge, g_A/g_V, the finite volume effect scales with a single dimensionless quantity, m_pi L, the product of the calculated pion mass and the spatial lattice extent. Our results indicate that for this quantity, m_pi L > 6 is required to ensure that finite volume effects are below 1%.

preprint2009arXiv

Singly and Doubly Charmed J=1/2 Baryon Spectrum from Lattice QCD

We compute the masses of the singly and doubly charmed baryons in full QCD using the relativistic Fermilab action for the charm quark. For the light quarks we use domain-wall fermions in the valence sector and improved Kogut-Susskind sea quarks. We use the low-lying charmonium spectrum to tune our heavy-quark action and as a guide to understanding the discretization errors associated with the heavy quark. Our results are in good agreement with experiment within our systematics. For the Ξ_{cc}, we find the isospin-averaged mass to be M_{Ξ_{cc}} = 3665 +/- 17 +/- 14 +0 -78 MeV; the three given uncertainties are statistical, systematic and an estimate of lattice discretization errors, respectively. In addition, we predict the mass splitting of the (isospin-averaged) spin-1/2 Ω_{cc} with the Ξ_{cc} to be M_{Ω_{cc}} - M_{Ξ_{cc}} = 98 +/- 9 +/- 22 +/- 13 MeV (in this mass splitting, the leading discretization errors are also suppressed by SU(3) symmetry). Combining this splitting with our determination of M_{Ξ_{cc}} leads to our prediction of the spin-1/2 Ω_{cc} mass, M_{Ω_{cc}} = 3763 +/- 19 +/- 26 +13 -79 MeV.

preprint2009arXiv

Strange Baryon Electromagnetic Form Factors and SU(3) Flavor Symmetry Breaking

We study the nucleon, Sigma and cascade octet baryon electromagnetic form factors and the effects of SU(3) flavor symmetry breaking from 2+1-flavor lattice calculations. We find that electric and magnetic radii are similar; the maximum discrepancy is about 10%. In the pion-mass region we explore, both the quark-component and full-baryon moments have small SU(3) symmetry breaking. We extrapolate the charge radii and the magnetic moments using three-flavor heavy-baryon chiral perturbation theory (HBXPT). The systematic errors due to chiral and continuum extrapolations remain significant, giving rise to charge radii for $p$ and $Σ^-$ that are 3--4 standard deviations away from the known experimental ones. Within these systematics the predicted $Σ^+$ and $Ξ^-$ radii are 0.67(5) and 0.306(15) fm$^2$ respectively. When the next-to-next-to-leading order of HBXPT is included, the extrapolated magnetic moments are less than 3 standard deviations away from PDG values, and the discrepancy is possibly due to remaining chiral and continuum extrapolation errors.

preprint2008arXiv

Cascade Baryon Spectrum from Lattice QCD

A comprehensive study of the cascade baryon spectrum using lattice QCD affords the prospect of predicting the masses of states not yet discovered experimentally, and determining the spin and parity of those states for which the quantum numbers are not yet known. The study of the cascades, containing two strange quarks, is particularly attractive for lattice QCD in that the chiral effects are reduced compared to states composed only of u/d quarks, and the states are typically narrow. We report preliminary results for the cascade spectrum obtained by using anisotropic N_f = 2 Wilson lattices with temporal lattice spacing inverse 5.56 GeV.

preprint2008arXiv

Excited nucleon spectrum with two flavors of dynamical fermions

We compute the spectrum of excited nucleons using the anisotropic Wilson lattice with two flavors of dynamical fermions. Using optimized sets of operators which transform irreducibly under the octahedral group, matrices of correlation functions are computed. We apply the variational method to these matrices to extract excited energy eigenstates. We obtain several states for each irrep and identify the continuum spin for the lowest-lying states, including a $J^P={5/2}^-$ state.

preprint2008arXiv

Light-Meson Two-Photon Decays in Full QCD

We present a study of two-photon decays of light mesons, focusing on the neutral pion decay. This important process highlights the effects of the axial anomaly in QCD but has been little studied on the lattice. By applying the Lehmann-Symanzik-Zimmermann (LSZ) reduction formula, we reconstruct the electromagnetic matrix elements from three-point vector-vector Green functions calculated on 2+1-flavor isotropic clover lattices.

preprint2007arXiv

Charm spectroscopy on dynamical 2+1 flavor domain wall fermion lattices with a relativistic heavy quark action

We present a preliminary calculation of the charmonium spectrum using the dynamical 2+1 flavor $24^3\times 64$ domain wall fermion lattice configurations generated by the RBC and UKQCD collaborations. We use the relativistic heavy quark action with 3 parameters non-perturbatively determined by matching to experimental quantities. Chiral extrapolation is done on four light sea quark masses from 0.005 to 0.03, with $m_s=0.04$ and $m_{res}=0.003$. We can either predict meson masses assuming the lattice spacing is known from other methods, or calculate the lattice spacing using those quantities.

preprint2007arXiv

Charmed spectroscopy from a nonperturbatively determined relativistic heavy quark action in full QCD

We present a preliminary calculation of the charmed meson spectrum using the 2+1 flavor domain wall fermion lattice configurations currently being generated by the RBC and UKQCD collaborations. The calculation is performed using the 3-parameter, relativistic heavy quark action with nonperturbatively determined coefficients. We will also demonstrate a step-scaling procedure for determining these coefficients nonperturbatively using a series of quenched, gauge field ensembles generated for three different lattice spacings.

preprint2007arXiv

Parameter Tuning of Three-Flavor Dynamical Anisotropic Clover Action

In this work, we perform parameter tuning with dynamical anisotropic clover lattices using the Schrödinger functional and stout-smearing in the fermion field. We find that $ξ_R/ξ_0$ is relatively close to 1 in our parameter search, which allows us to fix $ξ_0$ in our runs. We proposed to determine the gauge and fermion anisotropy in a Schrödinger-background small box using Wilson loop ratios and PCAC masses. We demonstrate that these ideas are equivalent to but more efficient than the conventional meson dispersion approach. The spatial and temporal clover coefficients are fixed to the tree-level tadpole-improved clover values, and we demonstrate that they satisfy the nonperturbative condition determined by Schrödinger functional method.

preprint2006arXiv

Nucleon structure in lattice QCD with dynamical domain-wall fermions quarks

We report RBC and RBC/UKQCD lattice QCD numerical calculations of nucleon electroweak matrix elements with dynamical domain-wall fermions (DWF) quarks. The first, RBC, set of dynamical DWF ensembles employs two degenerate flavors of DWF quarks and the DBW2 gauge action. Three sea quark mass values of 0.04, 0.03 and 0.02 in lattice units are used with about 200 gauge configurations each. The lattice cutoff is about 1.7 GeV and the spatial volume is about (1.9 fm)^3. Despite the small volume, the ratio of the isovector vector and axial charges g_A/g_V and that of structure function moments <x>_{u-d}/<x>_{Delta u - Delta d} are in agreement with experiment, and show only very mild quark mass dependence. The second, RBC/UK, set of ensembles employs one strange and two degenerate (up and down) dynamical DWF quarks and Iwasaki gauge action. The strange quark mass is set at 0.04, and three up/down mass values of 0.03, 0.02 and 0.01 in lattice units are used. The lattice cutoff is about 1.6 GeV and the spatial volume is about (3.0 fm)^3. Even with preliminary statistics of 25-30 gauge configurations, the ratios g_A/g_V and <x>_{u-d}/<x>_{Delta u - Delta d} are consistent with experiment and show only very mild quark mass dependence. Another structure function moment, d_1, though yet to be renormalized, appears small in both sets.