Showing 1–50 of 100 results for author: Liu, Z

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

Critical properties of metallic and deconfined quantum critical points in Dirac systems

Authors: Zi Hong Liu, Matthias Vojta, Fakher F. Assaad, Lukas Janssen

Abstract: We use large-scale fermion quantum Monte Carlo simulations to study metallic and deconfined quantum phase transitions in a bilayer honeycomb model, focusing on their quantum critical and finite-temperature properties. At weak interaction, a fully symmetric Dirac semimetal state is realized. At intermediate and strong interaction, respectively, two long-range-ordered phases that break different sym… ▽ More We use large-scale fermion quantum Monte Carlo simulations to study metallic and deconfined quantum phase transitions in a bilayer honeycomb model, focusing on their quantum critical and finite-temperature properties. At weak interaction, a fully symmetric Dirac semimetal state is realized. At intermediate and strong interaction, respectively, two long-range-ordered phases that break different symmetries are stabilized. The ordered phases feature partial and full, respectively, gap openings in the fermion spectrum. We clarify the symmetries of the different zero-temperature phases and the symmetry breaking patterns across the two quantum phase transitions between them. The first transition between the disordered and long-range-ordered semimetallic phases has previously been argued to be described by the $(2+1)$-dimensional Gross-Neveu-SO(3) field theory. By performing simulations with an improved symmetric Trotter decomposition, we further substantiate this claim by computing the critical exponents $1/ν$, $η_φ$, and $η_ψ$, which turn out to be consistent with the field-theoretical expectation within numerical and analytical uncertainties. The second transition between the two long-range-ordered phases has previously been proposed as a possible instance of a metallic deconfined quantum critical point. We further develop this scenario by analyzing the spectral functions in the single-particle, particle-hole, and particle-particle channels. Our results indicate gapless excitations with a unique velocity, supporting the emergence of Lorentz symmetry at criticality. We also determine the finite-temperature phase boundaries above the fully gapped state at large interaction, which smoothly vanish near the putative metallic deconfined quantum critical point, consistent with a continuous or weakly first-order transition. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 13 pages, 9 figures

Three ways to decipher the nature of exotic hadrons: multiplets, three-body hadronic molecules, and correlation functions

Authors: Ming-Zhu Liu, Ya-Wen Pan, Zhi-Wei Liu, Tian-Wei Wu, Jun-Xu Lu, Li-Sheng Geng

Abstract: In the past two decades, a plethora of hadronic states beyond the conventional quark model of $q\bar{q}$ mesons and $qqq$ baryons have been observed experimentally, which motivated extensive studies to understand their nature and the non-perturbative strong interaction. Since most of these exotic states are near the mass thresholds of a pair of conventional hadrons, the prevailing picture is that… ▽ More In the past two decades, a plethora of hadronic states beyond the conventional quark model of $q\bar{q}$ mesons and $qqq$ baryons have been observed experimentally, which motivated extensive studies to understand their nature and the non-perturbative strong interaction. Since most of these exotic states are near the mass thresholds of a pair of conventional hadrons, the prevailing picture is that they are primarily hadronic molecules. In principle, one can verify the molecular nature of these states by thoroughly comparing their masses, decay widths, and production rates in a particular picture with experimental data. However, this is difficult or impossible. First, quantum mechanics allows for the mixing of configurations allowed by symmetries and quantum numbers. Second, data are relatively scarce because of their small production rates and the many difficulties in the experimental measurements. As a result, other alternatives need to be explored. This review summarizes three such approaches that can help disentangle the nature of the many exotic hadrons discovered. In the first approach, based on the molecular interpretations for some exotic states, we study the likely existence of multiplets of hadronic molecules related by various symmetries, such as isospin symmetry, SU(3)-flavor symmetry, heavy quark spin/flavor symmetry, and heavy antiquark diquark symmetry. In the second approach, starting from some hadronic molecular candidates, one can derive the underlying hadron-hadron interactions. With these interactions, one can study related three-body systems and check whether three-body bound states/resonances exist. In the third approach, one can turn to the femtoscopy technique to derive the hadron-hadron interactions, hence inaccessible. This technique provided an unprecedented opportunity to understand the interactions between unstable hadrons. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 140 pages, 40 figures; any comments are appreciated!

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

Lattice QCD calculation of the $D_s^{*}$ radiative decay with (2+1)-flavor Wilson-clover ensembles

Authors: Yu Meng, Jin-Long Dang, Chuan Liu, Zhaofeng Liu, Tinghong Shen, Haobo Yan, Ke-Long Zhang

