Showing 1–50 of 129 results for author: Chen, Y

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

Fully heavy tetraquark resonant states with different flavors

Authors: Wei-Lin Wu, Yao Ma, Yan-Ke Chen, Lu Meng, Shi-Lin Zhu

Abstract: We use the quark potential model to calculate the mass spectrum of the S-wave fully heavy tetraquark systems with different flavors, including the $ bc\bar b\bar c, bb\bar c\bar c, cc\bar c\bar b $ and $ bb\bar b\bar c $ systems. We employ the Gaussian expansion method to solve the four-body Schrödinger equation, and the complex scaling method to identify resonant states. The… ▽ More We use the quark potential model to calculate the mass spectrum of the S-wave fully heavy tetraquark systems with different flavors, including the $ bc\bar b\bar c, bb\bar c\bar c, cc\bar c\bar b $ and $ bb\bar b\bar c $ systems. We employ the Gaussian expansion method to solve the four-body Schrödinger equation, and the complex scaling method to identify resonant states. The $ bc\bar b\bar c, bb\bar c\bar c, cc\bar c\bar b $ and $ bb\bar b\bar c $ resonant states are obtained in the mass regions of $ (13.2,13.5) $, $ (13.3,13.6) $, $ (10.0,10.3) $, $ (16.5,16.7) $ GeV, respectively. Among these states, the $ bc\bar b\bar c $ tetraquark states are the most promising ones to be discovered in the near future. We recommend the experimental exploration of the $ 1^{++} $ and $ 2^{++} $ $ bc\bar b\bar c $ states with masses near $ 13.3 $ GeV in the $ J/ψΥ$ channel. From the root-mean-square radii, we find that all the resonant states we have identified are compact tetraquark states. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 10 pages,7 figures,8 tables

Probing the pole origin of $X(3872)$ with the coupled-channel dynamics

Authors: Jun-Zhang Wang, Zi-Yang Lin, Yan-Ke Chen, Lu Meng, Shi-Lin Zhu

Abstract: The $X(3872)$, as the first and the most crucial member in the exotic charmoniumlike $XYZ$ family, has been studied for a long time. However, its dynamical origin, whether stemming from a $D\bar{D}^*$ hadronic molecule or the first excited $P$-wave charmonium $χ_{c1}(2P)$, remains controversial. In this Letter, we demonstrate that the $X(3872)$ definitely does not result from the mass shift of the… ▽ More The $X(3872)$, as the first and the most crucial member in the exotic charmoniumlike $XYZ$ family, has been studied for a long time. However, its dynamical origin, whether stemming from a $D\bar{D}^*$ hadronic molecule or the first excited $P$-wave charmonium $χ_{c1}(2P)$, remains controversial. In this Letter, we demonstrate that the $X(3872)$ definitely does not result from the mass shift of the higher bare $χ_{c1}(2P)$ resonance pole in the coupled-channel dynamics involving a short-distance $c\bar{c}$ core and the long-distance $D\bar{D}^*$ channels. Instead, it originates from either the $D\bar{D}^*$ molecular pole or the shadow pole associated with the $P$-wave charmonium, which depends on the concrete coupling mode between the $c\bar{c}$ and $D\bar{D}^*$. In order to further exploit the nature of $X(3872)$, we carefully investigate potential mechanisms that contribute to its pole width, which suggests that the coupled-channel dynamics plays a critical role in causing a noticeable discrepancy between the pole widths of $X(3872)$ and $T_{cc}^+$. Interestingly, we bridge the quantitative connection among the dynamics origin of $X(3872)$, its pole width and the properties of the predicted new resonance. The precise measurement of the pole width of $X(3872)$ and the search for the new charmoniumlike resonance become highly significant and can be anticipated in future LHCb, BESIII and Belle II experiments. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 11 pages, 5 figures

The radiative decay of scalar glueball from lattice QCD

Authors: Jintao Zou, Long-Cheng Gui, Ying Chen, Jian Liang, Xiangyu Jiang, Wen Qin, Yi-Bo Yang

Abstract: We perform the first lattice QCD study on the radiative decay of the scalar glueball to the vector meson $φ$ in the quenched approximation. The calculations are carried out on three gauge ensembles with different lattice spacings, which enable us to do the continuum extrapolation. We first revisit the radiative $J/ψ$ decay into the scalar glueball $G$ and obtain the partial decay width… ▽ More We perform the first lattice QCD study on the radiative decay of the scalar glueball to the vector meson $φ$ in the quenched approximation. The calculations are carried out on three gauge ensembles with different lattice spacings, which enable us to do the continuum extrapolation. We first revisit the radiative $J/ψ$ decay into the scalar glueball $G$ and obtain the partial decay width $Γ(J/ψ\to γG)=0.578(86)~\text{keV}$ and the branching fraction $\text{Br}(J/ψ\to γG) = 6.2(9)\times 10^{-3}$. We then extend the similar calculation to the process $G\to γφ$ and get the partial decay width $Γ(G \to γφ)= 0.074(47)~\text{keV}$, which implies that the combined branching fraction of $J/ψ\toγG\to γγφ$ is as small as $\mathcal{O}(10^{-9})$ such that this process is hardly detected by the BESIII experiment even with the large $J/ψ$ sample of $\mathcal{O}(10^{10})$. With the vector meson dominance model, the two-photon decay width of the scalar glueball is estimated to be $Γ(G\toγγ)=0.53(46)~\text{eV}$, which results in a large stickiness $S(G)\sim \mathcal{O}(10^4)$ of the scalar glueball by assuming the stickiness of $f_2(1270)$ to be one. △ Less

Submitted 15 July, 2024; v1 submitted 1 April, 2024; originally announced April 2024.

Comments: 13 pages,11 figures. This version includes updated numerical results and corrections to previous calculations

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

Benchmark calculations of fully heavy compact and molecular tetraquark states

