BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Keh-Fei Liu

Keh-Fei Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phnucl-thnucl-exhep-exhep-th

Catalog footprint

What is connected

51works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2023arXiv

Nucleon Electric Dipole Moment from the $θ$ Term with Lattice Chiral Fermions

We calculate the nucleon electric dipole moment (EDM) from the $θ$ term with overlap fermions on three domain wall lattices with different sea pion masses at lattice spacing 0.11 fm. Due to the chiral symmetry conserved by the overlap fermions, we have well defined topological charge and chiral limit for the EDM. Thus, the chiral extrapolation can be carried out reliably at nonzero lattice spacings. We use three to four different partially quenched valence pion masses for each sea pion mass and find that the EDM dependence on the valence and sea pion masses behaves oppositely, which can be described by partially quenched chiral perturbation theory. With the help of the cluster decomposition error reduction (CDER) technique, we determine the neutron and proton EDM at the physical pion mass to be $d_{n}=-0.00148\left(14\right)\left(31\right)\barθ$ e$\cdot$fm and $d_{p}=0.0038\left(11\right)\left(8\right)\barθ$ e$\cdot$fm. This work is a clear demonstration of the advantages of using chiral fermions in the nucleon EDM calculation and paves the road to future precise studies of the strong $CP$ violation effects.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Proton momentum and angular momentum decompositions with overlap fermions

We present a calculation of the proton momentum and angular momentum decompositions using overlap fermions on a $2+1$-flavor RBC/UKQCD domain-wall lattice at 0.143 fm with a pion mass of 171 MeV which is close to the physical one. A complete determination of the momentum and angular momentum fractions carried by up, down, strange and glue inside the proton has been done with valence pion masses varying from 171 to 391 MeV. We have utilized fast Fourier transform on the stochastic-sandwich method for connected-insertion parts and the cluster-decomposition error reduction technique for disconnected-insertion parts has been used to reduce statistical errors. The full nonperturbative renormalization and mixing between the quark and glue operators are carried out. The final results are normalized with the momentum and angular momentum sum rules and reported at the physical valence pion mass at ${\overline{\rm {MS}}}\, (μ= 2\ {\rm{GeV}})$. The renormalized momentum fractions for the quarks and glue are $\langle x \rangle^q = 0.491(20)(23)$ and $\langle x \rangle^g = 0.509(20)(23)$, respectively, and the renormalized total angular momentum fractions for quarks and glue are $2 J^q = 0.539(22)(44)$ and $2 J^g = 0.461(22)(44)$, respectively. The quark spin fraction is $Σ= 0.405(25)(37)$ from our previous work and the quark orbital angular momentum fraction is deduced from $2 L^q = 2 J^q - Σ$ to be $0.134(22)(44)$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

RI/MOM and RI/SMOM renormalization of quark bilinear operators using overlap fermions

We present the vector, scalar and tensor renormalization constants (RCs) using overlap fermions with either regularization independent momentum subtraction (RI/MOM) or symmetric momentum subtraction (RI/SMOM) as the intermediate scheme on the lattice with lattice spacings $a$ from 0.04 fm to 0.12 fm. Our gauge field configurations from the MILC and RBC/UKQCD collaborations include sea quarks using either the domain wall or the HISQ action, respectively. The results show that RI/MOM and RI/SMOM can provide consistent renormalization constants to the $\overline{\textrm{MS}}$ scheme, after proper $a^2p^2$ extrapolations. But at $p\sim 2$\,GeV, both RI/MOM and RI/SMOM suffer from nonperturbative effects which cannot be removed by the perturbative matching. The comparison between the results with different sea actions also suggests that the renormalization constant is discernibly sensitive to the lattice spacing but not to the bare gauge coupling in the gauge action.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Status on Lattice Calculations of the Proton Spin Decomposition