Abstract: We perform a lattice calculation on the radiative decay of $D_s^*$ using the (2+1)-flavor Wilson-clover gauge ensembles generated by CLQCD collaboration. A method allowing us to calculate the form factor with zero transfer momentum is proposed and applied to the radiative transition $D_s^*\rightarrow D_sγ$ and the Dalitz decay $D_s^*\rightarrow D_s e^+e^-$. After a continuum extrapolation using th… ▽ More We perform a lattice calculation on the radiative decay of $D_s^*$ using the (2+1)-flavor Wilson-clover gauge ensembles generated by CLQCD collaboration. A method allowing us to calculate the form factor with zero transfer momentum is proposed and applied to the radiative transition $D_s^*\rightarrow D_sγ$ and the Dalitz decay $D_s^*\rightarrow D_s e^+e^-$. After a continuum extrapolation using three lattice spacings, we obtain $Γ(D_s^*\rightarrow D_s γ)=0.0549(54)$ keV, where the error is purely statistical. The result is consistent with previous lattice calculations but with a error reduced to only a fifth of the before. The Dalitz decay rate is also calculated for the first time and the ratio with the radiative transition is found to be $R_{ee}=0.624(3)\%$. A total decay width of $D_s^*$ can then be determined as 0.0587(54) keV taking into account the experimental branching fraction. Combining with the most recent experimental measurement on the branching fraction of the purely leptonic decay $D_s^{+,*}\rightarrow e^+ν_e$, we obtain the quantity $f_{D_s^*}|V_{cs}|=(190.5^{+55.1}_{-41.7_{\textrm{stat.}}}\pm 12.6_{\textrm{syst.}})$ MeV, where the stat. is only the statistical error from the experiment, and syst. results from the experimental systematic uncertainty and the lattice statistical error. Our result leads to an improved systematic uncertainty compared to $42.7_{\textrm{syst.}}$ obtained using previous lattice prediction of total decay width $0.070(28)$ keV as the input. △ Less

Submitted 29 April, 2024; v1 submitted 24 January, 2024; originally announced January 2024.

Comments: 8 pages, 6 figures, published version

Journal ref: Physical Review D 109,074511(2024)

Understanding the nature of baryon resonances

Authors: Derek B. Leinweber, Curtis D. Abell, Liam C. Hockley, Waseem Kamleh, Zhan-Wei Liu, Finn M. Stokes, Anthony W. Thomas, Jia-Jun Wu

Abstract: This presentation opens with a brief review of lattice QCD calculations showing the $2s$ radial excitation of the nucleon sits at approximately 2 GeV, well above the Roper resonance position. We then proceed to reconcile this observation with experimental scattering data. While the idea of dressing quark-model states in a coupled-channel analysis to describe scattering data has been around for dec… ▽ More This presentation opens with a brief review of lattice QCD calculations showing the $2s$ radial excitation of the nucleon sits at approximately 2 GeV, well above the Roper resonance position. We then proceed to reconcile this observation with experimental scattering data. While the idea of dressing quark-model states in a coupled-channel analysis to describe scattering data has been around for decades, it's now possible to bring these descriptions to the finite-volume of lattice QCD for confrontation with lattice-QCD calculations. This combination of lattice QCD and experiment demands that we reconsider our preconceived notions about the quark-model and its excitation spectrum. We close with a discussion of an unanticipated resolution to the missing baryon resonances problem. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 8 pages, 2 figures. To appear as a plenary contribution to the proceedings of HADRON 2023, the 20th International Conference on Hadron Spectroscopy and Structure, 5th-9th June 2023, Genova, Italy

Report number: ADP-24-01/T1240

Structure of the $\mathbf{Λ(1670)}$ resonance

Authors: Jiong-Jiong Liu, Zhan-Wei Liu, Kan Chen, Dan Guo, Derek B. Leinweber, Xiang Liu, Anthony W. Thomas

Abstract: We examine the internal structure of the $Λ(1670)$ through an analysis of lattice QCD simulations and experimental data within Hamiltonian effective field theory. Two scenarios are presented. The first describes the $Λ(1670)$ as a bare three-quark basis state, which mixes with the $πΣ$, $\bar{K}N$, $ηΛ$ and $KΞ$ meson-baryon channels. In the second scenario, the $Λ(1670)$ is dynamically generated… ▽ More We examine the internal structure of the $Λ(1670)$ through an analysis of lattice QCD simulations and experimental data within Hamiltonian effective field theory. Two scenarios are presented. The first describes the $Λ(1670)$ as a bare three-quark basis state, which mixes with the $πΣ$, $\bar{K}N$, $ηΛ$ and $KΞ$ meson-baryon channels. In the second scenario, the $Λ(1670)$ is dynamically generated from these isospin-0 coupled channels. The $K^-p$ scattering data and the pole structures of the $Λ(1405)$ and the $Λ(1670)$ can be simultaneously described well in both scenarios. However, a comparison of the finite-volume spectra to lattice QCD calculations reveals significant differences between these scenarios, with a clear preference for the first case. Thus the lattice QCD results play a crucial role in allowing us to distinguish between these two scenarios for the internal structure of the $Λ(1670)$. △ Less

Submitted 7 March, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 13 pages, 5 figures, 1 table; accepted by Phys. Rev. D

Report number: ADP-23-29/T1238