Authors: Wei-Lin Wu, Yan-Ke Chen, Lu Meng, Shi-Lin Zhu

Abstract: We calculate the mass spectrum of the S-wave fully heavy tetraquark systems $ QQ\bar Q\bar Q~(Q=c,b) $ with both normal $ (J^{PC}=0^{++},1^{+-},2^{++}) $ and exotic $ (J^{PC}=0^{+-},1^{++},2^{+-}) $ C-parities using three different quark potential models (AL1, AP1, BGS). The exotic C-parity systems refer to the ones that cannot be composed of two S-wave ground heavy quarkonia. We incorporate the m… ▽ More We calculate the mass spectrum of the S-wave fully heavy tetraquark systems $ QQ\bar Q\bar Q~(Q=c,b) $ with both normal $ (J^{PC}=0^{++},1^{+-},2^{++}) $ and exotic $ (J^{PC}=0^{+-},1^{++},2^{+-}) $ C-parities using three different quark potential models (AL1, AP1, BGS). The exotic C-parity systems refer to the ones that cannot be composed of two S-wave ground heavy quarkonia. We incorporate the molecular dimeson and compact diquark-antidiquark spatial correlations simultaneously, thereby discerning the actual configurations of the states. We employ the Gaussian expansion method to solve the four-body Schrödinger equation, and the complex scaling method to identify the resonant states. The mass spectra in three different models qualitatively agree with each other. We obtain several resonant states with $ J^{PC} = 0^{++}, 1^{+-}, 2^{++}, 1^{++} $ in the mass region $(6.92,7.30)\, \mathrm{GeV}$, some of which are good candidates of the experimentally observed $X(6900)$ and $X(7200)$. We also obtain several exotic C-parity zero-width states with $ J^{PC}=0^{+-} $ and $ 2^{+-} $. These zero-width states have no corresponding S-wave diquarkonium threshold and can only decay strongly to final states with P-wave quarkonia. With the notation $T_{4Q,J(C)}(M)$, we deduce from the root mean square radii that the $ X(7200) $ candidates $ T_{4c,0(+)}(7173), T_{4c,2(+)}(7214) $ and the state $ T_{4c,1(-)}(7191) $ look like molecular states although most of the resonant and zero-width states are compact states. △ Less

Submitted 26 March, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 16 pages, 6 figures, 10 tables. Version accepted by PRD

Journal ref: Phys. Rev. D 109, 054034 (2024)

Soft modes in hot QCD matter

Authors: Jens Braun, Yong-rui Chen, Wei-jie Fu, Fei Gao, Chuang Huang, Friederike Ihssen, Jan M. Pawlowski, Fabian Rennecke, Franz R. Sattler, Yang-yang Tan, Rui Wen, Shi Yin

Abstract: The chiral crossover of QCD at finite temperature and vanishing baryon density turns into a second order phase transition if lighter than physical quark masses are considered. If this transition occurs sufficiently close to the physical point, its universal critical behaviour would largely control the physics of the QCD phase transition. We quantify the size of this region in QCD using functional… ▽ More The chiral crossover of QCD at finite temperature and vanishing baryon density turns into a second order phase transition if lighter than physical quark masses are considered. If this transition occurs sufficiently close to the physical point, its universal critical behaviour would largely control the physics of the QCD phase transition. We quantify the size of this region in QCD using functional approaches, both Dyson-Schwinger equations and the functional renormalisation group. The latter allows us to study both critical and non-critical effects on an equal footing, facilitating a precise determination of the scaling regime. We find that the physical point is far away from the critical region. Importantly, we show that the physics of the chiral crossover is dominated by soft modes even far beyond the critical region. While scaling functions determine all thermodynamic properties of the system in the critical region, the order parameter potential is the relevant quantity away from it. We compute this potential in QCD using the functional renormalisation group and Dyson-Schwinger equations and provide a simple parametrisation for phenomenological applications. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 7+8 pages, 5+4 figures

Unified description of the $Qs \bar q \bar q$ molecular bound states, molecular resonances and compact tetraquark states in the quark potential model

Authors: Yan-Ke Chen, Wei-Lin Wu, Lu Meng, Shi-Lin Zhu

Abstract: We calculate the mass spectrum of the $Qs\bar q \bar q$ $(Q=c, b)$ tetraquark states with $J^P=(0,1,2)^+$ using the AL1 quark potential model, which successfully describes the conventional hadron spectrum. We employ the Gaussian expansion method to solve the four-body Schrödinger equation, and use the complex scaling method to identify the resonances. With the notation… ▽ More We calculate the mass spectrum of the $Qs\bar q \bar q$ $(Q=c, b)$ tetraquark states with $J^P=(0,1,2)^+$ using the AL1 quark potential model, which successfully describes the conventional hadron spectrum. We employ the Gaussian expansion method to solve the four-body Schrödinger equation, and use the complex scaling method to identify the resonances. With the notation $T_{Q s, I(J)}^{\text {Theo. }}(M)$, we find several near-threshold bound states and resonances, including $T_{cs,0(0)}^{\mathrm{Theo.}}(2350)$, $T_{cs,0(0)}^{\mathrm{Theo.}}(2906)$, $T_{bs,0(0)}^{\mathrm{Theo.}}(5781)$, $T_{bs,0(1)}^{\mathrm{Theo.}}(5840)$, and $T_{bs,0(0)}^{\mathrm{Theo.}}(6240)$ which are close to the $D\bar{K}$, $D^*\bar{K}^*$, $\bar{B}\bar{K}$, $\bar{B}^*\bar{K}$ and $\bar{B}^*\bar{K}^*$ thresholds, respectively. Furthermore, their spatial structures clearly support their molecular natures. The resonance $T_{cs,0(0)}^{\mathrm{Theo.}}(2906)$ has a mass of $2906$ MeV, a width of $20$ MeV, and quantum numbers $I(J^P)=0(0^+)$, which may serve as a good candidate for the experimental $T_{cs0}(2900)$ state. We strongly urge the experimental search of the predicted states. △ Less