Lattice calculations of the proton spin components is reviewed. The lattice results of the quark spin from the axial-vector current matrix element at $\sim 0.3- 0.4$ is smaller than those from the constituent quark models. This is largely due to the fact that the vacuum polarization contribution from the disconnected insertion is negative. Its connection with the anomalous Ward identity is clarified and verified numerically. This resolves the contentious issue in the `proton spin crisis'. The glue spin and angular momentum are found to be large and there is notable contribution from the quark orbital angular momentum. Renormalization, mixing and normalization of the quark and glue angular momenta are discussed. With sufficient precision, they can be compared with more precise experimental measurements when the electron-ion collider facility is available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors $G^c_{E,M}(Q^2)$ in the momentum transfer range $0\leq Q^2 \leq 1.4$ $\rm GeV^2$. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment $μ^c_M=-0.00127(38)_{\rm stat}(5)_{\rm sys}$, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero $G^c_{E}(Q^2)$ indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the $[c(x)-\bar{c}(x)]$ distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PDFs and Neutrino-Nucleon Scattering from Hadronic Tensor

We review the Euclidean path-integral formulation of the nucleon hadronic tensor and classify the gauge invariant and topologically distinct insertions in terms of connected and disconnected insertions and also in terms of leading and higher-twist contributions in the DIS region. Converting the Euclidean hadronic tensor back to the Minkowski space requires solving an inverse problem of the Laplace transform. We have investigated several inverse algorithms and studied the pros and cons of each. We show a result with a relatively large momentum transfer ($Q^2 \sim 4\, {\rm GeV^2}$) to suppress the elastic scattering and reveal the contributions from the resonance and inelastic region of the neutrino-nucleon scattering. For elastic scattering, the hadronic tensor is the the product of the elastic form factors for the two corresponding currents. We checked numerically for the case of two charge vector currents ($V_4$) with the electric form factor calculated from the three-point function and found they agree within errors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ratio of strange to $u/d$ momentum fraction in disconnected insertions

The ratio of the strange quark momentum fraction $\langle x\rangle_{s+\bar{s}}$ to that of light quark $u$ or $d$ in disconnected insertions (DI) is calculated on the lattice with overlap fermions on four domain wall fermion ensembles. These ensembles cover three lattice spacings, three volumes and several pion masses including the physical one, from which a global fitting is carried out. A complete nonperturbative renormalization and the mixing between the quark and glue operators are taken into account. We find the ratio to be $\langle x\rangle_{s+\bar{s}}/\langle x\rangle_{u+\bar{u}} ({\rm DI})=0.795(79)(77)$ at $μ= 2$ GeV in the $\overline{\rm MS}$ scheme. This ratio can be used as a constraint to better determine the strange parton distribution especially in the small $x$ region in the global fittings of PDFs when the connected and disconnected sea are fitted and evolved separately, demonstrating a new way that connects lattice calculations with global analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Roper State from Overlap Fermions

The Roper state is extracted with valence overlap fermions on a $2+1$-flavor domain-wall fermion lattice (spacing $a = 0.114$ fm and $m_π = 330$ MeV) using both the Sequential Empirical Bayes (SEB) method and the variational method. The results are consistent, provided that a large smearing-size interpolation operator is included in the variational calculation to have better overlap with the lowest radial excitation. Similar calculations carried out for an anisotropic clover lattice with similar parameters find the Roper $\approx 280$ MeV higher than that of the overlap fermion. The fact that the prediction of the Roper state by overlap fermions is consistently lower than those of clover fermions, chirally improved fermions, and twisted-mass fermions over a wide range of pion masses has been dubbed a "Roper puzzle." To understand the origin of this difference, we study the hairpin $Z$-diagram in the isovector scalar meson ($a_0$) correlator in the quenched approximation. Comparing the $a_0$ correlators for clover and overlap fermions, at a pion mass of 290 MeV, we find that the spectral weight of the ghost state with clover fermions is smaller than that of the overlap at $a = 0.12$ fm and $0.09$ fm, whereas the whole $a_0$ correlators of clover and overlap at $a = 0.06$ fm coincide within errors. This suggests that chiral symmetry is restored for clover at $a \le 0.06$ fm and that the Roper should come down at and below this $a$. We conclude that this work supports a resolution of the "Roper puzzle" due to $Z$-graph type chiral dynamics. This entails coupling to higher components in the Fock space (e.g. $Nπ$, $Nππ$ states) to induce the effective flavor-spin interaction between quarks as prescribed in the chiral quark model, resulting in the parity-reversal pattern as observed in the experimental excited states of $N, Δ$ and $Λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The nucleon isovector tensor charge from lattice QCD using chiral fermions

In this work we present the isovector flavor combination for the nucleon tensor charge extracted from lattice QCD simulations using overlap fermions on $N_f=2+1$ domain-wall configurations. The pion mass dependence is studied using six valence quark masses, each reproducing a value for the pion mass in the valence sector between 147 and 330 MeV. We investigate and eliminate systematic uncertainties due to contamination by excited states, by employing several values for the source-sink separation that span from 1 fm to 1.6 fm. We apply a chiral extrapolation in the valence sector using a quadratic and a logarithmic term to fit the pion mass dependence, which describes well the lattice data. The lattice matrix element is renormalized non-perturbatively, and the final result is $g_T=1.096(30)$ in the $\overline{\rm MS}$ scheme at a renormalization scale of 2 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weak Neutral Current Axial Form Factor Using $(\barν)ν$-Nucleon Scattering and Lattice QCD Inputs

We present a determination of the neutral current weak axial charge $G^Z_A(0)=-0.654(3)_{\rm stat}(5)_{\rm sys}$ using the strange quark axial charge $G^s_A(0)$ calculated with lattice QCD. We then perform a phenomenological analysis, where we combine the strange quark electromagnetic form factor from lattice QCD with (anti)neutrino-nucleon scattering differential cross section from MiniBooNE experiments in a momentum transfer region $0.24\lesssim Q^2 \lesssim 0.71$ GeV$^2$ to determine the neutral current weak axial form factor $G^Z_A(Q^2)$ in the range of $0\lesssim Q^2\leq 1$ GeV$^2$. This yields a phenomenological value of $G^Z_A(0)=-0.687(89)_{\rm stat}(40)_{\rm sys}$. The value of $G^Z_A(0)$ constrained by the lattice QCD calculation of $G^s_A(0)$, when compared to its phenomenological determination, provides a significant improvement in precision and accuracy and can be used to provide a constraint on the fit to $G^Z_A(Q^2)$ for $Q^2>0$. This constrained fit leads to an unambiguous determination of (anti)neutrino-nucleon neutral current elastic scattering differential cross section near $Q^2=0$ and can play an important role in numerically isolating nuclear effects in this region. We show a consistent description of $G^Z_A(Q^2)$ obtained from the (anti)neutrino-nucleon scattering cross section data requires a nonzero contribution of the strange quark electromagnetic form factor. We demonstrate the robustness of our analysis by providing a post-diction of the BNL E734 experimental data.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Towards the nucleon hadronic tensor from lattice QCD

We present the first calculation of the hadronic tensor on the lattice for the nucleon. The hadronic tensor can be used to extract the structure functions in deep inelastic scatterings and also provide information for the neutrino-nucleon scattering which is crucial to the neutrino-nucleus scattering experiments at low energies. The most challenging part in the calculation is to solve an inverse problem. We have implemented and tested three algorithms using mock data, showing that the Bayesian Reconstruction method has the best resolution in extracting peak structures while the Backus-Gilbert and Maximum Entropy methods are somewhat more stable for the flat spectral function. Numerical results are presented for both the elastic case (clover fermions on domain wall configuration with $m_π\sim$ 370 MeV and $a\sim$ 0.06 fm) and a case (anisotropic clover lattice with $m_π\sim$ 380 MeV and $a_t\sim$ 0.035 fm) with large momentum transfer. For the former case, the reconstructed Minkowski hadronic tensor gives precisely the vector charge which proves the feasibility of the approach. While for the latter case, the nucleon resonances and possibly shallow inelastic scattering contributions around $ν=1$ GeV are clearly observed but no information is obtained for higher excited states with $ν>2$ GeV. A check of the effective masses of $ρ$ meson with different lattice setups indicates that, in order to reach higher energy transfers, using lattices with smaller lattice spacings is essential.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$π$N and strangeness sigma terms at the physical point with chiral fermions

Lattice QCD calculations with chiral fermions of the $π$N sigma term $σ_{πN}$ and strangeness sigma term $σ_{sN}$ including chiral interpolation with continuum and volume corrections are provided in this work, with the excited-state contaminations subtracted properly. We calculate the scalar matrix element for the light/strange quark directly and find $σ_{πN}=45.9(7.4)(2.8)$ MeV, with the disconnected insertion part contributing 20(12)(4)\%, and $σ_{sN}=40.2(11.7)(3.5)$ MeV, which is somewhat smaller than $σ_{πN}$. The ratio of the strange/light scalar matrix elements is $y$ = 0.09(3)(1).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Glue Spin of the Proton

We report the progress on the lattice QCD calculation of the glue spin contribution to proton spin. This calculation is carried out with valence overlap fermion on 2+1 flavor DWF gauge configurations at two lattice spacings with the momentum of the frame in the kinematic range $0\leq p^2 \leq 2$ GeV$^2$. A mild frame dependence is observed. The matching and mixing with large-momentum effective field theory are in progress. The unrenormalized result at $p^2=4$ GeV$^2$ with $O(a^2)$ correction gives $S_G$ = 0.13(3).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parton Distribution Function from the Hadronic Tensor on the Lattice

The path-integral formulation of the hadronic tensor W_{μν} of deep inelastic scattering is reviewed. It is shown that there are 3 gauge invariant and topologically distinct contributions. The separation of the connected sea partons from those of the disconnected sea can be achieved with a combination of the global fit of the parton distribution function (PDF), the semi-inclusive DIS data on the strange PDF and the lattice calculation of the ratio of the strange to $u/d$ momentum fraction in the disconnected insertion. We shall discuss numerical issues associated with lattice calculation of the hadronic tensor involving a four-point function, such as large hadron momenta and improved maximum entropy method to obtain the spectral density from the hadronic tensor in Euclidean time. We also draw a comparison between the large momentum approach to the parton distribution function (PDF) and the hadronic tensor approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic method with low mode substitution for nucleon isovector matrix elements

We introduce a stochastic sandwich method with low-mode substitution to evaluate the connected three-point functions. The isovector matrix elements of the nucleon for the axial-vector coupling $g_A^3$, scalar couplings $g_S^3$ and the quark momentum fraction $\langle x\rangle_{u -d}$ are calculated with overlap fermion on 2+1 flavor domain-wall configurations on a $24^3 \times 64$ lattice at $m_π = 330$ MeV with lattice spacing $a = 0.114$ fm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Charm and strange quark masses and $f_{D_s}$ from overlap fermions

We use overlap fermions as valence quarks to calculate meson masses in a wide quark mass range on the $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaborations. The well-defined quark masses in the overlap fermion formalism and the clear valence quark mass dependence of meson masses observed from the calculation facilitate a direct derivation of physical current quark masses through a global fit to the lattice data, which incorporates $O(a^2)$ and $O(m_c^4a^4)$ corrections, chiral extrapolation, and quark mass interpolation. Using the physical masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs, Sommer's scale parameter $r_0$ and the masses of charm quark and strange quark in the $\overline{\rm MS}$ scheme are determined to be $r_0=0.465(4)(9)$ fm, $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.118(6)(24)$ GeV (or $m_c^{\overline{\rm MS}}(m_c)=1.304(5)(20)$ GeV), and $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.101(3)(6)\,{\rm GeV}$, respectively. Furthermore, we observe that the mass difference of the vector meson and the pseudoscalar meson with the same valence quark content is proportional to the reciprocal of the square root of the valence quark masses. The hyperfine splitting of charmonium, $M_{J/ψ}-M_{η_c}$, is determined to be 119(2)(7) MeV, which is in good agreement with the experimental value. We also predict the decay constant of $D_s$ to be $f_{D_s}=254(2)(4)$ MeV. The masses of charmonium $P$-wave states $χ_{c0}, χ_{c1}$ and $h_c$ are also in good agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Glue Spin $S_G$ in The Longitudinally Polarized Nucleon

We present a lattice QCD calculation of the glue spin $S_G$ in the nucleon for the first time. It was recently shown that the first moment of the glue helicity distribution could be obtained through the cross-product of the the electric field $\vec{E}$ and the physical gauge field $\vec{A}_{phys}$ with the non-Abelian Coulomb gauge condition, i.e. $\int d^3 x\, \,\vec{E}(x) \times \vec{A}_{phys}(x)$ in the infinite momentum frame. We use the gauge field tensor from the overlap Dirac operator to check the frame dependence and calculate glue spin with several momenta. The calculation is carried out with valence overlap fermion on 2+1 flavor DWF gauge configurations on the $24^3 \times 64$ lattice with $a^{-1}=1.77$ GeV with the light sea quark mass corresponding to a pion mass of 330 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Meson Mass Decomposition from Lattice QCD

Hadron masses can be decomposed as a sum of quark and glue components which are defined through hadronic matrix elements of QCD operators. The components consist of the quark mass term, the quark energy term, the glue energy term, and the trace anomaly term. We calculate these components for mesons with lattice QCD for the first time. The calculation is carried out with overlap fermion on $2+1$ flavor domain-wall fermion gauge configurations. We confirm that $\sim 50\%$ of the light pion mass comes from the quark mass term and $\sim 10\%$ comes from the quark energy; whereas, while for the $ρ$ meson, the quark energy contributes roughly half of its mass but the quark mass term contributes little. The combined glue components contribute $\sim 40 - 50\%$ for both mesons. It is interesting to observe that the quark mass contribution to the mass of the vector meson is almost linear in quark mass over a large quark mass region below the charm quark mass. For heavy mesons, the quark mass term dominates the masses, while the contribution from the glue components is about $200$ MeV (a bare value around 2GeV) for the heavy pseudoscalar and vector mesons. The charmonium hyperfine splitting is found to be dominated by the quark energy term which is consistent with the picture of the quark potential model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orbital Angular Momentum and Generalized Transverse Momentum Distribution

We show that, when boosted to the infinite momentum frame, the quark and gluon orbital angular momentum operators defined in the nucleon spin sum rule of X. S. Chen et al. are the same as those derived from generalized transverse momentum distributions. This completes the connection between the infinite momentum limit of each term in that sum rule and experimentally measurable observables. We also show that these orbital angular momentum operators can be defined locally, and discuss the strategies of calculating them in lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quark and Glue Components of the Proton Spin from Lattice Calculation

The status of lattice calculations of the quark spin, the quark orbital angular momentum, the glue angular momentum and glue spin in the nucleon is summarized. The quark spin calculation is recently carried out from the anomalous Ward identity with chiral fermions and is found to be small mainly due to the large negative anomaly term which is believed to be the source of the `proton spin crisis'. We also present the first calculation of the glue spin at finite nucleon momenta.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quark Spin in Proton from Anomalous Ward Indentity

