Showing 1–50 of 60 results for author: Gong, M

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

Form factor for Dalitz decays from $J/ψ$ to light pseudoscalars

Authors: Chunjiang Shi, Ying Chen, Xiangyu Jiang, Ming Gong, Zhaofeng Liu, Wei Sun

Abstract: We calculate the form factor $M(q^2)$ for the Dalitz decay $J/ψ\to γ^*(q^2)η_{(N_f=1)}$ with $η_{(N_f)}$ being the SU($N_f$) flavor singlet pseudoscalar meson. The difference among the partial widths $Γ(J/ψ\to γη_{(N_f)})$ at different $N_f$ can be attributed in part to the $\mathbf{U}_A(1)$ anomaly that induces a $N_f$ scaling. $M(q^2)$'s in $N_f=1,2$ are both well described by the single pole mo… ▽ More We calculate the form factor $M(q^2)$ for the Dalitz decay $J/ψ\to γ^*(q^2)η_{(N_f=1)}$ with $η_{(N_f)}$ being the SU($N_f$) flavor singlet pseudoscalar meson. The difference among the partial widths $Γ(J/ψ\to γη_{(N_f)})$ at different $N_f$ can be attributed in part to the $\mathbf{U}_A(1)$ anomaly that induces a $N_f$ scaling. $M(q^2)$'s in $N_f=1,2$ are both well described by the single pole model $M(q^2)=M(0)/(1-q^2/Λ^2)$. Combined with the known experimental results of the Dalitz decays $J/ψ\to Pe^+e^-$, the pseudoscalar mass $m_P$ dependence of the pole parameter $Λ$ is approximated by $Λ(m_P^2)=Λ_1(1-m_P^2/Λ_2^2)$ with $Λ_1=2.64(4)~\mathrm{GeV}$ and $Λ_2=2.97(33)~\mathrm{GeV}$. These results provide inputs for future theoretical and experimental studies on the Dalitz decays $J/ψ\to Pe^+e^-$. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 9 pages, 5 figures

$X(3872)$ Relevant $D\bar{D}^*$ Scattering in $N_f=2$ Lattice QCD

Authors: Haozheng Li, Chunjiang Shi, Ying Chen, Ming Gong, Juzheng Liang, Zhaofeng Liu, Wei Sun

Abstract: We study the $S$-wave $D\bar{D}^*(I=0)$ scattering at four different pion masses $m_π$ ranging from 250 MeV to 417 MeV from $N_f=2$ lattice QCD. A bound state near the $D\bar{D}^*$ threshold likely exists even after considering the possible left-hand cut effect due to the one pion exchange interaction. At $m_π\approx 417$ MeV where the effective range expansion is valid, the compositeness of the b… ▽ More We study the $S$-wave $D\bar{D}^*(I=0)$ scattering at four different pion masses $m_π$ ranging from 250 MeV to 417 MeV from $N_f=2$ lattice QCD. A bound state near the $D\bar{D}^*$ threshold likely exists even after considering the possible left-hand cut effect due to the one pion exchange interaction. At $m_π\approx 417$ MeV where the effective range expansion is valid, the compositeness of the bound state is $X\approx 1$ and indicates a predominant $D\bar{D}^*$ component. This state may correspond to $X(3872)$. On the other hand, our results of the finite volume energies also hint at the existence of a $1^{++}$ resonance below 4.0 GeV with a width around 60 MeV. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 17 pages, 11 figures

Quark masses and low energy constants in the continuum from the tadpole improved clover ensembles

Authors: Zhi-Cheng Hu, Bo-Lun Hu, Ji-Hao Wang, Ming Gong, Liuming Liu, Peng Sun, Wei Sun, Wei Wang, Yi-Bo Yang, Dian-Jun Zhao

Abstract: We present the light-flavor quark masses and low energy constants using the 2+1 flavor full-QCD ensembles with stout smeared clover fermion action and Symanzik gauge actions. Both the fermion and gauge actions are tadpole improved self-consistently. The simulations are performed on 11 ensembles at 3 lattice spacings $a\in[0.05,0.11]$ fm, 4 spatial sizes $L\in[2.5, 5.1]$ fm, 7 pion masses… ▽ More We present the light-flavor quark masses and low energy constants using the 2+1 flavor full-QCD ensembles with stout smeared clover fermion action and Symanzik gauge actions. Both the fermion and gauge actions are tadpole improved self-consistently. The simulations are performed on 11 ensembles at 3 lattice spacings $a\in[0.05,0.11]$ fm, 4 spatial sizes $L\in[2.5, 5.1]$ fm, 7 pion masses $m_π\in[135,350]$ MeV, and several values of the strange quark mass. The quark mass is defined through the partially conserved axial current (PCAC) relation and renormalized to $\overline{\mathrm{MS}}$ 2 GeV through the intermediate regularization independent momentum subtraction (RI/MOM) scheme. The systematic uncertainty of using the symmetric momentum subtraction (SMOM) scheme is also included. Eventually, we predict $m_u=2.45(22)(20)$ MeV, $m_d=4.74(11)(09)$ MeV, and $m_s=98.8(2.9)(4.7)$ MeV with the systematic uncertainties from lattice spacing determination, continuum extrapolation and renormalization constant included. We also obtain the chiral condensate $Σ^{1/3}=268.6(3.6)(0.7)$ MeV and the pion decay constant $F=86.6(7)(1.4) $ MeV in the $N_f=2$ chiral limit, and the next-to-leading order low energy constants $\ell_3=2.43(54)(05)$ and $\ell_4=4.322(75)(96)$. △ Less

Submitted 7 January, 2024; v1 submitted 1 October, 2023; originally announced October 2023.

Comments: Version accepted by PRD. 7 pages, 4 figures, with more details in the appendix

Decays of $1^{-+}$ Charmoniumlike Hybrid

Authors: Chunjiang Shi, Ying Chen, Ming Gong, Xiangyu Jiang, Zhaofeng Liu, Wei Sun

Abstract: By extracting the transition amplitudes, we give the first lattice QCD prediction of the two-body decay partial widths of the $1^{-+}$ charmoniumlike hybrid $η_{c1}$. Given the calculated mass value $m_{η_{c1}}=4.329(36)$ GeV, the $η_{c1}$ decay is dominated by the open charm modes $D_1\bar{D}$, $D^*\bar{D}$ and $D^*\bar{D}^*$ with partial widths of $258(133)$ MeV, $88(18)$ MeV and $150(118)$ MeV,… ▽ More By extracting the transition amplitudes, we give the first lattice QCD prediction of the two-body decay partial widths of the $1^{-+}$ charmoniumlike hybrid $η_{c1}$. Given the calculated mass value $m_{η_{c1}}=4.329(36)$ GeV, the $η_{c1}$ decay is dominated by the open charm modes $D_1\bar{D}$, $D^*\bar{D}$ and $D^*\bar{D}^*$ with partial widths of $258(133)$ MeV, $88(18)$ MeV and $150(118)$ MeV, respectively. The coupling of $η_{c1}$ to $χ_{c1}$ plus a flavor singlet pseudoscalar is not small, but $χ_{c1}η$ decay is suppressed by the small $η-η'$ mixing angle. The partial width of $η_{c1}\to η_cη'$ is estimated to be around 1 MeV. We suggest experiments to search for $η_{c1}$ in the $P$-wave $D^*\bar{D}$ and $D^*\bar{D}^*$ systems. Especially, the polarization of $D^*\bar{D}^*$ can be used to distinguish the $1^{-+}$ product (total spin $S=1$) from $1^{--}$ products ($S=0$). △ Less

Submitted 19 March, 2024; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 17 pages, 7 figures. Revised to article format

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

$1^{-+}$ Hybrid in $J/ψ$ Radiative Decays from Lattice QCD

Authors: Feiyu Chen, Xiangyu Jiang, Ying Chen, Ming Gong, Zhaofeng Liu, Chunjiang Shi, Wei Sun