Journal ref: Phys. Rev. D 109, 054025 (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

Low-lying odd-parity nucleon resonances as quark-model like states

Authors: Curtis D. Abell, Derek B. Leinweber, Zhan-Wei Liu, Anthony W. Thomas, Jia-Jun Wu

Abstract: Recent lattice QCD results for the low-lying odd-parity excitations of the nucleon near the $N^{*}(1535)$ and $N^{*}(1650)$ resonance positions have revealed that the lattice QCD states have magnetic moments consistent with predictions from a constituent-quark-model. Using Hamiltonian Effective Field Theory (HEFT) to describe pion-nucleon scattering in the… ▽ More Recent lattice QCD results for the low-lying odd-parity excitations of the nucleon near the $N^{*}(1535)$ and $N^{*}(1650)$ resonance positions have revealed that the lattice QCD states have magnetic moments consistent with predictions from a constituent-quark-model. Using Hamiltonian Effective Field Theory (HEFT) to describe pion-nucleon scattering in the $I(J^{P}) = \frac{1}{2}(\frac{1}{2}^{-})$ channel, we represent these two quark-model like states as two single-particle bare basis states, dressed and mixed by meson-baryon scattering channels. By constraining the free parameters of the Hamiltonian with $S_{11}$ pion-nucleon scattering data, we perform the first calculation of the finite-volume spectrum using two bare-baryon basis states. By comparing this spectrum to contemporary lattice QCD results at three lattice volumes, we analyse the eigenvectors of the Hamiltonian to gain insight into the structure and composition of these two low-lying resonances. We find that an interpretation of the two low-lying nucleon resonances as quark-model like states dressed by meson-baryon interactions is consistent with both the $S_{11}$ scattering data and lattice QCD. We introduce a novel HEFT formalism for estimating scattering-state contaminations in lattice QCD correlation functions constructed with standard three-quark operators. Not only are historical lattice QCD results described with excellent accuracy, but correlation functions with large scattering-state contaminations are identified. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 17 pages, 12 figures

Report number: ADP-23-17-T1226

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

The Future of US Particle Physics -- The Snowmass 2021 Energy Frontier Report

Authors: Meenakshi Narain, Laura Reina, Alessandro Tricoli, Michael Begel, Alberto Belloni, Tulika Bose, Antonio Boveia, Sally Dawson, Caterina Doglioni, Ayres Freitas, James Hirschauer, Stefan Hoeche, Yen-Jie Lee, Huey-Wen Lin, Elliot Lipeles, Zhen Liu, Patrick Meade, Swagato Mukherjee, Pavel Nadolsky, Isobel Ojalvo, Simone Pagan Griso, Christophe Royon, Michael Schmitt, Reinhard Schwienhorst, Nausheen Shah , et al. (10 additional authors not shown)

Abstract: This report, as part of the 2021 Snowmass Process, summarizes the current status of collider physics at the Energy Frontier, the broad and exciting future prospects identified for the Energy Frontier, the challenges and needs of future experiments, and indicates high priority research areas. This report, as part of the 2021 Snowmass Process, summarizes the current status of collider physics at the Energy Frontier, the broad and exciting future prospects identified for the Energy Frontier, the challenges and needs of future experiments, and indicates high priority research areas. △ Less

Submitted 3 January, 2023; v1 submitted 20 November, 2022; originally announced November 2022.

Comments: 128 pages, 41 figures, 17 tables, contribution to Snowmass 2021

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.

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

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

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

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

Comments: v2: Published in prd

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

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

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

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

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

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

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

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

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

Comments: 6 pages, 4 figures, 2 tables

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

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

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

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

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

Comments: 13 pages, 12 figures

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

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

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

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

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

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

Strangeness $S = -2$ baryon-baryon interactions and femtoscopic correlation functions in covariant chiral effective field theory

Authors: Zhi-Wei Liu, Kai-Wen Li, Li-Sheng Geng

Abstract: We study the baryon-baryon interactions with strangeness $S = -2$ and corresponding momentum correlation functions in leading order covariant chiral effective field theory. The relevant low energy constants are determined by fitting to the latest HAL QCD simulations, taking into account all the coupled channels. Extrapolating the so-obtained strong interactions to the physical point and considerin… ▽ More We study the baryon-baryon interactions with strangeness $S = -2$ and corresponding momentum correlation functions in leading order covariant chiral effective field theory. The relevant low energy constants are determined by fitting to the latest HAL QCD simulations, taking into account all the coupled channels. Extrapolating the so-obtained strong interactions to the physical point and considering both quantum statistical effects and the Coulomb interaction, we calculate the $ΛΛ$ and $Ξ^-p$ correlation functions with a spherical Gaussian source and compare them with the recent experimental data. We find remarkable agreement between our predictions and the experimental measurements by using the source radius determined in proton-proton correlations, which demonstrates the consistency between theory, experiment, and lattice QCD simulations. Moreover, we predict the $Σ^+Σ^+$, $Σ^+Λ$, and $Σ^+Σ^-$ interactions and corresponding momentum correlation functions. We further investigate the influence of the source shape and size of the hadron pair on the correlation functions studied and show that the current data are not very sensitive to the source shape. Future experimental measurement of the predicted momentum correlation functions will provide a non-trivial test of not only SU(3) flavor symmetry and its breaking but also the baryon-baryon interactions derived in covariant chiral effective field theory. △ Less

Submitted 4 October, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: 21 pages, 17 figures; version accepted for publication in Chinese Physics C

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

Metallic and Deconfined Quantum Criticality in Dirac Systems

Authors: Zi Hong Liu, Matthias Vojta, Fakher F. Assaad, Lukas Janssen

Abstract: Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)$\times$U(1) symmetry. Using large-scale auxiliary-field quantum Monte Carlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuo… ▽ More Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)$\times$U(1) symmetry. Using large-scale auxiliary-field quantum Monte Carlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuous transition from the weakly-interacting semimetal to a different semimetallic phase in which the SO(3) symmetry is spontaneously broken and where two out of three Dirac cones acquire a mass gap. The associated quantum critical point can be understood in terms of a Gross-Neveu-SO(3) theory. Second, we subsequently observe a transition towards an insulating phase in which the SO(3) symmetry is restored and the U(1) symmetry is spontaneously broken. While strongly first order at the mean-field level, the QMC data is consistent with a direct and continuous transition. It is thus a candidate for a new type of deconfined quantum critical point that features gapless fermionic degrees of freedom. △ Less

Submitted 24 February, 2022; v1 submitted 13 August, 2021; originally announced August 2021.

Comments: 12 pages, 11 figures

