Showing 1–50 of 182 results for author: Liu, K

Charm physics with overlap fermions on 2+1-flavor domain wall fermion configurations

Authors: Donghao Li, Ying Chen, Ming Gong, Keh-Fei Liu, Zhaofeng Liu, Tingxiao Wang

Abstract: Decay constants of pseudoscalar mesons $D$, $D_s$, $η_c$ and vector mesons $D^*$, $D_s^*$, $J/ψ$ are determined from $N_f=2+1$ lattice QCD at a lattice spacing $a\sim0.08$ fm. For vector mesons, the decay constants defined by tensor currents are given in the $\overline{\rm MS}$ scheme at $2$ GeV. The calculation is performed on domain wall fermion configurations generated by the RBC-UKQCD Collabor… ▽ More Decay constants of pseudoscalar mesons $D$, $D_s$, $η_c$ and vector mesons $D^*$, $D_s^*$, $J/ψ$ are determined from $N_f=2+1$ lattice QCD at a lattice spacing $a\sim0.08$ fm. For vector mesons, the decay constants defined by tensor currents are given in the $\overline{\rm MS}$ scheme at $2$ GeV. The calculation is performed on domain wall fermion configurations generated by the RBC-UKQCD Collaborations and the overlap fermion action is used for the valence quarks. Comparing the current results with our previous ones at a coarser lattice spacing $a\sim0.11$ fm gives us a better understanding of the discretization error. We obtain $f_{D_s^*}^T(\overline{\rm MS},\text{2 GeV})/f_{D_s^*}=0.907(20)$ with a better precision than our previous result. Combining our $f_{D_s^*}=277(11)$ MeV with the total width of $D_s^*$ determined in a recent work gives a branching fraction $4.26(52)\times10^{-5}$ for $D_s^*$ leptonic decay. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, 12 tables

Trace anomaly form factors from lattice QCD

Authors: Bigeng Wang, Fangcheng He, Gen Wang, Terrence Draper, Jian Liang, Keh-Fei Liu, Yi-Bo Yang

Abstract: The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the… ▽ More The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to $Q^2\sim 4.3~\mathrm{GeV}^2$ and the nucleon up to $Q^2\sim 1~\mathrm{GeV}^2$. The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from $\sim 250$ to $\sim 540$ MeV, including the unitary point $\sim 340$ MeV. We calculate the radius of the glue trace anomaly for the pion and the nucleon from the $z$ expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion's form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. △ Less

Submitted 9 May, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: This is the version accepted for publication in Physical Review D

Elastic and resonance structures of the nucleon from hadronic tensor in lattice QCD: implications for neutrino-nucleon scattering and hadron physics

Authors: Jian Liang, Raza Sabbir Sufian, Bigeng Wang, Terrence Draper, Tanjib Khan, Keh-Fei Liu, Yi-Bo Yang

Abstract: Understanding the transitions of nucleons into various resonance structures through electromagnetic interactions plays a pivotal role in advancing our comprehension of the strong interactions within the domain of quark confinement. Furthermore, gaining precise insights into the elastic and resonance structures of nucleons is indispensable for deciphering the physics from neutrino-nucleus scatterin… ▽ More Understanding the transitions of nucleons into various resonance structures through electromagnetic interactions plays a pivotal role in advancing our comprehension of the strong interactions within the domain of quark confinement. Furthermore, gaining precise insights into the elastic and resonance structures of nucleons is indispensable for deciphering the physics from neutrino-nucleus scattering cross sections experimental data, which remain theoretically challenging, even in the context of neutrino-nucleon interactions whose profound understanding is imperative for the neutrino oscillation experiments. One promising avenue involves the direct evaluation of the lepton-nucleon scattering cross sections across quasi-elastic, resonance, shallow-inelastic, and deep inelastic regions, which can be achieved through the hadronic tensor formalism in lattice QCD. In this work, we present the determination of the nucleon's Sachs electric form factor using the hadronic tensor formalism and verify that it is consistent with that from the conventional three-point function calculation. We additionally obtain the transition form factor from the nucleon to its first radial excited state within a finite volume. Consequently, we identify the latter with the nucleon-to-Roper transition form factor $G_E^*(Q^2)$, determine the corresponding longitudinal helicity amplitude $S_{1/2}(Q^2)$ and compare our findings with experimental measurements, for the first time using the hadronic tensor formalism. The limitations and systematic improvements of the approach are also discussed. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 17 pages, 9 figures

Fast Fermion Smearing Scheme with Gaussian-like Profile

Authors: ChuanYang Li, Terrence Draper, Jun Hua, Jian Liang, Keh-Fei Liu, Jun Shi, Nan Wang, Yi-bo Yang

Abstract: We propose a novel smearing scheme which gives a Gaussian-like profile and is more efficient than the traditional Gaussian smearing in terms of computer time consumption. We also carry out a detailed analysis of the profiles, smearing sizes, and the behaviors of hadron effective masses of different smearing schemes, and point out that having a sufficient number of gauge paths in a smearing scheme… ▽ More We propose a novel smearing scheme which gives a Gaussian-like profile and is more efficient than the traditional Gaussian smearing in terms of computer time consumption. We also carry out a detailed analysis of the profiles, smearing sizes, and the behaviors of hadron effective masses of different smearing schemes, and point out that having a sufficient number of gauge paths in a smearing scheme is essential to produce strong smearing effects. For a moderate smearing size $\bar{r}\sim 10a$, the time cost for the novel smearing is less than $1/8$ of that for the traditional Gaussian smearing. In practical lattice calculations with larger smearing sizes or finer lattice spacings the improvement will be more substantial. △ Less

Submitted 27 February, 2024; v1 submitted 3 October, 2023; originally announced October 2023.

Lattice QCD Calculation of Electroweak Box Contributions to Superallowed Nuclear and Neutron Beta Decays

Authors: Peng-Xiang Ma, Xu Feng, Mikhail Gorchtein, Lu-Chang Jin, Keh-Fei Liu, Chien-Yeah Seng, Bi-Geng Wang, Zhao-Long Zhang

Abstract: We present the first lattice QCD calculation of the universal axial $γW$-box contribution $\square_{γW}^{VA}$ to both superallowed nuclear and neutron beta decays. This contribution emerges as a significant component within the theoretical uncertainties surrounding the extraction of $|V_{ud}|$ from superallowed decays. Our calculation is conducted using two domain wall fermion ensembles at the phy… ▽ More We present the first lattice QCD calculation of the universal axial $γW$-box contribution $\square_{γW}^{VA}$ to both superallowed nuclear and neutron beta decays. This contribution emerges as a significant component within the theoretical uncertainties surrounding the extraction of $|V_{ud}|$ from superallowed decays. Our calculation is conducted using two domain wall fermion ensembles at the physical pion mass. To construct the nucleon 4-point correlation functions, we employ the random sparsening field technique. Furthermore, we incorporate long-distance contributions to the hadronic function using the infinite-volume reconstruction method. Upon performing the continuum extrapolation, we arrive at $\square_{γW}^{VA}=3.65(8)_{\mathrm{lat}}(1)_{\mathrm{PT}}\times10^{-3}$. Consequently, this yields a slightly higher value of $|V_{ud}|=0.97386(11)_{\mathrm{exp.}}(9)_{\mathrm{RC}}(27)_{\mathrm{NS}}$, reducing the previous $2.1σ$ tension with the CKM unitarity to $1.8σ$. Additionally, we calculate the vector $γW$-box contribution to the axial charge $g_A$, denoted as $\square_{γW}^{VV}$, and explore its potential implications. △ Less

Submitted 18 April, 2024; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 11 pages, 13 figures

Separation of Infrared and Bulk in Thermal QCD

Authors: Xiao-Lan Meng, Peng Sun, Andrei Alexandru, Ivan Horváth, Keh-Fei Liu, Gen Wang, Yi-Bo Yang