Abstract: We present the first theoretical prediction of the production rate of $1^{-+}$ light hybrid meson $η_1$ in $J/ψ$ radiative decays. In the $N_f=2$ lattice QCD formalism with the pion mass $m_π\approx 350$ MeV, the related electromagnetic multipole form factors are extracted from the three-point functions that involve necessarily quark annihilation diagrams, which are calculated through the distilla… ▽ More We present the first theoretical prediction of the production rate of $1^{-+}$ light hybrid meson $η_1$ in $J/ψ$ radiative decays. In the $N_f=2$ lattice QCD formalism with the pion mass $m_π\approx 350$ MeV, the related electromagnetic multipole form factors are extracted from the three-point functions that involve necessarily quark annihilation diagrams, which are calculated through the distillation method. The partial width of $J/ψ\to γη_1$ is determined to be $2.29(77)~\mathrm{eV}$ at the $η_1$ mass $m_{η_1}=2.23(4)$ GeV. If $η_1$ corresponds to the recently observed $η_1(1855)$ in the process $J/ψ\to γη_1(1855)\to γηη'$ by BESIII, then the branching fraction $\mathrm{Br}(J/ψ\to γη_1(1855))$ is estimated to be $6.2(2.2)\times 10^{-5}$, which implies $\mathrm{Br}(η_1(1855)\to ηη')\sim 4.3\%$. △ Less

Submitted 2 April, 2023; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: v2: Published in prd

$T_{cc}^{+}(3875)$ relevant $DD^*$ scattering from $N_f=2$ lattice QCD

Authors: Siyang Chen, Chunjiang Shi, Ying Chen, Ming Gong, Zhaofeng Liu, Wei Sun, Renqiang Zhang

Abstract: The $S$-wave $DD^*$ scattering in the isospin $I=0,1$ channels is studied in $N_f=2$ lattice QCD at $m_π\approx 350$ MeV. It is observed that the $DD^*$ interaction is repulsive in the $I=1$ channel when the $DD^*$ energy is near the $DD^*$ threshold. In contrast, the $DD^*$ interaction in the $I=0$ channel is definitely attractive in a wide range of the $DD^*$ energy. This is consistent with the… ▽ More The $S$-wave $DD^*$ scattering in the isospin $I=0,1$ channels is studied in $N_f=2$ lattice QCD at $m_π\approx 350$ MeV. It is observed that the $DD^*$ interaction is repulsive in the $I=1$ channel when the $DD^*$ energy is near the $DD^*$ threshold. In contrast, the $DD^*$ interaction in the $I=0$ channel is definitely attractive in a wide range of the $DD^*$ energy. This is consistent with the isospin assignment $I=0$ for $T_{cc}^+(3875)$. By analyzing the components of the $DD^*$ correlation functions, it turns out that the quark diagram responsible for the different properties of $I=0,1$ $DD^*$ interactions can be understood as the charged $ρ$ meson exchange effect. This observation provides direct information on the internal dynamics of $T_{cc}^+(3875)$. △ Less

Submitted 17 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 12 pages, 14 figures, version to appear in Physics Letters B

Radiative Decay Width of $J/ψ\to γη_{(2)}$ from $N_f=2$ Lattice QCD

Authors: Xiangyu Jiang, Feiyu Chen, Ying Chen, Ming Gong, Ning Li, Zhaofeng Liu, Wei Sun, Renqiang Zhang

Abstract: The large radiative production rate for pseudoscalar mesons in the $J/ψ$ radiative decay remains elusive. We present the first lattice QCD calculation of partial decay width of $J/ψ$ radiatively decaying into $η_{(2)}$, the $\mathrm{SU(2)}$ flavor singlet pseudoscalar meson, which confirms QCD $\mathrm{U_A(1)}$ anomaly enhancement to the coupling of gluons with flavor singlet pseudoscalar mesons.… ▽ More The large radiative production rate for pseudoscalar mesons in the $J/ψ$ radiative decay remains elusive. We present the first lattice QCD calculation of partial decay width of $J/ψ$ radiatively decaying into $η_{(2)}$, the $\mathrm{SU(2)}$ flavor singlet pseudoscalar meson, which confirms QCD $\mathrm{U_A(1)}$ anomaly enhancement to the coupling of gluons with flavor singlet pseudoscalar mesons. The lattice simulation is carried out using $N_f=2$ lattice QCD gauge configurations at the pion mass $m_π \approx 350$ MeV. In particular, the distillation method has been utilized to calculate light quark loops. The results are reported here with the mass $m_{η_{(2)}}= 718(8)$ MeV and the decay width $Γ(J/ψ\toγη_{(2)})=0.385(45)$ keV. By assuming the dominance of $\mathrm{U_A(1)}$ anomaly and flavor singlet-octet mixing angle $θ=-24.5^\circ$, the production rates for the physical $η$ and $η'$ in $J/ψ$ radiative decay are predicted to be $1.15(14)\times 10^{-3}$ and $4.49(53)\times 10^{-3}$, respectively, which agree well with the experimental measurement data. Our study manifests the potential of lattice QCD studies on the light hadron production in $J/ψ$ radiative decays. △ Less

Submitted 11 February, 2023; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 6 pages, 4 figures, 2 tables

Journal ref: Phys. Rev. Lett. 130, 061901 (2023)

$η$-glueball mixing from $N_f=2$ lattice QCD

Authors: Xiangyu Jiang, Wei Sun, Feiyu Chen, Ying Chen, Ming Gong, Zhaofeng Liu, Renqiang Zhang

Abstract: We perform the first lattice study on the mixing of the isoscalar pseudoscalar meson $η$ and the pseudoscalar glueball $G$ in the $N_f=2$ QCD at the pion mass $m_π\approx 350$ MeV. The $η$ mass is determined to be $m_η=714(6)(16)$ MeV. Through the Witten-Veneziano relation, this value can be matched to a mass value of $\sim 981$ MeV for the $\mathrm{SU(3)}$ counterpart of $η$. Based on a large gau… ▽ More We perform the first lattice study on the mixing of the isoscalar pseudoscalar meson $η$ and the pseudoscalar glueball $G$ in the $N_f=2$ QCD at the pion mass $m_π\approx 350$ MeV. The $η$ mass is determined to be $m_η=714(6)(16)$ MeV. Through the Witten-Veneziano relation, this value can be matched to a mass value of $\sim 981$ MeV for the $\mathrm{SU(3)}$ counterpart of $η$. Based on a large gauge ensemble, the $η-G$ mixing energy and the mixing angle are determined to be $|x|=107(15)(2)$ MeV and $|θ|=3.46(46)^\circ$ from the $η-G$ correlators that are calculated using the distillation method. We conclude that the $η-G$ mixing is tiny and the topology induced interaction contributes most of $η$ mass owing to the QCD $\mathrm{U_A(1)}$ anomaly. △ Less

Submitted 25 May, 2023; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 13 pages, 12 figures

Journal ref: Phys. Rev. D 107, 094510 (2023)

Annihilation diagram contribution to charmonium masses

Authors: Renqiang Zhang, Wei Sun, Feiyu Chen, Ying Chen, Ming Gong, Xiangyu Jiang, Zhaofeng Liu

Abstract: In this work, we generate gauge configurations with $N_f=2$ dynamical charm quarks on anisotropic lattices. The mass shift of $1S$ and $1P$ charmonia owing to the charm quark annihilation effect can be investigated directly in a manner of unitary theory. The distillation method is adopted to treat the charm quark annihilation diagrams at a very precise level. For $1S$ charmonia, the charm quark an… ▽ More In this work, we generate gauge configurations with $N_f=2$ dynamical charm quarks on anisotropic lattices. The mass shift of $1S$ and $1P$ charmonia owing to the charm quark annihilation effect can be investigated directly in a manner of unitary theory. The distillation method is adopted to treat the charm quark annihilation diagrams at a very precise level. For $1S$ charmonia, the charm quark annihilation effect almost does not change the $J/ψ$ mass, but lifts the $η_c$ mass by approximately 3-4 MeV. For $1P$ charmonia, this effect results in positive mass shifts of approximately 1 MeV for $χ_{c1}$ and $h_c$, but decreases the $χ_{c2}$ mass by approximately 3 MeV. We have not obtain a reliable result for the mass shift of $χ_{c0}$. In addition, it is observed that the spin averaged mass of the spin-triplet $1P$ charmonia is in a good agreement with the $h_c$, as expected by the non-relativistic quark model and measured by experiments. △ Less