Journal ref: Phys. Rev. Lett. 128, 087201 (2022)

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

Test of the hyperon-nucleon interaction within leading order covariant chiral effective field theory

Authors: Jing Song, Zhi-Wei Liu, Kai-Wen Li, Li-Sheng Geng

Abstract: Motivated by the recent experimental measurements of differential cross sections of the $Σ^{-}p$ elastic scattering in the momentum range of $470$ to $850$ MeV$/c$ by the J-PARC E$40$ experiment, we extend our previous studies of $S=-1$ hyperon-nucleon interactions to relatively higher energies up to $900$ MeV$/c$ for both the coupled-channel $Λp\rightarrow(Λp, Σ^{+}n, Σ^{0}p)$,… ▽ More Motivated by the recent experimental measurements of differential cross sections of the $Σ^{-}p$ elastic scattering in the momentum range of $470$ to $850$ MeV$/c$ by the J-PARC E$40$ experiment, we extend our previous studies of $S=-1$ hyperon-nucleon interactions to relatively higher energies up to $900$ MeV$/c$ for both the coupled-channel $Λp\rightarrow(Λp, Σ^{+}n, Σ^{0}p)$, $Σ^{-}p\rightarrow(Λn, Σ^{0}n, Σ^{-}p)$ and single-channel $Σ^{+}p\rightarrowΣ^{+}p$ reactions. We show that although the leading order covariant chiral effective field theory is only constrained by the low energy data, it can describe the high energy data reasonably well, in particular, the J-PARC E40 differential cross sections. The predicted cusp structure close to the $ΣN$ threshold in the $Λp\to Λp$ reaction agrees with the latest ALICE observation as well as with the results of the next-to-leading order heavy baryon chiral effective theory. On the other hand, the comparison with the latest CLAS data on the $Λp$ cross sections between 0.9 and 2.0 GeV$/c$ clearly indicates the need of higher order chiral potentials for such high momenta. This is also the case for the latest J-PARC data on the $Σp \rightarrow Λn$ differential cross sections. Nevertheless, even for these cases, the predictions are in qualitative agreement with the data, albeit with large uncertainties, implying that the predicted total and differential cross sections are of relevance for ongoing and planned experiments. △ Less

Submitted 20 March, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures, comparison with the latest J-PARC and CLAS data made, uncertainty estimates updated, to appear in Physical Review C

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)

$^3S_1-{}^3D_1$ coupled channel $Λ_c N$ interactions: chiral effective field theory vs. lattice QCD

Authors: Jing Song, Yang Xiao, Zhi-Wei Liu, Kai-Wen Li, Li-Sheng Geng

Abstract: We study the lattice QCD $Λ_c N$ phase shifts for the $^3S_1-{}^3D_1$ coupled channel using both the leading order covariant chiral effective theory and the next-to-leading order non-relativistic chiral effective field theory. We show that although it is possible to describe simultaneously the $^3S_1$ and $^3D_1$ phase shifts and the inelasticity $η_1$, the fitted energy range is pretty small, onl… ▽ More We study the lattice QCD $Λ_c N$ phase shifts for the $^3S_1-{}^3D_1$ coupled channel using both the leading order covariant chiral effective theory and the next-to-leading order non-relativistic chiral effective field theory. We show that although it is possible to describe simultaneously the $^3S_1$ and $^3D_1$ phase shifts and the inelasticity $η_1$, the fitted energy range is pretty small, only up to $E_\mathrm{c.m.}=5$ MeV. This raises concerns regarding the consistency between leading/next-to-leading order chiral effective field theory and the lattice QCD simulations. △ Less

Submitted 24 October, 2022; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: 12 pages, 4 figures, to appear in Communications in Theoretical Physics

The ALF (Algorithms for Lattice Fermions) project release 2.4. Documentation for the auxiliary-field quantum Monte Carlo code

Authors: ALF Collaboration, F. F. Assaad, M. Bercx, F. Goth, A. Götz, J. S. Hofmann, E. Huffman, Z. Liu, F. Parisen Toldin, J. S. E. Portela, J. Schwab

Abstract: The Algorithms for Lattice Fermions package provides a general code for the finite-temperature and projective auxiliary-field quantum Monte Carlo algorithm. The code is engineered to be able to simulate any model that can be written in terms of sums of single-body operators, of squares of single-body operators and single-body operators coupled to a bosonic field with given dynamics. The package in… ▽ More The Algorithms for Lattice Fermions package provides a general code for the finite-temperature and projective auxiliary-field quantum Monte Carlo algorithm. The code is engineered to be able to simulate any model that can be written in terms of sums of single-body operators, of squares of single-body operators and single-body operators coupled to a bosonic field with given dynamics. The package includes five pre-defined model classes: SU(N) Kondo, SU(N) Hubbard, SU(N) t-V and SU(N) models with long range Coulomb repulsion on honeycomb, square and N-leg lattices, as well as $Z_2$ unconstrained lattice gauge theories coupled to fermionic and $Z_2$ matter. An implementation of the stochastic Maximum Entropy method is also provided. One can download the code from our Git instance at https://git.physik.uni-wuerzburg.de/ALF/ALF/-/tree/ALF-2.4 and sign in to file issues. △ Less

Submitted 28 December, 2022; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 126 pages, 12 figures. v7: update do ALF 2.4. Submission to SciPost. This article supersedes arXiv:1704.00131

Journal ref: SciPost Physics Codebases 1 (2022)

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

Trace anomaly and dynamical quark mass

