Source author record

Yi-Bo Yang

Yi-Bo Yang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

46 published item(s)

preprint2026arXiv

Total Gluon Helicity Contribution to the Proton Spin from Lattice QCD

We report a state-of-the-art lattice QCD calculation of the total gluon helicity contribution to the proton spin, $ΔG$. The calculation is done on ensembles with three different lattice spacings $a=\{0.08, 0.09, 0.11\}$ fm. By employing distillation and momentum smearing for proton external states, we extract the bare matrix elements of the topological current $K^μ$ using 5-HYP smeared Coulomb gauge fixing configurations. Furthermore, we apply a non-perturbative $\mathrm{RI/MOM}$ renormalization scheme augmented by the Cluster Decomposition Error Reduction (CDER) technique to determine the renormalization constants of $K^μ$. The results obtained from different components $K^{t,i}$ (with $i$ being the direction of proton momentum or polarization) are consistent with Lorentz covariance within uncertainties. After extrapolating to the continuum limit, $ΔG$ is found to be $ΔG = 0.231(17)^{\mathrm{sta.}}(44)^{\mathrm{sym.}}$ at the $\overline{\mathrm{MS}}$ scale $μ^2=10\ \mathrm{GeV}^2$, which constitutes approximately $46(9)\%$ of the proton spin.

preprint2024arXiv

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

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

preprint2023arXiv

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

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

preprint2022arXiv

Distribution Amplitudes of $K^*$ and $ϕ$ at Physical Pion Mass from Lattice QCD

We present the first lattice QCD calculation of the distribution amplitudes of longitudinally and transversely polarized vector mesons $K^*$ and $ϕ$ using large momentum effective theory. We use the clover fermion action on three ensembles with 2+1+1 flavors of highly improved staggered quarks (HISQ) action, generated by MILC collaboration, at physical pion mass and \{0.06, 0.09, 0.12\} fm lattice spacings, and choose three different hadron momenta $P_z=\{1.29, 1.72, 2.15\}$ GeV. The resulting lattice matrix elements are nonperturbatively renormalized in a hybrid scheme proposed recently. An extrapolation to the continuum and infinite momentum limit is carried out. We find that while the longitudinal distribution amplitudes tend to be close to the asymptotic form, the transverse ones deviate rather significantly from the asymptotic form. Our final results provide crucial {\it ab initio} theory inputs for analyzing pertinent exclusive processes.

preprint2022arXiv

First Lattice QCD determination of semileptonic decays of charmed-strange baryons $Ξ_c$

While the standard model is the most successfully theory to describe all interactions and constituents in elementary particle physics, it has been constantly examined for over four decades. Weak decays of charm quarks can measure the coupling strength of quarks in different families and serve as an ideal probe for CP violation. As the lowest charm-strange baryons with three different flavors, $Ξ_c$ baryons (made of $csu$ or $csd$) have been extensively studied in experiments at the large hadron collider and in electron-positron collision. However the lack of reliable knowledge in theory becomes the unavoidable obstacle in the way. In this work, we use the state-of-the-art Lattice QCD techniques, and generate 2+1 clover fermion ensembles with two lattice spacings, $a=(0.108{\rm fm},0.080{\rm fm})$. We then present the first {\it ab-initio} lattice QCD determination of form factors governing $Ξ_{c}\to Ξ\ell^+ν_{\ell}$, analogous with the notable $β$-decay of nuclei. Our theoretical results for decay widths are consistent with and about two times more precise than the latest measurements by ALICE and Belle collaborations. Together with experimental measurements, we independently determine the quark-mixing matrix element $|V_{cs}|$, which is found in good agreement with other determinations.

preprint2022arXiv

Low energy constant and mixed-action effect