We report a quark spin calculation from the anomalous Ward identity with overlap fermions on 2+1 flavor dynamical fermion configurations with light sea quark masses. Such a formulation decomposes the divergence of the flavor-singlet axial-vector current into a quark pseudoscalar term and a triangle anomaly term, flavor by flavor. A large negative contribution from the anomaly term is observed and it is canceled within errors by the contribution from the pseudoscalar term in the disconnected insertion in the heavy quark region. On the other hand, net negative contributions are obtained for the light and strange quarks in the disconnected insertion, since their quark pseudoscalar terms are smaller than that of the heavy quark. Our results are obtained from the 2+1 flavor domain wall fermion configurations on the 24^3*64 lattice with a-1 = 1.78(5) GeV and the light sea quark at m_π = 330 MeV. We use the overlap fermion for the valence and the quark loop so that the renormalization constants Z_m and Z_P cancel in the pseudoscalar operator 2mP. In addition, the overlap Dirac operator is used to calculate the local topological charge in the anomaly so that there is no renormalization for the anomaly term either. In this study, we find the total quark spin to be small mainlyly due to the large negative anomaly term which could be the source for the 'proton spin crisis'.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Revisiting Scalar Glueballs

It is commonly believed that the lowest-lying scalar glueball lies somewhere in the isosinglet scalar mesons $f_0(1370), f_0(1500)$ and $f_0(1710)$ denoted generically by $f_{0}$. In this work we consider lattice calculations and experimental data to infer the glue and $q\bar q$ components of $f_0$. These include the calculations of the scalar glueball masses in quenched and unquenched lattice QCD, measurements of the radiative decays $J/ψ\toγf_{0}$, the ratio of $f_{0}$ decays to $ππ$ and $K\bar K$, the ratio of $J/ψ$ decays to $f_0(1710)ω$ and $f_0(1710)ϕ$, the $f_0$ contributions to $B_s\to J/ψπ^+π^-$, and the near mass degeneracy of $a_0(1450)$ and $K_0^*(1430)$. All analyses suggest the prominent glueball nature of $f_0(1710)$ and the flavor octet structure of $f_0(1500)$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strange quark momentum fraction from overlap fermion