Abstract: New thermal phase of QCD, featuring scale invariance in the infrared (IR), was proposed to exist both in the pure-glue (N$_f$=0) and the ``real-world" (N$_f$=2+1) settings. Among key aspects of the proposal is that the system in this {\em IR phase} separates into two independent components: the scale-invariant IR part and the non-invariant bulk. Such scenario requires non-analyticities in the theo… ▽ More New thermal phase of QCD, featuring scale invariance in the infrared (IR), was proposed to exist both in the pure-glue (N$_f$=0) and the ``real-world" (N$_f$=2+1) settings. Among key aspects of the proposal is that the system in this {\em IR phase} separates into two independent components: the scale-invariant IR part and the non-invariant bulk. Such scenario requires non-analyticities in the theory and, in case of pure-glue QCD, they were found to arise via Anderson-like mobility edges in Dirac spectra ($λ_\fir \!=\! 0$, $\pm λ_\text{A} \!\neq\! 0$) manifested in dimension function $d_\fir(λ)$. Here we present first evidence that this mechanism is also at work in ``real-world QCD" (N$_f$=2+1 theory at physical quark masses and $a\!=\!0.105\,$fm), supporting the existence of the proposed IR regime. Dimensional jump between zero-modes and lowest near-zero modes very close to unity ($d_\fir\!=\!3$ vs $d_\fir \!\simeq \!2$) was found. △ Less

Submitted 26 May, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: 9 pages, 8 figures, references updated

TMD Handbook

Authors: Renaud Boussarie, Matthias Burkardt, Martha Constantinou, William Detmold, Markus Ebert, Michael Engelhardt, Sean Fleming, Leonard Gamberg, Xiangdong Ji, Zhong-Bo Kang, Christopher Lee, Keh-Fei Liu, Simonetta Liuti, Thomas Mehen, Andreas Metz, John Negele, Daniel Pitonyak, Alexei Prokudin, Jian-Wei Qiu, Abha Rajan, Marc Schlegel, Phiala Shanahan, Peter Schweitzer, Iain W. Stewart, Andrey Tarasov , et al. (4 additional authors not shown)

Abstract: This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum… ▽ More This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum and spin structure of hadrons, and are a key ingredient in the description of many collider physics cross sections. Understanding TMDs requires a combination of theoretical techniques from quantum field theory, nonperturbative calculations using lattice QCD, and phenomenological analysis of experimental data. The handbook covers a wide range of topics, from theoretical foundations to experimental analyses, as well as recent developments and future directions. It is intended to provide an essential reference for researchers and graduate students interested in understanding the structure of hadrons and the dynamics of partons in high energy collisions. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 471 pages, many figures

Report number: JLAB-THY-23-3780, LA-UR-21-20798, MIT-CTP/5386

Hadrons, Superconductor Vortices, and Cosmological Constant

Authors: Keh-Fei Liu

Abstract: We explore the roles of the trace anomaly in several hadron properties. We derive the scale invariant expression for the pressure from the gravitational form factors (GFF) of QCD which results in consistent results for the mass and rest energy from the GFF and those from the trace and the Hamiltonian of the energy-momentum tensor (EMT) operators. It is shown that the energy-equilibrium corresponde… ▽ More We explore the roles of the trace anomaly in several hadron properties. We derive the scale invariant expression for the pressure from the gravitational form factors (GFF) of QCD which results in consistent results for the mass and rest energy from the GFF and those from the trace and the Hamiltonian of the energy-momentum tensor (EMT) operators. It is shown that the energy-equilibrium correspondence of hadrons infers an equation of state where the trace anomaly matrix element, emerging from the glue condensate in the vacuum, gives a negative constant pressure that leads to confinement, much like the confinement mechanism for the vortices in type II superconductors where the negative constant pressure is due to the cost of depleting the superconducting condensate. We also note that both the trace anomaly in the QCD energy-momentum tensor and the cosmological constant in Einstein's equation are associated with the metric term which contributes to both energy and pressure. Their difference in terms of the role the pressure plays is discussed. Finally, we note that a lattice calculation of the trace anomaly distribution in the pion has addressed a question about the trace anomaly contribution to the pion mass and suggests that there might be a connection between the conformal symmetry breaking and chiral symmetry breaking in this case. △ Less

Submitted 20 December, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: Version to be published in Phys. Lett. B

RI/(S)MOM renormalizations of overlap quark bilinears with different levels of hypercubic smearing

Authors: Yujiang Bi, Ying Chen, Ming Gong, Fangcheng He, Keh-Fei Liu, Zhaofeng Liu, Yi-Bo Yang, Dian-Jun Zhao

Abstract: On configurations with 2+1-flavor dynamical domain-wall fermions, we calculate the RI/(S)MOM renormalization constants (RC) of overlap quark bilinears. Hypercubic (HYP) smearing is used to construct the overlap Dirac operator. We investigate the possible effects of the smearing on discretization errors in the RCs by varying the level of smearing from 0 to 1 and 2. The lattice is of size… ▽ More On configurations with 2+1-flavor dynamical domain-wall fermions, we calculate the RI/(S)MOM renormalization constants (RC) of overlap quark bilinears. Hypercubic (HYP) smearing is used to construct the overlap Dirac operator. We investigate the possible effects of the smearing on discretization errors in the RCs by varying the level of smearing from 0 to 1 and 2. The lattice is of size $32^3\times64$ and with lattice spacing $1/a=2.383(9)$ GeV. The RCs in the $\overline{\rm MS}$ scheme at 2 GeV are given at the end, with the uncertainty of $Z_T$ reaching $\le1$% for the tensor current. Results of the renormalized quark masses and hadron matrix elements show that the renormalization procedure suppresses the $\sim$ 30% difference of the bare quantities with or without HYP smearing into the 3%-5% level. △ Less

Submitted 1 October, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 33 pages, 14 figures. Statistics increased for computing $Z_A$. Two graphs added for illustration. Main results not changed. References and acknowledgements added. Match the version published on Phys. Rev. D

Journal ref: Phys.Rev.D 108 (2023) 5, 054506

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

Authors: Jian Liang, Andrei Alexandru, Terrence Draper, Keh-Fei Liu, Bigeng Wang, Gen Wang, Yi-Bo Yang

Abstract: 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 spacing… ▽ More 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. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 6 pages, 4 figures, and supplementary materials

Partial wave analysis of the charmed baryon hadronic decay $Λ_c^+\toΛπ^+π^0$

Authors: BESIII Collaboration, M. Ablikim, M. N. Achasov, P. Adlarson, M. Albrecht, R. Aliberti, A. Amoroso, M. R. An, Q. An, X. H. Bai, Y. Bai, O. Bakina, R. Baldini Ferroli, I. Balossino, Y. Ban, V. Batozskaya, D. Becker, K. Begzsuren, N. Berger, M. Bertani, D. Bettoni, F. Bianchi, J. Bloms, A. Bortone, I. Boyko , et al. (555 additional authors not shown)

Abstract: Based on $e^+e^-$ collision samples corresponding to an integrated luminosity of 4.4 $\mbox{fb$^{-1}$}$ collected with the BESIII detector at center-of-mass energies between $4.6\,\,\mathrm{GeV}$ and $4.7\,\,\mathrm{GeV}$, a partial wave analysis of the charmed baryon hadronic decay $Λ_c^+\toΛπ^+π^0$ is performed, and the decays $Λ_c^+\toΛρ(770)^{+}$ and $Λ_c^+\toΣ(1385)π$ are studied for the firs… ▽ More Based on $e^+e^-$ collision samples corresponding to an integrated luminosity of 4.4 $\mbox{fb$^{-1}$}$ collected with the BESIII detector at center-of-mass energies between $4.6\,\,\mathrm{GeV}$ and $4.7\,\,\mathrm{GeV}$, a partial wave analysis of the charmed baryon hadronic decay $Λ_c^+\toΛπ^+π^0$ is performed, and the decays $Λ_c^+\toΛρ(770)^{+}$ and $Λ_c^+\toΣ(1385)π$ are studied for the first time. Making use of the world-average branching fraction $\mathcal{B}(Λ_c^+\toΛπ^+π^0)$, their branching fractions are determined to be \begin{eqnarray*} \begin{aligned} \mathcal{B}(Λ_c^+\toΛρ(770)^+)=&(4.06\pm0.30\pm0.35\pm0.23)\times10^{-2},\\ \mathcal{B}(Λ_c^+\toΣ(1385)^+π^0)=&(5.86\pm0.49\pm0.52\pm0.35)\times10^{-3},\\ \mathcal{B}(Λ_c^+\toΣ(1385)^0π^+)=&(6.47\pm0.59\pm0.66\pm0.38)\times10^{-3},\\ \end{aligned} \end{eqnarray*} where the first uncertainties are statistical, the second are systematic, and the third are from the uncertainties of the branching fractions $\mathcal{B}(Λ_c^+\toΛπ^+π^0)$ and $\mathcal{B}(Σ(1385)\toΛπ)$. In addition, %according to amplitudes determined from the partial wave analysis, the decay asymmetry parameters are measured to be $α_{Λρ(770)^+}=-0.763\pm0.053\pm0.045$, $α_{Σ(1385)^{+}π^0}=-0.917\pm0.069\pm0.056$, and $α_{Σ(1385)^{0}π^+}=-0.789\pm0.098\pm0.056$. △ Less