Authors: Yi-Bo Yang, Jian Liang, Zhaofeng Liu, Peng Sun

Abstract: We investigated the origin of the RI'/MOM quark mass under the Landau gauge at the non-perturbative scale, using the chiral fermion with different quark masses and lattice spacings. Our result confirms that such a mass is non-vanishing based on the linear extrapolation to the chiral and continuum limit, and shows that such a mass comes from the spontaneous chiral symmetry breaking induced by the n… ▽ More We investigated the origin of the RI'/MOM quark mass under the Landau gauge at the non-perturbative scale, using the chiral fermion with different quark masses and lattice spacings. Our result confirms that such a mass is non-vanishing based on the linear extrapolation to the chiral and continuum limit, and shows that such a mass comes from the spontaneous chiral symmetry breaking induced by the near zero modes with the eigenvalue $λ<{\cal O}(5m_q)$, and is proportional to the quark matrix element of the trace anomaly at least down to $\sim $1.3 GeV. △ Less

Submitted 28 March, 2020; originally announced March 2020.

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

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

The nucleon isovector tensor charge from lattice QCD using chiral fermions

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

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

Submitted 16 February, 2020; originally announced February 2020.

Comments: 8 pages, 5 Figures, 1 Table

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

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

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)

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)

Magnetic moments of the spin-${3\over 2}$ singly heavy baryons

Authors: Lu Meng, Guang-Juan Wang, Chang-Zhi Leng, Zhan-Wei Liu, Shi-Lin Zhu

Abstract: We calculate the magnetic moments of spin-$3\over 2$ singly charmed baryons in the heavy baryon chiral perturbation theory (HBChPT). The analytical expressions are given up to $\mathcal{O}(p^3)$. The heavy quark symmetry is used to reduce the number of low energy constants (LECs). With the lattice QCD simulation data as the magnetic moments of the charmed baryons, the numerical results are given u… ▽ More We calculate the magnetic moments of spin-$3\over 2$ singly charmed baryons in the heavy baryon chiral perturbation theory (HBChPT). The analytical expressions are given up to $\mathcal{O}(p^3)$. The heavy quark symmetry is used to reduce the number of low energy constants (LECs). With the lattice QCD simulation data as the magnetic moments of the charmed baryons, the numerical results are given up to $\mathcal{O}(p^3)$ in three scenarios. In the first scenario, we use the results in the quark model as the leading order input. In the second scenario, we use the heavy quark symmetry and neglect the contribution of heavy quark. In the third scenario, the heavy quark contribution is considered on the basis of the scenario II and the magnetic moments of singly bottom baryons are given as a by-product. △ Less

Submitted 15 November, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 18 pages, 2 figures

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

Nucleon Excited States from Lattice QCD and Hamiltonian Effective Field Theory

Authors: Jia-jun Wu, Jonathan M. M. Hall, H. Kamano, Waseem Kamleh, T. -S. H. Lee, Derek B. Leinweber, Zhan-Wei Liu, Finn M. Stokes, Anthony W. Thomas

Abstract: An approach for relating the nucleon excited states extracted from lattice QCD and the nucleon resonances of experimental data has been developed using the Hamiltonian effective field theory (HEFT) method. By formulating HEFT in the finite volume of the lattice, the eigenstates of the Hamiltonian model can be related to the energy eigenstates observed in Lattice simulations. By taking the infinite… ▽ More An approach for relating the nucleon excited states extracted from lattice QCD and the nucleon resonances of experimental data has been developed using the Hamiltonian effective field theory (HEFT) method. By formulating HEFT in the finite volume of the lattice, the eigenstates of the Hamiltonian model can be related to the energy eigenstates observed in Lattice simulations. By taking the infinite-volume limit of HEFT, information from the lattice is linked to experiment. The approach opens a new window for the study of experimentally-observed resonances from the first principles of lattice QCD calculations. With the Hamiltonian approach, one not only describes the spectra of lattice-QCD eigenstates through the eigenvalues of the finite-volume Hamiltonian matrix, but one also learns the composition of the lattice-QCD eigenstates via the eigenvectors of the Hamiltonian matrix. One learns the composition of the states in terms of the meson-baryon basis states considered in formulating the effective field theory. One also learns the composition of the resonances observed in Nature. In this paper, we will focus on recent breakthroughs in our understanding of the structure of the $N^*(1535)$, $N^*(1440)$ and $Λ^*(1405)$ resonances using this method. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 20 pages, 9 figures, Proceedings of The 11th International Workshop on the Physics of Excited Nucleons, August 20 - 23, 2017, at the University of South Carolina, Columbia, SC

Magnetic moments of the spin-$1\over 2$ singly charmed baryons in chiral perturbation theory

Authors: Guang-Juan Wang, Lu Meng, Hao-Song Li, Zhan-Wei Liu, Shi-Lin Zhu

Abstract: We systematically derive the analytical expressions of the magnetic moments of the spin-$1\over 2$ singly charmed baryons to the next-to-next-to-leading order in the heavy baryon chiral perturbation theory (HBChPT). We discuss the analytical relations between the magnetic moments. We estimate the-low energy constants (LECs) in two scenarios. In the first scenario, we use the quark model and Lattic… ▽ More We systematically derive the analytical expressions of the magnetic moments of the spin-$1\over 2$ singly charmed baryons to the next-to-next-to-leading order in the heavy baryon chiral perturbation theory (HBChPT). We discuss the analytical relations between the magnetic moments. We estimate the-low energy constants (LECs) in two scenarios. In the first scenario, we use the quark model and Lattice QCD simulation results as input. In the second scenario, the heavy quark symmetry is adopted to reduce the number of the independent LECs, which are then fitted using the data from the Lattice QCD simulations. We give the numerical results to the next-to-leading order for the antitriplet charmed baryons and to the next-to-next-to-leading order for the sextet states. △ Less