Submitted 14 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, 3 tables. Version accepted by PRD

Journal ref: Phys.Rev.D 109, 014010 (2024)

Tetraquark bound states in constituent quark models: benchmark test calculations

Authors: Lu Meng, Yan-Ke Chen, Yao Ma, Shi-Lin Zhu

Abstract: We investigate the tetraquark bound states that are manifestly exotic using three distinct few-body methods: Gaussian Expansion Method (GEM), Resonating Group Method (RGM), and Diffusion Monte Carlo (DMC). We refer to manifestly exotic states that do not involve a mixture with the conventional mesons through the creation and annihilation of $n\bar{n}$, where $n=u, d$. Our calculations are conducte… ▽ More We investigate the tetraquark bound states that are manifestly exotic using three distinct few-body methods: Gaussian Expansion Method (GEM), Resonating Group Method (RGM), and Diffusion Monte Carlo (DMC). We refer to manifestly exotic states that do not involve a mixture with the conventional mesons through the creation and annihilation of $n\bar{n}$, where $n=u, d$. Our calculations are conducted with two types of quark models: the pure constituent quark model featuring one-gluon-exchange interactions and confinement interactions, and the chiral constituent quark model, supplemented by extra one-boson-exchange interactions. This study represents a comprehensive benchmark test of various few-body methods and quark models. Our findings reveal the superiority of GEM over RGM and DMC methods based on present implements for the tetraquark bound states. Additionally, we observe a tendency for the chiral quark model to overestimate the binding energies. We systematically explore the fully, triply, doubly, and singly heavy tetraquark states with $J^P=0^+,1^+,2^+$, encompassing over 150 states in total. We successfully identify several bound states, including $[cc\bar{n}\bar{n}]_{J^{P}=1^{+}}^{I=0}$, $[bb\bar{n}\bar{n}]_{J^{P}=1^{+}}^{I=0}$, $[bc\bar{n}\bar{n}]_{J^{P}=0^{+},1^{+},2^{+}}^{I=0}$, $[bs\bar{n}\bar{n}]_{J^{P}=0^{+},1^{+}}^{I=0}$, $[cs\bar{n}\bar{n}]_{J^{P}=0^{+}}^{I=0}$, and $[bb\bar{n}\bar{s}]_{J^{P}=1^{+}}$, all found to be bound states below the dimeson thresholds. △ Less

Submitted 25 November, 2023; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 15 pages, 7 figures and 5 Tables. Version accepted by PRD

Doubly heavy tetraquark states in the constituent quark model using diffusion Monte Carlo method

Authors: Yao Ma, Lu Meng, Yan-Ke Chen, Shi-Lin Zhu

Abstract: We use the diffusion Monte Carlo method to calculate the doubly heavy tetraquark $T_{cc}$ system in two kinds of constituent quark models, the pure constituent quark model AL1/AP1 and the chiral constituent quark model. When the discrete configurations are complete and no spatial clustering is preseted, the AL1/AP1 model gives an energy of $T_{cc}$ close to the $DD^*$ threshold, and the chiral con… ▽ More We use the diffusion Monte Carlo method to calculate the doubly heavy tetraquark $T_{cc}$ system in two kinds of constituent quark models, the pure constituent quark model AL1/AP1 and the chiral constituent quark model. When the discrete configurations are complete and no spatial clustering is preseted, the AL1/AP1 model gives an energy of $T_{cc}$ close to the $DD^*$ threshold, and the chiral constituent quark model yields a deeply bound state. We further calculate all doubly heavy tetraquark systems with $J^P=0^+,1^+,2^+$, and provide the binding energies of systems with bound states. The $I(J^P)=0(0^+)$ $bc\bar{n}\bar{n}$, $0(1^+)$ $bb\bar{n}\bar{n}$, $0(1^+)$ $bc\bar{n}\bar{n}$, $\frac{1}{2}(1^+)$ $bb\bar{s}\bar{n}$ systems have bound states in all three models. Since the DMC method has almost no restriction on the spatial part, the resulting bound states have greater binding energies than those obtained in previous works. △ Less

Submitted 2 April, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 9 pages, 1 figure. Comments are welcomed

Journal ref: Phys. Rev. D 109, 074001 (2024)

Radiative transition decay width of $ψ_2(3823)\rightarrowγχ_{c1}$ from lattice QCD

Authors: Ning Li, Yan Gao, Feiyu Chen, Ying Chen, Xiangyu Jiang, Chunjiang Shi, Wei Sun

Abstract: We present an exploratory $N_f=2$ lattice QCD study of $ψ_2(3823)\to γχ_{c1}$ at a pion mass $m_π\approx 350$~MeV. The related two-point and three-piont functions are calculated using the distillation method. The electromagnetic multipole form factor $\hat{V}(0)=2.083(11)$ for $J/ψ\toγη_c$ is consistent with previous lattice results, the form factors $\hat{E}_1(0)$, $\hat{M}_2(0)$ and… ▽ More We present an exploratory $N_f=2$ lattice QCD study of $ψ_2(3823)\to γχ_{c1}$ at a pion mass $m_π\approx 350$~MeV. The related two-point and three-piont functions are calculated using the distillation method. The electromagnetic multipole form factor $\hat{V}(0)=2.083(11)$ for $J/ψ\toγη_c$ is consistent with previous lattice results, the form factors $\hat{E}_1(0)$, $\hat{M}_2(0)$ and $\hat{E}_3(0)$ for $Γ(χ_{c2}\toγJ/ψ)$ have the same hierarchy as that derived from experiments and the predicted decay width $Γ(χ_{c2}\toγJ/ψ)=368(5)~\text{keV}$ is in excellent agreement with the PDG value $374(10)~\text{keV}$ and previous lattice QCD results in the quenched approximation. The same strategy is applied to the study of the process $ψ_2(3823)\to γχ_{c1}$ and the partial decay width is predicted to be $337(27)~\text{keV}$. According to the BESIII constraints on the $ψ_2(3823)$ decay channels and some phenomenological results, we estimate the total width $Γ(ψ_2(3823))=520(100)~\text{keV}$. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Virtual states in the coupled-channel problems with an improved complex scaling method