Submitted 13 December, 2022; v1 submitted 17 September, 2022; originally announced September 2022.

Lattice QCD and Particle Physics

Authors: Andreas S. Kronfeld, Tanmoy Bhattacharya, Thomas Blum, Norman H. Christ, Carleton DeTar, William Detmold, Robert Edwards, Anna Hasenfratz, Huey-Wen Lin, Swagato Mukherjee, Konstantinos Orginos, Richard Brower, Vincenzo Cirigliano, Zohreh Davoudi, Bálint Jóo, Chulwoo Jung, Christoph Lehner, Stefan Meinel, Ethan T. Neil, Peter Petreczky, David G. Richards, Alexei Bazavov, Simon Catterall, Jozef J. Dudek, Aida X. El-Khadra , et al. (57 additional authors not shown)

Abstract: Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). △ Less

Submitted 2 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

Report number: FERMILAB-CONF-22-531-T

Connected and Disconnected Sea Partons from CT18 Parametrization of PDFs

Authors: Tie-Jiun Hou, Mengshi Yan, Jian Liang, Keh-Fei Liu, C. -P. Yuan

Abstract: The separation of the connected and disconnected sea partons, which were uncovered in the Euclidean path-integral formulation of the hadronic tensor, is accommodated with an alternative parametrization of the non-perturbative parton distribution functions in the CT18 global analysis. This is achieved with the help of the distinct small $x$ behaviours of these two sea partons and the constraint fro… ▽ More The separation of the connected and disconnected sea partons, which were uncovered in the Euclidean path-integral formulation of the hadronic tensor, is accommodated with an alternative parametrization of the non-perturbative parton distribution functions in the CT18 global analysis. This is achieved with the help of the distinct small $x$ behaviours of these two sea partons and the constraint from the lattice calculation of the ratio of the strange momentum fraction to that of the $\bar u$ or $\bar d$ in the disconnected insertion. The whole dataset of CT18 is used in this CT18CS fit. The impact of the recent SeaQuest data on the $\bar{d}(x)-\bar{u}(x)$ distribution of CT18CS is also discussed. The separate momentum fractions for the valence, the connected sea and disconnected sea of $u$ and $d$, the strange and the gluon partons are presented at $μ=1.3$ GeV for the first time. They can be compared term-by-term with systematic error controlled lattice calculations. △ Less

Submitted 1 November, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: Revised manuscript accepted for publication in PRD

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

Authors: Fangcheng He, Yu-Jiang Bi, Terrence Draper, Keh-Fei Liu, Zhaofeng Liu, Yi-Bo Yang

Abstract: 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… ▽ More 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. △ Less

Submitted 13 November, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 25 pages, 20 figures; Version accepted for publication in PRD

Muon g-2 with overlap valence fermions

Authors: Gen Wang, Terrence Draper, Keh-Fei Liu, Yi-Bo Yang

Abstract: We present a lattice calculation of the leading order (LO) hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment for the connected light and strange quarks, $a^{\rm W}_{{\rm con}, l/s}$ in the widely used window $t_0=0.4~\mathrm{fm}$, $t_1=1.0~\mathrm{fm}$, $Δ=0.15~\mathrm{fm}$, and also of $a^{\rm S}_{{\rm con}, l/s}$ in the short distance region. We use overlap fe… ▽ More We present a lattice calculation of the leading order (LO) hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment for the connected light and strange quarks, $a^{\rm W}_{{\rm con}, l/s}$ in the widely used window $t_0=0.4~\mathrm{fm}$, $t_1=1.0~\mathrm{fm}$, $Δ=0.15~\mathrm{fm}$, and also of $a^{\rm S}_{{\rm con}, l/s}$ in the short distance region. We use overlap fermions on 4 physical-point ensembles. Two 2+1 flavor RBC/UKQCD ensembles use domain wall fermions (DWF) and Iwasaki gauge actions at $a = 0.084$ and 0.114 fm, and two 2+1+1 flavor MILC ensembles use the highly improved staggered quark (HISQ) and Symanzik gauge actions at $a = 0.088$ and 0.121 fm. We have incorporated infinite volume corrections from 3 additional DWF ensembles at ${\rm L}$ = 4.8, 6.4 and 9.6 fm and physical pion mass. For $a^{\rm W}_{{\rm con}, l}$, we find that our results on the two smaller lattice spacings are consistent with those using the unitary setup, but those at the two coarser lattice spacings are slightly different. Eventually, we predict $a^{\rm W}_{{\rm con}, l}=206.7(1.5)(1.0)$ and $a^{\rm W}_{{\rm con}, s}=26.8(0.1)(0.3)$, using linear extrapolation in $a^2$, with systematic uncertainties estimated from the difference of the central values from the RBC/UKQCD and MILC ensembles. △ Less

Submitted 14 February, 2023; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: matched to the version accepted by PRD

Prospects for precise predictions of $a_μ$ in the Standard Model

Authors: G. Colangelo, M. Davier, A. X. El-Khadra, M. Hoferichter, C. Lehner, L. Lellouch, T. Mibe, B. L. Roberts, T. Teubner, H. Wittig, B. Ananthanarayan, A. Bashir, J. Bijnens, T. Blum, P. Boyle, N. Bray-Ali, I. Caprini, C. M. Carloni Calame, O. Catà, M. Cè, J. Charles, N. H. Christ, F. Curciarello, I. Danilkin, D. Das , et al. (57 additional authors not shown)

Abstract: We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction. We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

Report number: FERMILAB-CONF-22-236-T, LTH 1303, MITP-22-030

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Authors: L. Alvarez Ruso, A. M. Ankowski, S. Bacca, A. B. Balantekin, J. Carlson, S. Gardiner, R. Gonzalez-Jimenez, R. Gupta, T. J. Hobbs, M. Hoferichter, J. Isaacson, N. Jachowicz, W. I. Jay, T. Katori, F. Kling, A. S. Kronfeld, S. W. Li, H. -W. Lin, K. -F. Liu, A. Lovato, K. Mahn, J. Menendez, A. S. Meyer, J. Morfin, S. Pastore , et al. (36 additional authors not shown)

Abstract: Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neut… ▽ More Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators. △ Less

Submitted 20 April, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 81 pages, contribution to Snowmass 2021

Report number: DESY-22-05, FERMILAB-FN-1161-T, MITP-22-027

Lattice QCD Calculations of Parton Physics

Authors: Martha Constantinou, Luigi Del Debbio, Xiangdong Ji, Huey-Wen Lin, Keh-Fei Liu, Christopher Monahan, Kostas Orginos, Peter Petreczky, Jian-Wei Qiu, David Richards, Nobuo Sato, Phiala Shanahan, C. -P. Yuan, Jian-Hui Zhang, Yong Zhao

Abstract: 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 c… ▽ More 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. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 14 pages

Report number: JLAB-THY-22-3564,MIT-CTP/5408,MSUHEP-22-004

Status on Lattice Calculations of the Proton Spin Decomposition

Authors: Keh-Fei Liu

Abstract: 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 clarifi… ▽ More 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. △ Less

Submitted 26 April, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 13 pages, 3 figures, version published in Association of Asia Pacific Physical Societies (AAPPS) Bulletin

Glueballs at Physical Pion Mass

Authors: Feiyu Chen, Xiangyu Jiang, Ying Chen, Keh-Fei Liu, Wei Sun, Yi-Bo Yang