Submitted 7 March, 2022; v1 submitted 4 October, 2021; originally announced October 2021.

Journal ref: Chinese Phys. C 46(2022) 043102

The Glueball content of $η_c$

Authors: Renqiang Zhang, Wei Sun, Ying Chen, Ming Gong, Long-Cheng Gui, Zhaofeng Liu

Abstract: We carry out the first lattice QCD derivation of the mixing energy and the mixing angle of the pseudoscalar charmonium and glueball on two gauge ensembles with $N_f=2$ degenerate dynamical charm quarks. The mixing energy is determined to be $49(6)$ MeV on the near physical charm ensemble, which seems insensitive to charm quark mass. By the assumption that $X(2370)$ is predominantly a pseudoscalar… ▽ More We carry out the first lattice QCD derivation of the mixing energy and the mixing angle of the pseudoscalar charmonium and glueball on two gauge ensembles with $N_f=2$ degenerate dynamical charm quarks. The mixing energy is determined to be $49(6)$ MeV on the near physical charm ensemble, which seems insensitive to charm quark mass. By the assumption that $X(2370)$ is predominantly a pseudoscalar glueball, the mixing angle is determined to be approximately $4.6(6)^\circ$, which results in a $+3.9(9)$ MeV mass shift of the ground state pseudoscalar charmonium. In the mean time, the mixing can raise the total width of the pseudoscalar charmonium by 7.2(8) MeV, which explains to some extent the relative large total width of the $η_c$ meson. As a result, the branching fraction of $η_c\to γγ$ can be understood in this $c\bar{c}$-glueball mixing framework. On the other hand, the possible discrepancy of the theoretical predictions and the experimental results of the partial width of $J/ψ\toγη_c$ cannot be alleviated by the $c\bar{c}$-glueball mixing picture yet, which demands future precise experimental measurements of this partial width. △ Less

Submitted 24 February, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 9 pages, 9 figures;v2: accepted by Physics Letter B

MSC Class: 81T25

Journal ref: Phys. Lett. B, 827(2022)136960

Charmed and $φ$ meson decay constants from 2+1-flavor lattice QCD

Authors: Ying Chen, Wei-Feng Chiu, Ming Gong, Zhaofeng Liu, Yunheng Ma

Abstract: On a lattice with 2+1-flavor dynamical domain-wall fermions at the physical pion mass, we calculate the decay constants of $D_{s}^{(*)}$, $D^{(*)}$ and $φ$. The lattice size is $48^3\times96$, which corresponds to a spatial extension of $\sim5.5$ fm with the lattice spacing $a\approx 0.114$ fm. For the valence light, strange and charm quarks, we use overlap fermions at several mass points close to… ▽ More On a lattice with 2+1-flavor dynamical domain-wall fermions at the physical pion mass, we calculate the decay constants of $D_{s}^{(*)}$, $D^{(*)}$ and $φ$. The lattice size is $48^3\times96$, which corresponds to a spatial extension of $\sim5.5$ fm with the lattice spacing $a\approx 0.114$ fm. For the valence light, strange and charm quarks, we use overlap fermions at several mass points close to their physical values. Our results at the physical point are $f_D=213(5)$ MeV, $f_{D_s}=249(7)$ MeV, $f_{D^*}=234(6)$ MeV, $f_{D_s^*}=274(7)$ MeV, and $f_φ=241(9)$ MeV. The couplings of $D^*$ and $D_s^*$ to the tensor current ($f_V^T$) can be derived, respectively, from the ratios $f_{D^*}^T/f_{D^*}=0.91(4)$ and $f_{D_s^*}^T/f_{D_s^*}=0.92(4)$, which are the first lattice QCD results. We also obtain the ratios $f_{D^*}/f_D=1.10(3)$ and $f_{D_s^*}/f_{D_s}=1.10(4)$, which reflect the size of heavy quark symmetry breaking in charmed mesons. The ratios $f_{D_s}/f_{D}=1.16(3)$ and $f_{D_s^*}/f_{D^*}=1.17(3)$ can be taken as a measure of SU(3) flavor symmetry breaking. △ Less

Submitted 25 October, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 20 pages, 5 figures, 10 tables; references and acknowledgements added; minor changes, version to be published in Chinese Physics C

Journal ref: Chin. Phys. C45, No. 2 (2021) 023109

Strangeonium-like hybrids on the lattice

Authors: Yunheng Ma, Ying Chen, Ming Gong, Zhaofeng Liu

Abstract: The strangeonium-like $s\bar{s}g$ hybrids are investigated from lattice QCD in the quenched approximation. In the Coulomb gauge, spatially extended operators are constructed for $1^{--}$ and $(0,1,2)^{-+}$ states with the color octet $s\bar{s}$ component being separated from the chromomagnetic field strength by spatial distances $r$, whose matrix elements between the vacuum and the corresponding s… ▽ More The strangeonium-like $s\bar{s}g$ hybrids are investigated from lattice QCD in the quenched approximation. In the Coulomb gauge, spatially extended operators are constructed for $1^{--}$ and $(0,1,2)^{-+}$ states with the color octet $s\bar{s}$ component being separated from the chromomagnetic field strength by spatial distances $r$, whose matrix elements between the vacuum and the corresponding states are interpreted as Bethe-Salpeter (BS) wave functions. In each of the $(1,2)^{-+}$ channels, the masses and the BS wave functions are reliably derived. The $1^{-+}$ ground state mass is around 2.1-2.2 GeV, and that of $2^{-+}$ is around 2.3-2.4 GeV, while the masses of the first excited states are roughly 1.4 GeV higher. This mass splitting is much larger than the expectation of the phenomenological flux-tube model or constituent gluon model for hybrids, which is usually a few hundred MeV. The BS wave functions with respect to $r$ show clear radial nodal structures of non-relativistic two-body system, which imply that $r$ is a meaningful dynamical variable for these hybrids and motivate a color halo picture of hybrids that the color octet $s\bar{s}$ is surrounded by gluonic degrees of freedom. In the $1^{--}$ channel, the properties of the lowest two states comply with those of $φ(1020)$ and $φ(1680)$. We have not obtained convincing information relevant to $φ(2170)$ yet, however, we argue that whether $φ(2170)$ is a conventional $s\bar{s}$ meson or a $s\bar{s}g$ hybrid within the color halo scenario, the ratio of partial decay widths $Γ(φη)$ and $Γ(φη')$ observed by BESIII can be understood by the mechanism of hadronic transition of a strangeonium-like meson along with the $η-η'$ mixing. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 10 pages, 6 figures, to be submitted to Chin. Phys. C

A Lattice Study of the Two-photon Decay Widths for Scalar and Pseudo-scalar Charmonium

Authors: Ying Chen, Ming Gong, Ning Li, Chuan Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Yu Meng, Chao Xiong, Ke-Long Zhang

Abstract: In this exploratory study, two photon decay widths of pseudo-scalar ($η_c$) and scalar ($χ_{c0}$) charmonium are computed using two ensembles of $N_f=2$ twisted mass lattice QCD gauge configurations. The simulation is performed two lattice ensembles with lattice spacings $a=0.067$ fm with size $32^3\times{64}$ and $a=0.085$ fm with size $24^3\times{48}$, respectively. The results for the decay wid… ▽ More In this exploratory study, two photon decay widths of pseudo-scalar ($η_c$) and scalar ($χ_{c0}$) charmonium are computed using two ensembles of $N_f=2$ twisted mass lattice QCD gauge configurations. The simulation is performed two lattice ensembles with lattice spacings $a=0.067$ fm with size $32^3\times{64}$ and $a=0.085$ fm with size $24^3\times{48}$, respectively. The results for the decay widths for the two charmonia are obtained which are in the right ballpark however smaller than the experimental ones. Possible reasons for these discrepancies are discussed. △ Less

Submitted 30 May, 2020; v1 submitted 22 March, 2020; originally announced March 2020.

Comments: version accepted by Chinese Physics C

Lattice QCD package GWU-code and QUDA with HIP

Authors: Yu-Jiang Bi, Yi Xiao, Ming Gong, Wei-Yi Guo, Peng Sun, Shun Xu, Yi-Bo Yang

Abstract: The open source HIP platform for GPU computing provides an uniform framework to support both the NVIDIA and AMD GPUs, and also the possibility to porting the CUDA code to the HIP- compatible one. We present the porting progress on the Overlap fermion inverter (GWU-code) and also the general Lattice QCD inverter package - QUDA. The manual of using QUDA on HIP and also the tips of porting general CU… ▽ More The open source HIP platform for GPU computing provides an uniform framework to support both the NVIDIA and AMD GPUs, and also the possibility to porting the CUDA code to the HIP- compatible one. We present the porting progress on the Overlap fermion inverter (GWU-code) and also the general Lattice QCD inverter package - QUDA. The manual of using QUDA on HIP and also the tips of porting general CUDA code into the HIP framework are also provided. △ Less

Submitted 25 January, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: 6 pages, 1 figures; presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China. Author list and typo fixed

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)