We present a calculation of $< x >_s$ for the strange quark in the nucleon. We also report the ratio of the strange $< x >$ to that of $u/d$ in the disconnected insertion which will be useful in constraining the global fit of parton distribution functions at small $x$. We adopt overlap fermion action on $2 + 1$ flavor domain-wall fermion configurations on the $24^3 \times 64$ lattice with a light sea quark mass which corresponds to $m_π=330$ MeV. Smeared grid $Z_3$ sources are deployed to calculate the nucleon propagator with low-mode substitution. Even-odd grid sources and time-dilution technique with stochastic noises are used to calculate the high mode contribution to the quark loop. Low mode averaging (LMA) for the quark loop is applied to reduce the statistical error of the disconnected insertion calculation. We find the ratio $< x >_s/< x >_{u/d}^{\mathrm{DI}}= 0.78(3)$ in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Parton Orbital Angular Momentum: Status and Prospects

Theoretical progress on the formulation and classification of the quark and gluon orbital angular momenta (OAM) is reviewed. Their relation to parton distributions and open questions and puzzles are discussed. We give a status report on the lattice calculation of the parton kinetic and canonical OAM and point out several strategies to calculate the quark and gluon canonical OAM on the lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From Nuclear Structure to Nucleon Structure

Similarities between nuclear structure study with many-body theory approach and nucleon structure calculations with lattice QCD are pointed out. We will give an example of how to obtain the connected sea partons from a combination of the experimental data, a global fit of parton distribution functions and a lattice calculation. We also present a complete calculation of the quark and glue decomposition of the proton momentum and angular momentum in the quenched approximation. It is found that the quark orbital angular momentum constitutes about 50% of the proton spin.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Meson Mass Decomposition