Submitted 25 September, 2018; v1 submitted 1 March, 2018; originally announced March 2018.

Comments: Published version in Phys. Rev. D

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

EMUS-QMC: Elective Momentum Ultra-Size Quantum Monte Carlo Method

Authors: Zi Hong Liu, Xiao Yan Xu, Yang Qi, Kai Sun, Zi Yang Meng

Abstract: One bottleneck of quantum Monte Carlo (QMC) simulation of strongly correlated electron systems lies at the scaling relation of computational complexity with respect to the system sizes. For generic lattice models of interacting fermions, the best methodology at hand still scales with $βN^3$ where $β$ is the inverse temperature and $N$ is the system size. Such scaling behavior has greatly hampered… ▽ More One bottleneck of quantum Monte Carlo (QMC) simulation of strongly correlated electron systems lies at the scaling relation of computational complexity with respect to the system sizes. For generic lattice models of interacting fermions, the best methodology at hand still scales with $βN^3$ where $β$ is the inverse temperature and $N$ is the system size. Such scaling behavior has greatly hampered the accessibility of the universal infrared (IR) physics of many interesting correlated electron models at (2+1)D, let alone (3+1)D. To reduce the computational complexity, we develop a new QMC method with inhomogeneous momentum-space mesh, dubbed elective momentum ultra-size quantum Monte Carlo (EQMC) method. Instead of treating all fermionic excitations on an equal footing as in conventional QMC methods, by converting the fermion determinant into the momentum space, our method focuses on fermion modes that are directly associated with low-energy (IR) physics in the vicinity of the so-called hot-spots, while other fermion modes irrelevant for universal properties are ignored. As shown in the manuscript, for any cutoff-independent quantities, e.g. scaling exponents, this method can achieve the same level of accuracy with orders of magnitude increase in computational efficiency. We demonstrate this method with a model of antiferromagnetic itinerant quantum critical point, realized via coupling itinerant fermions with a frustrated transverse-field Ising model on a triangle lattice. The system size of $48 \times 48 \times 32$ ($L\times L\timesβ$, almost 3 times of previous investigations) are comfortably accessed with EQMC. With much larger system sizes, the scaling exponents are unveiled with unprecedentedly high accuracy, and this result sheds new light on the open debate about the nature and the universality class of itinerant quantum critical points. △ Less

Submitted 2 December, 2018; v1 submitted 30 December, 2017; originally announced January 2018.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. B 99, 085114 (2019)

Structure of the Nucleon and its Excitations

Authors: Waseem Kamleh, Derek Leinweber, Zhan-wei Liu, Finn Stokes, Anthony Thomas, Samuel Thomas, Jia-jun Wu

Abstract: The structure of the ground state nucleon and its finite-volume excitations are examined from three different perspectives. Using new techniques to extract the relativistic components of the nucleon wave function, the node structure of both the upper and lower components of the nucleon wave function are illustrated. A non-trivial role for gluonic components is manifest. In the second approach, the… ▽ More The structure of the ground state nucleon and its finite-volume excitations are examined from three different perspectives. Using new techniques to extract the relativistic components of the nucleon wave function, the node structure of both the upper and lower components of the nucleon wave function are illustrated. A non-trivial role for gluonic components is manifest. In the second approach, the parity-expanded variational analysis (PEVA) technique is utilised to isolate states at finite momenta, enabling a novel examination of the electric and magnetic form factors of nucleon excitations. Here the magnetic form factors of low-lying odd-parity nucleons are particularly interesting. Finally, the structure of the nucleon spectrum is examined in a Hamiltonian effective field theory analysis incorporating recent lattice-QCD determinations of low-lying two-particle scattering-state energies in the finite volume. The Roper resonance of Nature is observed to originate from multi-particle coupled-channel interactions while the first radial excitation of the nucleon sits much higher at approximately 1.9 GeV. △ Less

Submitted 15 November, 2017; originally announced November 2017.

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

Report number: ADP-17-38/T1044

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

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

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

Submitted 2 November, 2017; originally announced November 2017.

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

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

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

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

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

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

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

Report number: MSUHEP-17-019

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

Magnetic moments of the spin-${3\over 2}$ doubly heavy baryons

Authors: Lu Meng, Hao-Song Li, Zhan-Wei Liu, Shi-Lin Zhu