Color halo scenario of charmonium-like hybrids

Authors: Yunheng Ma, Wei Sun, Ying Chen, Ming Gong, Zhaofeng Liu

Abstract: The internal structures of $J^{PC}=1^{--}, (0,1,2)^{-+}$ charmonium-like hybrids are investigated under lattice QCD in the quenched approximation. We define the Bethe-Salpeter wave function $Φ_n(r)$ in the Coulomb gauge as the matrix element of a spatially extended hybrid-like operator $\bar{c}{c}g$ between the vacuum and $n$-th state for each $J^{PC}$ with $r$ being the spatial separation between… ▽ More The internal structures of $J^{PC}=1^{--}, (0,1,2)^{-+}$ charmonium-like hybrids are investigated under lattice QCD in the quenched approximation. We define the Bethe-Salpeter wave function $Φ_n(r)$ in the Coulomb gauge as the matrix element of a spatially extended hybrid-like operator $\bar{c}{c}g$ between the vacuum and $n$-th state for each $J^{PC}$ with $r$ being the spatial separation between a localized $\bar{c}c$ component and the chromomagnetic strength tensor. These wave functions exhibit some similarities for states with the aforementioned different quantum numbers, and their $r$-behaviors (no node for the ground states and one node for the first excited states) imply that $r$ can be a meaningful dynamical variable for these states. Additionally, the mass splittings of the ground states and first excited states of charmonium-like hybrids in these channels are obtained for the first time to be approximately 1.2-1.4 GeV. These results do not support the flux-tube description of heavy-quarkonium-like hybrids in the Born-Oppenheimer approximation. In contrast, a charmonium-like hybrid can be viewed as a "color halo" charmonium for which a relatively localized color octet $\bar{c}c$ is surrounded by gluonic degrees of freedom, which can readily decay into a charmonium state along with one or more light hadrons. The color halo picture is compatible with the decay properties of $Y(4260)$ and suggests LHCb and BelleII to search for $(0,1,2)^{-+}$ charmonium-like hybrids in $χ_{c0,1,2}η$ and $J/ψω(φ)$ final states. △ Less

Submitted 19 August, 2021; v1 submitted 22 October, 2019; originally announced October 2019.

Comments: 17 pages, 9 figures. The contents are considerably enriched, more references are added. Match the publication version in Chin. Phys. C (in press)

Journal ref: Chinese Physics C 45, No. 9 (2021) 093111

A coupled-channel lattice study on the resonance-like structure $Z_c(3900)$

Authors: Ting Chen, Ying Chen, Ming Gong, Chuan Liu, Liuming Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Markus Werner, Jian-Bo Zhang

Abstract: In this exploratory study, near-threshold scattering of $D$ and $\bar{D}^*$ meson is investigated using lattice QCD with $N_f=2+1+1$ twisted mass fermion configurations. The calculation is performed within the coupled-channel Lüscher's finite-size formalism. The study focuses on the channel with $I^G(J^{PC})=1^+(1^{+-})$ where the resonance-like structure $Z_c(3900)$ was discovered. We first ident… ▽ More In this exploratory study, near-threshold scattering of $D$ and $\bar{D}^*$ meson is investigated using lattice QCD with $N_f=2+1+1$ twisted mass fermion configurations. The calculation is performed within the coupled-channel Lüscher's finite-size formalism. The study focuses on the channel with $I^G(J^{PC})=1^+(1^{+-})$ where the resonance-like structure $Z_c(3900)$ was discovered. We first identify the most relevant two channels of the problem and the lattice study is performed within the two-channel scattering model. Combined with a two-channel Ross-Shaw theory, scattering parameters are extracted from the energy levels by solving the generalized eigenvalue problem. Our results on the scattering length parameters suggest that, at the particular lattice parameters that we studied, the best fitted parameters do not correspond to a peak behavior in the elastic scattering cross section near the threshold. Furthermore, within the zero-range Ross-Shaw theory, the scenario of a narrow resonance close to the threshold is disfavored beyond $3σ$ level. △ Less

Submitted 11 August, 2019; v1 submitted 7 July, 2019; originally announced July 2019.

Comments: 15 pages, 6 figures, minor changes. Version accepted by Chinese Physics C

Journal ref: Chin.Phys. C43 (2019) no.10, 103103

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)

Glueball relevant study on isoscalars from $N_f=2$ lattice QCD

Authors: Wei Sun, Long-cheng Gui, Ying Chen, Ming Gong, Zhaofeng Liu

Abstract: We perform a glueball-relevant study on isoscalars based on anisotropic $N_f=2$ lattice QCD gauge configurations. In the scalar channel, we identify the ground state obtained through gluonic operators to be a single-particle state through its dispersion relation. When $q\bar{q}$ operator is included, we find the mass of this state does not change, and the $q\bar{q}$ operator couples very weakly to… ▽ More We perform a glueball-relevant study on isoscalars based on anisotropic $N_f=2$ lattice QCD gauge configurations. In the scalar channel, we identify the ground state obtained through gluonic operators to be a single-particle state through its dispersion relation. When $q\bar{q}$ operator is included, we find the mass of this state does not change, and the $q\bar{q}$ operator couples very weakly to this state. So this state is most likely a glueball state. For pseudoscalars, along with the exiting lattice results, our study implies that both the conventional $q\bar{q}$ state $η_2$ (or $η'$ in flavor $SU(3)$) and a heavier glueball-like state with a mass of roughly 2.6 GeV exist in the spectrum of lattice QCD with dynamical quarks. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: 8 pages, 3 figures, 3 tables, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

Journal ref: EPJ Web of Conferences 175, 05016 (2018)

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)

Glueball spectrum from $N_f=2$ lattice QCD study on anisotropic lattices

Authors: Wei Sun, Long-Cheng Gui, Ying Chen, Ming Gong, Chuan Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Jian-Bo Zhang

