Processing math: 100%

Investigating S-wave bound states composed of two pseudoscalar mesons

  • In this study, we systematically investigated two-pseudoscalar meson systems with the Bethe-Salpeter equation in the ladder and instantaneous approximations. By solving the Bethe-Salpeter equation numerically with the kernel containing the one-particle exchange diagrams, we found that the KˉK, DK, BˉK, DˉD, BˉB, BD, DˉK, BK, and BˉD systems with I=0 can exist as bound states. We also studied the contributions from heavy meson (J/ψ and Υ) exchanges and found that the contributions from heavy meson exchanges cannot be ignored.
  • The European Muon Collaboration (EMC) effect, a type of nuclear medium modification in the valence quark regime of 0.3x0.75, refers to the noticeable deviation from unity of the structure-function ratio between the heavy nucleus (A>2) and the deuteron [19]. The deuteron is considered here as the reference nucleus, which is approximately regarded as a system of a free proton plus a free neutron. The EMC effect was first discovered in the muon-induced deep inelastic scattering (DIS) off the iron nucleus [1, 2]. Soon after this discovery, the EMC effect was confirmed by the electron-nucleus DIS data at SLAC [3, 4]. To date, many experimental measurements of the EMC effect have been reported on various nuclear targets.

    The EMC effect is surprising and attracts a signifiant amount of interest from theorists. The nuclear structure function is measured with the hard probe of momentum above GeV, while the per-nucleon binding energy inside the nucleus is approximately in the MeV scale. At the fundamental level of quarks and gluons, it is not clear how the relatively long-range nuclear force modifies the short-distance structure inside the nucleon. In the quark-parton model, the structure function is the incoherent summation of the quark distributions inside the nucleon. The discovery of the EMC effect implies that the quark distribution is evidently modified by the nuclear medium. If the quarks are completely confined inside the nucleon and do not play a role in the emergence of nuclear force, then the quark distribution should not be modified with the presence of surrounding nucleons. After decades of studies, many models have been constructed that properly describe the EMC effect, such as the off-shell correction [10, 11], x-rescaling model [12, 13], nucleon swelling and dynamical rescaling model [1419], cluster model [2027], point-like configuration suppression model [28, 29], and statistical model [30]. To differentiate the various models, more experiments and new observables beyond the F2 ratio are expected.

    It is speculated that the strength of the EMC effect depends on the local density instead of the global average density of the nucleus, according to the measurements of the very light nuclei 3He and 9Be [31]. This interesting finding stimulates the physicists to imagine that the nuclear EMC effect emerges mainly from the local cluster structures inside the nucleus. The other hypothesis is that the EMC effect primarily results from the high-virtuality nucleons, which belong to the short-distance configurations of nucleons. This is supported by the unexpected linear correlation between the EMC effect and nucleon-nucleon short-range correlations (N-N SRC) [32, 33].

    The N-N SRC pairs are the temporary close-proximity fluctuations of two strongly interacting nucleons [3437]. Experimentally, the N-N SRC pairs are identified as the nucleon pairs of high relative momentum between nucleons and small center-of-mass momentum of the pair [3842]. The nucleons in SRC exhibit a much higher momenta than the nuclear Fermi momentum kF. The abundance of N-N SRC pairs can be simply characterized by the probability of finding high-momentum nucleons [4345]. These nucleons can form a close-proximity configuration and are also sensitive to the repulsive core of nucleon-nucleon interaction [46]. It is widely accepted that the intermediate-distance tensor force is the primary source of the formation of N-N SRC [4649].

    Inspired by the observed linear correlation between the EMC effect and N-N SRC [32, 33], some nuclear physicists suggest that SRC pairs may be the underlying source of the EMC effect. This assumption is close to the traditional cluster model for the EMC effect [2027]. The difference is that the cluster model is at the parton level with the six-quark bag picture, whereas the SRC explanation is based on the nucleon degrees of freedom, from which the properties of SRC nucleon are greatly modified. In a recent theoretical work, it was argued that the linear correlation between the EMC effect and N-N SRC is the natural result of the scale separation of the nucleon structure part (Λ-independent) and twist-four part (nuclear modification, Λ-dependent) of nuclear matrix elements [50]. It was demonstrated that the linear correlation between the EMC effect and N-N SRC can be derived in the effective field theory.

    One key and intriguing question is whether there is causality between the EMC effect and SRC. Recently, the CLAS collaboration tested the SRC-driven EMC model with the simultaneous measurements of DIS and quasi-elastic inclusive process on the deuteron and heavier nuclei [45]. They extracted the modification function of the nucleon structure in SRC pairs and found that this modification function is nucleus-independent [45]. They showed that the EMC effect in all measured nuclei is consistent with the universal modification function of SRC pairs, and that the magnitude of the EMC effect in the nucleus can be described by the number of SRC pairs. In their view, the EMC effect is not the traditional static modification on all the independent nucleons but a strong dynamical effect for short time intervals of two strongly interacting nucleons fluctuating into a temporary high-local-density SRC pair [45]. The universal modification function of SRC was also carefully studied by J. Arrington and N. Fomin [51]. They found that there is almost no A-dependence of the universal modification function extracted with the Local-Density model, while there is weak A-dependence with the High-Virtuality model. The universal modification function resulting from data of various nuclei is consistent with a truly universal function, and the Local-Density hypothesis is favored [51].

    A different approach has also been investigated and reported. Recently, in Ref. [52], the relationship between SRC and the EMC effect was further examined by incorporating the nuclear binding and nucleon off-shell effects. The authors argued that their analysis does not support the hypothesis that there is a causal connection between nucleons residing in SRCs and the EMC effect [52]. The EMC effect of the low-momentum and high-momentum nuclei are separately studied. The authors found that the Fermi motion effect is much stronger than the off-shell effect for SRC nucleons, with three different models for the off-shell effect [52]. Hence, they concluded that the SRC nucleons do not cause the dominant EMC effect [52]. This conclusion is contrary to what one expected in the past. It is worth further examining the relationship between the EMC effect and SRC from different viewpoints or theoretical models.

    Currently, the approaches for describing the EMC effect can be classified into the following three categories: (i) all the nucleons are slightly modified when embedded in the nuclear medium; (ii) nucleons are unmodified most of the time but greatly modified when they fluctuate into N-N SRC; (iii) mean-field uncorrelated nucleons are slightly modified, and the nucleons are substantially modified for a short-time interval in the temporary SRC state. It is a prominent research topic to study whether the nuclear EMC effect entirely comes from the N-N SRC. Hence, in this work, we focus on the last two approaches to explain the EMC effect. In Sec. II, we illustrate the models we use to calculate the nuclear EMC effect. In Sec. III, we present the results of the EMC effect from SRC and mean-field nucleons. Finally, a short summary is provided in Sec. IV.

    It is known that the traditional nuclear structure is almost irrelevant to the nuclear EMC effect. The nucleon momentum distribution leads to the Fermi motion effect around x=1. The per-nucleon nuclear binding energy is much smaller than the high-momentum virtual photon probe or nucleon mass. Nevertheless, the identity of nucleons inside nucleus is well established and it is the core of the traditional nuclear physics. The nucleon structure should be connected to the properties of the nucleon.

    In this study, we applied the x-rescaling model to evaluate the EMC effect, which is analyzed on the view that the nucleus is a sum of quasi-particles (bound nucleons). The mass is a fundamental property of the nucleon, and the change in the nucleon mass inside the nucleus should be taken into account for the nuclear medium effect. The nucleon effective mass in nucleus has been successfully used to describe the nuclear EMC effect [12, 13]. The Bjorken scaling variable is defined in terms of the free nucleon mass m as x=Q2/(2mν). However, the struck nucleon in lepton-nucleus DIS could be far off-shell. The true scaling variable for nuclear DIS should be assumed to be x=Q2/(2mν)=xm/m=xη, where m is the effective mass of the bound nucleon. Here, η=m/m is the rescaling factor of x, and the per-nucleon nuclear structure function FA2 is given by

    FA2(x,Q2)=FN2(xη,Q2),

    (1)

    where FN2 is the free nucleon structure function. The rescaling of x is taken into account for the off-shell correction of the bound nucleon [12, 13, 53]. It was also pointed out that the exchanged virtual meson would take away a fraction of the nucleon momentum, thus resulting in the x-rescaling of the nuclear structure function [54].

    The nucleon effective mass in the x-scaling model is used to describe the off-shellness of the nucleon, with E2=p2+m2 [12, 13, 53]. The nucleon effective mass was also defined by Brueckner in 1950s within a non-relativistic many-body theory to account for the momentum-dependence of potential energy of a single particle, with E(k)=k2/2m+V(0)+bk2=k2/2m+V(0), V(k)=V(0)+bk2+..., and m=m/(1+2bm) [55]. Therefore, the effective mass of the SRC nucleon in this study was different from Brueckner's definition. Brueckner's nucleon effective mass reflects leading effects of the space-time non-locality of the underlying nuclear interactions [56, 57], while the effective mass of SRC nucleon arises from the local interactions at short distance. The relations between these two effective masses should be investigated in the future.

    The intriguing question we aimed to answer through the present study is whether the N-N SRCs are totally responsible for the nuclear EMC effect. Therefore, we consider the first model, referred as model-A, for the convenience of discussion, in which only the short-range correlated nucleons are substantially modified while the uncorrelated nucleons are nearly unmodified. This model strongly relies on the causality between the SRC and the EMC effect, i.e., the N-N SRC is the primary source of the EMC effect. For model-A, the nuclear structure function FA2 is decomposed as

    FA2=[nASRCFpinSRC2+nASRCFninSRC2+(ZnASRC)Fp2+(AZnASRC)Fn2]/A,

    (2)

    where nASRC is the number of proton-neutron SRC pairs in nucleus A, FpinSRC2 and FninSRC2 are the modified nucleon structure functions in the SRC pair, and Fp2 and Fn2 are free nucleon structure functions. In Eq. (2), Z, N, and A are respectively the proton number, neutron number, and mass number. Here, the number of SRC pairs should be viewed as the time-averaged value for the dynamical system. Given that the deuteron is in the SRC configuration occasionally, the time-averaged number of SRC pairs in the deuteron is less than one, that is, ndSRC<1.

    The SRC universality and isophobic property of N-N SRC pairs are the other two foundations of model-A. The universality of SRC can be described by a similar form of nuclear wave function at high nucleon momentum, which is confirmed by the experimental observations of the x-independence and the weak Q2-dependence of the cross section ratio between two different nuclei in the region of 1.4x2 [4345]. Different experiments have revealed that most of the SRC pairs are the proton-neutron pairs [37, 39, 40, 42, 58, 59]. This isophobic property supports the point that the immediate tensor force is the primary source for the formation of N-N SRC pairs [4649].

    For model-A, the number of SRC pairs in nucleus A and modified nucleon structure functions in SRC pair are key inputs. The number of SRC pairs in nucleus A is closely related to the measured SRC scaling ratio a2(nucleus A over the deuteron) and number of SRC pairs in the deuteron, which is written as

    nASRC=[A×a2(A)×ndSRC]/2.

    (3)

    Note that the above relation (Eq. (3)) is a simplified assumption. The SRC scaling ratio a2 is measured using the high-energy electron inclusive scattering process off the nuclear targets [4345], and the number of SRC pairs in the deuteron was determined in a previous analysis [60]. The free nucleon structure functions can be calculated with the parton distribution functions fi(x,Q2), as FN2(x,Q2)=ie2ixfi(x,Q2). In this study, the proton parton distribution functions were extracted from global analyses such as CT14 [61] and CJ15 [62]. The parton distributions of the free neutron are easily given by the parton distributions of the proton under the assumption of isospin symmetry, i.e., un=dp and dn=up. By using the x-rescaling model, the structure function of the SRC nucleon is connected with the free nucleon structure function, which is expressed as

    FpinSRC2(x,Q2)=Fp2(xηSRC,Q2),FninSRC2(x,Q2)=Fn2(xηSRC,Q2),

    (4)

    in which ηSRCis the rescaling factor for the SRC nucleon; ηSRC is directly connected with the effective mass of SRC nucleon as ηSRC=m/mSRC, which is a universal factor among different nuclei. Given that the effective mass of SRC nucleon mSRC was extracted from a correlation analysis between the nuclear mass and SRC scaling ratio a2, the rescaling factor for SRC nucleon was computed to be ηSRC=1.10 [60].

    According to the nuclear shell model, the nucleons move independently in the mutual potential created by all the nucleons, which is usually approximated with the mean field. These mean-field nucleons are mainly governed by the long-range nuclear force. As revealed by a high energy electron probe, we know that nucleon-nucleon short-range correlations exist and they constitute an important microscopic structure of the nucleus. Although short-range correlated nucleons interact intensively, they are minorities in a temporary state. A more general hypothesis is that the structure function of mean-field nucleons is slightly modified, whereas the structure function of SRC nucleons is strongly modified. In other words, the N-N SRC may not generate the enough EMC effect.

    For the second model, referred as model-B for the convenience of discussion, we propose that both the mean-field nucleons and SRC nucleons are modified by the nuclear medium or correlated partner nucleon. The nuclear structure function in model-B is decomposed as

    FA2=[nASRCFpinSRC2+nASRCFninSRC2+(ZnASRC)Fp2+(AZnASRC)Fn2]/A,

    (5)

    where FpinSRC2 and FninSRC2 denote the structure functions of SRC nucleons, and Fp2 and Fn2 denote the structure functions of mean-field nucleons. Here, the number of SRC pairs nASRC and structure functions of SRC nucleons are assumed to be the same as those of model-A. In model-B, the structure functions of mean-field nucleons are also calculated with the x-rescaling model, which is expressed as

    Fp2(x,Q2)=Fp2(xηMF,Q2),Fn2(x,Q2)=Fn2(xηMF,Q2).

    (6)

    Different from the situation for SRC nucleons, we assume that the rescaling factor ηMF for the mean-field nucleon is nucleus-dependent, given that the effective mass of the mean-field nucleon depends on the nucleus. The nucleon densities of different nuclei are different. In this analysis, we let ηMF be a free parameter for each nucleus. Note that the rescaling factor ηMF for mean-field nucleons should be smaller than the rescaling factor ηSRC for SRC nucleons.

    Figures 1 and 2 depict recent experimental measurements of the nuclear EMC effects in light and heavy nuclei, respectively. The predictions of model-A and model-B are also shown in the figures for the sake of comparison. The experimental data are extracted from the recent high-precision measurements by CLAS at JLab [45]. Note that the experimental data points are distributed in the valence quark region of x smaller than 0.6. Given that the data are far from the Fermi motion region near x1, the Fermi motion correction was neglected in this work.

    Figure 1

    Figure 1.  (color online) The predicted EMC ratios from the x-rescaling models are shown along with the experimental data (light nuclei). See the main text for details of the models. The experimental data are extracted from JLab Hall C [31]. Q2 was set as 5.3 GeV2 in the model calculations to be consistent with the experiment.

    Figure 2

    Figure 2.  (color online) The predicted EMC ratios from the x-rescaling models are shown with the experimental data (heavy nuclei). See the main text for details of the models. The experimental data are extracted from CLAS at JLab [45]. Q2 was set to 2 GeV2 in the model calculations to be consistent with the experiment.

    We found that the EMC effect from model-A is much weaker than the experimental observations. For the calculations of structure function ratios in model-A, we used the parton distribution functions of CT14 and CJ15. Note that the dependence on the data set of parton distribution functions is weak. In conclusion, based on the x-rescaling model, only the nuclear modifications from short-range correlated nucleons are not enough to reproduce the nuclear EMC effect in experiments. We speculate that the valence distribution of mean-field nucleons is also modified, the modification of valence distribution in N-N SRC is not universal in different nuclei, or some other short-distance structures beyond N-N SRC exist with strong modifications on the inner nucleon structure, such as 3N-SRC and α clusters.

    Note also that the number of proton-neutron SRC pairs in the deuteron is estimated to be ndSRC=0.041 by K. S. Egiyan et al. [43]; this value is much larger than the value from a previous analysis of ours [60]. In their analysis, the number of nucleons in N-N SRC pairs was defined as the number of nucleons of high momenta k>kF275 MeV/c [39, 43]. With this definition, a small fraction of mean-field nucleons may be misidentified as SRC nucleons, resulting in more SRC pairs than in our previous analysis. Nevertheless, Fig. 3 shows the predicted EMC ratios from model-A, taking the SRC a2 data averaged from experiments [33, 45] and ndSRC=0.041 [43]. The predicted EMC slopes are still smaller than the data by CLAS collaboration. Therefore, based on either our estimation on SRC numbers or the estimation by K. S. Egiyan et al., only the modifications on SRC nucleons are not enough to interpret the nuclear EMC effect. Then, we let the number of SRC pairs inside the deuteron be a free parameter. We found that, as the number of SRC pairs inside the deuteron increases by approximately 10%, the EMC effect can be explained only with the SRC nucleons. However, this high number of SRC pairs inside the deuteron is contradictory (much higher) with respect to the analysis based on the experimental data.

    Figure 3

    Figure 3.  (color online) The predicted EMC ratios from a x-rescaling model (model-A) are shown along with the experimental data (heavy nuclei), with different input values for the parameter ndSRC. See the main text for details of the model. The experimental data are extracted from CLAS at JLab [45]. Q2 was set to 2 GeV2 in the model calculations to be consistent with the experiment.

    In model-B, the mean-field nucleons are also modified, in addition to the SRC nucleons. We also assume in model-B that the rescaling factor ηMF is a free parameter and it depends on the nuclear medium. Thus, we performed the least square fit of model-B to the EMC ratio data in the range of 0.35<xB<0.65 to find the optimal parameter ηMF for each measured nucleus. The resulting values of ηMF are listed in Table 1. In model-B, ηMF of the deuteron is simply one; ηMF of 208Pb was determined to be 1.022 ± 0.002, which is a much smaller value than the rescaling factor for the SRC nucleon. Nevertheless, the mean-field nucleons in 208Pb are evidently modified, judged by the obtained rescaling factor ηMF. By introducing the EMC effect of the mean-field nucleon, the model-B successfully explains the nuclear EMC effect.

    Table 1

    Table 1.  The values of the fitted rescaling factor ηMF for the mean-field nucleon are listed under the framework of model-B. In this model, modifications on both the SRC and mean-field nucleons lead to the observed nuclear EMC effect. The errors come only from the fits to the EMC effect data. The uncertainties of the parameters ndSRC and a2 are not included.
    nucleus ηMF nucleus ηMF
    4He1.008 ± 0.001 9Be1.005 ± 0.002
    12C1.016 ± 0.00227Al1.021 ± 0.002
    56Fe1.027 ± 0.001208Pb1.022 ± 0.002
    DownLoad: CSV
    Show Table

    Furthermore, let us analyze the nuclear dependence of the rescaling factor ηMF for the mean-field nucleon in model-B. The correlations between ηMF and ln(A), ηMF and the nucleon density, and ηMF and the proton density are shown in Fig. 4. The nucleon density and proton density are calculated using A/(43πR3) and Z/(43πR3), respectively, in which R is the charge radius of a nucleus. The data of nuclear charge radii are extracted from Ref. [63]. Given that the radius of the neutron distribution in the nucleus may not be the same as the charge radius, we also plot the correlation between ηMF and the proton density of the nucleus. Although the linear correlation is not perfect, the rescaling factor ηMF of the mean-field nucleon is more or less correlated with the nucleon density. The obtained rescaling factor of the mean-field nucleon is proportional to the average nuclear density.

    Figure 4

    Figure 4.  (color online) Correlations among the rescaling factor ηMF of mean-field nucleon, ln(A), the average nucleon density, and the average proton density. See the main text for how these densities are calculated.

    In the x-rescaling model, we tested the idea that N-N SRC is the dominant source for the nuclear EMC effect. The nuclear EMC effects of some nuclei were calculated within the x-re scaling model under the assumptions that the SRC nucleon is universal among different nuclei and that only the inner structure of short-range correlated nucleons are modified. The input mass of the N-N SRC pair and number of SRC pairs inside the deuteron are extracted from a previous analysis of the a2 data and nuclear mass [60]. We found that the nuclear medium correction on N-N SRC is not enough to explain the EMC effect observed in experiments, if the model applied in this study is correct. This conclusion is consistent with results analyzed with the off-shellness correction [52].

    If we assume that the rescaling factor ηSRC is A-dependent, Model-A can effectively describe the experimental data. However, this assumption breaks the universality of N-N SRC, which is basically supported by experimental observations [4345] and some theoretical predictions [64, 65]. For Model-B, we found that the rescaling factor for uncorrelated nucleon is approximately linearly correlated with the nuclear density. This linear relation can be tested with further experimental measurements on more nuclear targets of different densities.

    We speculate that more origins of nucleon structure modifications beyond the short-distance configurations are needed, such as 3N SRC and α clusters. Other possible interpretations are that either the mean-field nucleon is significantly modified, or the N-N SRC pairs in different nuclei have different nuclear medium modifications. For a preliminary exploration, we show that the EMC effect can be explained if we just assume that the mean-field nucleon is also modified. Moreover, the nuclear modification on the mean-field nucleon scales with the density of the nucleus. In summary, we conclude that the SRC universality is wrong, the mean-field nucleon is also slightly modified, there are other sources beyond N-N SRC for the EMC effect such as the α cluster, or the applied x-rescaling model needs improvement.

    Strong evidence of 3N SRC has not been found in the inclusive4He/3He cross section ratio at JLab, and it is shown that isolating 3N SRC is much more challenging compared to isolating 2N SRC [66]. However, the theorists suggest that the scaling phenomenon from inclusive scattering on 3N SRC requires a high Q23 GeV2 and the current experimental situation should be improved [67]. Within Model-A, ndSRC needs to increase from 0.021 to 0.1 to explain the experimental data of the EMC effect. Considering that the 3N SRC and α cluster also contribute to the EMC effect, the numbers of 3N SRC pairs and α clusters should be of the same order as N-N SRC if the structure-function modifications of N-N SRC, 3N SRC and α are similar. Given that there is no 3N SRC and α cluster in the deuteron and that the nuclear modifications inside 3N SRC and α could be stronger than that inside N-N SRC, the numbers of 3N SRC pairs and α clusters in heavy nuclei could be smaller than that of N-N SRC pairs, inferred from the current data of the EMC effect. More experimental measurements are needed for searching other short-distance structures beyond N-N SRC.

    [1] M. Y. Han and Y. Nambu, Phys. Rev. 139, B1006-B1010 (1965) doi: 10.1103/PhysRev.139.B1006
    [2] P. A. Zyla et al. (Particle Data Group), PTEP 2020, 083C01 (2020) doi: 10.1093/ptep/ptaa104
    [3] A. Astier, L. Montanet, M. Baubillier et al., Phys. Lett. B 25, 294-297 (1967) doi: 10.1016/0370-2693(67)90020-2
    [4] R. Ammar, R. Davis, W. Kropac et al., Phys. Rev. Lett. 21, 1832-1835 (1968) doi: 10.1103/PhysRevLett.21.1832
    [5] B. Aubert et al. (BaBar), Phys. Rev. Lett. 90, 242001 (2003) doi: 10.1103/PhysRevLett.90.242001
    [6] S. Godfrey and K. Moats, Phys. Rev. D 93, 034035 (2016) doi: 10.1103/PhysRevD.93.034035
    [7] S. L. Olsen, T. Skwarnicki, and D. Zieminska, Rev. Mod. Phys. 90, 015003 (2018) doi: 10.1103/RevModPhys.90.015003
    [8] J. D. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 (1982) doi: 10.1103/PhysRevLett.48.659
    [9] J. D. Weinstein and N. Isgur, Phys. Rev. D 27, 588 (1983)
    [10] J. D. Weinstein and N. Isgur, Phys. Rev. D 41, 2236 (1990) doi: 10.1103/PhysRevD.41.2236
    [11] J. A. Oller, Nucl. Phys. A 714, 161-182 (2003) doi: 10.1016/S0375-9474(02)01360-X
    [12] V. Baru, J. Haidenbauer, C. Hanhart et al., Phys. Lett. B 586, 53-61 (2004) doi: 10.1016/j.physletb.2004.01.088
    [13] L. Y. Dai and M. R. Pennington, Phys. Rev. D 90, 036004 (2014) doi: 10.1103/PhysRevD.90.036004
    [14] H. A. Ahmed and C. W. Xiao, Phys. Rev. D 101, 094034 (2020) doi: 10.1103/PhysRevD.101.094034
    [15] Z. Q. Wang, X. W. Kang, J. A. Oller et al., Phys. Rev. D 105, 074016 (2022) doi: 10.1103/PhysRevD.105.074016
    [16] A. P. Szczepaniak, Phys. Lett. B 567, 23-26 (2003) doi: 10.1016/S0370-2693(03)00865-7
    [17] J. Hofmann and M. F. M. Lutz, Nucl. Phys. A 733, 142-152 (2004) doi: 10.1016/j.nuclphysa.2003.12.013
    [18] F. K. Guo, P. N. Shen, H. C. Chiang et al., Phys. Lett. B 641, 278-285 (2006) doi: 10.1016/j.physletb.2006.08.064
    [19] J. M. Flynn and J. Nieves, Phys. Rev. D 75, 074024 (2007) doi: 10.1103/PhysRevD.75.074024
    [20] Z. X. Xie, G. Q. Feng, and X. H. Guo, Phys. Rev. D 81, 036014 (2010) doi: 10.1103/PhysRevD.81.036014
    [21] Z. H. Guo, U. G. Meißner, and D. L. Yao, Phys. Rev. D 92, 094008 (2015) doi: 10.1103/PhysRevD.92.094008
    [22] M. L. Du, F. K. Guo, U. G. Meißner et al., Eur. Phys. J. C 77, 728 (2017) doi: 10.1140/epjc/s10052-017-5287-6
    [23] T. W. Wu, M. Z. Liu, L. S. Geng et al., Phys. Rev. D 100, 034029 (2019) doi: 10.1103/PhysRevD.100.034029
    [24] S. Y. Kong, J. T. Zhu, D. Song et al., Phys. Rev. D 104, 094012 (2021) doi: 10.1103/PhysRevD.104.094012
    [25] B. L. Huang, Z. Y. Lin, and S. L. Zhu, Phys. Rev. D 105, 036016 (2022) doi: 10.1103/PhysRevD.105.036016
    [26] Y. J. Zhang, H. C. Chiang, P. N. Shen et al., Phys. Rev. D 74, 014013 (2006) doi: 10.1103/PhysRevD.74.014013
    [27] A. Faessler, T. Gutsche, V. E. Lyubovitskij et al., Phys. Rev. D 77, 114013 (2008) doi: 10.1103/PhysRevD.77.114013
    [28] G. Q. Feng, Z. X. Xie, and X. H. Guo, Phys. Rev. D 83, 016003 (2011) doi: 10.1103/PhysRevD.83.016003
    [29] V. M. Abazov et al. (D0), Phys. Rev. Lett. 117, 022003 (2016) doi: 10.1103/PhysRevLett.117.022003
    [30] C. J. Xiao and D. Y. Chen, Eur. Phys. J. A 53, 127 (2017) doi: 10.1140/epja/i2017-12310-x
    [31] S. S. Agaev, K. Azizi, and H. Sundu, Eur. Phys. J. Plus 131, 351 (2016) doi: 10.1140/epjp/i2016-16351-8
    [32] X. Chen and J. Ping, Eur. Phys. J. C 76, 351 (2016) doi: 10.1140/epjc/s10052-016-4210-x
    [33] R. Chen and X. Liu, Phys. Rev. D 94, 034006 (2016) doi: 10.1103/PhysRevD.94.034006
    [34] J. X. Lu, X. L. Ren, and L. S. Geng, Eur. Phys. J. C 77, 94 (2017) doi: 10.1140/epjc/s10052-017-4660-9
    [35] Z. Y. Wang, J. J. Qi, and X. H. Guo, Adv. High Energy Phys. 2019, 7576254 (2019) doi: 10.1155/2019/7576254
    [36] Y. R. Liu and Z. Y. Zhang, Phys. Rev. C 80, 015208 (2009) doi: 10.1103/PhysRevC.80.015208
    [37] S. Ohkoda, Y. Yamaguchi, S. Yasui et al., Phys. Rev. D 86, 034019 (2012) doi: 10.1103/PhysRevD.86.034019
    [38] Z. M. Ding, H. Y. Jiang, and J. He, Eur. Phys. J. C 80, 1179 (2020) doi: 10.1140/epjc/s10052-020-08754-6
    [39] D. Gamermann, E. Oset, D. Strottman et al., Phys. Rev. D 76, 074016 (2007) doi: 10.1103/PhysRevD.76.074016
    [40] L. R. Dai, J. J. Xie, and E. Oset, Eur. Phys. J. C 76, 121 (2016) doi: 10.1140/epjc/s10052-016-3946-7
    [41] L. Dai, G. Toledo, and E. Oset, Eur. Phys. J. C 80, 510 (2020) doi: 10.1140/epjc/s10052-020-8058-8
    [42] E. Wang, H. S. Li, W. H. Liang et al., Phys. Rev. D 103, 054008 (2021) doi: 10.1103/PhysRevD.103.054008
    [43] O. Deineka, I. Danilkin, and M. Vanderhaeghen, Phys. Lett. B 827, 136982 (2022) doi: 10.1016/j.physletb.2022.136982
    [44] F. Aceti, M. Bayar, J. M. Dias et al., Eur. Phys. J. A 50, 103 (2014) doi: 10.1140/epja/i2014-14103-1
    [45] F. Aceti, M. Bayar, E. Oset et al., Phys. Rev. D 90, 016003 (2014) doi: 10.1103/PhysRevD.90.016003
    [46] Z. M. Ding, H. Y. Jiang, D. Song et al., Eur. Phys. J. C 81, 732 (2021) doi: 10.1140/epjc/s10052-021-09534-6
    [47] T. Branz, T. Gutsche, and V. E. Lyubovitskij, Phys. Rev. D 78, 114004 (2008) doi: 10.1103/PhysRevD.78.114004
    [48] G. Li, F. l. Shao, C. W. Zhao et al., Phys. Rev. D 87, 034020 (2013) doi: 10.1103/PhysRevD.87.034020
    [49] M. J. Zhao, Z. Y. Wang, C. Wang et al., Phys. Rev. D 105, 096016 (2022) doi: 10.1103/PhysRevD.105.096016
    [50] Z. Y. Wang, J. J. Qi, J. Xu et al., Phys. Rev. D 102, 036008 (2020) doi: 10.1103/PhysRevD.102.036008
    [51] G. J. Ding, Phys. Rev. D 79, 014001 (2009) doi: 10.1103/PhysRevD.79.014001
    [52] F. K. Guo, C. Hanhart, U. G. Meißner et al., Rev. Mod. Phys. 90, 015004 (2018) doi: 10.1103/RevModPhys.90.015004
    [53] W. A. Bardeen, E. J. Eichten, and C. T. Hill, Phys. Rev. D 68, 054024 (2003) doi: 10.1103/PhysRevD.68.054024
    [54] C. W. Xiao and E. Oset, Eur. Phys. J. A 49, 52 (2013) doi: 10.1140/epja/i2013-13052-5
    [55] D. Gamermann, E. Oset, and B. S. Zou, Eur. Phys. J. A 41, 85-91 (2009) doi: 10.1140/epja/i2009-10767-8
    [56] L. L. Wei, H. S. Li, E. Wang et al., Phys. Rev. D 103, 114013 (2021) doi: 10.1103/PhysRevD.103.114013
    [57] S. Prelovsek, S. Collins, D. Mohler et al., JHEP 06, 035 (2021) doi: 10.1007/JHEP06(2021)035
    [58] X. Liu, Z. G. Luo, Y. R. Liu et al., Eur. Phys. J. C 61, 411-428 (2009) doi: 10.1140/epjc/s10052-009-1020-4
    [59] M. Z. Liu, D. J. Jia, and D. Y. Chen, Chin. Phys. C 41, 053105 (2017) doi: 10.1088/1674-1137/41/5/053105
    [60] S. Sakai, L. Roca, and E. Oset, Phys. Rev. D 96, 054023 (2017) doi: 10.1103/PhysRevD.96.054023
    [61] R. Aaij et al. (LHCb), Nature Commun. 13, 3351 (2022) doi: 10.1038/s41467-022-30206-w
    [62] N. Li, Z. F. Sun, X. Liu et al., Phys. Rev. D 88, 114008 (2013) doi: 10.1103/PhysRevD.88.114008
    [63] L. M. Abreu, arXiv: 2206.01166[hep-ph]
    [64] S. Ohkoda, Y. Yamaguchi, S. Yasui et al., Phys. Rev. D 86, 014004 (2012) doi: 10.1103/PhysRevD.86.014004
    [65] L. R. Dai, E. Oset, A. Feijoo et al., Phys. Rev. D 105, 074017 (2022) doi: 10.1103/PhysRevD.105.074017
    [66] J. He, Phys. Rev. D 92, 034004 (2015) doi: 10.1103/PhysRevD.92.034004
    [67] https://home.web.cern.ch/news/news/physics/lhcb-discovers-three-new-exotic-particles.
  • [1] M. Y. Han and Y. Nambu, Phys. Rev. 139, B1006-B1010 (1965) doi: 10.1103/PhysRev.139.B1006
    [2] P. A. Zyla et al. (Particle Data Group), PTEP 2020, 083C01 (2020) doi: 10.1093/ptep/ptaa104
    [3] A. Astier, L. Montanet, M. Baubillier et al., Phys. Lett. B 25, 294-297 (1967) doi: 10.1016/0370-2693(67)90020-2
    [4] R. Ammar, R. Davis, W. Kropac et al., Phys. Rev. Lett. 21, 1832-1835 (1968) doi: 10.1103/PhysRevLett.21.1832
    [5] B. Aubert et al. (BaBar), Phys. Rev. Lett. 90, 242001 (2003) doi: 10.1103/PhysRevLett.90.242001
    [6] S. Godfrey and K. Moats, Phys. Rev. D 93, 034035 (2016) doi: 10.1103/PhysRevD.93.034035
    [7] S. L. Olsen, T. Skwarnicki, and D. Zieminska, Rev. Mod. Phys. 90, 015003 (2018) doi: 10.1103/RevModPhys.90.015003
    [8] J. D. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 (1982) doi: 10.1103/PhysRevLett.48.659
    [9] J. D. Weinstein and N. Isgur, Phys. Rev. D 27, 588 (1983)
    [10] J. D. Weinstein and N. Isgur, Phys. Rev. D 41, 2236 (1990) doi: 10.1103/PhysRevD.41.2236
    [11] J. A. Oller, Nucl. Phys. A 714, 161-182 (2003) doi: 10.1016/S0375-9474(02)01360-X
    [12] V. Baru, J. Haidenbauer, C. Hanhart et al., Phys. Lett. B 586, 53-61 (2004) doi: 10.1016/j.physletb.2004.01.088
    [13] L. Y. Dai and M. R. Pennington, Phys. Rev. D 90, 036004 (2014) doi: 10.1103/PhysRevD.90.036004
    [14] H. A. Ahmed and C. W. Xiao, Phys. Rev. D 101, 094034 (2020) doi: 10.1103/PhysRevD.101.094034
    [15] Z. Q. Wang, X. W. Kang, J. A. Oller et al., Phys. Rev. D 105, 074016 (2022) doi: 10.1103/PhysRevD.105.074016
    [16] A. P. Szczepaniak, Phys. Lett. B 567, 23-26 (2003) doi: 10.1016/S0370-2693(03)00865-7
    [17] J. Hofmann and M. F. M. Lutz, Nucl. Phys. A 733, 142-152 (2004) doi: 10.1016/j.nuclphysa.2003.12.013
    [18] F. K. Guo, P. N. Shen, H. C. Chiang et al., Phys. Lett. B 641, 278-285 (2006) doi: 10.1016/j.physletb.2006.08.064
    [19] J. M. Flynn and J. Nieves, Phys. Rev. D 75, 074024 (2007) doi: 10.1103/PhysRevD.75.074024
    [20] Z. X. Xie, G. Q. Feng, and X. H. Guo, Phys. Rev. D 81, 036014 (2010) doi: 10.1103/PhysRevD.81.036014
    [21] Z. H. Guo, U. G. Meißner, and D. L. Yao, Phys. Rev. D 92, 094008 (2015) doi: 10.1103/PhysRevD.92.094008
    [22] M. L. Du, F. K. Guo, U. G. Meißner et al., Eur. Phys. J. C 77, 728 (2017) doi: 10.1140/epjc/s10052-017-5287-6
    [23] T. W. Wu, M. Z. Liu, L. S. Geng et al., Phys. Rev. D 100, 034029 (2019) doi: 10.1103/PhysRevD.100.034029
    [24] S. Y. Kong, J. T. Zhu, D. Song et al., Phys. Rev. D 104, 094012 (2021) doi: 10.1103/PhysRevD.104.094012
    [25] B. L. Huang, Z. Y. Lin, and S. L. Zhu, Phys. Rev. D 105, 036016 (2022) doi: 10.1103/PhysRevD.105.036016
    [26] Y. J. Zhang, H. C. Chiang, P. N. Shen et al., Phys. Rev. D 74, 014013 (2006) doi: 10.1103/PhysRevD.74.014013
    [27] A. Faessler, T. Gutsche, V. E. Lyubovitskij et al., Phys. Rev. D 77, 114013 (2008) doi: 10.1103/PhysRevD.77.114013
    [28] G. Q. Feng, Z. X. Xie, and X. H. Guo, Phys. Rev. D 83, 016003 (2011) doi: 10.1103/PhysRevD.83.016003
    [29] V. M. Abazov et al. (D0), Phys. Rev. Lett. 117, 022003 (2016) doi: 10.1103/PhysRevLett.117.022003
    [30] C. J. Xiao and D. Y. Chen, Eur. Phys. J. A 53, 127 (2017) doi: 10.1140/epja/i2017-12310-x
    [31] S. S. Agaev, K. Azizi, and H. Sundu, Eur. Phys. J. Plus 131, 351 (2016) doi: 10.1140/epjp/i2016-16351-8
    [32] X. Chen and J. Ping, Eur. Phys. J. C 76, 351 (2016) doi: 10.1140/epjc/s10052-016-4210-x
    [33] R. Chen and X. Liu, Phys. Rev. D 94, 034006 (2016) doi: 10.1103/PhysRevD.94.034006
    [34] J. X. Lu, X. L. Ren, and L. S. Geng, Eur. Phys. J. C 77, 94 (2017) doi: 10.1140/epjc/s10052-017-4660-9
    [35] Z. Y. Wang, J. J. Qi, and X. H. Guo, Adv. High Energy Phys. 2019, 7576254 (2019) doi: 10.1155/2019/7576254
    [36] Y. R. Liu and Z. Y. Zhang, Phys. Rev. C 80, 015208 (2009) doi: 10.1103/PhysRevC.80.015208
    [37] S. Ohkoda, Y. Yamaguchi, S. Yasui et al., Phys. Rev. D 86, 034019 (2012) doi: 10.1103/PhysRevD.86.034019
    [38] Z. M. Ding, H. Y. Jiang, and J. He, Eur. Phys. J. C 80, 1179 (2020) doi: 10.1140/epjc/s10052-020-08754-6
    [39] D. Gamermann, E. Oset, D. Strottman et al., Phys. Rev. D 76, 074016 (2007) doi: 10.1103/PhysRevD.76.074016
    [40] L. R. Dai, J. J. Xie, and E. Oset, Eur. Phys. J. C 76, 121 (2016) doi: 10.1140/epjc/s10052-016-3946-7
    [41] L. Dai, G. Toledo, and E. Oset, Eur. Phys. J. C 80, 510 (2020) doi: 10.1140/epjc/s10052-020-8058-8
    [42] E. Wang, H. S. Li, W. H. Liang et al., Phys. Rev. D 103, 054008 (2021) doi: 10.1103/PhysRevD.103.054008
    [43] O. Deineka, I. Danilkin, and M. Vanderhaeghen, Phys. Lett. B 827, 136982 (2022) doi: 10.1016/j.physletb.2022.136982
    [44] F. Aceti, M. Bayar, J. M. Dias et al., Eur. Phys. J. A 50, 103 (2014) doi: 10.1140/epja/i2014-14103-1
    [45] F. Aceti, M. Bayar, E. Oset et al., Phys. Rev. D 90, 016003 (2014) doi: 10.1103/PhysRevD.90.016003
    [46] Z. M. Ding, H. Y. Jiang, D. Song et al., Eur. Phys. J. C 81, 732 (2021) doi: 10.1140/epjc/s10052-021-09534-6
    [47] T. Branz, T. Gutsche, and V. E. Lyubovitskij, Phys. Rev. D 78, 114004 (2008) doi: 10.1103/PhysRevD.78.114004
    [48] G. Li, F. l. Shao, C. W. Zhao et al., Phys. Rev. D 87, 034020 (2013) doi: 10.1103/PhysRevD.87.034020
    [49] M. J. Zhao, Z. Y. Wang, C. Wang et al., Phys. Rev. D 105, 096016 (2022) doi: 10.1103/PhysRevD.105.096016
    [50] Z. Y. Wang, J. J. Qi, J. Xu et al., Phys. Rev. D 102, 036008 (2020) doi: 10.1103/PhysRevD.102.036008
    [51] G. J. Ding, Phys. Rev. D 79, 014001 (2009) doi: 10.1103/PhysRevD.79.014001
    [52] F. K. Guo, C. Hanhart, U. G. Meißner et al., Rev. Mod. Phys. 90, 015004 (2018) doi: 10.1103/RevModPhys.90.015004
    [53] W. A. Bardeen, E. J. Eichten, and C. T. Hill, Phys. Rev. D 68, 054024 (2003) doi: 10.1103/PhysRevD.68.054024
    [54] C. W. Xiao and E. Oset, Eur. Phys. J. A 49, 52 (2013) doi: 10.1140/epja/i2013-13052-5
    [55] D. Gamermann, E. Oset, and B. S. Zou, Eur. Phys. J. A 41, 85-91 (2009) doi: 10.1140/epja/i2009-10767-8
    [56] L. L. Wei, H. S. Li, E. Wang et al., Phys. Rev. D 103, 114013 (2021) doi: 10.1103/PhysRevD.103.114013
    [57] S. Prelovsek, S. Collins, D. Mohler et al., JHEP 06, 035 (2021) doi: 10.1007/JHEP06(2021)035
    [58] X. Liu, Z. G. Luo, Y. R. Liu et al., Eur. Phys. J. C 61, 411-428 (2009) doi: 10.1140/epjc/s10052-009-1020-4
    [59] M. Z. Liu, D. J. Jia, and D. Y. Chen, Chin. Phys. C 41, 053105 (2017) doi: 10.1088/1674-1137/41/5/053105
    [60] S. Sakai, L. Roca, and E. Oset, Phys. Rev. D 96, 054023 (2017) doi: 10.1103/PhysRevD.96.054023
    [61] R. Aaij et al. (LHCb), Nature Commun. 13, 3351 (2022) doi: 10.1038/s41467-022-30206-w
    [62] N. Li, Z. F. Sun, X. Liu et al., Phys. Rev. D 88, 114008 (2013) doi: 10.1103/PhysRevD.88.114008
    [63] L. M. Abreu, arXiv: 2206.01166[hep-ph]
    [64] S. Ohkoda, Y. Yamaguchi, S. Yasui et al., Phys. Rev. D 86, 014004 (2012) doi: 10.1103/PhysRevD.86.014004
    [65] L. R. Dai, E. Oset, A. Feijoo et al., Phys. Rev. D 105, 074017 (2022) doi: 10.1103/PhysRevD.105.074017
    [66] J. He, Phys. Rev. D 92, 034004 (2015) doi: 10.1103/PhysRevD.92.034004
    [67] https://home.web.cern.ch/news/news/physics/lhcb-discovers-three-new-exotic-particles.
  • 加载中