Abstract: We study glueballs on two $N_f=2+1$ RBC/UKQCD gauge ensembles with physical quark masses at two lattice spacings. The statistical uncertainties of the glueball correlation functions are considerably reduced through the cluster decomposition error reduction (CDER) method. The Bethe-Salpeter wave functions are obtained for the scalar, tensor and pseudoscalar glueballs by using spatially extended glu… ▽ More We study glueballs on two $N_f=2+1$ RBC/UKQCD gauge ensembles with physical quark masses at two lattice spacings. The statistical uncertainties of the glueball correlation functions are considerably reduced through the cluster decomposition error reduction (CDER) method. The Bethe-Salpeter wave functions are obtained for the scalar, tensor and pseudoscalar glueballs by using spatially extended glueball operators defined through the gauge potential $A_μ(x)$ in the Coulomb gauge. These wave functions show similar features of non-relativistic two-gluon systems, and they are used to optimize the signals of the related correlation functions at the early time regions. Consequently, the ground state masses can be extracted precisely. To the extent that the excited state contamination is not important, our calculation gives glueball masses at the physical pion mass for the first time. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Proton momentum and angular momentum decompositions with overlap fermions

Authors: Gen Wang, Yi-Bo Yang, Jian Liang, Terrence Draper, Keh-Fei Liu

Abstract: 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 va… ▽ More 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)$. △ Less

Submitted 5 August, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: matched to the version accepted by PRD

Study of the decay $D^+\to K^*(892)^+ K_S^0$ in $D^+\to K^+ K_S^0 π^0$

Authors: BESIII Collaboration, M. Ablikim, M. N. Achasov, P. Adlarson, S. Ahmed, M. Albrecht, R. Aliberti, A. Amoroso, M. R. An, Q. An, X. H. Bai, Y. Bai, O. Bakina, R. Baldini Ferroli, I. Balossino, Y. Ban, K. Begzsuren, N. Berger, M. Bertani, D. Bettoni, F. Bianchi, J. Bloms, A. Bortone, I. Boyko, R. A. Briere , et al. (492 additional authors not shown)

Abstract: Based on an $e^{+}e^{-}$ collision data sample corresponding to an integrated luminosity of 2.93 $\mathrm{fb}^{-1}$ collected with the BESIII detector at $\sqrt{s}=3.773 \mathrm{GeV}$, the first amplitude analysis of the singly Cabibbo-suppressed decay $D^{+}\to K^+ K_S^0 π^0$ is performed. From the amplitude analysis, the $K^*(892)^+ K_S^0$ component is found to be dominant with a fraction of… ▽ More Based on an $e^{+}e^{-}$ collision data sample corresponding to an integrated luminosity of 2.93 $\mathrm{fb}^{-1}$ collected with the BESIII detector at $\sqrt{s}=3.773 \mathrm{GeV}$, the first amplitude analysis of the singly Cabibbo-suppressed decay $D^{+}\to K^+ K_S^0 π^0$ is performed. From the amplitude analysis, the $K^*(892)^+ K_S^0$ component is found to be dominant with a fraction of $(57.1\pm2.6\pm4.2)\%$, where the first uncertainty is statistical and the second systematic. In combination with the absolute branching fraction $\mathcal{B}(D^+\to K^+ K_S^0 π^0)$ measured by BESIII, we obtain $\mathcal{B}(D^+\to K^*(892)^+ K_S^0)=(8.69\pm0.40\pm0.64\pm0.51)\times10^{-3}$, where the third uncertainty is due to the branching fraction $\mathcal{B}(D^+\to K^+ K_S^0 π^0)$. The precision of this result is significantly improved compared to the previous measurement. △ Less

Submitted 16 July, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 12 pages, 15 figures

Journal ref: Phys. Rev. D 104, 012006 (2021)

Proton Mass Decomposition and Hadron Cosmological Constant

Authors: Keh-Fei Liu

Abstract: Lattice results on sigma terms and global analysis of parton momentum fractions are used to give the quark and glue fractions of the proton mass and rest energy. The mass decomposition in terms of the trace of the energy-momentum tensor is renormalization group invariant. The decomposition of the rest energy from the Hamiltonian and the gravitational form factors are scheme and scale dependent. Th… ▽ More Lattice results on sigma terms and global analysis of parton momentum fractions are used to give the quark and glue fractions of the proton mass and rest energy. The mass decomposition in terms of the trace of the energy-momentum tensor is renormalization group invariant. The decomposition of the rest energy from the Hamiltonian and the gravitational form factors are scheme and scale dependent. The separation of the energy-momentum tensor into the traceless part which is composed of the quark and glue parton momentum fractions and the trace part has the minimum scheme dependence. We identify the glue part of the trace anomaly $\langle H_β\rangle $ as the vacuum energy from the glue condensate in the vacuum. From the metric term of the gravitational form factors, which is the stress part of the stress-energy-momentum tensor, we find that the trace part of the rest energy, dominated by $\langle H_β\rangle$, gives a {\it constant} restoring pressure which balances that from the traceless part of the Hamiltonian to confine the hadron, much like the cosmological constant Einstein introduced for a static universe. From a lattice calculation of $\langle H_β\rangle$ in the charmonium, we deduce the associated string tension which turns out to be in good agreement with that from a Cornell potential which fits the charmonium spectrum. △ Less

Submitted 29 September, 2021; v1 submitted 29 March, 2021; originally announced March 2021.

Comments: Version to be published in Phys. Rev. D

Nucleon isovector scalar charge from overlap fermions

Authors: Liuming Liu, Ting Chen, Terrence Draper, Jian Liang, Keh-Fei Liu, Gen Wang, Yi-Bo Yang

Abstract: We calculate the nucleon isovector scalar charge in lattice QCD using overlap fermions on five ensembles of gauge configurations generated by the RBC/UKQCD collaboration using the domain-wall quark action with $2+1$ dynamical flavors. The five ensembles cover five pion masses, $m_π\approx$ 139, 171, 302, 337 and 371 MeV, and four lattice spacings, $a \approx $ 0.06, 0.08, 0.11 and 0.14 fm. Three t… ▽ More We calculate the nucleon isovector scalar charge in lattice QCD using overlap fermions on five ensembles of gauge configurations generated by the RBC/UKQCD collaboration using the domain-wall quark action with $2+1$ dynamical flavors. The five ensembles cover five pion masses, $m_π\approx$ 139, 171, 302, 337 and 371 MeV, and four lattice spacings, $a \approx $ 0.06, 0.08, 0.11 and 0.14 fm. Three to six valence quark masses are computed on each ensemble to investigate the pion mass dependence. The extrapolation to the physical pion mass, continuum and infinite volume limits is obtained by a global fit of all data to a formula originated from partially quenched chiral perturbation theory. The excited-states contamination is carefully analyzed with 3--5 sink-source separations and multi-state fits. Our final result, in the $\overline{\text{MS}}$ scheme at 2 GeV, is $g_{S}^{u-d}= 0.94 (10)_{stat}(8)_{sys}$, where the first error is the statistical error and the second is the systematic error. △ Less

Submitted 12 November, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 18 pages, 7 figures, 4 tables

Journal ref: Phys. Rev. D 104, 094503(2021)

Detecting flavor content of the vacuum using the Dirac operator spectrum

Authors: Jian Liang, Andrei Alexandru, Yu-Jiang Bi, Terrence Draper, Keh-Fei Liu, Yi-Bo Yang

Abstract: We compute the overlap Dirac spectrum on three ensembles generated using 2+1 flavor domain wall fermions. The spectral density is determined up to $λ\sim$100 MeV with sub-percentage statistical uncertainty. The three ensembles have different lattice spacings and two of them have quark masses tuned to the physical point. We show that we can resolve the flavor content of the sea quarks and constrain… ▽ More We compute the overlap Dirac spectrum on three ensembles generated using 2+1 flavor domain wall fermions. The spectral density is determined up to $λ\sim$100 MeV with sub-percentage statistical uncertainty. The three ensembles have different lattice spacings and two of them have quark masses tuned to the physical point. We show that we can resolve the flavor content of the sea quarks and constrain their masses using the Dirac spectral density. We find that the density is close to a constant below $λ\le$ 20 MeV (but 10\% higher than that in the 2-flavor chiral limit) as predicted by chiral perturbative theory ($χ$PT), and then increases linearly due to the strange quark mass. Using the next to leading order $χ$PT, one can extract the light and strange quark masses with $\sim$20\% uncertainties. Using the non-perturbative RI/MOM renormalization, we obtain the chiral condensates at $\overline{\textrm{MS}}$ 2 GeV as $Σ=(260.3(0.7)(1.3)(0.7)(0.8)\ \textrm{MeV})^3$ in the $N_f=2$ (keeping the strange quark mass at the physical point) chiral limit and $Σ_0=(232.6(0.9)(1.2)(0.7)(0.8)\ \textrm{MeV})^3$ in the $N_f=3$ chiral limit, where the four uncertainties come from the statistical fluctuation, renormalization constant, continuum extrapolation and lattice spacing determination. Note that {$Σ/Σ_0=1.40(2)(2)$ is much larger than 1} due to the strange quark mass effect. △ Less