Abstract: The lowest-lying glueballs are investigated in lattice QCD using $N_f=2$ clover Wilson fermion on anisotropic lattices. We simulate at two different and relatively heavy quark masses, corresponding to physical pion mass of $m_π\sim 938$ MeV and $650$ MeV. The quark mass dependence of the glueball masses have not been investigated in the present study. Only the gluonic operators built from Wilson l… ▽ More The lowest-lying glueballs are investigated in lattice QCD using $N_f=2$ clover Wilson fermion on anisotropic lattices. We simulate at two different and relatively heavy quark masses, corresponding to physical pion mass of $m_π\sim 938$ MeV and $650$ MeV. The quark mass dependence of the glueball masses have not been investigated in the present study. Only the gluonic operators built from Wilson loops are utilized in calculating the corresponding correlation functions. In the tensor channel, we obtain the ground state mass to be 2.363(39) GeV and 2.384(67) GeV at $m_π\sim 938$ MeV and $650$ MeV, respectively. In the pseudoscalar channel, when using the gluonic operator whose continuum limit has the form of $ε_{ijk}TrB_iD_jB_k$, we obtain the ground state mass to be 2.573(55) GeV and 2.585(65) GeV at the two pion masses. These results are compatible with the corresponding results in the quenched approximation. In contrast, if we use the topological charge density as field operators for the pseudoscalar, the masses of the lowest state are much lighter (around 1GeV) and compatible with the expected masses of the flavor singlet $q\bar{q}$ meson. This indicates that the operator $ε_{ijk}TrB_iD_jB_k$ and the topological charge density couple rather differently to the glueball states and $q\bar{q}$ mesons. The observation of the light flavor singlet pseudoscalar meson can be viewed as the manifestation of effects of dynamical quarks. In the scalar channel, the ground state masses extracted from the correlation functions of gluonic operators are determined to be around 1.4-1.5 GeV, which is close to the ground state masses from the correlation functions of the quark bilinear operators. In all cases, the mixing between glueballs and conventional mesons remains to be further clarified in the future. △ Less

Submitted 7 September, 2018; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: 14 pages, 8 figures, 9 tables; typos fixed, minor modifications, match the published version

Journal ref: Chin.Phys. C42 (2018) no.9, 093103

Quark chiral condensate from the overlap quark propagator

Authors: Chao Wang, Yujiang Bi, Hao Cai, Ying Chen, Ming Gong, Zhaofeng Liu

Abstract: From the overlap lattice quark propagator calculated in the Landau gauge, we determine the quark chiral condensate by fitting operator product expansion formulas to the lattice data. The quark propagators are computed on domain wall fermion configurations generated by the RBC-UKQCD Collaborations with $N_f=2+1$ flavors. Three ensembles with different light sea quark masses are used at one lattice… ▽ More From the overlap lattice quark propagator calculated in the Landau gauge, we determine the quark chiral condensate by fitting operator product expansion formulas to the lattice data. The quark propagators are computed on domain wall fermion configurations generated by the RBC-UKQCD Collaborations with $N_f=2+1$ flavors. Three ensembles with different light sea quark masses are used at one lattice spacing $1/a=1.75(4)$ GeV. We obtain $\langle\barψψ\rangle^{\overline{\rm MS}}(2\mbox{ GeV})=(-305(15)(21)\mbox{ MeV})^3$ in the SU(2) chiral limit. △ Less

Submitted 14 December, 2016; originally announced December 2016.

Comments: 17 pages, 6 figures

Journal ref: Chin. Phys. C41 (2017) no.5, 053102

Exotic vector charmonium and its leptonic decay width

Authors: Ying Chen, Wei-Feng Chiu, Ming Gong, Long-Cheng Gui, Zhaofeng Liu

Abstract: We propose a novel type of interpolating field operators, which manifests the hybrid-like configuration that the charm quark-antiquark pair recoils against gluonic degrees of freedom. A heavy vector charmonium-like state with a mass of $4.33(2)\,{\rm GeV}$ is disentangled from the conventional charmonium states in the quenched approximation. This state has affinity for the hybrid-like operators bu… ▽ More We propose a novel type of interpolating field operators, which manifests the hybrid-like configuration that the charm quark-antiquark pair recoils against gluonic degrees of freedom. A heavy vector charmonium-like state with a mass of $4.33(2)\,{\rm GeV}$ is disentangled from the conventional charmonium states in the quenched approximation. This state has affinity for the hybrid-like operators but couples less to the relevant quark bilinear operator. We also try to extract its leptonic decay constant and give a tentative upper limit that it is less than one tenth of that of $J/ψ$, which corresponds to a leptonic decay width about dozens of eV. The connection of this state with $X(4260)$ is also discussed. △ Less

Submitted 12 April, 2016; originally announced April 2016.

Comments: 9 pages, 7 figures, to be submitted to Chinese Physics C

Journal ref: Chinese Physics C 40, 081002 (2016)

A Lattice Study of $(\bar{D}_1 D^{*})^\pm$ Near-threshold Scattering

Authors: Ting Chen, Ying Chen, Ming Gong, Yu-Hong Lei, Ning Li, Chuan Liu, Yu-Bin Liu, Zhao-Feng Liu, Jian-Ping Ma, Zhan-Lin Wang, Jian-Bo Zhang

Abstract: In this exploratory lattice study, low-energy near threshold scattering of the $(\bar{D}_1 D^{*})^\pm$ meson system is analyzed using lattice QCD with $N_f=2$ twisted mass fermion configurations. Both s-wave ($J^P=0^-$) and p-wave ($J^P=1^+$) channels are investigated. It is found that the interaction between the two charmed mesons is attractive near the threshold in both channels. This calculatio… ▽ More In this exploratory lattice study, low-energy near threshold scattering of the $(\bar{D}_1 D^{*})^\pm$ meson system is analyzed using lattice QCD with $N_f=2$ twisted mass fermion configurations. Both s-wave ($J^P=0^-$) and p-wave ($J^P=1^+$) channels are investigated. It is found that the interaction between the two charmed mesons is attractive near the threshold in both channels. This calculation provides some hints in the searching of resonances or bound states around the threshold of $(\bar{D}_1 D^{*})^\pm$ system. △ Less

Submitted 28 June, 2016; v1 submitted 31 January, 2016; originally announced February 2016.

Comments: 20 pages, 15 figures, matches the version on PRD

Journal ref: Phys. Rev. D 93, 114501 (2016)

Two Photon Decays of $η_c$ from Lattice QCD

Authors: Ting Chen, Ying Chen, Ming Gong, Yu-Hong Lei, Ning Li, Chuan Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Wei-Feng Qiu, Zhan-Lin Wang, Jian-Bo Zhang

Abstract: We present an exploratory lattice study for the two-photon decay of $η_c$ using $N_f=2$ twisted mass lattice QCD gauge configurations generated by the European Twisted Mass Collaboration. Two different lattice spacings of $a=0.067$fm and $a=0.085$fm are used in the study, both of which are of physical size of 2$fm$. The decay widths are found to be $1.025(5)$KeV for the coarser lattice and… ▽ More We present an exploratory lattice study for the two-photon decay of $η_c$ using $N_f=2$ twisted mass lattice QCD gauge configurations generated by the European Twisted Mass Collaboration. Two different lattice spacings of $a=0.067$fm and $a=0.085$fm are used in the study, both of which are of physical size of 2$fm$. The decay widths are found to be $1.025(5)$KeV for the coarser lattice and $1.062(5)$KeV for the finer lattice respectively where the errors are purely statistical. A naive extrapolation towards the continuum limit yields $Γ\simeq 1.122(14)$KeV which is smaller than the previous quenched result and most of the current experimental results. Possible reasons are discussed. △ Less

Submitted 28 June, 2016; v1 submitted 30 January, 2016; originally announced February 2016.

Comments: 13 pages, 7 figures; matches the published version

Spectrum and Bethe-Salpeter amplitudes of $Ω$ baryons from lattice QCD

Authors: Jian Liang, Wei Sun, Ying Chen, Wei-Feng Chiu, Ming Gong, Chuan Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Jian-Bo Zhang