Hadron masses can be decomposed as a sum of components which are defined through hadronic matrix elements of QCD operators. The components consist of the quark mass term, the quark energy term, the glue energy term and the trace anomaly term. We calculate these components of mesons with lattice QCD for the first time. The calculation is carried out with overlap fermion on $2+1$ flavor domain-wall fermion gauge configurations. We confirm that $\sim 50\%$ of the light pion mass comes from the quark mass and $\sim 10\%$ comes from the quark energy, whereas, the contributions are found to be the other way around for the $ρ$ mass. The combined glue components contribute $\sim 40 - 50\%$ for both mesons. It is interesting to observe that the quark mass contribution to the mass of the vector meson is almost linear in quark mass over a large quark mass region below the charm quark mass. For heavy mesons, the quark mass term dominates the masses, while the contribution from the glue components is about $400\sim500$ MeV for the heavy pseudoscalar and vector mesons. The charmonium hyperfine splitting is found to be dominated by the quark energy term which is consistent with the quark potential model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-perturbative renormalization of overlap quark bilinears on 2+1-flavor domain wall fermion configurations

We present renormalization constants of overlap quark bilinear operators on 2+1-flavor domain wall fermion configurations. This setup is being used by the chiQCD collaboration in calculations of physical quantities such as strangeness in the nucleon and the strange and charm quark masses. The scale independent renormalization constant for the axial vector current is computed using the Ward Identity. The renormalization constants for scalar, pseudoscalar and vector current are calculated in the RI-MOM scheme. Results in the MS-bar scheme are also given. The step scaling function of quark masses in the RI-MOM scheme is computed as well. The analysis uses, in total, six different ensembles of three sea quarks each on two lattices with sizes 24^3x64 and 32^3x64 at spacings a=(1.73 GeV)^{-1} and (2.28 GeV)^{-1}, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Momentum Dependence of the Flavor Structure of the Nucleon Sea