Abstract: In this work, we investigate the chiral corrections to the magnetic moments of the spin-$3\over 2$ doubly charmed baryons systematically up to next-to-next-to-leading order with the heavy baryon chiral perturbation theory. The numerical results are given up to next-to-leading order: $μ_{Ξ^{*++}_{cc}}=1.72μ_{N}$, $μ_{Ξ^{*+}_{cc}}=-0.09μ_{N}$, $μ_{Ω^{*+}_{cc}}=0.99μ_{N}$. As a by-product, we have al… ▽ More In this work, we investigate the chiral corrections to the magnetic moments of the spin-$3\over 2$ doubly charmed baryons systematically up to next-to-next-to-leading order with the heavy baryon chiral perturbation theory. The numerical results are given up to next-to-leading order: $μ_{Ξ^{*++}_{cc}}=1.72μ_{N}$, $μ_{Ξ^{*+}_{cc}}=-0.09μ_{N}$, $μ_{Ω^{*+}_{cc}}=0.99μ_{N}$. As a by-product, we have also calculated the magnetic moments of the spin-$3\over 2$ doubly bottom baryons and charmed bottom baryons: $μ_{Ξ^{*0}_{bb}}=0.63μ_{N}$, $μ_{Ξ^{*-}_{bb}}=-0.79μ_{N}$, $μ_{Ω^{*-}_{bb}}=0.12μ_{N}$, $μ_{Ξ^{*+}_{bc}}=1.12μ_{N}$, $μ_{Ξ^{*0}_{bc}}=-0.40μ_{N}$, $μ_{Ω^{*0}_{bc}}=0.56μ_{N}$. △ Less

Submitted 7 June, 2018; v1 submitted 19 October, 2017; originally announced October 2017.

Comments: 10 pages,2 figures. arXiv admin note: text overlap with arXiv:1707.02765. Replace the published version

$D D^{*}$ potentials in chiral perturbation theory and possible molecular states

Authors: Hao Xu, Bo Wang, Zhan-Wei Liu, Xiang Liu

Abstract: The $DD^{*}$ potentials are studied within the framework of heavy meson chiral effective field theory. We obtain the effective potentials of the $DD^{*}$ system up to $O(ε^2)$ at the one-loop level. In addition to the one-pion exchange contribution, the contact and two-pion exchange interactions are also investigated in detail. Furthermore, we search for the possible molecular states by solving th… ▽ More The $DD^{*}$ potentials are studied within the framework of heavy meson chiral effective field theory. We obtain the effective potentials of the $DD^{*}$ system up to $O(ε^2)$ at the one-loop level. In addition to the one-pion exchange contribution, the contact and two-pion exchange interactions are also investigated in detail. Furthermore, we search for the possible molecular states by solving the Schrödinger equation with the potentials. We notice that the contact and two-pion exchange potentials are numerically non-negligible and important for the existence of a bound state. In our results, no bound state is found in the $I=1$ channel within a wide range of cutoff parameter, while there exists a bound state in the $I=0$ channel as the cutoff is near $m_ρ$ in our approach. △ Less

Submitted 5 December, 2021; v1 submitted 23 August, 2017; originally announced August 2017.

Comments: 18 pages, 12 figures. Published in Phys. Rev. D. We correct three typos in which $I=1$ and $I=0$ were mistaken for each other

Journal ref: Phys. Rev. D 99, 014027 (2019)

Radiative decays of the doubly charmed baryons in chiral perturbation theory

Authors: Hao-Song Li, Lu Meng, Zhan-Wei Liu, Shi-Lin Zhu

Abstract: We have systematically investigated the spin-$\frac{3}{2}$ to spin-$\frac{1}{2}$ doubly charmed baryon transition magnetic moments to the next-to-next-to-leading order in the heavy baryon chiral perturbation theory (HBChPT). Numerical results of transition magnetic moments and decay widths are presented to the next-to-leading order: $μ_{Ξ_{cc}^{*++}\rightarrowΞ_{cc}^{++}}=-2.35μ_{N}$,… ▽ More We have systematically investigated the spin-$\frac{3}{2}$ to spin-$\frac{1}{2}$ doubly charmed baryon transition magnetic moments to the next-to-next-to-leading order in the heavy baryon chiral perturbation theory (HBChPT). Numerical results of transition magnetic moments and decay widths are presented to the next-to-leading order: $μ_{Ξ_{cc}^{*++}\rightarrowΞ_{cc}^{++}}=-2.35μ_{N}$, $μ_{Ξ_{cc}^{*+}\rightarrowΞ_{cc}^{+}}=1.55μ_{N}$, $μ_{Ω_{cc}^{*+}\rightarrowΩ_{cc}^{+}}=1.54μ_{N}$, $Γ_{Ξ_{cc}^{*++}\rightarrowΞ_{cc}^{++}}=22.0$ keV, $Γ_{Ξ_{cc}^{*+}\rightarrowΞ_{cc}^{+}}=9.57$ keV, $Γ_{Ω_{cc}^{*+}\rightarrowΩ_{cc}^{+}}=9.45$ keV. △ Less

Submitted 26 September, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: arXiv admin note: text overlap with arXiv:1707.02765, arXiv:1706.06458

Magnetic moments of the doubly charmed and bottomed baryons

Authors: Hao-Song Li, Lu Meng, Zhan-Wei Liu, Shi-Lin Zhu