Abstract: The $Ω$ baryons with $J^P=3/2^\pm, 1/2^\pm$ are studied on the lattice in the quenched approximation. Their mass levels are ordered as $M_{3/2^+}<M_{3/2^-}\approx M_{1/2^-}<M_{1/2^+}$, as is expected from the constituent quark model. The mass values are also close to those of the four $Ω$ states observed in experiments, respectively. We calculate the Bethe-Salpeter amplitudes of $Ω(3/2^+)$ and… ▽ More The $Ω$ baryons with $J^P=3/2^\pm, 1/2^\pm$ are studied on the lattice in the quenched approximation. Their mass levels are ordered as $M_{3/2^+}<M_{3/2^-}\approx M_{1/2^-}<M_{1/2^+}$, as is expected from the constituent quark model. The mass values are also close to those of the four $Ω$ states observed in experiments, respectively. We calculate the Bethe-Salpeter amplitudes of $Ω(3/2^+)$ and $Ω(1/2^+)$ and find there is a radial node for the $Ω(1/2^+)$ Bethe-Salpeter amplitude, which may imply that $Ω(1/2^+)$ is an orbital excitation of $Ω$ baryons as a member of the $(D,L_N^P)=(70,0_2^+)$ supermultiplet in the $SU(6)\bigotimes O(3)$ quark model description. Our results are helpful for identifying the quantum number of experimentally observed $Ω$ states. △ Less

Submitted 12 April, 2016; v1 submitted 13 November, 2015; originally announced November 2015.

Comments: 7 pages, 5 figures, submitted to Chinese Physics C

Journal ref: Chin. Phys. C 40 (2016) 041001

Strange and Charm Quark Spins from Anomalous Ward Identity

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

Abstract: We present a calculation of the strange and charm quark contributions to the nucleon spin from the anomalous Ward identity (AWI). It is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations on a $24^3 \times 64$ lattice with the light sea mass at $m_π = 330$ MeV. To satisfy the AWI, the overlap fermion for the pseudoscalar density and the overlap Dirac operat… ▽ More We present a calculation of the strange and charm quark contributions to the nucleon spin from the anomalous Ward identity (AWI). It is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations on a $24^3 \times 64$ lattice with the light sea mass at $m_π = 330$ MeV. To satisfy the AWI, the overlap fermion for the pseudoscalar density and the overlap Dirac operator for the topological density, which do not have multiplicative renormalization, are used to normalize the form factor of the local axial-vector current at finite $q^2$. For the charm quark, we find that the negative pseudoscalar term almost cancels the positive topological term. For the strange quark, the pseudoscalar term is less negative than that of the charm. By imposing the AWI, the strange $g_A(q^2)$ at $q^2 =0$ is obtained by a global fit of the pseudoscalar and the topological form factors, together with $g_A(q^2)$ and the induced pseudoscalar form factor $h_A(q^2)$ at finite $q^2$. The chiral extrapolation to the physical pion mass gives $Δs + Δ{\bar{s}} = -0.0403(44)(78)$. △ Less

Submitted 24 April, 2017; v1 submitted 11 November, 2015; originally announced November 2015.

Comments: 8 pages, 9 figures. Updated version where a sign error is corrected

Journal ref: Phys. Rev. D 95, 114509 (2017)

Diquark mass differences from unquenched lattice QCD

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

Abstract: We calculate diquark correlation functions in the Landau gauge on the lattice using overlap valence quarks and 2+1-flavor domain wall fermion configurations. Quark masses are extracted from the scalar part of quark propagators in the Landau gauge. Scalar diquark quark mass difference and axial vector scalar diquark mass difference are obtained for diquarks composed of two light quarks and of a str… ▽ More We calculate diquark correlation functions in the Landau gauge on the lattice using overlap valence quarks and 2+1-flavor domain wall fermion configurations. Quark masses are extracted from the scalar part of quark propagators in the Landau gauge. Scalar diquark quark mass difference and axial vector scalar diquark mass difference are obtained for diquarks composed of two light quarks and of a strange and a light quark. Light sea quark mass dependence of the results is examined. Two lattice spacings are used to check the discretization effects. The coarse and fine lattices are of sizes $24^3\times64$ and $32^3\times64$ with inverse spacings $1/a=1.75(4) {\rm~GeV}$ and $2.33(5) {\rm~GeV}$, respectively. △ Less

Submitted 25 October, 2015; originally announced October 2015.

Comments: 9 figures

Journal ref: Chin.Phys. C40 (2016) no.7, 073106

Stochastic method with low mode substitution for nucleon isovector matrix elements

Authors: Yi-Bo Yang, Andrei Alexandru, Terrence Draper, Ming Gong, Keh-Fei Liu

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

Submitted 23 January, 2016; v1 submitted 15 September, 2015; originally announced September 2015.

Comments: 15 pages, 13 figures, the version accepted by PRD

Journal ref: Phys. Rev. D 93, 034503 (2016)

Quark Spin in Proton from Anomalous Ward Indentity

Authors: Yi-Bo Yang, Ming Gong, Keh-Fei Liu, Mingyang Sun

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

Submitted 15 April, 2015; originally announced April 2015.

Comments: Talk presented at the 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University New York, NY. 7 pages, 4 figures

Journal ref: PoS(LATTICE2014)138

Phase transition in finite density and temperature lattice QCD

Authors: Rui Wang, Ying Chen, Ming Gong, Chuan Liu, Yu-Bin Liu, Zhao-Feng Liu, Jian-Ping Ma, Xiang-Fei Meng, Jian-Bo Zhang

Abstract: We investigate the behavior of the chiral condensate in lattice QCD at finite temperature and finite chemical potential. The study was done using two flavors of light quarks and with a series of $β$ and $ma$ at the lattice size $24\times12^{2}\times6$. The calculation was done in the Taylar expansion formalism. We are able to calculate the first and second order derivatives of… ▽ More We investigate the behavior of the chiral condensate in lattice QCD at finite temperature and finite chemical potential. The study was done using two flavors of light quarks and with a series of $β$ and $ma$ at the lattice size $24\times12^{2}\times6$. The calculation was done in the Taylar expansion formalism. We are able to calculate the first and second order derivatives of $\langle\barψψ\rangle$ in both isoscalar and isovector channels. With the first derivatives being small, we find that the second derivatives are sizable close to the phase transition and the magnitude of $\barψψ$ decreases under the influence of finite chemical potential in both channels. △ Less

Submitted 8 April, 2015; originally announced April 2015.

Comments: 4 pages, 4 figures, to be submitted to Chinese Physics C

Low-energy Scattering of $(D^{*}\bar{D}^{*})^\pm$ System and the Resonance-like Structure $Z_c(4025)$

Authors: Ying Chen, Ming Gong, Yu-Hong Lei, Ning Li, Jian Liang, Chuan Liu, Jin-Long Liu, Yong-Fu Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Zhan-Lin Wang, Jian-Bo Zhang

Abstract: In this paper, low-energy scattering of the $(D^{*}\bar{D}^{*})^\pm$ meson system is studied within Lüscher's finite-size formalism using $N_{f}=2$ twisted mass gauge field configurations. With three different pion mass values, the $s$-wave threshold scattering parameters, namely the scattering length $a_0$ and the effective range $r_0$, are extracted in $J^P=1^+$ channel. Our results indicate tha… ▽ More In this paper, low-energy scattering of the $(D^{*}\bar{D}^{*})^\pm$ meson system is studied within Lüscher's finite-size formalism using $N_{f}=2$ twisted mass gauge field configurations. With three different pion mass values, the $s$-wave threshold scattering parameters, namely the scattering length $a_0$ and the effective range $r_0$, are extracted in $J^P=1^+$ channel. Our results indicate that, in this particular channel, the interaction between the two vector charmed mesons is weakly repulsive in nature hence do not support the possibility of a shallow bound state for the two mesons, at least for the pion mass values being studied. This study provides some useful information on the nature of the newly discovered resonance-like structure $Z_c(4025)$ observed in various experiments. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: 11 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:1403.1318