Authors: Yan-Ke Chen, Lu Meng, Zi-Yang Lin, Shi-Lin Zhu

Abstract: We improve the complex scaling method (CSM) to obtain virtual states, which were previously challenging in the conventional CSM. Our approach solves the Schrödinger equation in the momentum space as an eigenvalue problem by choosing the flexible contours. It proves to be highly effective in identifying the poles across the different Riemann sheets in the multichannel scatterings. It is more straig… ▽ More We improve the complex scaling method (CSM) to obtain virtual states, which were previously challenging in the conventional CSM. Our approach solves the Schrödinger equation in the momentum space as an eigenvalue problem by choosing the flexible contours. It proves to be highly effective in identifying the poles across the different Riemann sheets in the multichannel scatterings. It is more straightforward and efficient than searching for the zeros of the Fredholm determinant of the Lippmann-Schwinger equation using the root-finding algorithms. This advancement significantly extends the capabilities of the CSM in accurately characterizing the resonances and virtual states in quantum systems. △ Less

Submitted 20 February, 2024; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 12 pages, 14 figures, 2 tables. Version accepted by PRD

Journal ref: Phys.Rev.D 109, 034006 (2024)

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

Ground state baryons in the flux-tube three-body confinement model using Diffusion Monte Carlo

Authors: Yao Ma, Lu Meng, Yan-Ke Chen, Shi-Lin Zhu

Abstract: We make a systematical diffusion Monte Carlo (DMC) calculation for all ground state baryons in two confinement scenarios, the pairwise confinement and the three-body flux-tube confinement. With the baryons as an example, we illustrate a feasible procedure to investigate the few-quark states with possible few-body confinement mechanisms, which can be extended to the multiquark states easily. For ea… ▽ More We make a systematical diffusion Monte Carlo (DMC) calculation for all ground state baryons in two confinement scenarios, the pairwise confinement and the three-body flux-tube confinement. With the baryons as an example, we illustrate a feasible procedure to investigate the few-quark states with possible few-body confinement mechanisms, which can be extended to the multiquark states easily. For each baryon, we extract the mass, mean-square radius, charge radius, and the quark distributions. We use the Jackknife resampling method to estimate the statistical uncertainties of masses to be less than 1 MeV. To determine the baryon charge radii, we include the constituent quark size effect, which is fixed by the experimental and lattice QCD results. Our results show that both two-body and three-body confinement mechanisms can give a good description of the experimental data if the parameters are chosen properly. In the flux-tube confinement, introducing different tension parameters for the baryons and mesons are necessary, specifically, $σ_Y= 0.9204 σ_{Q\bar{Q}}$. The lesson from the calculation of the nucleon mass with the DMC method is that the improper pre-assignment of the channels may prevent us from obtaining the real ground state. With this experience, we obtain the real ground state (the $η_c η_c$ threshold with the di-meson configuration) of the $cc\bar{c}\bar{c}$ system with $J^{PC}=0^{++}$ starting from the diquark-antidiquark spin-color channels alone, which is hard to achieve in the variational method and was not obtained in the previous DMC calculations. △ Less

Submitted 20 February, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 24 pages, 18 figures. The Supplement material is attached in the source code of LaTeX. Comments are welcomed

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

Rediscovery of Numerical Lüscher's Formula from the Neural Network

Authors: Yu Lu, Yi-Jia Wang, Ying Chen, Jia-Jun Wu

Abstract: We present that by predicting the spectrum in discrete space from the phase shift in continuous space, the neural network can remarkably reproduce the numerical Lüscher's formula to a high precision. The model-independent property of the Lüscher's formula is naturally realized by the generalizability of the neural network. This exhibits the great potential of the neural network to extract model-in… ▽ More We present that by predicting the spectrum in discrete space from the phase shift in continuous space, the neural network can remarkably reproduce the numerical Lüscher's formula to a high precision. The model-independent property of the Lüscher's formula is naturally realized by the generalizability of the neural network. This exhibits the great potential of the neural network to extract model-independent relation between model-dependent quantities, and this data-driven approach could greatly facilitate the discovery of the physical principles underneath the intricate data. △ Less

Submitted 8 April, 2024; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 7 figures, accepted by Chinese Physics C

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.

Fundamental Physics in Small Experiments

Authors: T. Blum, P. Winter, T. Bhattacharya, T. Y. Chen, V. Cirigliano, D. DeMille, A. Gerarci, N. R. Hutzler, T. M. Ito, O. Kim, R. Lehnert, W. M. Morse, Y. K. Semertzidis