Difference between the $\bar u$ and $\bar d$ sea quark distributions in the proton was first observed in the violation of the Gottfried sum rule in deep-inelastic scattering (DIS) experiments. The parton momentum fraction $x$ dependence of this difference has been measured over the region $0.02 < x < 0.35$ from Drell-Yan and semi-inclusive DIS experiments. The Drell-Yan data suggested a possible sign-change for $\bar d(x)-\bar u(x)$ near $x \sim 0.3$, which has not yet been explained by existing theoretical models. We present an independent evidence for the $\bar d(x)-\bar u(x)$ sign-change at $x \sim 0.3$ from an analysis of the DIS data. We further discuss the $x$-dependence of $\bar d(x)-\bar u(x)$ in the context of meson cloud model and the lattice QCD formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Oscillatory behavior of the domain wall fermions revisited

In the generic domain wall fermion formulation of chiral fermions on the lattice, the zero modes of the four-dimensional Wilson fermion operator with the negative mass parameter $-M_5$ introduce unphysical massive modes propagating in the four-dimensional spacetime. In the free fermion case, the pole mass of this kind of unphysical modes is given by $\tilde{E}=\ln(1-M_5)$, which acquires an imaginary part, $iπ$, when $M_5>1$ and results in an oscillatory behavior of the domain wall fermion propagator in time. The existence of the unphysical modes in the presence of gauge fields is investigated in the mean field approximation, and their physical consequences are discussed. In addition, we also give a semiquantitative criterion for tuning $M_5$ in the realistic numerical study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Roper Puzzle

We carried out a calculation of the Roper state with the Sequential Empirical Bayesian (SEB) method with overlap valence fermion on 2+1-flavor domain-wall fermion configurations on the 24^3 x 64 lattice with a^{-1} = 1.73 GeV. The light sea quark mass corresponds to a pion mass of 330 MeV. The mass of the Roper, chirally extrapolated to the physical pion mass, is 1404(112) MeV which is consistent with the experimental value at 1440 MeV. When compared to the Roper state calculation with variational method for Clover and twisted mass fermions, it is found that the Roper states from SEB with overlap fermion are systematically lower by 400 - 800 MeV for all the quark masses ranging from light to the strange mass region. We study the origin of the difference by exploring the size of the interpolation field in relation to the radial wavefunction of the Roper and also the dynamical influence of the higher Fock space.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Is $1^-+$ Meson a Hybrid?

We calculate the vacuum to meson matrix elements of the dimension-4 operator $\barψγ_4\nblr_i ψ$ and dimension-5 operator $\barψ\epsγ_jψB_k$ of the $1^{-+}$ meson on the lattice and compare them to the corresponding matrix elements of the ordinary mesons to discern if it is a hybrid. For the charmoniums and strange quarkoniums, we find that the matrix elements of $1^{-+}$ are comparable in size as compared to other known $q\bar{q}$ mesons. They are particularly similar to those of the $2^{++}$ meson, since their dimension-4 operators are in the same Lorentz multiplet. Based on these observations, we find no evidence to support the notion that the lowest $1^{-+}$ mesons in the $c\bar{c}$ and $s\bar{s}$ regions are hybrids. As for the exotic quantum number is concerned, the non-relativistic reduction reveals that the leading terms in the dimension-4 and dimension-5 operators of $1^{-+}$ are identical up to a proportional constant and it involves a center-of-mass momentum operator of the quark-antiquark pair. This explains why $1^{-+}$ is an exotic quantum number in the constituent quark model where the center of mass of the $q\bar{q}$ is not a dynamical degree of freedom. Since QCD has gluon fields in the context of the flux-tube which is appropriate for heavy quarkoniums to allow the valence $q\bar{q}$ to recoil against them, it can accommodate such states as $1^{-+}$. By the same token, hadronic models with additional constituents besides the quarks can also accommodate the $q\bar{q}$ center-of-mass motion. To account for the quantum numbers of these $q\bar{q}$ mesons in QCD and hadron models in the non-relativistic case, the parity and total angular momentum should be modified to $P = (-)^{L + l +1}$ and $\vec{J} = \vec{L} + \vec{l} + \vec{S}$, where $L$ is the orbital angular momentum of the $q\bar{q}$ pair in the meson.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-perturbative renormalization of overlap quark bilinears on domain wall fermion configurations