We present the pion mass and decay constant using the overlap fermion valence on domain wall (DW) fermion sea at several lattice spacings. The mixed action effect in the lattice calculation is also studied, and the result suggests that the mixed action effect with the overlap valence on DW sea would be proportional to the fourth power of the lattice spacing. We obtain the pion decay constant at the physical pion mass and $N_f=2$ chiral limit to be $92.4(3)(2)~{\rm MeV}$ and $87.0(5)(7)~{\rm MeV}$ respectively; and the physical u/d averaged quark mass at $\overline{\textrm{MS}}$ 2 GeV is $3.74(4)(5)(5)(3)~{\rm MeV}$ with the linear ${\cal O}(a^2)$ continuum extrapolation. Using the FLAG value of the NLO low energy constant, we obtain the $N_f=2$ chiral condensate to be $Σ^{\overline{\textrm{MS}}(2~{\rm GeV})}=\big(266(2)(1)~{\rm MeV}\big)^3$.

preprint2022arXiv

Nonperturbative Determination of Collins-Soper Kernel from Quasi Transverse-Momentum Dependent Wave Functions

In the framework of large-momentum effective theory at one-loop matching accuracy, we perform a lattice calculation of the Collins-Soper kernel which governs the rapidity evolution of transverse-momentum-dependent (TMD) distributions. We first obtain the quasi TMD wave functions at three different meson momenta on a lattice with valence clover quarks on a dynamical HISQ sea and lattice spacing $a=0.12$~fm from MILC, and renormalize the pertinent linear divergences using Wilson loops. Through one-loop matching to the light-cone wave functions, we determine the Collins-Soper kernel with transverse separation up to 0.6~fm. We study the systematic uncertainties from operator mixing and scale dependence, as well as the impact from higher power corrections. Our results potentially allow for a determination of the soft function and other transverse-momentum dependent quantities at one-loop accuracy.

preprint2022arXiv

Proton momentum and angular momentum decompositions with overlap fermions

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

preprint2022arXiv

Renormalization of transverse-momentum-dependent parton distribution on the lattice

To calculate the transverse-momentum-dependent parton distribution functions (TMDPDFs) from lattice QCD, an important goal yet to be realized, it is crucial to establish a viable non-perturbative renormalization approach for linear divergences in the corresponding Euclidean quasi-TMDPDF correlators in large-momentum effective theory. We perform a first systematic study of the renormalization property of the quasi-TMDPDFs by calculating the relevant matrix elements in a pion state at 5 lattice spacings ranging from 0.03 fm to 0.12 fm. We demonstrate that the square root of the Wilson loop combined with the short distance hadron matrix element provides a successful method to remove all ultraviolet divergences of the quasi-TMD operator, and thus provide the necessary justification to perform a continuum limit calculation of TMDPDFs. In contrast, the popular RI/MOM renormalization scheme fails to eliminate all linear divergences.

preprint2022arXiv

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

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

preprint2021arXiv

A Hybrid Renormalization Scheme for Quasi Light-Front Correlations in Large-Momentum Effective Theory

In large-momentum effective theory (LaMET), calculating parton physics starts from calculating coordinate-space-$z$ correlation functions $\tilde h(z, a,P^z)$ in a hadron of momentum $P^z$ in lattice QCD. Such correlation functions involve both linear and logarithmic divergences in lattice spacing $a$, and thus need to be properly renormalized. We introduce a hybrid renormalization procedure to match these lattice correlations to those in the continuum $\overline{\rm MS}$ scheme, without introducing extra non-perturbative effects at large $z$. We analyze the effect of ${\cal O}(Λ_{\rm QCD})$ ambiguity in the Wilson line self-energy subtraction involved in this hybrid scheme. To obtain the momentum-space distributions, we recommend to extrapolate the lattice data to the asymptotic $z$-region using the generic properties of the coordinate space correlations at moderate and large $P^z$, respectively.

preprint2020arXiv

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

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

preprint2020arXiv

Lattice QCD package GWU-code and QUDA with HIP

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

preprint2020arXiv

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

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

preprint2020arXiv

Roper State from Overlap Fermions

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

preprint2020arXiv

The nucleon isovector tensor charge from lattice QCD using chiral fermions

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