Abstract: High energy physics aims to understand the fundamental laws of particles and their interactions at both the largest and smallest scales of the universe. This typically means probing very high energies or large distances or using high-intensity beams, which often requires large-scale experiments. A complementary approach is offered through high-precision measurements in small- and mid-scale size ex… ▽ More High energy physics aims to understand the fundamental laws of particles and their interactions at both the largest and smallest scales of the universe. This typically means probing very high energies or large distances or using high-intensity beams, which often requires large-scale experiments. A complementary approach is offered through high-precision measurements in small- and mid-scale size experiments, often at lower energies. The field of such high-precision experiments has seen tremendous progress and importance for particle physics for at least two reasons. First, they exploit synergies to adjacent areas of particle physics and benefit by many recent advances in experimental techniques. Together with intensified phenomenological explorations, these advances led to the realization that challenges associated with weak couplings or the expected suppression factors from the mass scale of new physics can be overcome with such methods. Second, many of these measurements add a new set of particle physics phenomena and observables that can be reached compared to the more conventional methodologies using high energies. Combining high-precision, smaller-scale measurements with the large-scale efforts therefore casts a wider and tighter net for possible effects originating from physics beyond the Standard Model. This report presents a broad set of small-scale research projects that could provide key new precision measurements in the areas of electric dipole moments, magnetic dipole moments, fermion flavor violation, tests of spacetime symmetries, and tests with gravity. The growing impact of these high-precision studies in high energy physics and the complementary input they provide compared to large-scale efforts warrants strong support over the next decades. In particular, EDM searches are expected to improve sensitivities by four or more orders of magnitude in the next decade or two. △ Less

Submitted 27 October, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Snowmass 2021 Community Study on the Future of Particle Physics, Rare Processes and Precision Measurements Frontier, Topical Group RF3 Report v2: 3 additional references and one co-author added

$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

Variational Tensor Network Operator

Authors: Yu-Hsueh Chen, Ke Hsu, Wei-Lin Tu, Hyun-Yong Lee, Ying-Jer Kao

Abstract: We propose a simple and generic construction of the variational tensor network operators to study the quantum spin systems by the synergy of ideas from the imaginary-time evolution and variational optimization of trial wave functions. By applying these operators to simple initial states, accurate variational ground state wave functions with extremely few parameters can be obtained. Furthermore, th… ▽ More We propose a simple and generic construction of the variational tensor network operators to study the quantum spin systems by the synergy of ideas from the imaginary-time evolution and variational optimization of trial wave functions. By applying these operators to simple initial states, accurate variational ground state wave functions with extremely few parameters can be obtained. Furthermore, the framework can be applied to study spontaneously symmetry breaking, symmetry protected topological, and intrinsic topologically ordered phases, and we show that symmetries of the local tensors associated with these phases can emerge directly after the optimization without any gauge fixing. This provides a universal way to identify quantum phase transitions without prior knowledge of the system. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 16 pages, 12 figures

Journal ref: Phys. Rev. Research 4, 043153 (2022)

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

Topological susceptibility in finite temperature QCD with physical $(u/d, s, c)$ domain-wall quarks

Authors: Yu-Chih Chen, Ting-Wai Chiu, Tung-Han Hsieh

Abstract: We perform hybrid Monte-Carlo (HMC) simulation of lattice QCD with $N_f=2+1+1$ domain-wall quarks at the physical point, on the $64^3 \times (64,20,16,12,10,8,6)$ lattices, each with three lattice spacings. The lattice spacings and the bare quark masses are determined on the $64^4$ lattices. The resulting gauge ensembles provide a basis for studying finite temperature QCD with $N_f=2+1+1 $ domain-… ▽ More We perform hybrid Monte-Carlo (HMC) simulation of lattice QCD with $N_f=2+1+1$ domain-wall quarks at the physical point, on the $64^3 \times (64,20,16,12,10,8,6)$ lattices, each with three lattice spacings. The lattice spacings and the bare quark masses are determined on the $64^4$ lattices. The resulting gauge ensembles provide a basis for studying finite temperature QCD with $N_f=2+1+1 $ domain-wall quarks at the physical point. In this paper, we determine the topological susceptibility of the QCD vacuum for $T > T_c \sim 150 $ MeV. The topological charge of each gauge configuration is measured by the clover charge in the Wilson flow at the same flow time in physical units, and the topological susceptibility $ χ_t(a,T) $ is determined for each ensemble with lattice spacing $a$ and temperature $T$. Using the topological susceptibility $χ_t(a,T) $ of 15 gauge ensembles with three lattice spacings and different temperatures in the range $T \sim 155-516 $ MeV, we extract the topological susceptibility $χ_t(T)$ in the continuum limit. To compare our results with others, we survey the continuum extrapolated $χ_t(T)$ in lattice QCD with $N_f=2+1(+1)$ dynamical quarks at/near the physical point, and discuss their discrepancies. Moreover, a detailed discussion on the reweighting method for domain-wall fermion is presented. △ Less

Submitted 20 September, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 47 pages, 7 figures, v2: add a section on comparison with other lattice results and an appendix on renormalized chiral condensate, accepted for publication in PRD

Report number: NTUTH-22-505A

Journal ref: Phys. Rev. D 106, 074501 (2022)

Finite temperature QCD with physical $(u/d, s, c)$ domain-wall quarks

Authors: Yu-Chih Chen, Ting-Wai Chiu, Tung-Han Hsieh

Abstract: In order to understand the role of QCD in the early universe, we perform hybrid Monte-Carlo simulation of lattice QCD with $N_f=2+1+1$ optimal domain-wall quarks at the physical point, on the $64^3 \times (6,8,10,12,16,20,64)$ lattices, each with three lattice spacings. The lattice spacings and the bare quark masses are determined on the $64^4$ lattices. The resulting gauge ensembles provide a bas… ▽ More In order to understand the role of QCD in the early universe, we perform hybrid Monte-Carlo simulation of lattice QCD with $N_f=2+1+1$ optimal domain-wall quarks at the physical point, on the $64^3 \times (6,8,10,12,16,20,64)$ lattices, each with three lattice spacings. The lattice spacings and the bare quark masses are determined on the $64^4$ lattices. The resulting gauge ensembles provide a basis for studying finite temperature QCD with $N_f=2+1+1 $ domain-wall quarks at the physical point. In this Proceeding, we present our first result on the topological susceptibility of the QCD vacuum. The topological charge of each gauge configuration is measured by the clover charge in the Wilson flow at the same flow time in physical units, and the topological susceptibility $ χ_t(a,T) $ is determined for each ensemble with lattice spacing $a$ and temperature $T$. Using the topological susceptibility $χ_t(a,T) $ of 15 gauge ensembles with three lattice spacings and different temperatures in the range $T \sim 155-516 $~MeV, we extract the topological susceptibility $χ_t(T)$ in the continuum limit. △ Less