We present renormalization constants of overlap quark bilinear operators on 2+1-flavor domain wall fermion configurations. Both overlap and domain wall fermions have chiral symmetry on the lattice. The scale independent renormalization constant for the local axial vector current is computed using a Ward Identity. The renormalization constants for the scalar, pseudoscalar and vector current are calculated in the RI-MOM scheme. Results in the MS-bar scheme are obtained by using perturbative conversion ratios. The analysis uses in total six ensembles with lattice sizes 24^3x64 and 32^3x64.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Connected-Sea Partons

According to the path-integral formalism of the hadronic tensor, the nucleon sea contains two distinct components called connected sea (CS) and disconnected sea (DS). We discuss how the CS and DS are accessed in the lattice QCD calculation of the moments of the parton distributions. We show that the CS and DS components for $\bar u(x) + \bar d(x)$ can be extracted by using recent data on the strangeness parton distribution, the CT10 global fit, and the lattice result of the ratio of the strange to $u(d)$ moments in the disconnected insertion. The extracted CS and DS for $\bar u(x) + \bar d(x)$ have distinct Bjorken $x$ dependence in qualitative agreement with expectation. The analysis also shows that the momentum fraction of the $\bar u(x) + \bar d(x)$ is about equally divided between CS and DS at $Q^2 = 2.5 {\rm GeV}^2$. Implications on future global analysis for parton distributions are presented.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Critical point of $N_f = 3$ QCD from lattice simulations in the canonical ensemble

A canonical ensemble algorithm is employed to study the phase diagram of $N_f = 3$ QCD using lattice simulations. We lock in the desired quark number sector using an exact Fourier transform of the fermion determinant. We scan the phase space below $T_c$ and look for an S-shape structure in the chemical potential, which signals the coexistence phase of a first order phase transition in finite volume. Applying Maxwell construction, we determine the boundaries of the coexistence phase at three temperatures and extrapolate them to locate the critical point. Using an improved gauge action and improved Wilson fermions on lattices with a spatial extent of $1.8 \fm$ and quark masses close to that of the strange, we find the critical point at $T_E = 0.925(5) T_c$ and baryon chemical potential $μ_B^E = 2.60(8) T_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Finite density phase transition of QCD with $N_f=4$ and $N_f=2$ using canonical ensemble method

In a progress toward searching for the QCD critical point, we study the finite density phase transition of $N_f = 4$ and 2 lattice QCD at finite temperature with the canonical ensemble approach. We develop a winding number expansion method to accurately project out the particle number from the fermion determinant which greatly extends the applicable range of baryon number sectors to make the study feasible. Our lattice simulation was carried out with the clover fermions and improved gauge action. For a given temperature, we calculate the baryon chemical potential from the canonical approach to look for the mixed phase as a signal for the first order phase transition. In the case of $N_f=4$, we observe an "S-shape" structure in the chemical potential-density plane due to the surface tension of the mixed phase in a finite volume which is a signal for the first order phase transition. We use the Maxwell construction to determine the phase boundaries for three temperatures below $T_c$. The intersecting point of the two extrapolated boundaries turns out to be at the expected first order transition point at $T_c$ with $μ= 0$. This serves as a check for our method of identifying the critical point. We also studied the $N_f =2$ case, but do not see a signal of the mixed phase for temperature as low as 0.83 $T_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nucleon strangeness form factors and moments of PDF

The calculation of the nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD is presented. Disconnected insertions are evaluated using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. We find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), which is consistent with experimental values, and has an order of magnitude smaller error. Preliminary results for the strangeness contribution to the second moment of the parton distribution function are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Neutron Electric Dipole Moment at Fixed Topology

We describe the finite volume effects of CP-odd quantities, such as the neutron electric dipole moment and the anapole moment in the $θ$-vacuum, under different topological sectors. We evaluate the three-point Green's functions for the electromagnetic current in a fixed non-trivial topological sector in order to extract these CP-odd observables. We discuss the role of zero modes in the CP-odd Green's function and show that, in the quenched approximation, there is a power divergence in the quark mass for CP-odd quantities at finite volume.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD

We present the N_f=2+1 clover fermion lattice QCD calculation of the nucleon strangeness form factors. We evaluate disconnected insertions using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. We find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), where the first error is statistical, and the second is the uncertainties in Q^2 and chiral extrapolations. This is consistent with experimental values, and has an order of magnitude smaller error.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The calculation of nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD

We study the strangeness electromagnetic form factors of the nucleon from the N_f=2+1 clover fermion lattice QCD calculation. The disconnected insertions are evaluated using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. In addition to increasing the number of Z(4) noises, we find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), where the first error is statistical, and the second is the uncertainties in Q^2 and chiral extrapolations. This is consistent with experimental values, and has an order of magnitude smaller error. We also study the strangeness second moment of the partion distribution function of the nucleon, <x^2>_{s-\bar{s}}.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

2+1 flavor QCD calculation of <x> and <x^2>

We calculate the connected insertions of the nucleon three-point function to study the first few moments of the unpolarized structure functions of the nucleon. (The disconnected insertions are discussed elsewhere in these proceedings). The calculation employs the CP-PACS/JLQCD 2+1 dynamical clover fermions on a 16^3x32 lattice with lattice spacing a=0.1219 fm. The sequential source technique, using non-zero and zero momentum point nucleon field as the secondary source, is applied enabling a study of different currents at various momentum transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Light scalar mesons in 2+1 flavor full QCD

We study the a_0 and K_0^* light scalar mesons in 2+1 flavor full QCD. Particular attention is paid to fitting excited states, with an eye toward determining whether scattering states are revealed. An ultimate goal will be to see how dynamical quarks affect the picture outlined with an earlier quenched study using overlap fermions, namely, that it is the a_0(1450), not the a_0(980), which is the lowest $\bar{q}q$ isovector scalar state.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Strangeness and glue in the nucleon from lattice QCD

We study the strangeness contribution to nucleon matrix elements using Nf=2+1 dynamical clover fermion configurations generated by the CP-PACS/JLQCD collaboration. In order to evaluate the disconnected insertion (DI), we use the Z(4) stochastic method, along with unbiased subtraction from the hopping parameter expansion which reduces the off-diagonal noises in the stochastic method. Furthermore, we find that using many nucleon sources for each configuration is effective in improving the signal. Our results for the quark contribution to the first moment <x>_q in the DI, and the strangeness magnetic moment show that the statistical errors are under control with these techniques. We also study the gluonic contribution to the nucleon using the overlap operator to construct the gauge field tensor, F_{mu,nu}. The application to the calculation of first moment, <x>_G, gives a good signal in quenched lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Winding number expansion for the canonical approach to finite density simulations

The canonical partition function approach was designed to avoid the overlap problem that affects the lattice simulations of nuclear matter at high density. The method employs the projections of the quark determinant on a fix quark number sector. When the quark number is large, the evaluation of the projected determinant becomes numerically unstable. In this paper a different evaluation method based on expanding the determinant in terms of loops winding around the lattice is studied. We show that this method is stable and significantly faster than our original algorithm. This greatly expands the range of quark numbers that we can simulate effectively.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Lattice QCD at finite density via a new canonical approach

We carry out a finite density calculation based on a canonical approach which is designed to address the overlap problem. Two degenerate flavor simulations are performed using Wilson gauge action and Wilson fermions on $4^4$ lattices, at temperatures close to the critical temperature $T_c\approx 170\MeV$ and large densities (5 to 20 times nuclear matter density). In this region, we find that the algorithm works well. We compare our results with those from other approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Progress on a canonical finite density algorithm

We test the finite density algorithm in the canonical ensemble which combines the HMC update with the accept/reject step according to the ratio of the fermion number projected determinant to the unprojected one as a way of avoiding the determinant fluctuation problem. We report our preliminary results on the Polyakov loop in different baryon number sectors which exhibit deconfinement transitions on small lattices. The largest density we obtain around $T_c$ is an order of magnitude larger than that of nuclear matter. From the conserved vector current, we calculate the quark number and verify that the mixing of different baryon sectors is small.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Origin of Difference between $\overline{d}$ and $\overline{u}$ Partons in the Nucleon

Using the Euclidean path-integral formulation for the hadronic tensor, we show that the violation of the Gottfried sum rule does not come from the disconnected quark-loop insertion. Rather, it comes from the connected (quark line) insertion involving quarks propagating in the backward time direction. We demonstrate this by studying sum rules in terms of the scalar and axial- vector matrix elements in lattice gauge calculations. The effects of eliminating backward time propagation are presented.

0 citations0 reviewsNo code linkedOpen access