preprint2020arXiv

Trace anomaly and dynamical quark mass

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.

preprint2019arXiv

Towards the nucleon hadronic tensor from lattice QCD

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

preprint2019arXiv

Unpolarized isovector quark distribution function from Lattice QCD: A systematic analysis of renormalization and matching

We present a detailed Lattice QCD study of the unpolarized isovector quark Parton Distribution Function (PDF) using large-momentum effective theory framework. We choose a quasi-PDF defined by a spatial correlator which is free from mixing with other operators of the same dimension. In the lattice simulation, we use a Gaussian-momentum-smeared source at $M_π=356$ MeV and $P_z \in \{1.8,2.3\}$ GeV. To control the systematics associated with the excited states, we explore {five different source-sink separations}. The nonperturbative renormalization is conducted in a regularization-independent momentum subtraction scheme, and the matching between the renormalized quasi-PDF and $\bar{\rm MS}$ PDF is calculated based on perturbative QCD up to one-loop order. Systematic errors due to renormalization and perturbative matching are also analyzed in detail. Our results for lightcone PDF are in reasonable agreement with the latest phenomenological analysis.

preprint2016arXiv

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

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

preprint2016arXiv

Glue Spin of the Proton

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

preprint2016arXiv

Stochastic method with low mode substitution for nucleon isovector matrix elements

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

preprint2016arXiv

Unphysical Poles of Domain Wall Fermions at finite $L_s$

We investigate the origin and behavior of oscillations observed in the hadron correlators constructed from the domain wall fermion (DWF) for different parameters involved in lattice QCD simulations. This oscillatory behavior at early time slices hinders the extraction of excited states in hadron spectroscopy. Furthermore, the deviation from exponential decay may have a significant impact on fermion loop calculations performed on the lattice. We present results for several well-known implementations of the DWF actions. We extend the study of Shamir DWF action to include Boriçi and Möbius DWF actions by analyzing the poles of 4D quark propagator. For each action considered, we find an unphysical mode when analyzing the pole structure of the free DWF propagator for a finite extent of the $5^{\text{th}}$ dimension $L_s$, and we show that this mode is responsible for the oscillatory behavior observed in hadron correlators. We have performed numerical checks on these results and have found that the presence of oscillatory behavior is sensitive to the DWF parameters $a_5$, $b_5$, $c_5$ and DW height $M$. To minimize oscillations, our results suggest that one should choose $M a_5<1$ for the Shamir and Boriçi DWFs, and $M(b_5-c_5)<1$ for Möbius DWF when $M>1$. For each calculation considered, the Borçi DWF displayed the smallest magnitude of oscillation when compared to the Shamir and Möbius DWF actions using the same input parameters.

preprint2015arXiv

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

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

preprint2015arXiv

Diquark mass differences from unquenched lattice QCD

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

preprint2015arXiv

Glue Spin $S_G$ in The Longitudinally Polarized Nucleon

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

preprint2015arXiv

Meson Mass Decomposition from Lattice QCD

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

preprint2015arXiv

Quark Spin in Proton from Anomalous Ward Indentity

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

preprint2015arXiv

Strange quark momentum fraction from overlap fermion

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

preprint2014arXiv

Glueballs in charmonia radiative decays

Scalar \cite{scalar_paper} and tensor \cite{tensor_paper} glueballs created in $J/ψ$ radiative decays are studied in quenched lattice QCD. Using two anisotropic lattices to approach the continuum limit, we compute the relevant form factors responsible for the decay rates for $J/ψ\rightarrowγG_{0^{++}}$ and $J/ψ\rightarrowγG_{2^{++}}$. Comparing with the existing experimental data, it is argued that $f_0(1710)$ is a favorable candidate for scalar glueball. The decay rate for $J/ψ\rightarrowγG_{2^{++}}$ is found to be quite substantial. A comprehensive search in the tensor channel on BESIII is therefore suggested.

preprint2014arXiv

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

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