Submitted 5 April, 2022; v1 submitted 4 December, 2021; originally announced December 2021.

Comments: 9 pages, 1 figure, Proceedings of the 38th International Symposium on Lattice Field Theory, LATTICE 2021, 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Report number: NTUTH-21-505A

Journal ref: PoS LATTICE2021 (2022) 574

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.

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

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)

Gluons in charmoniumlike states

Authors: Wei Sun, Ying Chen, Peng Sun, Yi-Bo Yang

Abstract: The mass components of charmoniumlike states are investigated through the decomposition of QCD energy-momentum tensor (EMT) on lattice. The quark mass contribution $\langle H_m\rangle$ and the momentum fraction $\langle x\rangle$ of valence charm quark and antiquark are calculated for conventional $1S,1P,1D$ charmonia and the exotic $1^{-+}$ charmoniumlike state, based on the $N_f=2+1$ gauge confi… ▽ More The mass components of charmoniumlike states are investigated through the decomposition of QCD energy-momentum tensor (EMT) on lattice. The quark mass contribution $\langle H_m\rangle$ and the momentum fraction $\langle x\rangle$ of valence charm quark and antiquark are calculated for conventional $1S,1P,1D$ charmonia and the exotic $1^{-+}$ charmoniumlike state, based on the $N_f=2+1$ gauge configurations generated by the RBC/UKQCD collaboration. It is found that $\langle H_m\rangle$ is close to each other and around 2.0 to 2.2 GeV for these states, which implies that the mass splittings among these states come almost from the gluon contribution of QCD trace anomaly. The $\langle x\rangle$ of the $1^{-+}$ state is only around 0.55, while that in conventional charmonia is around 0.7 to 0.8. This difference manifests that the proportion of light quarks and gluons in the $1^{-+}$ charmoniumlike state is significantly larger than conventional states. △ Less

Submitted 10 May, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 9 pages, 6 figures; v2: match the published version

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

Confinement and the Global $SU(3)$ Color Symmetry

Authors: Ying Chen

Abstract: The global $SU(3)$ color symmetry and its physical consequences are discussed. The Nöther current is actually governed by the conserved matter current of color charges if the color field generated by this charge is properly polarized. The color field strength of a charge can have a uniform part due to the nontrivial QCD vacuum field and the nonzero gluon condensate, which implies that the self-ene… ▽ More The global $SU(3)$ color symmetry and its physical consequences are discussed. The Nöther current is actually governed by the conserved matter current of color charges if the color field generated by this charge is properly polarized. The color field strength of a charge can have a uniform part due to the nontrivial QCD vacuum field and the nonzero gluon condensate, which implies that the self-energy of a system with a net color charge is infinite and thereby cannot exist as a free state. This is precisely what the color confinement means. Accordingly, the Cornell type potential with the feature of the Casimir scaling is derived for a color singlet system composed of a static color charge and an anti-charge. The uniform color field also implies that a hadron has a minimal size and a minimal energy. Furthermore, the global $SU(3)$ color symmetry requires that the minimal irreducible color singlet systems can only be $q\bar{q}$, $qqq$, $gg$, $ggg$, $q\bar{q}g$, $qqqg$ and $\bar{q}\bar{q}\bar{q}g$, etc., as such a multi-quark systems can only exist as a molecular configurations if there are no other binding mechanisms. △ Less

Submitted 16 October, 2020; v1 submitted 29 September, 2020; originally announced October 2020.

Comments: 5 pages, 3 figures. Comments are welcome. Typos corrected. Submitted to Chin. Phys. C

Diquark Correlations in Hadron Physics: Origin, Impact and Evidence

Authors: M. Yu. Barabanov, M. A. Bedolla, W. K. Brooks, G. D. Cates, C. Chen, Y. Chen, E. Cisbani, M. Ding, G. Eichmann, R. Ent, J. Ferretti, R. W. Gothe, T. Horn, S. Liuti, C. Mezrag, A. Pilloni, A. J. R. Puckett, C. D. Roberts, P. Rossi, G. Salme, E. Santopinto, J. Segovia, S. N. Syritsyn, M. Takizawa, E. Tomasi-Gustafsson , et al. (2 additional authors not shown)

Abstract: The last decade has seen a marked shift in how the internal structure of hadrons is understood. Modern experimental facilities, new theoretical techniques for the continuum bound-state problem and progress with lattice-regularised QCD have provided strong indications that soft quark+quark (diquark) correlations play a crucial role in hadron physics. For example, theory indicates that the appearanc… ▽ More The last decade has seen a marked shift in how the internal structure of hadrons is understood. Modern experimental facilities, new theoretical techniques for the continuum bound-state problem and progress with lattice-regularised QCD have provided strong indications that soft quark+quark (diquark) correlations play a crucial role in hadron physics. For example, theory indicates that the appearance of such correlations is a necessary consequence of dynamical chiral symmetry breaking, viz. a corollary of emergent hadronic mass that is responsible for almost all visible mass in the universe; experiment has uncovered signals for such correlations in the flavour-separation of the proton's electromagnetic form factors; and phenomenology suggests that diquark correlations might be critical to the formation of exotic tetra- and penta-quark hadrons. A broad spectrum of such information is evaluated herein, with a view to consolidating the facts and therefrom moving toward a coherent, unified picture of hadron structure and the role that diquark correlations might play. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: 113 pages, 41 figures, 8 tables

Report number: NJU-INP 024/20

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

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

New $N_f=2$ Pseudofermion Action for Monte-Carlo Simulation of Lattice Field Theory with Domain-Wall Fermions