Abstract: The chiral corrections to the magnetic moments of the spin-$\frac{1}{2}$ doubly charmed baryons are systematically investigated up to next-to-next-to-leading order with heavy baryon chiral perturbation theory (HBChPT). The numerical results are calculated up to next-to-leading order: $μ_{Ξ^{++}_{cc}}=-0.25μ_{N}$, $μ_{Ξ^{+}_{cc}}=0.85μ_{N}$, $μ_{��^{+}_{cc}}=0.78μ_{N}$. We also calculate the magneti… ▽ More The chiral corrections to the magnetic moments of the spin-$\frac{1}{2}$ doubly charmed baryons are systematically investigated up to next-to-next-to-leading order with heavy baryon chiral perturbation theory (HBChPT). The numerical results are calculated up to next-to-leading order: $μ_{Ξ^{++}_{cc}}=-0.25μ_{N}$, $μ_{Ξ^{+}_{cc}}=0.85μ_{N}$, $μ_{Ω^{+}_{cc}}=0.78μ_{N}$. We also calculate the magnetic moments of the other doubly heavy baryons, including the doubly bottomed baryons (bbq) and the doubly heavy baryons containing a light quark, a charm quark and a bottom quark ($\{bc\}q$ and $[bc]q$): $μ_{Ξ^{0}_{bb}}=-0.84μ_{N}$, $μ_{Ξ^{-}_{bb}}=0.26μ_{N}$, $μ_{Ω^{-}_{bb}}=0.19μ_{N}$, $μ_{Ξ^{+}_{\{bc\}q}}=-0.54μ_{N}$, $μ_{Ξ^{0}_{\{bc\}q}}=0.56μ_{N}$, $μ_{Ω^{0}_{\{bc\}q}}=0.49μ_{N}$, $μ_{Ξ^{+}_{[bc]q}}=0.69μ_{N}$, $μ_{Ξ^{0}_{[bc]q}}=-0.59μ_{N}$, $μ_{Ω^{0}_{[bc]q}}=0.24μ_{N}$. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 11 pages, 2 figures

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

Decuplet to octet baryon transitions in chiral perturbation theory

Authors: Hao-Song Li, Zhan-Wei Liu, Xiao-Lin Chen, Wei-Zhen Deng, Shi-Lin Zhu

Abstract: We have systematically investigated the decuplet (T) to octet (B) baryon ($T\rightarrow Bγ$) transition magnetic moments to the next-to-next-to-leading order and electric quadruple moments to the next-to-leading order in the framework of the heavy baryon chiral perturbation theory. Our calculation includes the contributions from both the intermediate decuplet and octet baryon states in the loops.… ▽ More We have systematically investigated the decuplet (T) to octet (B) baryon ($T\rightarrow Bγ$) transition magnetic moments to the next-to-next-to-leading order and electric quadruple moments to the next-to-leading order in the framework of the heavy baryon chiral perturbation theory. Our calculation includes the contributions from both the intermediate decuplet and octet baryon states in the loops. Our results show reasonably good convergence of the chiral expansion and agreement with the experimental data. The analytical expressions may be useful to the chiral extrapolation of the lattice simulations of the decuplet electromagnetic properties. △ Less

Submitted 18 June, 2017; originally announced June 2017.

Comments: 20 pages, 2 figures. arXiv admin note: substantial text overlap with arXiv:1608.04617

Spectroscopy with Local Multi-hadron Interpolators in Lattice QCD

Authors: Adrian L. Kiratidis, Waseem Kamleh, Derek B. Leinweber, Zhan-Wei Liu, Finn M. Stokes, Anthony W. Thomas

Abstract: The positive-parity nucleon spectrum is studied in 2 + 1 flavour lattice QCD in an attempt to discover novel low-lying energy eigenstates in the region of the Roper resonance. In this work, we employ standard three-quark interpolating fields and introduce new local five-quark meson-baryon operators that hold the possibility of revealing new states that have been missed in previous studies. Motivat… ▽ More The positive-parity nucleon spectrum is studied in 2 + 1 flavour lattice QCD in an attempt to discover novel low-lying energy eigenstates in the region of the Roper resonance. In this work, we employ standard three-quark interpolating fields and introduce new local five-quark meson-baryon operators that hold the possibility of revealing new states that have been missed in previous studies. Motivated by phenomenological arguments, five-quark interpolators based on the $σ{N}$, $π{N}$ and $a_0{N}$ channels are constructed. Despite the introduction of qualitatively different operators, no novel energy levels are extracted near the regime of the Roper resonance. △ Less

Submitted 28 April, 2017; originally announced April 2017.

Comments: INPC 2016 Conference proceedings; 7 pages, 2 Figures, 1 Table

Report number: ADP-17-15/T1021

Structure of the Roper Resonance from Lattice QCD Constraints

Authors: Jia-jun Wu, Derek B. Leinweber, Zhan-wei Liu, Anthony W. Thomas

Abstract: Two different descriptions of the existing pion-nucleon scattering data in the region of the Roper resonance are constructed. Both descriptions fit the experimental data very well. In one scenario the resonance is the result of strong rescattering between coupled meson-baryon channels, while in the other scenario, the resonance has a large bare-baryon (or quark-model like) component. The predictio… ▽ More Two different descriptions of the existing pion-nucleon scattering data in the region of the Roper resonance are constructed. Both descriptions fit the experimental data very well. In one scenario the resonance is the result of strong rescattering between coupled meson-baryon channels, while in the other scenario, the resonance has a large bare-baryon (or quark-model like) component. The predictions of these two scenarios are compared with the latest lattice QCD simulation results in this channel. Consideration of the finite volume spectra, the manner in which the states are excited from the vacuum in lattice QCD and the composition of the states in Hamiltonian effective field theory enable a discrimination of these two different descriptions. We find the second scenario is not consistent with lattice QCD results whereas the first agrees with the lattice QCD constraints. In this scenario, the mass of the quark-model like state is approximately 2 GeV and in the finite volume of the lattice is dressed to place the first radial excitation of the nucleon at 1.9 GeV. Within this description, the infinite-volume Roper resonance is best described as a resonance generated dynamically through strongly coupled meson-baryon channels. △ Less

Submitted 14 May, 2018; v1 submitted 30 March, 2017; originally announced March 2017.

Comments: 5 pages, 5 figures

Report number: ADP-17-13/T1019

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

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

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

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

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

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

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