preprint2014arXiv

Meson Mass Decomposition

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

preprint2014arXiv

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

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

preprint2014arXiv

Oscillatory behavior of the domain wall fermions revisited

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

preprint2013arXiv

Is $1^-+$ Meson a Hybrid?

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

preprint2013arXiv

Lattice Study of Radiative $J/ψ$ Decay to a Tensor Glueball

The radiative decay of $J/ψ$ into a pure gauge tensor glueball is studied in the quenched lattice QCD formalism. With two anisotropic lattices, the mutlipole amplitudes E_1(0), M_2(0) and E_3(0) are obtained to be 0.114(12)(6)GeV, -0.011(5)(1)GeV, and 0.023(8)(1)GeV, respectively. The first error comes from the statistics, the Q^2 interpolation, and the continuum extrapolation, while the second is due to the uncertainty of the scale parameter r_0^{-1}=410(20) MeV. Thus the partial decay width $Γ(J/ψ\rightarrow γG_{2^{++}})$ is estimated to be 1.01(22)(10) keV which corresponds to a large branch ratio 1.1(2)(1)x10^{-2}. The phenomenological implication of this result is also discussed.

preprint2013arXiv

Lattice study on $η_{c2}$ and X(3872)

Properties of $2^{-+}$ charmonium $η_{c2}$ are investigated in quenched lattice QCD. The mass of $η_{c2}$ is determined to be 3.80(3) GeV, which is close to the mass of $D$-wave charmonium $ψ(3770)$ and in agreement with quark model predictions. The transition width of $η_{c2}\to γJ/ψ$ is also obtained with a value $Γ=3.8(9)$ keV. Since the possible $2^{-+}$ assignment to X(3872) has not been ruled out by experiments, our results help to clarify the nature of X(3872).

preprint2013arXiv

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

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

preprint2013arXiv

Scalar Glueball in Radiative $J/ψ$ Decay on Lattice

The form factors in the radiative decay of $J/ψ$ to a scalar glueball are studied within quenched lattice QCD on anisotropic lattices. The continuum extrapolation is carried out by using two different lattice spacings. With the results of these form factors, the partial width of $J/ψ$ radiatively decaying into the pure gauge scalar glueball is predicted to be 0.35(8) keV, which corresponds to a branching ratio of 3.8(9)x10^{-3}. By comparing with the experiments, out results indicate that f_0(1710) has a larger overlap with the pure gauge glueball than other related scalar mesons.

preprint2012arXiv

Charmless $B_s\to PP, PV, VV$ Decays Based on the six-quark Effective Hamiltonian with Strong Phase Effects II

We provide a systematic study of charmless $B_s \to PP, PV, VV$ decays ($P$ and $V$ denote pseudoscalar and vector mesons, respectively) based on an approximate six-quark operator effective Hamiltonian from QCD. The calculation of the relevant hard-scattering kernels is carried out, the resulting transition form factors are consistent with the results of QCD sum rule calculations. By taking into account important classes of power corrections involving "chirally-enhanced" terms and the vertex corrections as well as weak annihilation contributions with non-trivial strong phase, we present predictions for the branching ratios and CP asymmetries of $B_s$ decays into PP, PV and VV final states, and also for the corresponding polarization observables in VV final states. It is found that the weak annihilation contributions with non-trivial strong phase have remarkable effects on the observables in the color-suppressed and penguin-dominated decay modes. In addition, we discuss the SU(3) flavor symmetry and show that the symmetry relations are generally respected.

preprint2011arXiv

Predictive AdS/QCD Model for Mass Spectra of Mesons with Three Flavors