Authors: Yu-Chih Chen, Ting-Wai Chiu

Abstract: We construct a novel $ N_f = 2 $ pseudofermion action for Monte-Carlo simulation of lattice gauge theory with domain-wall fermions (DWF), of which the effective four-dimensional lattice Dirac operator is equal to the overlap-Dirac operator with the argument of the sign function equal to $ H = c γ_5 D_w (1 + d D_w)^{-1} $, where $ c $ and $ d $ are parameters, and $D_w$ is the standard Wilson-Dirac… ▽ More We construct a novel $ N_f = 2 $ pseudofermion action for Monte-Carlo simulation of lattice gauge theory with domain-wall fermions (DWF), of which the effective four-dimensional lattice Dirac operator is equal to the overlap-Dirac operator with the argument of the sign function equal to $ H = c γ_5 D_w (1 + d D_w)^{-1} $, where $ c $ and $ d $ are parameters, and $D_w$ is the standard Wilson-Dirac operator plus a negative parameter $-m_0 \; (0 < m_0 < 2)$. This new action is particularly useful for the challenging simulations of lattice gauge theories with large $N_f = 2n $ DWF, on the large lattices, and in the strong-coupling regime. △ Less

Submitted 6 July, 2019; originally announced July 2019.

Comments: 14 pages, 3 EPS figures

Report number: NTUTH-19-505A

Journal ref: Phys. Rev. D 100, 054513 (2019)

Chromopolarizabilities of bottomonia from the $Υ(2S,3S,4S) \to Υ(1S,2S)ππ$ transitions

Authors: Yun-Hua Chen, Feng-Kun Guo

Abstract: The dipion transitions $Υ(2S,3S,4S) \to Υ(1S,2S)ππ$ are systematically studied by considering the mechanisms of the hadronization of soft gluons, exchanging the bottomoniumlike $Z_b$ states, and the bottom-meson loops. The strong pion-pion final-state interaction, especially including the channel coupling to $K\bar{K}$ in the $S$-wave, is taken into account in a model-independent way using the dis… ▽ More The dipion transitions $Υ(2S,3S,4S) \to Υ(1S,2S)ππ$ are systematically studied by considering the mechanisms of the hadronization of soft gluons, exchanging the bottomoniumlike $Z_b$ states, and the bottom-meson loops. The strong pion-pion final-state interaction, especially including the channel coupling to $K\bar{K}$ in the $S$-wave, is taken into account in a model-independent way using the dispersion theory. Through fitting to the available experimental data, we extract values of the transition chromopolarizabilities $|α_{Υ(mS)Υ(nS)}|$, which measure the chromoelectric couplings of the bottomonia with soft gluons. It is found that the $Z_b$ exchange has a slight impact on the extracted chromopolarizablity values, and the obtained $|α_{Υ(2S)Υ(1S)}|$ considering the $Z_b$ exchange is $(0.29\pm 0.20)~\text{GeV}^{-3}$. Our results could be useful in studying the interactions of bottomonium with light hadrons. △ Less

Submitted 26 September, 2019; v1 submitted 13 June, 2019; originally announced June 2019.

Comments: 14 pages, 3 figures, more discussions added

Journal ref: Phys. Rev. D 100, 054035 (2019)

Study of the pseudoscalar glueball in $J/ψ$ radiative decays

Authors: Long-Cheng Gui, Jia-Mei Dong, Ying Chen, Yi-Bo Yang

Abstract: We aim to explore the production rate of the pseudoscalar glueball in $J/ψ$ radiative decay by lattice QCD in quenched approximation. The calculation is performed on three anisotropic lattices with the spatial lattice spacing ranging from 0.222(2) fm to 0.110(1) fm. As a calibration of some systematical uncertainties, we first extract the $M1$ form factor $\hat{V}(0)$ of the process… ▽ More We aim to explore the production rate of the pseudoscalar glueball in $J/ψ$ radiative decay by lattice QCD in quenched approximation. The calculation is performed on three anisotropic lattices with the spatial lattice spacing ranging from 0.222(2) fm to 0.110(1) fm. As a calibration of some systematical uncertainties, we first extract the $M1$ form factor $\hat{V}(0)$ of the process $J/ψ\toγη_{c}$ and get the result $\hat{V}(0)=1.933(41)$ in the continuum limit, which gives the partial width $Γ(J/ψ\toγη_{c})=2.47(11)$ keV. These results are in agreement with that of previous lattice studies. As for the pseudoscalar glueball $G_{0^{-+}}$, its mass is derived to be $2.395(14)$ GeV, and the form factor $\hat{V}(0)$ of the process $J/ψ\toγG_{0^{-+}}$ is determined to be $\hat{V}(0)=0.0246(43)$ after continuum extrapolation. Finally, the production rate of the pseudoscalar glueball is predicted to be $2.31(90)\times10^{-4}$, which is much smaller than that of conventional light $q\bar{q}$ $η$ states. After the subtraction of the phase space factor, the couplings of $J/ψXγ$ are similar where $X$ stands for $η$ states and the pseudoscalar glueball. Possibly, the $U_{A}(1)$ anomaly plays an important role for the large couplings of gluons to the flavor singlet $η$ states in $J/ψ$ radiative decays. △ Less

Submitted 29 September, 2019; v1 submitted 9 June, 2019; originally announced June 2019.

Comments: 11 pages, 16 figures, 7 tables, revised version. To appear in PRD

Journal ref: Phys. Rev. D 100, 054511 (2019)

Chromopolarizability of charmonium and $ππ$ final state interaction revisited

Authors: Yun-Hua Chen