Submitted 8 July, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 10 pages, 6 figures

PDFs and Neutrino-Nucleon Scattering from Hadronic Tensor

Authors: Jian Liang, Keh-Fei Liu

Abstract: 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… ▽ More 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. △ Less

Submitted 27 August, 2020; originally announced August 2020.

Comments: 9 pages, 7 figures, proceedings for the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

PDF in PDFs from Hadronic Tensor and LaMET

Authors: Keh-Fei Liu

Abstract: We point out a problem of the phenomenological definition of the valence partons as the difference between the partons and antipartons in the context of the NNLO evolution equations. After demonstrating that the classification of the parton degrees of freedom (PDF) of the parton distribution functions (PDFs) are the same in the QCD path-intergral formulations of the hadronic tensor and the quasi-P… ▽ More We point out a problem of the phenomenological definition of the valence partons as the difference between the partons and antipartons in the context of the NNLO evolution equations. After demonstrating that the classification of the parton degrees of freedom (PDF) of the parton distribution functions (PDFs) are the same in the QCD path-intergral formulations of the hadronic tensor and the quasi-PDF with LaMET, we resolve the problem by showing that the proper definition of the valence should be in terms of those in the connected insertions only. We also prove that the strange partons appear as the disconnected sea in the nucleon. △ Less

Submitted 17 September, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: References added, version to be published in PRD

Journal ref: Phys. Rev. D 102, 074502 (2020)

Parton distributions and lattice QCD calculations: toward 3D structure

Authors: Martha Constantinou, Aurore Courtoy, Markus A. Ebert, Michael Engelhardt, Tommaso Giani, Tim Hobbs, Tie-Jiun Hou, Aleksander Kusina, Krzysztof Kutak, Jian Liang, Huey-Wen Lin, Keh-Fei Liu, Simonetta Liuti, Cédric Mezrag, Pavel Nadolsky, Emanuele R. Nocera, Fred Olness, Jian-Wei Qiu, Marco Radici, Anatoly Radyushkin, Abha Rajan, Ted Rogers, Juan Rojo, Gerrit Schierholz, C. -P. Yuan , et al. (2 additional authors not shown)

Abstract: The strong force which binds hadrons is described by the theory of Quantum Chromodynamics (QCD). Determining the character and manifestations of QCD is one of the most important and challenging outstanding issues necessary for a comprehensive understanding of the structure of hadrons. Within the context of the QCD parton picture, the Parton Distribution Functions (PDFs) have been remarkably succes… ▽ More The strong force which binds hadrons is described by the theory of Quantum Chromodynamics (QCD). Determining the character and manifestations of QCD is one of the most important and challenging outstanding issues necessary for a comprehensive understanding of the structure of hadrons. Within the context of the QCD parton picture, the Parton Distribution Functions (PDFs) have been remarkably successful in describing a wide variety of processes. However, these PDFs have generally been confined to the description of collinear partons within the hadron. New experiments and facilities provide the opportunity to additionally explore the transverse structure of hadrons which is described by Generalized Parton Distributions (GPDs) and Transverse Momentum Dependent Parton Distribution Functions (TMD PDFs). In our previous review, we compared and contrasted the two main approaches used to determine the collinear PDFs: the first based on perturbative QCD factorization theorems, and the second based on lattice QCD calculations. In the present report, we provide an update of recent progress on the collinear PDFs, and also expand the scope to encompass the generalized PDFs (GPDs and TMD PDFs). We review the current state of the various calculations, and consider what new data might be available in the near future. We also examine how a shared effort can foster dialog between the PDF and Lattice QCD communities, and yield improvements for these generalized PDFs. △ Less

Submitted 7 July, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 2020 PDFLattice Report, 64 pages, 34 figures

Report number: Nikhef 2020-018; MIT-CTP/5213; MSUHEP-20-012; IFJPAN-IV-2020-3; SMU-HEP-20-03

Lattice Calculation of Pion Form Factor with Overlap Fermions

Authors: Gen Wang, Jian Liang, Terrence Draper, Keh-Fei Liu, Yi-Bo Yang

Abstract: We present a precise calculation of the pion form factor using overlap fermions on seven ensembles of 2+1-flavor domain-wall configurations with pion masses varying from 139 to 340 MeV. Taking advantage of the fast Fourier transform and other techniques to access many combinations of source and sink momenta, we find the pion mean square charge radius to be… ▽ More We present a precise calculation of the pion form factor using overlap fermions on seven ensembles of 2+1-flavor domain-wall configurations with pion masses varying from 139 to 340 MeV. Taking advantage of the fast Fourier transform and other techniques to access many combinations of source and sink momenta, we find the pion mean square charge radius to be $\langle {r_π^2} \rangle= 0.430(5)(13)\ {\rm{fm^2}}$, which agrees well with the experimental result, and includes the systematic uncertainties from chiral extrapolation, lattice spacing and finite-volume dependence. We also find that $\langle {r_π^2} \rangle$ depends on both the valence and sea quark masses strongly and predict the pion form factor up to $Q^2 = 1.0 \ {\rm{GeV^2}}$ which agrees with experiments very well. △ Less

Submitted 31 August, 2021; v1 submitted 9 June, 2020; originally announced June 2020.

Comments: 15 pages, 17 figures, matched to the version accepted by PRD

Journal ref: Phys. Rev. D 104, 074502 (2021)

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

Authors: Raza Sabbir Sufian, Tianbo Liu, Andrei Alexandru, Stanley J. Brodsky, Guy F. de Téramond, Hans Günter Dosch, Terrence Draper, Keh-Fei Liu, Yi-Bo Yang

Abstract: 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 m… ▽ More 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. △ Less

Submitted 20 September, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: Published version, Physics Letters B

Report number: JLAB-THY-20-3155, SLAC-PUB-17515

Journal ref: Phys. Lett. B 808 (2020), 135633

The nucleon isovector tensor charge from lattice QCD using chiral fermions

Authors: Derek Horkel, Yujiang Bi, Martha Constantinou, Terrence Draper, Jian Liang, Keh-Fei Liu, Zhaofeng Liu, Yi-Bo Yang

Abstract: 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 uncertaintie… ▽ More 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. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: 8 pages, 5 Figures, 1 Table

Journal ref: Phys. Rev. D 101, 094501 (2020)

Roper State from Overlap Fermions

Authors: Mingyang Sun, Ying Chen, Gen Wang, Andrei Alexandru, Shao-Jing Dong, Terrence Draper, Jacob Fallica, Ming Gong, Frank X. Lee, Anyi Li, Jian Liang, Keh-Fei Liu, Nilmani Mathur, Yi-Bo Yang

Abstract: 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 lowes… ▽ More 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 $Λ$. △ Less

Submitted 30 March, 2020; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: 15 pages, 16 figures, revised manuscript accepted for publication in Phys. Rev. D

Report number: INT-PUB-19-005

Journal ref: Phys. Rev. D 101, 054511 (2020)

Towards the nucleon hadronic tensor from lattice QCD

Authors: Jian Liang, Terrence Draper, Keh-Fei Liu, Alexander Rothkopf, Yi-Bo Yang

Abstract: 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 p… ▽ More 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. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 22 pages, 8 figures

Journal ref: Phys. Rev. D 101, 114503 (2020)

Lattice QCD and Neutrino-Nucleus Scattering

Authors: Andreas S. Kronfeld, David G. Richards, William Detmold, Rajan Gupta, Huey-Wen Lin, Keh-Fei Liu, Aaron S. Meyer, Raza Sufian, Sergey Syritsin

Abstract: 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 calcul… ▽ More 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. △ Less

Submitted 22 April, 2019; originally announced April 2019.