Journal ref: Phys. Rev. D 92, 054507 (2015)

Strange quark momentum fraction from overlap fermion

Authors: Mingyang Sun, Yi-Bo Yang, Keh-Fei Liu, Ming Gong

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

Submitted 9 March, 2015; v1 submitted 19 February, 2015; originally announced February 2015.

Comments: The 32nd International Symposium on Lattice Field Theory; Fix spelling of a coauthor's name (add dash)

Wave functions of $SU(3)$ pure gauge glueballs on the lattice

Authors: Jian Liang, Ying Chen, Wei-Feng Chiu, Long-Cheng Gui, Ming Gong, Zhaofeng Liu

Abstract: The Bethe-Salpeter wave functions of $SU(3)$ pure gauge glueballs are revisited in this study. The ground and the first excited states of scalar and tensor glueballs are identified unambiguously by using the variational method on the basis of large operator sets. We calculate their wave functions in the Coulomb gauge and use two lattices with different lattice spacings to check the discretization… ▽ More The Bethe-Salpeter wave functions of $SU(3)$ pure gauge glueballs are revisited in this study. The ground and the first excited states of scalar and tensor glueballs are identified unambiguously by using the variational method on the basis of large operator sets. We calculate their wave functions in the Coulomb gauge and use two lattices with different lattice spacings to check the discretization artifacts. For ground states, the radial wave functions are approximately Gaussian and the size of the tensor is twice as large as that of the scalar. For the first excited states, the radial nodes are clearly observed for both the scalar and the tensor glueballs, such that they can be interpreted as the first radial excitations. These observations may shed light on the theoretical understanding of the inner structure of glueballs. △ Less

Submitted 30 June, 2015; v1 submitted 4 November, 2014; originally announced November 2014.

Comments: published in PRD

Journal ref: Phys. Rev. D 91, 054513 (2015)

Meson Mass Decomposition

Authors: Yi-Bo Yang, Ying Chen, Terrence Draper, Ming Gong, Keh-Fei Liu, Zhaofeng Liu, Jian-Ping Ma

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

Submitted 4 November, 2014; originally announced November 2014.

Comments: 7 Pages, 4 figures, contribution to the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June 2014, Columbia University, New York, NY, USA

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

Authors: Yi-Bo Yang, Ying Chen, Andrei Alexandru, Shao-Jing Dong, Terrence Draper, Ming Gong, Frank X. Lee, Anyi Li, Keh-Fei Liu, Zhaofeng Liu, Michael Lujan

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

Submitted 30 July, 2015; v1 submitted 13 October, 2014; originally announced October 2014.

Comments: 19 pages, 15 figures, revised version accepted for publication in PRD

Journal ref: Phys. Rev. D 92, 034517 (2015)

Meson Mass Decomposition from Lattice QCD

Authors: Yi-Bo Yang, Ying Chen, Terrence Draper, Ming Gong, Keh-Fei Liu, Zhaofeng Liu, Jian-Ping Ma

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

Submitted 23 April, 2015; v1 submitted 17 May, 2014; originally announced May 2014.

Comments: 8 pages, 12 figures, revised version accept to PRD. Typo corrected

Journal ref: Phys. Rev. D 91, 074516 (2015)

The Roper Puzzle

Authors: Keh-Fei Liu, Ying Chen, Ming Gong, Raza Sufian, Mingyang Sun, Anyi Li

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

Submitted 26 March, 2014; originally announced March 2014.

Comments: 7 pages, 6 figures, presented at the 31st International Symposium on Lattice Field Theory (LATTICE 2013), 29 July - 3 August 2013, Mainz, Germany

Low-energy Scattering of $(D\bar{D}^{*})^\pm$ System And the Resonance-like Structure $Z_c(3900)$

Authors: Ying Chen, Ming Gong, Yu-Hong Lei, Ning Li, Jian Liang, Chuan Liu, Hang Liu, Jin-Long Liu, Liuming Liu, Yong-Fu Liu, Yu-Bin Liu, Zhaofeng Liu, Jian-Ping Ma, Zhan-Lin Wang, Yi-Bo Yang, Jian-Bo Zhang

Abstract: In this exploratory lattice study, low-energy scattering of the $(D\bar{D}^{*})^\pm$ meson system is analyzed using lattice QCD with $N_f=2$ twisted mass fermion configurations with three pion mass values. The calculation is performed within single-channel Lüscher's finite-size formalism. The threshold scattering parameters, namely the scattering length $a_0$ and the effective range $r_0$, for the… ▽ More In this exploratory lattice study, low-energy scattering of the $(D\bar{D}^{*})^\pm$ meson system is analyzed using lattice QCD with $N_f=2$ twisted mass fermion configurations with three pion mass values. The calculation is performed within single-channel Lüscher's finite-size formalism. The threshold scattering parameters, namely the scattering length $a_0$ and the effective range $r_0$, for the $s$-wave scattering in $J^P=1^+$ channel are extracted. For the cases in our study, the interaction between the two charmed mesons is weakly repulsive. Our lattice results therefore do not support the possibility of a shallow bound state for the two mesons for the pion mass values we studied. This calculation provides some useful information on the nature of the newly discovered resonance-like structure $Z_c(3900)$ by various experimental groups. △ Less

Submitted 29 April, 2014; v1 submitted 5 March, 2014; originally announced March 2014.

Comments: 14 pages, 6 figures. Version accepted by PRD

Journal ref: Phys. Rev. D 89, 094506 (2014)

Charmonium, $D_s$ and $D_s^*$ from overlap fermion on DWF configurations

Authors: Y. B. Yang, Y. Chen, A. Alexandru, S. J. Dong, T. Draper, M. Gong, F. X. Lee, A. Li, K. F. Liu, Z. Liu, M. Lujan, N. Mathur

Abstract: We take a new approach to determine the scale parameter $r_0$, the physical masses of strange and charm quarks through a global fit which incorporates continuum extrapolation, chiral extrapolation and quark mass interpolation to the lattice data. The charmonium and charm-strange meson spectrum are calculated with overlap valence quarks on $2+1$-flavor domain-wall fermion gauge configurations gener… ▽ More We take a new approach to determine the scale parameter $r_0$, the physical masses of strange and charm quarks through a global fit which incorporates continuum extrapolation, chiral extrapolation and quark mass interpolation to the lattice data. The charmonium and charm-strange meson spectrum are calculated with overlap valence quarks on $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaboration. We use the masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs and obtain $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.110(24)\,{\rm GeV}$, $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.104(9)\,{\rm GeV}$ and $r_0=0.458(11)\,{\rm fm}$. Subsequently, the hyperfine-splitting of charmonium and $f_{D_s}$ are predicted to be $112(5)\,{\rm MeV}$ and $254(5)\,{\rm MeV}$, respectively. △ Less

Submitted 7 January, 2014; originally announced January 2014.

Comments: 7 pages, 6 figures, POS for lattice 2013

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

Authors: Zhaofeng Liu, Ying Chen, Shao-Jing Dong, Michael Glatzmaier, Ming Gong, Anyi Li, Keh-Fei Liu, Yi-Bo Yang, Jian-Bo Zhang

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

Submitted 19 August, 2014; v1 submitted 30 December, 2013; originally announced December 2013.

Comments: 26 pages, 17 figures. More discussions on O(4) breaking effects, and on the perturbative running and a^2p^2 extrapolation of Zs. A subsection for the calculation of the step scaling function of quark mass is added. References added. Version to appear in PRD

Journal ref: Phys. Rev. D 90, 034505 (2014)

A Lattice Study of Quark and Glue Momenta and Angular Momenta in the Nucleon

Authors: M. Deka, T. Doi, Y. B. Yang, B. Chakraborty, S. J. Dong, T. Draper, M. Glatzmaier, M. Gong, H. W. Lin, K. F. Liu, D. Mankame, N. Mathur, T. Streuer