Figures(3) / Tables(1)

Get Citation
Zi-Long Li, Xing-Yue Tan, Zhu-Feng Zhang, Zhen-Yang Wang and Xin-Heng Guo. Investigating S-wave bound states composed of two pseudoscalar mesons[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac93ed
Zi-Long Li, Xing-Yue Tan, Zhu-Feng Zhang, Zhen-Yang Wang and Xin-Heng Guo. Investigating S-wave bound states composed of two pseudoscalar mesons[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac93ed shu
Milestone
Received: 2022-08-23
Article Metric

Article Views(1316)
PDF Downloads(37)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Investigating S-wave bound states composed of two pseudoscalar mesons

    Corresponding author: Zhen-Yang Wang, wangzhenyang@nbu.edu.cn
    Corresponding author: Xin-Heng Guo, xhguo@bnu.edu.cn
  • 1. Physics Department, Ningbo University, Zhejiang 315211, China
  • 2. College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China

Abstract: In this study, we systematically investigated two-pseudoscalar meson systems with the Bethe-Salpeter equation in the ladder and instantaneous approximations. By solving the Bethe-Salpeter equation numerically with the kernel containing the one-particle exchange diagrams, we found that the KˉK, DK, BˉK, DˉD, BˉB, BD, DˉK, BK, and BˉD systems with I=0 can exist as bound states. We also studied the contributions from heavy meson (J/ψ and Υ) exchanges and found that the contributions from heavy meson exchanges cannot be ignored.

    HTML

    I.   INTRODUCTION
    • Quantum chromodynamics (QCD) is the theory of the strong interaction between quarks mediated by gluons, which are color-charged [1]. In principle, QCD allows complex quark and gluon compositions of hadrons, such as multiquark hadrons, hadronic molecules, hybrid hadrons, and glueballs, which are nonstandard hadronic particles. Most of these nonstandard hadrons have unusual masses, decay widths, etc., which cannot be given a satisfactory explanation by the traditional quark model. To date, more than thirty non-qˉq state candidates in light and heavy sectors have been reported experimentally [2]. These resonances are crucial for the deep understanding of hadron spectroscopy and the nonperturbative nature and spontaneous chiral symmetry breaking of QCD.

      In the low-lying scalar meson sector, f0(980) [3], a0(980) [4], and Ds0(2317) [5] have nonexotic JP(=0+) quantum numbers. However, their masses are much lower than the quark model expectation for the corresponding P-wave qˉq states [6, 7]. Their natures are still under debate in spite of the efforts during the past several decades. Since their masses are near the threshold of the constituent particles and have spin-parity quantum numbers corresponding to the S-wave combinations of the constituent particles, one would naturally identify them as hadronic molecules, which are analogs of nuclei. f0(980) and/or a0(980) could be KˉK molecules [815] and Ds0(2317) could be a DK molecule [1626]. Such a picture leads to results consistent with the experiments. In addition to these particles, the possible S-wave bound states of the ˉBK, DˉD, BˉB, BD, ˉBˉK, DD, ˉBˉB, and ˉBD systems have not been observed experimentally.

      On the theoretical aspect, the authors in Ref. [26] systematically studied the possible S-wave bound state of two pseudoscalar mesons by the nonrelativistic Schr¨odinger (NRS) equation. Reference [18] predicted the existence of a BˉK bound state Bs0 with a mass of 5.725 ± 0.039 GeV based on the heavy chiral unitary approach. Subsequently, Refs. [27] and [28] confirmed the existence of Bs0 in the BˉK bound state scenario and further studied the decay widths of its possible decay channels. Recently, Kong et al. [24] systematically investigated DK/ˉBK and ˉDK/BK systems in a quasipotential Bethe-Salpter equation (qBSE) approach by considering the light meson exchange potential and found that only the isoscalar systems can exist as molecular states. However, the mass of X(5568) reported by the D0 collaboration [29] is considerably below the BK threshold, and hence, it cannot be aBK molecule [3035]. In Ref. [36], the authors reported their findings on the S-wave DˉD, BD, and BˉB systems in the chiral SU(3) quark model (QM); their calculations favor the existence of the isoscalar BˉB molecule, but the existence of isovector DˉD and BD molecules is disfavored. The authors of Ref. [37] studiedD()D() and B()B() molecular states by solving the coupled channel Schr¨odinger (CCS) equations; onlyI(JP)=1(0+)BBcan be a bound state in the PP (P=D,B) system because the kinetic term is suppressed in the bottom sector, and the effect of channel couplings becomes more important. With the qBSE approach [38], the existence of DˉD and BˉB molecular states with I(JP)=0(0+) were predicted, yet no bound state was produced from the DD and ˉBˉB interaction. In Ref. [39], a new hidden charm resonance with a mass of 3.7 GeV was predicted within the coupled channel unitary approach. Later, the DˉD bound state was searched in several processes, such as BDˉDK [40], ψ(3770)γD0ˉD0 [41], and γγDˉD [42, 43]. There are some differences in the results of different methods. Therefore, more efforts are needed to investigate the possible S-wave bound state composed of two pseudoscalar mesons.

      In the present study, we systematically investigate whether the S-wave bound states of two-pseudoscalar meson systems exist in the Bethe-Salpeter (BS) approach (in the ladder approximation and the instantaneous approximation for the kernel). For doubly heavy pseudoscalar meson systems, we consider not only the interaction through exchanged light mesons (ρ, ω and σ) but also the contribution of heavy vector mesons (J/ψ or Υ). As reported in Refs. [44, 45], in spite of the large mass ofJ/ψ, which suppresses the propagator of the exchanged J/ψ, it was found that the interaction could bind the DˉD and DˉD systems. Similarly, Refs. [38, 46] also reported that the contribution from a heavy meson exchange (J/ψ or Υ) is very important to form a molecular state, especially in systems in which the contributions from ρ and ω cancel each other out.

      The remainder of this paper is organized as follows. In Sec. II, we discuss and establish the BS equation for two-pseudoscalar meson systems. This equation is solved numerically, and the numerical results of the two-pseudoscalar meson systems are presented in Sec. III. In the last section, we summarize the study.

    II.   THE BETHE-SALPETER FORMALISM FOR THE TWO-PSEUDOSCALAR MESON SYSTEM
    • The BS wave function for the bound state |P composed of two pseudoscalar mesons have the following form:

      χ(x1,x2,P)=0|TP1(x1)P2(x2)|P,

      (1)

      where P1(x1) and P2(x2) are the field operators of the two constituent particles at space coordinates x1 and x2, respectively. The BS wave function in momentum space is defined as

      χP(x1,x2,P)=eiPXd4p(2π)4eipxχP(p),

      (2)

      where p represents the relative momentum of the two constituent particles, and p=λ2p1λ1p2 (p1=λ1P+p, p2=λ2Pp) with λ1=m1/(m1+m2) and λ2=m2/(m1+m2). p1(2) and m1(2) represent the momentum and mass of the constituent particle, respectively.

      The BS wave function χP(p) satisfies the following BS equation:

      χP(p)=SP1(p1)d4q(2π)4K(P,p,q)χP(q)SP2(p2),

      (3)

      where SP1(p1) and SP2(p2) are the propagators of constituent particles, and K(P,p,q) is the kernel, which is defined as the sum of all two-particle irreducible diagrams.

      In the following, we use the variables pl(=pv) and pt(=pplv) as the longitudinal and transverse projections of the relative momentum (p) along the bound state velocity (v), respectively. Then, the propagators of the constituent mesons can be expressed as

      SP1(λ1P+p)=i(λ1M+pl)2ω21+iϵ,

      (4)

      and

      SˉP2(λ2Pp)=i(λ2Mpl)2ω22+iϵ,

      (5)

      where ω1(2)=m21(2)+p2t (we have defined p2t=ptpt).

      To obtain the interaction kernel of the two-pseudoscalar meson system through exchanging light and heavy vector mesons and the light scalar meson, the following effective Lagrangians, as in Refs. [38, 47, 48], are needed:

      LKKV=igKKρρμ(KτμKμKτK)+i(gKKωωμ+gKKωϕμ)(KμKμKK),LDDV=igDDV(DbαDaDaαDb)Vαba+igDDJ/ψ(DαDDαD)J/ψα,LBBV=igBBV(BbαBaBaαBb)Vαba+igBBΥ(BαBBαB)Υα,LDDσ=gDDσDaDaσ,LBBσ=gBBσBaBaσ,

      (6)

      where J/ψα, Υα, and σ represent the J/ψ, Υ, and σ field operators, respectively, and the nonet vector meson matrix reads as

      V=(ρ02+ω2ρ+K+ρρ02+ω2K0KˉK0ϕ).

      (7)

      The coupling constants involved in Eq. (6) are taken as gKKρ=gKKω=gKKϕ=3, gDDV=gBBV=βgv/2 with gv=5.8, β=0.9, gDDJ/ϕ=mJ/ψ/fJ/ψ with fJ/ψ=405 MeV, and gBBΥ=mΥ/fΥ with fΥ=715.2 MeV.

      In the so-called ladder approximation, the interaction kernel K(P,p,q) can be derived in the lowest-order form as follows:

      K(p1,p2;q1,q1,mV)=(2π)2δ4(q1+q2p1p2)CIgPPVgPPV(p1+q1)μ(p2+q2)νΔμν(k,mV),K(p1,p2;q1,q1,mσ)=(2π)2δ4(q1+q2p1p2)CIg2PPσΔσ(k,mσ),

      (8)

      where mV represent the masses of the exchanged light and heavy vector mesons (ρ, ω, ψ, J/ψ, and Υ). Δμν(k,mV) and Δσ(k,mσ) represent the propagators for the vector and scalar mesons, respectively, and have the following forms:

      Δμν(k,mV)=ik2m2V(gμνkμkνm2V),Δσ(k,mσ)=ik2m2σ.

      (9)

      CI in Eq. (8) is the isospin coefficient for I=0 and I=1. For the KˉK, DK, ˉBK, DˉD, BˉB, and BD systems,

      C0={3/2forρ1/2forω1forϕ1forJ/ψ1forΥ1forσ,C1={1/2forρ1/2forω1forϕ1forJ/ψ1forΥ1forσ.

      (10)

      For the ˉKˉK, DˉK, ˉBˉK, DD, ˉBˉB, and ˉBD systems,

      C0={3/2forρ1/2forω1forϕ1forJ/ψ1forΥ1forσ,C1={1/2forρ1/2forω1forϕ1forJ/ψ1forΥ1forσ.

      (11)

      In Eqs. (10) and (11), the exchanged mesons of ϕ, J/ψ, and Υ only appear for the KˉK/ˉKˉK, DˉD/DD, and BˉB/ˉBˉB systems, and σ is only considered in the doubly heavy pseudoscalar meson systems.

      In order to manipulate the off shell effect of the exchanged mesons and finite size effect of the interacting hadrons, we introduce a form factor F(k2) at each vertex. Generally, the form factor has the following form:

      FM(k2)=Λ2m2Λ2k2,

      (12)

      where Λ, m, and k represent the cutoff parameter, mass, and momentum of the exchanged meson, respectively. This form factor is normalized at the on shell momentum of k2=m2. In contrast, if k2is taken to be infinitely large (), the form factor, which can be expressed as the overlap integral of the wave functions of the hadrons at the vertex, would approach zero. Considering the difference in the wave functions and masses of the light and heavy mesons, and ensuring a positive form factor, different magnitudes of cutoff Λ will be chosen for the heavy and light mesons.

      The propagators (4) and (5), interaction kernel (8), and form factor (12) are substituted into the BS equation, i.e., Eq. (3), and the instantaneous approximation (pl=ql, where the energy exchanged between the constituent particles of the binding system is neglected) in the kernel is considered. Consequently, when a vector meson and a scalar meson in the center-of-mass frame of the bound state (P=0) are exchanged, Eq. (3) becomes

      χP(pl,pt)=i[(λ1M+pl)2ω21+iϵ][(λ2Mpl)2ω22+iϵ]dql2πd3qt(2π)3

      ×{CIgPPVgPPV4(λ1M+pl)(λ2pl)+(pt+qt)2+(p2tq2t)/m2V(ptqt)m2V(Λ2m2V)2[Λ2+(pt+qt)2]2+CIgPPσgPPσ1(ptqt)m2σ(Λ2m2σ)2[Λ2+(pt+qt)2]2}χP(ql,qt).

      (13)

      In the above equation, there are poles in the pl plane at λ1Mω1+iϵ, λ1M+ω1iϵ, λ2M+ω2iϵ, and λ2Mω2+iϵ. After integrating pl on both sides of Eq. (13) by selecting the proper contour, we can obtain the three-dimensional integral equation for ˜χP(pt) (˜χP(pt)=dplχP(pl,pt)), which only depends on the the three momentum, pt. By completing the azimuthal integration, the three-dimensional BS equation becomes a one-dimensional integral equation as

      ˜χP(|pt|)=d|pt|A(|pt|,|qt|)˜χP(|qt|),

      (14)

      where the propagators and kernels after one-dimensional simplification are included in A(|pt|,|qt|). The numerical solutions for ˜χP(|pt|) can be obtained by discretizing the integration region into n pieces (with n sufficiently large). In this way, the integral equation becomes a matrix equation, and the BS scalar function ˜χP(|pt|) becomes an n dimensional vector.

    III.   NUMERICAL RESULTS
    • In this section, we will solve the BS equation numerically and study whether the S-wave bound states composed of two pseudoscalar mesons exist. In our model, there is only one parameter, the cutoff Λ, which comes from the form factor. The binding energy Eb is defined as Eb=Mm1m2 in the rest frame of the bound state. We take the averaged masses of the pseudoscalar mesons and the exchanged light and heavy mesons from PDG [2], mK = 494.988 MeV, mD = 1868.04 MeV, mB = 5279.44 MeV, mρ = 775.26 MeV, mω = 782.65 MeV, mϕ = 1019.461 MeV, mJ/ψ = 3096.9 MeV, and mΥ = 9460.3 MeV.

      In Fig. 1, we present some possible bound states composed of two pseudoscalar mesons when only the light meson (ρ, ω, ϕ, and σ) exchange contributions are considered. Here, we vary the binding energy from 0 to –50 MeV and the cutoff in a wide range of 0.8–5 GeV. We find that only the KˉK, DK, ˉBK, DˉD, BˉB, BD, DˉK, ˉBˉK, and ˉBD systems with I=0 can exist as bound states. The ˉKˉK, DD, and ˉBˉB systems with I=0 are forbidden because of Bose symmetry, and the interactions in I=1 systems are repulsive; hence, no bound states exist in the ˉKˉK, DD, and ˉBˉB systems. Furthermore, we cannot predict with certainty the masses of bound states that will be measured experimentally because our results are dependent on the cutoff Λ. The contribution of the σ exchange is included in our work, despite the large uncertainties in its mass and structure. In our previous works [49, 50] and Ref. [51], it was found that the contribution of the σ exchange is too small to form bound states, and the same result is found in our current work.

      Figure 1.  (color online) The numerical results for KˉK, DK, ˉBK, DˉD, BˉB, and BD systems with I=0 (a) and DˉK, ˉBˉK, and ˉBD systems with I=0 (b).

      The systems that may exist as bound states are presented in Fig. 1. It is noted that, in the hidden bottom system, the cutoff is the smallest in Fig. 1(a). This is because the B meson has the largest mass, which requires a smaller cutoff value as compared with those in the other systems, and the cutoff is determined by the overlap integrals of the wave functions and hadrons at the vertices. The size of the B meson is the smallest among the constituent particles. There is no bound state for the system with I=1. This is because the isospin coefficients of ρ and ω are –1/2 and 1/2, respectively, as shown in Eq. (10), and the masses of ρ and ω are almost equal. Therefore, contributions from the ρ and ω exchanges almost cancel each other out. Among these possible bound states, the KˉK and DK bound states can be related to the experimentally observed f0(980) and Ds0(2317), respectively [52]. Based on the heavy chiral unitary approach [18] and linear chiral symmetry [53], the authors predicted the existence of a b-partner state Bs0 of Ds0(2317) as the BˉK bound state, which can also be confirmed in our model with the cutoff Λ = 2436 MeV. The experimentally observed X(5568) cannot be a BˉK bound state in our model [35]. In Ref. [39], a new hidden charm resonance with a mass of 3.7 GeV (named as X(3700)) was predicted corresponding mostly to a DˉD state. Later, it has been searched in BDˉDK [40], e+eJ/ψDˉD [54], ψ(3770)X(3700)γ [55], γγDˉD [42, 43], ΛbΛDˉD [56], etc. Recently, lattice QCD also found a DˉD bound state just below the threshold with the binding energy Eb=4.0+3.75.0 MeV [57]. The existence of the BˉB bound state was also confirmed by the effective potential model [58], the heavy quark effective theory [59], the chiral SU(3) QM [36], and the qBSE [38]. The BD system, analogous to the DK system, can also be a bound state in the local hidden gauge symmetry (HGS) approach [60].

      Recently, the Tcc with the quantum number I(JP)=0(1+) and quark content ccˉuˉd was reported by the LHCb Collaboration [61], which is the first experimentally discovered open charmed tetraquark state. The Tcc has a mass just below the threshold of DD and can be an ideal candidate for the DD bound state. In fact, this inspired us to investigate the possibility of two-pseudoscalar meson systems as bound states with open flavor. In Ref. [62], the authors systematically investigated possible deuteron-like molecular states with two heavy quarks using the one-boson-exchange (OBE) model. According to their results, the I=1 DD system might not be a molecule; the I=1 ˉBˉB, I=0, and I=1 DˉB systems might be molecule candidates, but the results are sensitive to the cutoff. Based on the Heavy-Meson Effective Theory, the systems of DD with I=1, ˉBˉB with I=1, and DˉB with I=0 and I=1 can exist as shallow bound states [63].

      According to the results of the two-pseudoscalar meson systems, which are presented in Fig. 1(b), only the isoscalar system can exist as bound states. This is because, for the isovector systems, the isospin coefficients corresponding to the ρ and ω exchanges are –1/2, as shown in Eq. (11); hence, the total interaction is repulsive in the isovector systems. Comparing the results in Fig. 1 (a) and Fig. 1(b), we can see that the cutoff Λ is larger in Fig. 1(b), which is caused by the difference in the isospin coefficients, i.e., 1/2 and –1/2 in the systems with I=0 due to the ω exchange in Eq. (10) and Eq. (11), respectively. From Fig. 1(a) and Fig. 1(b), we can also find that, for the constituent particles with same masses, the larger the mass of the constituent particle, the smaller the cutoff Λ; moreover, the larger the difference in masses of the constituent particles, the larger the cutoff Λ. To easily compare the results of different theoretical models, we list the results of some models and our results in Table 1.

      KˉKDKˉBKDˉDBˉBBDˉKˉKDˉKˉBˉKDDˉBˉBˉBD
      I010101010101010101010101
      qBSE [24, 38]××××××××××
      HGS [60, 64]××××
      OBE [58, 62]××××××
      CCS [37, 65]××××
      QM [36]×××
      NRS [26]×××××××××××××××××
      Our results×××××××××××××××

      Table 1.  The results for different theoretical models. and × indicate whether the corresponding system is a bound state or not, respectively. indicates the corresponding system without any study.

      In Ref. [45], the authors systematically studied the interaction of DˉD in the isospin I=0 channel. In their work, it is shown that the exchange of a light qˉq is OZI forbidden in the I=1 channel. As a consequence, only the J/ψ exchange is allowed in the case of I=1, and the simultaneous two pion exchange, which was evaluated in Ref. [45,44], was found to be weaker than the exchange of the vector meson. In spite of the large mass of the J/ψ, which suppresses the propagator of the exchanged J/ψ, it was found in Refs. [45,44] that the interaction could bind the DˉD and DˉD systems with I=1 weakly. Subsequently, in Ref. [66], it was also found that the bound state DˉD with IG(JP)=1+(1+) would disappear if the J/ψ exchange was removed, which means that the J/ψ exchange is important to provide an attractive interaction and produce the pole in the isovector system. In the present work, we also consider the effects of exchanged heavy mesons. Considering the differences in the wave functions and masses of the light and heavy mesons, and ensuring a positive from factor, we choose different magnitudes of cutoffs ΛL and ΛH for the exchanged light and heavy mesons, respectively.

      We vary the cutoff ΛL of the exchanged light mesons in the range of 800–1500 MeV to find the cutoff ΛH of the exchanged heavy mesons that can form bound states. The results for some possible bound states of DˉD and BˉB are presented in Figs. 2 and 3. From these results in Figs. 2 and 3, we can see that the effect of the exchange of a heavy meson cannot be ignored. As the contributions of the exchanged ρ and ω almost cancel each other out in the I=1 DˉD and BˉB systems, the main contribution comes from the heavy meson exchange. It can be seen that the results for ΛH are almost twice the mass of the exchanged heavy meson in Figs. 2(b) and 3(b). This is similar to the case of considering only the light meson exchange, in which the value of the cutoff Λ is also approximately twice the mass of the exchange meson, as in Fig. 1(a). The possibility of bound states being formed when only considering the contribution of the heavy meson exchange, along with the existence of ˉDDand BˉB bound states with I=1, is yet to be confirmed experimentally.

      Figure 2.  (color online) The numerical results for I=0 (a) and I=1 (b) DˉD systems with J/ψ meson exchange included.

      Figure 3.  (color online) The numerical results for I=0 (a) and I=1 (b) BˉB systems with Υ meson exchange included.

    IV.   SUMMARY
    • In this study, we derived the BS equation for the S-wave KˉK, DK, BˉK, DˉD, BˉB, BD, KK, DˉK, ˉBˉK, DD, ˉBˉB, and ˉBD systems and systematically investigated the possible bound states of these systems with the ladder approximation and instantaneous approximation for the kernel. In our model, the kernel containing one-particle-exchange diagrams was induced by the light meson (ρ, ω, ϕ, and σ) and heavy meson (J/ψ and Υ) exchanges. To investigate the bound states, we have numerically solved the BS equations for S-wave systems composed of two pseudoscalar mesons. The possible S-wave bound states investigated in our study are helpful in explaining the structures of experimentally discovered exotic states and predicting unobserved exotic states.

      We found that the KˉK, DK, ˉBK, DˉD, BˉB, BD, DˉK, ˉBˉK, and ˉBDsystems with I=0 can exist as bound states. The ˉKˉK, DD, and ˉBˉB systems with I=0 are forbidden because of Bose symmetry, and the interactions in the I=1 systems are repulsive; hence, no bound states exist in the ˉKˉK, DD, and ˉBˉB systems. We also found that, for the constituent particles with the same mass, the larger the mass of the constituent particle, the smaller the cutoff Λ; moreover, the larger the difference in the masses of the constituent particles, the larger the cutoff Λ. The contribution of the σ exchange is too small to form bound states.

      In the calculation, we considered the heavy meson exchanges in the kernel. We found that the effect of the heavy meson exchange cannot be neglected for the DˉD and BˉB systems. Since the contributions from the ρ and ω exchanges almost cancel each other out in the I=1 DˉD and BˉB systems, the main contribution comes from the heavy meson exchanges, and the I=1 DˉD and BˉB systems can exist as bound states. However, since the cutoff ΛH for the heavy meson exchanges is very large, the possibility of bound states being formed when only considering the contribution of the heavy meson exchange, along with the existence of ˉDDand BˉB bound states with I=1, is yet to be confirmed experimentally.

      With the restarted LHC and other experiments, more experimental studies of exotic hadrons will be performed in the near future. Recently, the LHCb collaboration observed three never-before-seen particles: a new kind of pentaquark and the first-ever pair of tetraquarks, which includes a new type of tetraquark [67]. These will help physicists better understand how quarks bind together into exotic particles. The theoretical explanation of the structures of experimentally observed exotic hadrons and the existence of possible molecular states predicted theoretically remain controversial. Therefore, more precise experimental studies of the exotic states will be needed to test the results of theoretical studies and improve theoretical models.

    ACKNOWLEDGMENTS
    • One of the authors (Z.-Y. Wang) thanks Professor Jia-Jun Wu and Dr. Rui-Cheng Li for helpful discussions and useful suggestions.

Reference (67)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return