Comments: 33 pp. One of seven whitepapers from the USQCD Collaboration

Report number: FERMILAB-PUB-19-172-T

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

Authors: Jian Liang, Mingyang Sun, Yi-Bo Yang, Terrence Draper, Keh-Fei Liu

Abstract: 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 comple… ▽ More 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. △ Less

Submitted 29 August, 2020; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 11 pages, 9 figures, matched to the version accepted by PRD

Journal ref: Phys. Rev. D 102, 034514 (2020)

Coupling to Multihadron States with Chiral Fermions

Authors: Jacob Fallica, Keh-Fei Liu, Jian Liang, Gen Wang, Yi-Bo Yang

Abstract: Chiral symmerty is presumed to be a crucial component in the strong interaction and QCD, but its role in spectroscopy, especially for baryons, has not been fully explored. Compounding this, chiral fermions are uncommon in lattice calculations due to their expensive nature. We calculate $ηπ$, $Kπ$ and $Nπ$ states with $q\bar{q}$ and $qqq$ interpolation fields at $a=0.114\,\mathrm{fm}$ on a… ▽ More Chiral symmerty is presumed to be a crucial component in the strong interaction and QCD, but its role in spectroscopy, especially for baryons, has not been fully explored. Compounding this, chiral fermions are uncommon in lattice calculations due to their expensive nature. We calculate $ηπ$, $Kπ$ and $Nπ$ states with $q\bar{q}$ and $qqq$ interpolation fields at $a=0.114\,\mathrm{fm}$ on a $48^3\times 96$ mixed-action lattice at the physical pion mass, with domain-wall sea quarks and overlap valence quarks. We study the spectral weights of these states as a function of the valence pion mass, which ranges from $m_π=115-665\,\mathrm{MeV}$, to be compared with the results from non-chiral clover valence quarks on the same domain-wall lattice in order to examine their non-chiral effects, which are expected to decrease with the lattice spacing. △ Less

Submitted 5 December, 2018; originally announced December 2018.

Comments: The 36th Annual International Symposium on Lattice Field Theory, 22-28 July, 2018, Michigan State University. 7 pages, 8 figures

Pion Form Factor with Overlap Fermion

Authors: Gen Wang, Jian Liang, Terrence Draper, Keh-Fei Liu, Yi-Bo Yang

Abstract: We present a calculation of the pion form factor using overlap fermions on 2+1-flavor domain-wall configurations on a $24^3\times 64$ lattice with $a=0.11 \, {\rm{fm}}$ and on a $32^3 \times 64$ lattice with $a=0.143 \, {\rm{fm}}$ generated by the RBC/UKQCD collaboration. Using the multi-mass algorithm, a simulation has been done with various valence quark masses with a range of space-like $Q^2$ f… ▽ More We present a calculation of the pion form factor using overlap fermions on 2+1-flavor domain-wall configurations on a $24^3\times 64$ lattice with $a=0.11 \, {\rm{fm}}$ and on a $32^3 \times 64$ lattice with $a=0.143 \, {\rm{fm}}$ generated by the RBC/UKQCD collaboration. Using the multi-mass algorithm, a simulation has been done with various valence quark masses with a range of space-like $Q^2$ from 0.0 to 0.6 ${\rm{GeV^2}}$. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Comments: 7 pages, 5 figures, 2 table, Proceedings for the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

Observation of $D^+ \to f_0(500) e^+ν_e$ and Improved Measurements of $D \toρe^+ν_e$

Authors: M. Ablikim, M. N. Achasov, S. Ahmed, M. Albrecht, M. Alekseev, A. Amoroso, F. F. An, Q. An, Y. Bai, O. Bakina, R. Baldini Ferroli, Y. Ban, K. Begzsuren, D. W. Bennett, J. V. Bennett, N. Berger, M. Bertani, D. Bettoni, F. Bianchi, E. Boger, I. Boyko, R. A. Briere, H. Cai, X. Cai, A. Calcaterra , et al. (438 additional authors not shown)

Abstract: Using a data sample corresponding to an integrated luminosity of 2.93~fb$^{-1}$ recorded by the BESIII detector at a center-of-mass energy of $3.773$ GeV, we present an analysis of the decays $\bar{D}^0\toπ^+π^0 e^-\barν_e$ and $D^+\toπ^-π^+ e^+ν_e$. By performing a partial wave analysis, the $π^+π^-$ $S$-wave contribution to $D^+\toπ^-π^+ e^+ν_e$ is observed to be $(25.7\pm1.6\pm1.1)$% with a sta… ▽ More Using a data sample corresponding to an integrated luminosity of 2.93~fb$^{-1}$ recorded by the BESIII detector at a center-of-mass energy of $3.773$ GeV, we present an analysis of the decays $\bar{D}^0\toπ^+π^0 e^-\barν_e$ and $D^+\toπ^-π^+ e^+ν_e$. By performing a partial wave analysis, the $π^+π^-$ $S$-wave contribution to $D^+\toπ^-π^+ e^+ν_e$ is observed to be $(25.7\pm1.6\pm1.1)$% with a statistical significance greater than 10$σ$, besides the dominant $P$-wave contribution. This is the first observation of the $S$-wave contribution. We measure the branching fractions $\mathcal{B}(D^{0} \to ρ^- e^+ ν_e) = (1.445\pm 0.058 \pm 0.039) \times10^{-3}$, $\mathcal{B}(D^{+} \to ρ^0 e^+ ν_e) = (1.860\pm 0.070 \pm 0.061) \times10^{-3}$, and $\mathcal{B}(D^{+} \to f_0(500) e^+ ν_e, f_0(500)\toπ^+π^-) = (6.30\pm 0.43 \pm 0.32) \times10^{-4}$. An upper limit of $\mathcal{B}(D^{+} \to f_0(980) e^+ ν_e, f_0(980)\toπ^+π^-) < 2.8 \times10^{-5}$ is set at the 90% confidence level. We also obtain the hadronic form factor ratios of $D\to ρe^+ν_e$ at $q^{2}=0$ assuming the single-pole dominance parameterization: $r_{V}=\frac{V(0)}{A_{1}(0)}=1.695\pm0.083\pm0.051$, $r_{2}=\frac{A_{2}(0)}{A_{1}(0)}=0.845\pm0.056\pm0.039$. △ Less

Submitted 19 February, 2019; v1 submitted 17 September, 2018; originally announced September 2018.

Journal ref: Phys. Rev. Lett. 122, 062001 (2019)

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

Authors: Raza Sabbir Sufian, Keh-Fei Liu, David G. Richards

Abstract: 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 expe… ▽ More 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. △ Less

Submitted 1 February, 2020; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: Published Version (20 pages, 4 figures)

Report number: JLAB-THY-18-2792

Journal ref: JHEP 2001 (2020) 136

Proton Mass Decomposition from the QCD Energy Momentum Tensor

Authors: Yi-Bo Yang, Jian Liang, Yu-Jiang Bi, Ying Chen, Terrence Draper, Keh-Fei Liu, Zhaofeng Liu