Abstract: The chromopolarizability of a quarknonium describes the quarknonium's interaction with soft gluonic fields and can be measured in the heavy quarkonium decays. Within the framework of dispersion theory which consider the $ππ$ final state interaction (FSI) model-independently, we analyze the transition $ψ^\prime\to J/ψπ^+π^-$ and obtain the chromopolarizability $α_{ψ^\prime ψ}$ and the parameter… ▽ More The chromopolarizability of a quarknonium describes the quarknonium's interaction with soft gluonic fields and can be measured in the heavy quarkonium decays. Within the framework of dispersion theory which consider the $ππ$ final state interaction (FSI) model-independently, we analyze the transition $ψ^\prime\to J/ψπ^+π^-$ and obtain the chromopolarizability $α_{ψ^\prime ψ}$ and the parameter $κ$. It is found that the $ππ$ FSI plays an important role in extracting the chromopolarizability from the experimental data. The obtained chromopolarizability with the FSI is reduced to about 1/2 of that without the FSI. With the FSI, we determine the chromopolarizability $|α_{ψ^\primeψ}|=(1.44\pm 0.02)$ GeV$^{-3}$ and the parameter $κ=0.139\pm 0.005.$ Our results could be useful in studying the interactions of charmonium with light hadrons. △ Less

Submitted 2 May, 2019; v1 submitted 13 January, 2019; originally announced January 2019.

Comments: 8 pages, 1 figure

Journal ref: Advances in High Energy Physics, Volume 2019, Article ID 7650678

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)

The Belle II Physics Book

Authors: E. Kou, P. Urquijo, W. Altmannshofer, F. Beaujean, G. Bell, M. Beneke, I. I. Bigi, F. Bishara M. Blanke, C. Bobeth, M. Bona, N. Brambilla, V. M. Braun, J. Brod, A. J. Buras, H. Y. Cheng, C. W. Chiang, G. Colangelo, H. Czyz, A. Datta, F. De Fazio, T. Deppisch, M. J. Dolan, S. Fajfer, T. Feldmann, S. Godfrey , et al. (504 additional authors not shown)

Abstract: We present the physics program of the Belle II experiment, located on the intensity frontier SuperKEKB $e^+e^-$ collider. Belle II collected its first collisions in 2018, and is expected to operate for the next decade. It is anticipated to collect 50/ab of collision data over its lifetime. This book is the outcome of a joint effort of Belle II collaborators and theorists through the Belle II theor… ▽ More We present the physics program of the Belle II experiment, located on the intensity frontier SuperKEKB $e^+e^-$ collider. Belle II collected its first collisions in 2018, and is expected to operate for the next decade. It is anticipated to collect 50/ab of collision data over its lifetime. This book is the outcome of a joint effort of Belle II collaborators and theorists through the Belle II theory interface platform (B2TiP), an effort that commenced in 2014. The aim of B2TiP was to elucidate the potential impacts of the Belle II program, which includes a wide scope of physics topics: B physics, charm, tau, quarkonium, electroweak precision measurements and dark sector searches. It is composed of nine working groups (WGs), which are coordinated by teams of theorist and experimentalists conveners: Semileptonic and leptonic B decays, Radiative and Electroweak penguins, phi_1 and phi_2 (time-dependent CP violation) measurements, phi_3 measurements, Charmless hadronic B decay, Charm, Quarkonium(like), tau and low-multiplicity processes, new physics and global fit analyses. This book highlights "golden- and silver-channels", i.e. those that would have the highest potential impact in the field. Theorists scrutinised the role of those measurements and estimated the respective theoretical uncertainties, achievable now as well as prospects for the future. Experimentalists investigated the expected improvements with the large dataset expected from Belle II, taking into account improved performance from the upgraded detector. △ Less

Submitted 2 September, 2019; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 689 pages

Report number: KEK Preprint 2018-27, BELLE2-PUB-PH-2018-001, FERMILAB-PUB-18-398-T, JLAB-THY-18-2780, INT-PUB-18-047, UWThPh 2018-26

Journal ref: Prog Theor Exp Phys (2019)

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)

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)

Mass Preconditioning for the Exact One-Flavor Action in Lattice QCD with Domain-Wall Fermion

Authors: Yu-Chih Chen, Ting-Wai Chiu

Abstract: The mass-preconditioning (MP) technique has become a standard tool to enhance the efficiency of the hybrid Monte-Carlo simulation (HMC) of lattice QCD with dynamical quarks, for 2-flavors QCD with degenerate quark masses, as well as its extension to the case of one-flavor by taking the square-root of the fermion determinant of 2-flavors with degenerate masses. However, for lattice QCD with domain-… ▽ More The mass-preconditioning (MP) technique has become a standard tool to enhance the efficiency of the hybrid Monte-Carlo simulation (HMC) of lattice QCD with dynamical quarks, for 2-flavors QCD with degenerate quark masses, as well as its extension to the case of one-flavor by taking the square-root of the fermion determinant of 2-flavors with degenerate masses. However, for lattice QCD with domain-wall fermion, the fermion determinant of any single fermion flavor can be expressed as a functional integral with an exact pseudofermion action $ φ^\dagger H^{-1} φ$, where $ H^{-1} $ is a positive-definite Hermitian operator without taking square-root, and with the chiral structure \cite{Chen:2014hyy}. Consequently, the mass-preconditioning for the exact one-flavor action (EOFA) does not necessarily follow the conventional (old) MP pattern. In this paper, we present a new mass-preconditioning for the EOFA, which is more efficient than the old MP which we have used in Refs. \cite{Chen:2014hyy,Chen:2014bbc}. We perform numerical tests in lattice QCD with $ N_f = 1 $ and $ N_f = 1+1+1+1 $ optimal domain-wall quarks, with one mass-preconditioner applied to one of the exact one-flavor actions, and we find that the efficiency of the new MP is more than 20\% higher than that of the old MP. △ Less

Submitted 26 October, 2017; originally announced October 2017.

Comments: 17 pages, 5 figures

Report number: NTUTH-17-505C

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)