Abstract: We report a complete calculation of the quark and glue momenta and angular momenta in the proton. These include the quark contributions from both the connected and disconnected insertions. The quark disconnected insertion loops are computed with $Z_4$ noise, and the signal-to-noise is improved with unbiased subtractions. The glue operator is comprised of gauge-field tensors constructed from the ov… ▽ More We report a complete calculation of the quark and glue momenta and angular momenta in the proton. These include the quark contributions from both the connected and disconnected insertions. The quark disconnected insertion loops are computed with $Z_4$ noise, and the signal-to-noise is improved with unbiased subtractions. The glue operator is comprised of gauge-field tensors constructed from the overlap operator. The calculation is carried out on a $16^3 \times 24$ quenched lattice at $β= 6.0$ for Wilson fermions with $κ=0.154, 0.155$, and $0.1555$ which correspond to pion masses at $650, 538$, and $478$~MeV, respectively. The chirally extrapolated $u$ and $d$ quark momentum/angular momentum fraction is found to be $0.64(5)/0.70(5)$, the strange momentum/angular momentum fraction is $0.024(6)/0.023(7)$, and that of the glue is $0.33(6)/0.28(8)$. The previous study of quark spin on the same lattice revealed that it carries a fraction of $0.25(12)$ of proton spin. The orbital angular momenta of the quarks are then obtained from subtracting the spin from their corresponding angular momentum components. We find that the quark orbital angular momentum constitutes $0.47(13)$ of the proton spin with almost all of it coming from the disconnected insertions. △ Less

Submitted 12 February, 2015; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: Renormalization section is expanded to include more details. There are slight changes in the final numbers. A few modification and corrections are made in the rest of the text

Report number: RIKEN-QHP-109

Journal ref: Phys.Rev. D91 (2015) 1, 014505

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

Authors: Zhaofeng Liu, Ying Chen, Yi-Bo Yang, Shao-Jing Dong, Michael Glatzmaier, Ming Gong, Keh-Fei Liu, Anyi Li, Jian-Bo Zhang

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

Submitted 2 December, 2013; originally announced December 2013.

Comments: 7 pages, 10 figures, presented at the 31st International Symposium on Lattice Field Theory (LATTICE 2013), 29 July - 3 August 2013, Mainz, Germany

Journal ref: PoS(LATTICE 2013)307

Oscillatory behavior of the domain wall fermions revisited

Authors: Jian Liang, Ying Chen, Ming Gong, Long-Cheng Gui, Keh-Fei Liu, Zhaofeng Liu, Yi-Bo Yang

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

Submitted 8 June, 2014; v1 submitted 13 October, 2013; originally announced October 2013.

Comments: 7 pages, no figures. Minor corrections to match the published version

Journal ref: Physical Review D 89, 094507 (2014)

Strangeness and charmness content of nucleon from overlap fermions on 2+1-flavor domain-wall fermion configurations

Authors: M. Gong, A. Alexandru, Y. Chen, T. Doi, S. J. Dong, T. Draper, W. Freeman, M. Glatzmaier, A. Li, K. F. Liu, Z. Liu

Abstract: We present a calculation of the strangeness and charmness contents <N|\bar{s}s|N> and <N|\bar{c}c|N> of the nucleon from dynamical lattice QCD with 2+1 flavors. The calculation is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations. The configurations are generated by the RBC collaboration on a 24^3*64 lattice with sea quark mass am_l=0.005, am_s=0.04, and… ▽ More We present a calculation of the strangeness and charmness contents <N|\bar{s}s|N> and <N|\bar{c}c|N> of the nucleon from dynamical lattice QCD with 2+1 flavors. The calculation is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations. The configurations are generated by the RBC collaboration on a 24^3*64 lattice with sea quark mass am_l=0.005, am_s=0.04, and inverse lattice spacing a^{-1}=1.73GeV. Both actions have chiral symmetry which is essential in avoiding contamination due to the operator mixing with other flavors. Nucleon propagator and the quark loops are both computed with stochastic grid sources, while low-mode substitution and low-mode averaging methods are used respectively which substantially improve the signal to noise ratio. We obtain the strangeness matrix element f_{T_{s}} = m_s <N|\bar{s}s|N> / M_N = 0.0334(62), and the charmness content f_{T_{c}} = m_c <N|\bar{c}c|N> / M_N = 0.094(31) which is resolved from zero by 3σprecision for the first time. △ Less

Submitted 17 August, 2013; v1 submitted 3 April, 2013; originally announced April 2013.

Comments: 34 pages, 18 figures; correct typos and minor updates; Figures updated, final version published in Phys. Rev. D 88, 014503 (2013)

Journal ref: Phys. Rev. D 88, 014503 (2013)

The $Δ_{mix}$ parameter in the overlap on domain-wall mixed action

Authors: M. Lujan, A. Alexandru, T. Draper, W. Freeman, M. Gong, F. X. Lee, A. Li, K. F. Liu, N. Mathur

Abstract: A direct calculation of the mixed-action parameter $Δ_{mix}$ with valence overlap fermions on a domain-wall fermion sea is presented. The calculation is performed on four ensembles of the 2+1-flavor domain-wall gauge configurations: $24^3 \times 64$ ($a m_l= 0.005$, $a=0.114\fm$) and $32^3 \times 64$ ($a m_l = 0.004, 0.006, 0.008$, $a=0.085\fm$). For pion masses close to $300\MeV$ we find \hbox{… ▽ More A direct calculation of the mixed-action parameter $Δ_{mix}$ with valence overlap fermions on a domain-wall fermion sea is presented. The calculation is performed on four ensembles of the 2+1-flavor domain-wall gauge configurations: $24^3 \times 64$ ($a m_l= 0.005$, $a=0.114\fm$) and $32^3 \times 64$ ($a m_l = 0.004, 0.006, 0.008$, $a=0.085\fm$). For pion masses close to $300\MeV$ we find \hbox{$Δ_{mix}=0.030(6)\GeV^4$} at $a=0.114\fm$ and $Δ_{mix}=0.033(12)\GeV^4$ at $a=0.085\fm$. The results are quite independent of the lattice spacing and they are significantly smaller than the results for valence domain-wall fermions on Asqtad sea or those of valence overlap fermions on clover sea. Combining the results extracted from these two ensembles, we get $Δ_{mix}=0.030(6)(5)\GeV^4$, where the first error is statistical and the second is the systematic error associated with the fitting method. △ Less

Submitted 27 April, 2012; originally announced April 2012.

Comments: 6 pages, 6 figures

The Strangeness and Charmness of Nucleon from Overlap Fermions

Authors: χQCD Collaboration, :, M. Gong, A. Li, A. Alexandru, T. Draper, K. F. Liu

Abstract: The calculation of the strangeness and charmness of the nucleon is presented with overlap fermion action on 2+1 flavor domain wall fermion configurations. We adopt stochastic grid sources and the low mode substitution technique to improve the signals of nucleon correlation functions and the loops. The calculation is done on a $24^3\times 64$ lattice with $m_l=0.005$, $m_h=0.04$, and… ▽ More The calculation of the strangeness and charmness of the nucleon is presented with overlap fermion action on 2+1 flavor domain wall fermion configurations. We adopt stochastic grid sources and the low mode substitution technique to improve the signals of nucleon correlation functions and the loops. The calculation is done on a $24^3\times 64$ lattice with $m_l=0.005$, $m_h=0.04$, and $a^{-1}=1.73\,{\rm GeV}$. We find $ f_{T_{s}} = 0.048(15)$ and $f_{T_{c}} = 0.029(43)$. △ Less

Submitted 3 April, 2012; originally announced April 2012.

Comments: Talk presented at the XXIX International Symposium on Lattice Field Theory, Squaw Valley, Lake Tahoe, California, USA, July 10-16, 2011. 7 pages, 5 figures

Journal ref: PoS(Lattice 2011)156