Abstract: We report results on the proton mass decomposition and also on related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of $N_f = 2+1$ DWF configurations with three lattice spacings and three volumes, and several pion masses including the physical pion mass. With fully non-perturbative renormalization (and universal normalization on both quark… ▽ More We report results on the proton mass decomposition and also on related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of $N_f = 2+1$ DWF configurations with three lattice spacings and three volumes, and several pion masses including the physical pion mass. With fully non-perturbative renormalization (and universal normalization on both quark and gluon), we find that the quark energy and glue field energy contribute 33(4)(4)\% and 37(5)(4)\% respectively in the $\overline{MS}$ scheme at $μ= 2$ GeV. A quarter of the trace anomaly gives a 23(1)(1)\% contribution to the proton mass based on the sum rule, given 9(2)(1)\% contribution from the $u, d,$ and $s$ quark scalar condensates. The $u,d,s$ and glue momentum fractions in the $\overline{MS}$ scheme are in good agreement with global analyses at $μ= 2$ GeV. △ Less

Submitted 6 September, 2018; v1 submitted 26 August, 2018; originally announced August 2018.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. Lett. 121, 212001 (2018)

Gluon Quasi-PDF From Lattice QCD

Authors: Zhou-You Fan, Yi-Bo Yang, Adam Anthony, Huey-Wen Lin, Keh-Fei Liu

Abstract: We present the first attempt to access the $x$-dependence of the gluon unpolarized parton distribution function (PDF), based on lattice simulations using the large-momentum effective theory (LaMET) approach. The lattice calculation is carried out with pion masses of 340 and 678 MeV on a 2+1-flavor DWF configuration with lattice spacing $a=0.111$~fm, for the gluon quasi-PDF matrix element with the… ▽ More We present the first attempt to access the $x$-dependence of the gluon unpolarized parton distribution function (PDF), based on lattice simulations using the large-momentum effective theory (LaMET) approach. The lattice calculation is carried out with pion masses of 340 and 678 MeV on a 2+1-flavor DWF configuration with lattice spacing $a=0.111$~fm, for the gluon quasi-PDF matrix element with the nucleon momentum up to 0.93 GeV. Taking the normalization from similar matrix elements in the rest frame of the nucleon and pion, our results for these matrix elements are consistent with the Fourier transform of the global fit CT14 and PDF4LHC15 NNLO of the gluon PDF, within statistical uncertainty and the systematic one up to power corrections, perturbative ${\cal O}(α_s)$ matching and the mixing from the quark PDFs. △ Less

Submitted 5 November, 2018; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 5 pages, 5 figures, the version accepted by PRL

Journal ref: Phys. Rev. Lett. 121, 242001 (2018)

Quark spins and Anomalous Ward Identity

Authors: Jian Liang, Yi-Bo Yang, Terrence Draper, Ming Gong, Keh-Fei Liu

Abstract: We calculate the intrinsic quark spin contribution to the total proton spin using overlap valence quarks on three ensembles of $2+1$-flavor RBC/UKQCD domain-wall configurations with different lattice spacings. The lowest pion mass of the ensembles is around 171 MeV which is close to the physical point. With overlap fermions and topological charge derived from the overlap operator, we verify the an… ▽ More We calculate the intrinsic quark spin contribution to the total proton spin using overlap valence quarks on three ensembles of $2+1$-flavor RBC/UKQCD domain-wall configurations with different lattice spacings. The lowest pion mass of the ensembles is around 171 MeV which is close to the physical point. With overlap fermions and topological charge derived from the overlap operator, we verify the anomalous Ward identity between nucleon states with momentum transfer. Both the connected and disconnected insertions of the axial-vector current are calculated. For the disconnected-insertion part, the cluster-decomposition error reduction (CDER) technique is utilized for the lattice with the largest volume and the error can be reduced by $10\%\sim40\%$. Nonperturbative renormalization is carried out and the final results are all reported in the $\overline{\rm MS}$ scheme at 2 GeV. We determine the total quark spin contribution to the nucleon spin to be $ΔΣ=0.401(25)(37)$, which is consistent with the recent global fitting result of experimental data. The isovector axial coupling we obtain in this study is $g_A^3=1.256(16)(30)$, which agrees well with the experimental value of 1.2723(23). △ Less

Submitted 2 November, 2018; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: 37 pages, 16 figures, matched to the version accepted by PRD

Journal ref: Phys. Rev. D 98, 074505 (2018)

Nonperturbatively-renormalized glue momentum fraction at physical pion mass from Lattice QCD

Authors: Yi-Bo Yang, Ming Gong, Jian Liang, Huey-Wen Lin, Keh-Fei Liu, Dimitra Pefkou, Phiala Shanahan

Abstract: We present the first nonperturbatively-renormalized determination of the glue momentum fraction $\langle x \rangle_g$ in the nucleon, based on lattice-QCD simulations at physical pion mass using the cluster-decomposition error reduction (CDER) technique. We provide the first practical strategy to renormalize the glue energy-momentum tensor (EMT) nonperturbatively in the RI/MOM scheme, and convert… ▽ More We present the first nonperturbatively-renormalized determination of the glue momentum fraction $\langle x \rangle_g$ in the nucleon, based on lattice-QCD simulations at physical pion mass using the cluster-decomposition error reduction (CDER) technique. We provide the first practical strategy to renormalize the glue energy-momentum tensor (EMT) nonperturbatively in the RI/MOM scheme, and convert the results to the $\overline{\textrm{MS}}$ scheme with 1-loop matching. The simulation results show that the CDER technique can reduce the statistical uncertainty of its renormalization constant by a factor of ${\cal O}$(300) in calculations using typical state-of-the-art lattice volume, and the nonperturbatively-renormalized $\langle x \rangle_g$ is shown to be independent of the lattice definitions of the glue EMT up to discretization errors. We determine the renormalized $\langle x \rangle_g^{\overline{\textrm{MS}}}(2\textrm{ GeV})$ to be 0.47(4)(11) at physical pion mass, which is consistent with the experimentally-determined value. △ Less

Submitted 28 September, 2018; v1 submitted 1 May, 2018; originally announced May 2018.

Comments: 9 pages, 11 figures, the version accepted by PRD

Journal ref: Phys. Rev. D 98, 074506 (2018)

Comment on `Nucleon Structure Functions from Operator Product Expansion on the Lattice'

Authors: Keh-Fei Liu

Abstract: It is suggested in the paper by A.J. Chambers {\it et al.} (Phys. Rev. Lett. 118, 242001 (2017), arXiv:1703.01153) that the time-ordered current-curent correlator in the nucleon calculated on the lattice is to be identified as the forward Compton amplitude so that it is related to the sum of the even moments of the structure function as in the Minkowski space in the continuum. We point out two pro… ▽ More It is suggested in the paper by A.J. Chambers {\it et al.} (Phys. Rev. Lett. 118, 242001 (2017), arXiv:1703.01153) that the time-ordered current-curent correlator in the nucleon calculated on the lattice is to be identified as the forward Compton amplitude so that it is related to the sum of the even moments of the structure function as in the Minkowski space in the continuum. We point out two problems with this identification. First of all, the current-current correlator defined in the Euclidean space is not analytic everywhere on the rest of the complex $ν$ or $ω$ plane, besides the cuts on the real axis. As such, there is no dispersion relation to relate it to its imaginary part and hence the moments of the structure function. On the lattice, there is an additional difficulty in that the higher dimensional local operators from the operator production expansion (OPE) of the current-current product can mix with lower dimensional higher-twist operators which leads to divergences in the powers of inverse lattice spacing. This mixing needs to be removed before their matrix elements can be identified as the moments of the structure function. △ Less

Submitted 2 December, 2017; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: This article has been withdrawn by the author. The comment on the analyticity of the Euclidean current-current correlator is not relevant and actually validates the approach of the cited work in the stated kinematic region

Parton distributions and lattice QCD calculations: a community white paper

Authors: Huey-Wen Lin, Emanuele R. Nocera, Fred Olness, Kostas Orginos, Juan Rojo, Alberto Accardi, Constantia Alexandrou, Alessandro Bacchetta, Giuseppe Bozzi, Jiunn-Wei Chen, Sara Collins, Amanda Cooper-Sarkar, Martha Constantinou, Luigi Del Debbio, Michael Engelhardt, Jeremy Green, Rajan Gupta, Lucian A. Harland-Lang, Tomomi Ishikawa, Aleksander Kusina, Keh-Fei Liu, Simonetta Liuti, Christopher Monahan, Pavel Nadolsky, Jian-Wei Qiu , et al. (7 additional authors not shown)

Abstract: In the framework of quantum chromodynamics (QCD), parton distribution functions (PDFs) quantify how the momentum and spin of a hadron are divided among its quark and gluon constituents. Two main approaches exist to determine PDFs. The first approach, based on QCD factorization theorems, realizes a QCD analysis of a suitable set of hard-scattering measurements, often using a variety of hadronic obs… ▽ More In the framework of quantum chromodynamics (QCD), parton distribution functions (PDFs) quantify how the momentum and spin of a hadron are divided among its quark and gluon constituents. Two main approaches exist to determine PDFs. The first approach, based on QCD factorization theorems, realizes a QCD analysis of a suitable set of hard-scattering measurements, often using a variety of hadronic observables. The second approach, based on first-principle operator definitions of PDFs, uses lattice QCD to compute directly some PDF-related quantities, such as their moments. Motivated by recent progress in both approaches, in this document we present an overview of lattice-QCD and global-analysis techniques used to determine unpolarized and polarized proton PDFs and their moments. We provide benchmark numbers to validate present and future lattice-QCD calculations and we illustrate how they could be used to reduce the PDF uncertainties in current unpolarized and polarized global analyses. This document represents a first step towards establishing a common language between the two communities, to foster dialogue and to further improve our knowledge of PDFs. △ Less

Submitted 8 February, 2018; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 80 pages, 15 figures, 28 tables, minor typos corrected, published in Progress in Particle and Nuclear Physics

Report number: DESY 17-185, IFJPAN-IV-2017-19, INT-PUB-17-042, MSUHEP-17-017, Nikhef-2017-047, OUTP-17-15P, SMU-HEP-17-08

Lattice calculation of hadronic tensor of the nucleon

Authors: Jian Liang, Keh-Fei Liu, Yi-Bo Yang

Abstract: We report an attempt to calculate the deep inelastic scattering structure functions from the hadronic tensor calculated on the lattice. We used the Backus-Gilbert reconstruction method to address the inverse Laplace transformation for the analytic continuation from the Euclidean to the Minkowski space. We report an attempt to calculate the deep inelastic scattering structure functions from the hadronic tensor calculated on the lattice. We used the Backus-Gilbert reconstruction method to address the inverse Laplace transformation for the analytic continuation from the Euclidean to the Minkowski space. △ Less

Submitted 30 October, 2017; originally announced October 2017.

Comments: 8 pages, 5 figures; Proceedings of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

Proton mass decomposition

Authors: Yi-Bo Yang, Ying Chen, Terrence Draper, Jian Liang, Keh-Fei Liu

Abstract: We report the results on the proton mass decomposition and also on the related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of $N_f = 2+1$ DWF configurations with three lattice spacings and volumes, and several pion masses including the physical pion mass. With 1-loop perturbative calculation and proper normalization of the glue operator, w… ▽ More We report the results on the proton mass decomposition and also on the related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of $N_f = 2+1$ DWF configurations with three lattice spacings and volumes, and several pion masses including the physical pion mass. With 1-loop perturbative calculation and proper normalization of the glue operator, we find that the $u, d,$ and $s$ quark masses contribute 9(2)\% to the proton mass. The quark energy and glue field energy contribute 31(5)\% and 37(5)\% respectively in the $\overline{MS}$ scheme at $μ= 2$ GeV. The trace anomaly gives the remaining 23(1)\% contribution. The $u,d,s$ and glue momentum fractions in the $\overline{MS}$ scheme are consistent with the global analysis at $μ= 2$ GeV. △ Less

Submitted 24 October, 2017; originally announced October 2017.

Comments: 6 pages. Proceedings of the 35th International Symposium on Lattice Field Theory (Lattice2017), Granada, Spain

RI/MOM and RI/SMOM renormalization of overlap quark bilinears on domain wall fermion configurations

Authors: Yujiang Bi, Hao Cai, Ying Chen, Ming Gong, Keh-Fei Liu, Zhaofeng Liu, Yi-Bo Yang

Abstract: Renormalization constants (RCs) of overlap quark bilinear operators on 2+1-flavor domain wall fermion configurations are calculated by using the RI/MOM and RI/SMOM schemes. The scale independent RC for the axial vector current is computed by using a Ward identity. Then the RCs for the quark field and the vector, tensor, scalar and pseudoscalar operators are calculated in both the RI/MOM and RI/SMO… ▽ More Renormalization constants (RCs) of overlap quark bilinear operators on 2+1-flavor domain wall fermion configurations are calculated by using the RI/MOM and RI/SMOM schemes. The scale independent RC for the axial vector current is computed by using a Ward identity. Then the RCs for the quark field and the vector, tensor, scalar and pseudoscalar operators are calculated in both the RI/MOM and RI/SMOM schemes. The RCs are converted to the $\overline{\rm MS}$ scheme and we compare the numerical results from using the two intermediate schemes. The lattice size is $48^3\times96$ and the inverse spacing $1/a = 1.730(4) {\rm~GeV}$. △ Less

Submitted 4 November, 2017; v1 submitted 24 October, 2017; originally announced October 2017.

Comments: Minor changes and updates of Figure 10 and 15 to be more clear

Report number: MSUHEP-17-019

Journal ref: Phys. Rev. D 97, 094501 (2018)

Variance Reduction and Cluster Decomposition

Authors: Keh-Fei Liu, Jian Liang, Yi-Bo Yang

Abstract: It is a common problem in lattice QCD calculation of the mass of the hadron with an annihilation channel that the signal falls off in time while the noise remains constant. In addition, the disconnected insertion calculation of the three-point function and the calculation of the neutron electric dipole moment with the $θ$ term suffer from a noise problem due to the $\sqrt{V}$ fluctuation. We ident… ▽ More It is a common problem in lattice QCD calculation of the mass of the hadron with an annihilation channel that the signal falls off in time while the noise remains constant. In addition, the disconnected insertion calculation of the three-point function and the calculation of the neutron electric dipole moment with the $θ$ term suffer from a noise problem due to the $\sqrt{V}$ fluctuation. We identify these problems to have the same origin and the $\sqrt{V}$ problem can be overcome by utilizing the cluster decomposition principle. We demonstrate this by considering the calculations of the glueball mass, the strangeness content in the nucleon, and the CP violation angle in the nucleon due to the $θ$ term. It is found that for lattices with physical sizes of 4.5 - 5.5 fm, the statistical errors of these quantities can be reduced by a factor of 3 to 4. The systematic errors can be estimated from the Akaike information criterion. For the strangeness content, we find that the systematic error is of the same size as that of the statistical one when the cluster decomposition principle is utilized. This results in a 2 to 3 times reduction in the overall error. △ Less

Submitted 10 July, 2018; v1 submitted 17 May, 2017; originally announced May 2017.

Comments: 7 pages, 5 figures, appendix added to address the systematic error

Report number: MSUHEP-17-006

Journal ref: Phys. Rev. D 97, 034507 (2018)

Sea Quarks Contribution to the Nucleon Magnetic Moment and Charge Radius at the Physical Point

Authors: Raza Sabbir Sufian, Yi-Bo Yang, Jian Liang, Terrence Draper, Keh-Fei Liu

Abstract: We report a comprehensive analysis of the light and strange disconnected-sea quarks contribution to the nucleon magnetic moment, charge radius, and the electric and magnetic form factors. The lattice QCD calculation includes ensembles across several lattice volumes and lattice spacings with one of the ensembles at the physical pion mass. We adopt a model-independent extrapolation of the nucleon ma… ▽ More We report a comprehensive analysis of the light and strange disconnected-sea quarks contribution to the nucleon magnetic moment, charge radius, and the electric and magnetic form factors. The lattice QCD calculation includes ensembles across several lattice volumes and lattice spacings with one of the ensembles at the physical pion mass. We adopt a model-independent extrapolation of the nucleon magnetic moment and the charge radius. We have performed a simultaneous chiral, infinite volume, and continuum extrapolation in a global fit to calculate results in the continuum limit. We find that the combined light and strange disconnected-sea quarks contribution to the nucleon magnetic moment is $μ_M\,(\text{DI})=-0.022(11)(09)\,μ_N$ and to the nucleon mean square charge radius is $\langle r^2\rangle_E\,\text{(DI)}=-0.019(05)(05)$ fm$^2$ which is about $1/3$ of the difference between the $\langle r_p^2\rangle_E$ of electron-proton scattering and that of muonic atom and so cannot be ignored in obtaining the proton charge radius in the lattice QCD calculation. The most important outcome of this lattice QCD calculation is that while the combined light-sea and strange quarks contribution to the nucleon magnetic moment is small at about $1\%$, a negative $2.5(9)\%$ contribution to the proton mean square charge radius and a relatively larger positive $16.3(6.1)\%$ contribution to the neutron mean square charge radius come from the sea quarks in the nucleon. For the first time, by performing global fits, we also give predictions of the light and strange disconnected-sea quarks contributions to the nucleon electric and magnetic form factors at the physical point and in the continuum and infinite volume limits in the momentum transfer range of $0\leq Q^2\leq 0.5$ GeV$^2$. △ Less

Submitted 7 December, 2017; v1 submitted 16 May, 2017; originally announced May 2017.

Comments: Published Version, 26 pages, 8 figures

Journal ref: Phys. Rev. D 96, 114504 (2017)