The predictive soft-wall AdS/QCD model with a modified 5D metric at the infrared region is constructed to obtain a non-trivial dilaton solution for three flavor quarks $u$, $d$ and $s$. Such a model is shown to incorporate both the chiral symmetry breaking and linear confinement. After considering some high-order terms including the $U(1)_L\times U(1)_R$ chiral symmetry breaking term, we find that the resulting predictions for the SU(3) octet and singlet resonance states of pseudoscalar, scalar, vector and axial-vector mesons agree well with the experimentally confirmed resonance states. Contributions from the instanton effects given by the determinant term are also discussed. It is observed that the chiral symmetry breaking phenomena of $SU(3)_L\times SU(3)_R$ and $U(1)_L\times U(1)_R$ can be well described in this model, while the SU(3) flavor symmetry breaking effect due to quark mass difference in the source term is not enough to explain all of the current experimental data.

preprint2010arXiv

Charmless $B\to PP, PV, VV$ Decays Based on the six-quark Effective Hamiltonian with Strong Phase Effects

Based on an approximate six-quark operator effective Hamiltonian from perturbative QCD, we present a systematical study of charmless $B \to PP, PV, VV$ decays ($P$ and $V$ denoting pseudoscalar and vector mesons, respectively). The calculation of the relevant hard-scattering kernels is completed, the resulting transition form factors are consistent with the results of QCD sum rule calculation. Important classes of power corrections include "chirally-enhanced" terms, vertex corrections and weak annihilation contributions with non-trivial strong phase. With these considerations, predictions are presented for the branching ratios and CP asymmetries of B-meson decays into PP, PV and VV final states, and also for the corresponding polarization observables in VV final states. Several decay modes and observables, which are of particular interest phenomenologically, are discussed in detail, including the effect of annihilation amplitude with strong phase, the $ππ$, $πK$ and $πρ$ systems, the longitudinal polarization fraction $f_L$ in $ρK^*$ and $ϕK^*$ systems and so on. It is observed that predictions in our framework generally agree with the current experimental data.

preprint2010arXiv

Prediction for the Mass Spectra of Resonance Mesons in the Soft-Wall AdS/QCD with a Modified 5D Metric

A soft-wall anti-de Sitter/QCD model with a modified five-dimensional metric at the infrared region is constructed to obtain a nontrivial dilaton solution, which incorporates the chiral symmetry breaking and linear confinement. By taking the pion mass and decay constant as two input mass scales, the resulting predictions for the resonance states of pseudoscalar, scalar, vector and axial-vector mesons agree remarkably with the experimentally confirmed resonance states. The effects of the quartic interaction term are investigated by taking an appropriate sign and magnitude for maintaining the stability of the bulk scalar potential. It is shown that such a simply modified soft-wall anti-de Sitter/QCD model can lead to a consistent prediction for the mass spectra of resonance states in the pseudoscalar, scalar, vector and axial-vector mesons; the agreement with the experimental data is found to be better than 10% for the excited meson states. The resulting pion form factor also agrees well with the experimental data.

preprint2008arXiv

QCD Factorization Based on Six-Quark Operator Effective Hamiltonian from Perturbative QCD and Charmless Bottom Meson Decays $B_{(s)}\to ππ,πK, KK$

The charmless bottom meson decays are systematically investigated based on an approximate six quark operator effective Hamiltonian from perturbative QCD. It is shown that within this framework the naive QCD factorization method provides a simple way to evaluate the hadronic matrix elements of two body mesonic decays. The singularities caused by on mass-shell quark propagator and gluon exchanging interaction are appropriately treated. Such a simple framework allows us to make theoretical predictions for the decay amplitudes with reasonable input parameters. The resulting theoretical predictions for all the branching ratios and CP asymmetries in the charmless $B^0, B^+, B_s\to ππ, πK, KK$ decays are found to be consistent with the current experimental data except for a few decay modes. The observed large branching ratio in $B\to π^0π^0$ decay remains a puzzle though the predicted branching ratio may be significantly improved by considering the large vertex corrections in the effective Wilson coefficients. More precise measurements of charmless bottom meson decays, especially on CP-violations in $B\to K K$ and $B_s\to ππ, πK, KK$ decay modes, will provide a useful test and guide us to a better understanding on perturbative and nonperturbative QCD.