Correlation function and the inverse problem in the BD interaction

  • We study the correlation functions of the B0D+,B+D0 system, which develops a bound state of approximately 40MeV, using inputs consistent with the Tcc(3875) state. Then, we address the inverse problem starting from these correlation functions to determine the scattering observables related to the system, including the existence of the bound state and its molecular nature. The important output of the approach is the uncertainty with which these observables can be obtained, considering errors in the B0D+,B+D0 correlation functions typical of current values in correlation functions. We find that it is possible to obtain scattering lengths and effective ranges with relatively high precision and the existence of a bound state. Although the pole position is obtained with errors of the order of 50% of the binding energy, the molecular probability of the state is obtained with a very small error of the order of 6%. All these findings serve as motivation to perform such measurements in future runs of high energy hadron collisions.
  • Femtoscopic correlation functions are emerging as a tool for understanding the interactions of hadrons at small relative momenta. Considerable experimental research has already been conducted within the strange sector [114] (see also review paper [15]). Incursions in the D sector have been investigated [16], and plans to extend studies to this sector are expected in the near future [17]. In future runs of the LHC, the ALICE collaboration will likely also access the bottom sector.

    The theoretical community is also devoting considerable research to the subject [1833], and a model independent analysis of the correlation functions was very recently proposed [3438], where instead of contrasting theoretical models with experimental data, the inverse path was followed. The data were used to determine the scattering observables of coupled channels, explore the possibility of having several bound states, and eventually determine the nature of these bound states as molecular states or otherwise.

    In the present study, we address the BD interaction, where a bound state was predicted for the BD system with isospin I=0 in Ref. [39]. The observation of the Tcc state [40, 41] and subsequent theoretical research support this state as a bound state of the D0D+, D+D0 channels in I=0 [4253]. Two recent studies investigating Tcc from a different perspective also concluded that the Tcc state is indeed a molecular state [54, 55]. In the first study, the scattering length and effective range of the D0D+, D+D0 channels, as well as the shape of the D0D0π+ mass distributions, were theoretically studied [54], concluding that Tcc is a molecular state of the D0D+, D+D0 components, essentially constructing an I=0 state. In the second study, a different path was followed, assuming that Tcc could correspond to a non molecular, genuine, or preexisting state and be dressed by the meson-meson components where it is observed [55]. It was found that, in principle, it would be possible to have a Tcc state with negligible molecular probability, but at the heavy cost of a small D0D+ scattering length and large effective range, which are far from the observed experimental values.

    If Tcc is bound, using arguments of heavy quark symmetry, the ˉBD state should also be bound, even more than Tcc because there is a general rule that the heavier the quarks, the stronger the interaction and binding energies [5658]. This rule is also satisfied when the extension of the local hidden gauge approach [5962] is used as a source of interaction, exchanging vector mesons between the heavy meson components. In Ref. [39], different B()D(), BˉD() pairs were found to be bound using the same regulator for the loops as in Ref. [45], particularly the BD state in I=0, which was found to be bound by approximately 1530MeV. The existence of one bound state in the BD system is also supported by the phase moment obained in Ref. [63].

    The finding of the Tcc state has been extremely useful for placing constraints on the meson-meson interaction and its range, as reflected in the cutoff used to regularize the loops [64, 65]. The information and arguments on heavy quark symmetry used in Ref. [39] consolidated their results. With confidence in the predictions of Ref. [39] and their inputs consistent with the information obtained from Tcc, we present here one work which should stimulate experimental measurements to corroborate these findings. First, we reproduce the results of Ref. [39] and evaluate the correlation functions of D0B+ and D+B0. Next, assuming that these correlation functions correspond to actual data, we address the inverse problem of determining from them the value of the scattering observables, scattering length, and effective range for these two channels; the existence of a bound state below the threshold; the molecular probability of each of the D0B+ and D+B0 components; and most importantly, the precision with which we can determine these magnitudes, assuming reasonable error bars for the correlation function data. This information is important to obtain an idea of what to expect given the experimental constraints when such experiments are performed.

    The details can be found in Ref. [39]; hence, we do not cover them here. Light vector mesons are exchanged between the B and D components, as shown in Fig. 1. We can see that heavy quarks are spectators in the exchange of light vector mesons. Furthermore, assuming mρ=mω, the interaction for the B+D0(1), B0D+(2) channels is given by

    Figure 1

    Figure 1.  Vector exchange between the D0B+ and D+B0 components (in brackets, the four momentum of each particle).

    Vij=14f2Cij(p1+p2)(p3+p4);f=93 MeV,

    (1)

    with Cij the element of matrix C,

    Cij=(1111),

    (2)

    and projected over S-wave

    (p1+p2)(p3+p4)12[3s2(m2B+m2D)(m2Bm2D)2s].

    (3)

    The scattering matrix is then given by

    T=[1VG]1V,

    (4)

    where G=diag(G1,G2), with Gi being the BD meson loop function regularized with a cutoff

    Gi(s)=|q|<qmaxd3q(2π)3ω1(q)+ω2(q)2ω1(q)ω2(q)×1s[ω1(q)+ω2(q)]2+iε,

    (5)

    where ω1(q)=q2+m2i, ω2(q)=q2+M2i, and mi,Mi are the masses of the D and B mesons in channel i. With the isospin states (B+,B0), (D+,D0), the I=0 combination is given by

    |BD,I=0=12|B+D0+B0D+,

    (6)

    with which we can obtain

    BD,I=0|V|BD,I=0=12(V11+V22+2V12)=2V11,

    (7)

    which is the result obtained in Ref. [39], with Vij given by Eqs. (1), (2), and (3). The I=1 combination is given by

    |BD,I=1,I3=0=12|B+D0B0D+.

    (8)

    We reproduce the results of Ref. [39] using the cutoff regularization with an extra form factor (Eq. (27) from Ref. [39]) stemming from the S-wave projection of a vector meson exchange. However, practically identical results are obtained by ignoring this form factor and decreasing qmax. To match the formalism of the correlation functions of Ref. [30], we ignore this form factor and take qmax of the order of 420MeV, as used in Ref. [45], to obtain the binding of the related Tcc state.

    We can evaluate the scattering length a, r0 for the B+D0 and B0D+ channels by recalling the relationship between our T matrix and that used in quantum mechanics [64]

    T=8πsfQM8πs11a+12r0k2ik,

    (9)

    with

    k=λ1/2(s,m21,m22)2s,

    (10)

    from which we easily find

    1ai=(8πsT1ii)|sth,i,

    (11)

    r0,i=[2sμis(8πsT1ii+iki)]sth,i,

    (12)

    where μi is the reduced mass in channel i, and sth,i and ki are the square of the threshold mass and the center of mass momenta of the mesons for channel i, respectively.

    We find that there is a pole, which is below the threshold of the two channels and hence corresponds to a bound state. The couplings are obtained from the T matrix in the vicinity of the pole,

    Tij=gigjssp,

    (13)

    where sp is the square of the mass at the pole. Thus,

    g21=limssp(ssp)T11,

    (14)

    g1gj=limssp(ssp)T1j,

    (15)

    which determine the relative sign of g2 with respect to g1. Once the couplings are evaluated, we calculate the molecular probabilities of the B+D0 and B0D+ channels, as in Refs. [64, 66, 67].

    Pi=g2iGis|s=sp.

    (16)

    Another magnitude of relevance is the wave function at the origin in coordinate space, given by [64]

    ψi(r=0)=giGi|sth,i.

    (17)

    We follow the formalism of Ref. [30] and write the correlation functions for the two channels as

    CB+D0(pD0)=1+4πθ(qmaxpD0)drr2S12(r)×{|j0(pD0r)+TB+D0,B+D0(E)˜G(B+D0)(r;E)|2+|TB0D+,B+D0(E)˜G(B0D+)(r;E)|2j20(pD0r)},

    (18)

    CB0D+(pD+)=1+4πθ(qmaxpD+)drr2S12(r)×{|j0(pD+r)+TB0D+,B0D+(E)˜G(B0D+)(r;E)|2+|TB+D0,B0D+(E)˜G(B+D0)(r;E)|2j20(pD+r)},

    (19)

    where pi is the momentum of the particles in the rest frame of the pair,

    pi=λ1/2(s,m2i,M2i)2s,

    (20)

    S12(r) is the source function, parameterized as a Gaussian normalized to 1,

    S12(r)=1(4πR)3e(r2/4R2),

    (21)

    and the ˜G(i)(r;E) function is defined as

    ˜G(i)(r;E)=d3q(2π)3ω1(q)+ω2(q)2ω1(q)ω2(q)j0(qr)s[ω1(q)+ω2(q)]2+iε,

    (22)

    where j0(qr) is the spherical Bessel function, and E=s=m2i+p2i+M2i+p2i.

    Here, we assume that the correlation functions have already been measured and attempt to extract the maximum information available using a general framework in which no model assumptions are made. To perform the test, we use the correlation function with the model described in the previous sections, assuming errors at the order of ±0.02, which are slightly larger than those obtained in current measurements of correlation functions.

    We begin by assuming that there is an interaction between the coupled channels, given by

    V=(V11V12V12V22),

    (23)

    where Vij are unknown potentials to be determined, and the T matrix is given by Eq. (4), using the G function of Eq. (5) with an unknown qmax. We make no assumption on the isospin structure of a possible bound state but assume that the interaction is isospin symmetric, which implies, according to Eqs. (6), (7), and (8), that

    BD,I=0|V|BD,I=1=12(V11V22)=0,

    (24)

    V22=V11.

    (25)

    In addition, to consider possible sources of interaction originating from channels neglected in our approach, we introduce several energy dependent terms, as discussed in Refs. [68, 69] and studied in [3436]:

    V11=V11+αm2V(ssth,1),

    (26)

    V12=V12+βm2V(ssth,1),

    (27)

    where the factor m2V (with mV=800MeV) is introduced to make α,β dimensionless. Then, we have 4 free parameters for the interaction, as well as qmax and R, which is a total of 6 parameters to fit the two correlation functions. We must be aware that there are strong correlations between these parameters because the input used to obtain the correlation functions corresponds to an interaction at I=0; hence, what matters is the combination V11+V12. This means that the values we obtain for the parameters in fits to the pseudodata are not meaningful, and only the values of the observables obtained from them are significant. To manage with these correlations, we use the bootstrap or resampling method [7072], generating random centroids of the data with a Gaussian distribution weight and performing a large number of fits to the data with the new centroids and same errors. After each fit, the values of the observables are evaluated, and the average and dispersion for each are calculated.

    Next, we use the model in Section II with a cutoff regularization of qmax=420MeV , as in the study of the Tcc state in Ref. [45]. We obtain a pole at s=7110.41 MeV and the couplings given in Table 1. Similarly, the probabilities obtained and the wave functions at the origin are given in Table 2, and the scattering length and effective range are shown in Table 3. As shown, the probability obtained for the sum of the two B+D0,B0D+ channels is of the order of 96%. The small deviation from unity is due to the energy dependence of the original potential of Eqs. (1) and (3). We also observe that while a1,r0,1 are real, a2,r0,2 are complex because the B+D0 channel is open at the threshold of the B0D+ channel. The couplings are very similar, and so are the wave functions at the origin, indicating an I=0 state, according to Eq. (6). The results for the correlation functions of the two channels are shown in Fig. 2, calculated with R=1fm.

    Table 1

    Table 1.  Pole position and couplings with qmax=420MeV. [in units of MeV]
    sp g1 g2
    (7110.41+0i) 31636.8 31631.0
    DownLoad: CSV
    Show Table

    Table 2

    Table 2.  Probability Pi and wave function at the origin ψi(r=0) for channel i.
    P1 P2 ψ1(r=0) ψ2(r=0)
    0.52 0.44 14.75 13.61
    DownLoad: CSV
    Show Table

    Table 3

    Table 3.  Scattering length ai and effective range ri for channel i. [in units of fm]
    a1 a2 r1 r2
    0.71 0.500.16i 0.61 1.221.77i
    DownLoad: CSV
    Show Table

    Figure 2

    Figure 2.  (color online) Correlation functions for the B+D0 and B0D+ channels.

    Next, we discuss the results obtained from the resampling method fits to the data. The data with the assumed errors are shown in Fig. 3. A warning should be given: as discussed in detail in Ref. [36], the sign of V12 is undefined in the present procedure; however, we rely on arguments of heavy quark flavor symmetry to choose solutions with V11 and V12 of the same sign.

    Figure 3

    Figure 3.  (color online) Correlation functions for the B+D0 (left) and B0D+ (right) channels, with 26 points in each curve with an error of ±0.02. The centroids of the red data follow the theoretical curve. In blue, we plot the centroids obtained in one of the resampling runs, with a random Gaussian generation of the centroids of each point.

    In Table 4, we show the average values and dispersion of the parameters obtained. As we discussed above, they are not significative owing to the existing correlations. One indication of these correlations is the relatively large parameter errors. Nevertheless, the important aspect is the value of the observables. These can be found in Tables 5, 6, and 7.

    Table 4

    Table 4.  Values obtained for the parameters V11,V12, α,β, qmax , and R.
    V11 V12 α β qmax/MeV R/fm
    1537.54±918.93 1512.32±913.55 31.89±10.00 165.16±133.31 407.9±60.8 1.01±0.04
    DownLoad: CSV
    Show Table

    Table 5

    Table 5.  Average values and dispersion of the pole position and couplings. [in units of MeV]
    sp g1 g2
    7107.84±17.79 34623.08±14300.15 34506.64±14304.57
    DownLoad: CSV
    Show Table

    Table 6

    Table 6.  Average value and dispersion of the probability Pi and wave function at the origin ψi(r=0) for channel i.
    P1 P2 ψ1(r=0) ψ2(r=0)
    0.49±0.03 0.42±0.03 13.61±2.65 12.61±2.423
    DownLoad: CSV
    Show Table

    Table 7

    Table 7.  Average value and dispersion of the scattering length ai and effective range ri for channel i. [in units of fm]
    a1 a2 r1 r2
    0.72±0.03 (0.51±0.02)(0.17±0.01)i 0.61±0.19 (1.41±0.28)(1.65±0.07)i
    DownLoad: CSV
    Show Table

    In Table 5, we show the value of the energy at which the bound state is found, together with the values of the couplings. As shown, a bound state is found around 7108MeV, compatible with the bound state obtained with the original model within uncertainties. Interestingly, the error obtained is of the order of 18MeV. This is not a small error for a binding energy of 39MeV; however, this is what can be achieved with the assumed precision of the correlation data. More positively, using the data of the correlations at the BD threshold, we are still able to predict that there is a bound state with an approximately 40MeV binding.

    The couplings g1,g2 obtained are also compatible with the original ones, and the errors are also not small. Yet, the approximately equal values of the couplings suggest that we are dealing with an I=0 state.

    It is interesting to analyze our obtained probabilities of the states, which are shown in Table 6. We again obtain numbers for the probabilities of the two channels compatible with those obtained from the original model, but once again, it is the precision by which they can be obtained that is important. The uncertainties are very small, that is, of the order of 6%. This may be surprising in view of the formula used to obtain Pi (Eq. (16)), which is proportional to g2i, and gi has large errors according to Table 5. If g2i is bigger in a fit because the state is more bound, Gi/s also decreases in strength and the product becomes more stable. This is an interesting and fortunate result of our analysis, which allows us to conclude that the application of the inverse method from the femtoscopic correlation functions would allow us to determine the nature of the bound state obtained with a high accuracy.

    As shown in Table 7, it is also rewarding to find that we can determine the scattering lengths with good precision and the effective ranges with smaller precision but significant values.

    We should also stress that the inverse method allows us to obtain the size of the source function with a relative accuracy of approximately 4% (see Table 4).

    We address the problem of evaluating the correlation functions for the BD system using inputs extracted from a successful study of the Tcc(3875) state. We take two channels, B0D+ and B+D0, the small mass differences between which induce visible differences in the correlation functions. The system also develops a bound state of approximately 40 MeV. Once this is achieved, we address the inverse problem of obtaining the observables associated with this system, starting from the correlation functions, assuming errors as in current measurements. Although we obtain results compatible with those obtained from the original model, an important new result is the uncertainty by which we can obtain the observables of the system. We obtain the size of the source with a precision of approximately 4%. We also determine that there is a bound state using the data of the correlation functions above the threshold of the channels, although with an uncertainty of approximately 50% of the binding energy. However, remarkably, we can determine the molecular nature of the obtained state with a good precision of approximately 6%. The scattering lengths of the two channels are obtained with good precision and significant values, and the effective ranges are also obtained but with smaller precision. All these results indicate that the measurement of these correlation functions in the future will allow us to obtain valuable information on the BD interaction and the bound states associated with this interaction. This study could serve as motivation to conduct such measurements in the future.

    [1] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 114, 022301 (2015), arXiv:1408.4360[nucl-ex] doi: 10.1103/PhysRevLett.114.022301
    [2] S. Acharya et al. (ALICE), Phys. Lett. B 774, 64 (2017), arXiv:1705.04929[nucl-ex] doi: 10.1016/j.physletb.2017.09.009
    [3] S. Acharya et al. (ALICE), Phys. Rev. C 99, 024001 (2019), arXiv:1805.12455[nucl-ex] doi: 10.1103/PhysRevC.99.024001
    [4] J. Adam et al. (STAR), Phys. Lett. B 790, 490 (2019), arXiv:1808.02511[hep-ex] doi: 10.1016/j.physletb.2019.01.055
    [5] S. Acharya et al. (ALICE), Phys. Rev. Lett. 123, 112002 (2019), arXiv:1904.12198[nucl-ex] doi: 10.1103/PhysRevLett.123.112002
    [6] S. Acharya et al. (ALICE), Phys. Lett. B 797, 134822 (2019), arXiv:1905.07209[nucl-ex] doi: 10.1016/j.physletb.2019.134822
    [7] S. Acharya et al. (ALICE), Phys. Rev. Lett. 124, 092301 (2020), arXiv:1905.13470[nucl-ex] doi: 10.1103/PhysRevLett.124.092301
    [8] S. Acharya et al. (ALICE), Phys. Lett. B 805, 135419 (2020), arXiv:1910.14407[nucl-ex] doi: 10.1016/j.physletb.2020.135419
    [9] A. Collaboration et al. (ALICE), Nature 588, 232 (2020), arXiv:2005.11495[nucl-ex] doi: 10.1038/s41586-020-3001-6
    [10] S. Acharya et al. (ALICE), Phys. Rev. Lett. 127, 172301 (2021), arXiv:2105.05578[nucl-ex] doi: 10.1103/PhysRevLett.127.172301
    [11] S. Acharya et al. (ALICE), Phys. Lett. B 822, 136708 (2021), arXiv:2105.05683[nucl-ex] doi: 10.1016/j.physletb.2021.136708
    [12] S. Acharya et al. (ALICE), Phys. Lett. B 790, 22 (2019), arXiv:1809.07899[nucl-ex] doi: 10.1016/j.physletb.2018.12.033
    [13] S. Acharya et al. (ALICE), Phys. Lett. B 833, 137335 (2022), arXiv:2111.06611[nucl-ex] doi: 10.1016/j.physletb.2022.137335
    [14] S. Acharya et al. (ALICE), Eur. Phys. J. C 83, 340 (2023), arXiv:2205.15176[nucl-ex] doi: 10.1140/epjc/s10052-023-11476-0
    [15] L. Fabbietti, V. Mantovani Sarti, and O. Vazquez Doce, Ann. Rev. Nucl. Part. Sci. 71, 377 (2021), arXiv:2012.09806[nucl-ex] doi: 10.1146/annurev-nucl-102419-034438
    [16] S. Acharya et al. (ALICE), Phys. Rev. D 106, 052010 (2022), arXiv:2201.05352[nucl-ex] doi: 10.1103/PhysRevD.106.052010
    [17] S. Acharya et al. (ALICE), (2022), arXiv: 2211.02491[physics.ins-det]
    [18] K. Morita, T. Furumoto, and A. Ohnishi, Phys. Rev. C 91, 024916 (2015), arXiv:1408.6682[nucl-th] doi: 10.1103/PhysRevC.91.024916
    [19] A. Ohnishi, K. Morita, K. Miyahara et al., Nucl. Phys. A 954, 294 (2016), arXiv:1603.05761[nucl-th] doi: 10.1016/j.nuclphysa.2016.05.010
    [20] V. M. Sarti, A. Feijoo, I. Vidaña et al., (2023), arXiv: 2309.08756[hep-ph]
    [21] K. Morita, A. Ohnishi, F. Etminan et al., Phys. Rev. C 94, 031901 (2016), arXiv:1605.06765[hep-ph] doi: 10.1103/PhysRevC.94.031901
    [22] T. Hatsuda, K. Morita, A. Ohnishi et al., Nucl. Phys. A 967, 856 (2017), arXiv:1704.05225[nucl-th] doi: 10.1016/j.nuclphysa.2017.04.041
    [23] D. L. Mihaylov, V. Mantovani Sarti, O. W. Arnold et al., Eur. Phys. J. C 78, 394 (2018), arXiv:1802.08481[hep-ph] doi: 10.1140/epjc/s10052-018-5859-0
    [24] J. Haidenbauer, Nucl. Phys. A 981, 1 (2019), arXiv:1808.05049[hep-ph] doi: 10.1016/j.nuclphysa.2018.10.090
    [25] K. Morita, S. Gongyo, T. Hatsuda et al., Phys. Rev. C 101, 015201 (2020), arXiv:1908.05414[nucl-th] doi: 10.1103/PhysRevC.101.015201
    [26] Y. Kamiya, T. Hyodo, K. Morita et al., Phys. Rev. Lett. 124, 132501 (2020), arXiv:1911.01041[nucl-th] doi: 10.1103/PhysRevLett.124.132501
    [27] Y. Kamiya, K. Sasaki, T. Fukui et al., Phys. Rev. C 105, 014915 (2022), arXiv:2108.09644[hep-ph] doi: 10.1103/PhysRevC.105.014915
    [28] Y. Kamiya, T. Hyodo, and A. Ohnishi, Eur. Phys. J. A 58, 131 (2022), arXiv:2203.13814[hep-ph] doi: 10.1140/epja/s10050-022-00782-y
    [29] Z.-W. Liu, J.-X. Lu, and L.-S. Geng, Phys. Rev. D 107, 074019 (2023), arXiv:2302.01046[hep-ph] doi: 10.1103/PhysRevD.107.074019
    [30] I. Vidana, A. Feijoo, M. Albaladejo et al., Phys. Lett. B 846, 138201 (2023), arXiv:2303.06079[hep-ph] doi: 10.1016/j.physletb.2023.138201
    [31] M. Albaladejo, J. Nieves, and E. Ruiz-Arriola, Phys. Rev. D 108, 014020 (2023), arXiv:2304.03107[hep-ph] doi: 10.1103/PhysRevD.108.014020
    [32] Z.-W. Liu, J.-X. Lu, M.-Z. Liu et al., Phys. Rev. D 108, L031503 (2023), arXiv:2305.19048[hep-ph] doi: 10.1103/PhysRevD.108.L031503
    [33] Z.-W. Liu, K.-W. Li, and L.-S. Geng, Chin. Phys. C 47, 024108 (2023), arXiv:2201.04997[hep-ph] doi: 10.1088/1674-1137/ac988a
    [34] N. Ikeno, G. Toledo, and E. Oset, Phys. Lett. B 847, 138281 (2023), arXiv:2305.16431[hep-ph] doi: 10.1016/j.physletb.2023.138281
    [35] M. Albaladejo, A. Feijoo, I. Vidaña et al., (2023), arXiv: 2307.09873[hep-ph]
    [36] A. Feijoo, L. R. Dai, L. M. Abreu et al., Phys. Rev. D 109, 016014 (2024), arXiv:2309.00444 [hep-ph] doi: 10.1103/PhysRevD.109.016014
    [37] R. Molina, C.-W. Xiao, W.-H. Liang et al., Phys. Rev. D 109, 054002 (2024), arXiv:2310.12593 [hep-ph] doi: 10.1103/PhysRevD.109.054002
    [38] R. Molina, Z.-W. Liu, L.-S. Geng et al., (2023), arXiv: 2312.11993[hep-ph]
    [39] S. Sakai, L. Roca, and E. Oset, Phys. Rev. D 96, 054023 (2017), arXiv:1704.02196[hep-ph] doi: 10.1103/PhysRevD.96.054023
    [40] R. Aaij et al. (LHCb), Nature Phys. 18, 751 (2022), arXiv:2109.01038[hep-ex] doi: 10.1038/s41567-022-01614-y
    [41] R. Aaij et al. (LHCb), Nature Commun. 13, 3351 (2022), arXiv:2109.01056[hep-ex] doi: 10.1038/s41467-022-30206-w
    [42] N. Li, Z.-F. Sun, X. Liu et al., Phys. Rev. D 88, 114008 (2013), arXiv:1211.5007[hep-ph] doi: 10.1103/PhysRevD.88.114008
    [43] M.-Z. Liu, T.-W. Wu, M. Pavon Valderrama et al., Phys. Rev. D 99, 094018 (2019), arXiv:1902.03044[hep-ph] doi: 10.1103/PhysRevD.99.094018
    [44] Z.-M. Ding, H.-Y. Jiang, and J. He, Eur. Phys. J. C 80, 1179 (2020), arXiv:2011.04980[hep-ph] doi: 10.1140/epjc/s10052-020-08754-6
    [45] A. Feijoo, W. H. Liang, and E. Oset, Phys. Rev. D 104, 114015 (2021), arXiv:2108.02730[hep-ph] doi: 10.1103/PhysRevD.104.114015
    [46] X.-K. Dong, F.-K. Guo, and B.-S. Zou, Commun. Theor. Phys. 73, 125201 (2021), arXiv:2108.02673[hep-ph] doi: 10.1088/1572-9494/ac27a2
    [47] K. Chen, R. Chen, L. Meng et al., Eur. Phys. J. C 82, 581 (2022), arXiv:2109.13057[hep-ph] doi: 10.1140/epjc/s10052-022-10540-5
    [48] M. Albaladejo, Phys. Lett. B 829, 137052 (2022), arXiv:2110.02944[hep-ph] doi: 10.1016/j.physletb.2022.137052
    [49] M.-L. Du, V. Baru, X.-K. Dong et al., Phys. Rev. D 105, 014024 (2022), arXiv:2110.13765[hep-ph] doi: 10.1103/PhysRevD.105.014024
    [50] S. S. Agaev, K. Azizi, and H. Sundu, JHEP 06, 057 (2022), arXiv:2201.02788[hep-ph] doi: 10.1007/JHEP06(2022)057
    [51] H.-X. Chen, W. Chen, X. Liu et al., Rept. Prog. Phys. 86, 026201 (2023), arXiv:2204.02649[hep-ph] doi: 10.1088/1361-6633/aca3b6
    [52] G.-J. Wang, Z. Yang, J.-J. Wu et al., (2023), arXiv: 2306.12406[hep-ph]
    [53] C. Chen, C. Meng, Z.-G. et al., Chin. Phys. C 83, 043102 (2024), arXiv:2307.12069 [hep-ph] doi: 10.1088/1674-1137/ad1a0b
    [54] L. R. Dai, L. M. Abreu, A. Feijoo et al., Eur. Phys. J. C 83, 983 (2023), arXiv:2304.01870[hep-ph] doi: 10.1140/epjc/s10052-023-12159-6
    [55] L. R. Dai, J. Song, and E. Oset, Phys. Lett. B 846, 138200 (2023), arXiv:2306.01607[hep-ph] doi: 10.1016/j.physletb.2023.138200
    [56] J. P. Ader, J. M. Richard, and P. Taxil, Phys. Rev. D 25, 2370 (1982) doi: 10.1103/PhysRevD.25.2370
    [57] S. Zouzou, B. Silvestre-Brac, C. Gignoux et al., Z. Phys. C 30, 457 (1986) doi: 10.1007/BF01557611
    [58] J. Carlson, L. Heller, and J. A. Tjon, Phys. Rev. D 37, 744 (1988) doi: 10.1103/PhysRevD.37.744
    [59] M. Bando, T. Kugo, S. Uehara et al., Phys. Rev. Lett. 54, 1215 (1985) doi: 10.1103/PhysRevLett.54.1215
    [60] S. Furui, R. Kobayashi, and M. Nakagawa, Nuovo Cim. A 108, 241 (1995) doi: 10.1007/BF02816745
    [61] U. G. Meissner, Phys. Rept. 161, 213 (1988) doi: 10.1016/0370-1573(88)90090-7
    [62] H. Nagahiro, L. Roca, A. Hosaka et al., Phys. Rev. D 79, 014015 (2009), arXiv:0809.0943[hep-ph] doi: 10.1103/PhysRevD.79.014015
    [63] D.-L. Yao, P. Fernandez-Soler, F.-K. Guo et al., Phys. Rev. D 101, 034014 (2020), arXiv:1906.00727[hep-ph] doi: 10.1103/PhysRevD.101.034014
    [64] D. Gamermann, J. Nieves, E. Oset et al., Phys. Rev. D 81, 014029 (2010), arXiv:0911.4407[hep-ph] doi: 10.1103/PhysRevD.81.014029
    [65] J. Song, L. R. Dai, and E. Oset, Eur. Phys. J. A 58, 133 (2022), arXiv:2201.04414[hep-ph] doi: 10.1140/epja/s10050-022-00753-3
    [66] T. Hyodo, D. Jido, and A. Hosaka, Phys. Rev. C 85, 015201 (2012), arXiv:1108.5524[nucl-th] doi: 10.1103/PhysRevC.85.015201
    [67] T. Sekihara, T. Hyodo, and D. Jido, PTEP 2015, 063D04 (2015), arXiv:1411.2308[hep-ph] doi: 10.1093/ptep/ptv081
    [68] F. Aceti, L. R. Dai, L. S. Geng et al., Eur. Phys. J. A 50, 57 (2014), arXiv:1301.2554[hep-ph] doi: 10.1140/epja/i2014-14057-2
    [69] T. Hyodo, Int. J. Mod. Phys. A 28, 1330045 (2013), arXiv:1310.1176[hep-ph] doi: 10.1142/S0217751X13300457
    [70] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, (1992)
    [71] B. Efron and R. Tibshirani, Statist. Sci. 57, 54 (1986) doi: 10.1214/ss/1177013815
    [72] M. Albaladejo, D Jido, J. Nieves et al., Eur. Phys. J. C 76, 300 (2016), arXiv:1604.01193[hep-ph] doi: 10.1140/epjc/s10052-016-4144-3
  • [1] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 114, 022301 (2015), arXiv:1408.4360[nucl-ex] doi: 10.1103/PhysRevLett.114.022301
    [2] S. Acharya et al. (ALICE), Phys. Lett. B 774, 64 (2017), arXiv:1705.04929[nucl-ex] doi: 10.1016/j.physletb.2017.09.009
    [3] S. Acharya et al. (ALICE), Phys. Rev. C 99, 024001 (2019), arXiv:1805.12455[nucl-ex] doi: 10.1103/PhysRevC.99.024001
    [4] J. Adam et al. (STAR), Phys. Lett. B 790, 490 (2019), arXiv:1808.02511[hep-ex] doi: 10.1016/j.physletb.2019.01.055
    [5] S. Acharya et al. (ALICE), Phys. Rev. Lett. 123, 112002 (2019), arXiv:1904.12198[nucl-ex] doi: 10.1103/PhysRevLett.123.112002
    [6] S. Acharya et al. (ALICE), Phys. Lett. B 797, 134822 (2019), arXiv:1905.07209[nucl-ex] doi: 10.1016/j.physletb.2019.134822
    [7] S. Acharya et al. (ALICE), Phys. Rev. Lett. 124, 092301 (2020), arXiv:1905.13470[nucl-ex] doi: 10.1103/PhysRevLett.124.092301
    [8] S. Acharya et al. (ALICE), Phys. Lett. B 805, 135419 (2020), arXiv:1910.14407[nucl-ex] doi: 10.1016/j.physletb.2020.135419
    [9] A. Collaboration et al. (ALICE), Nature 588, 232 (2020), arXiv:2005.11495[nucl-ex] doi: 10.1038/s41586-020-3001-6
    [10] S. Acharya et al. (ALICE), Phys. Rev. Lett. 127, 172301 (2021), arXiv:2105.05578[nucl-ex] doi: 10.1103/PhysRevLett.127.172301
    [11] S. Acharya et al. (ALICE), Phys. Lett. B 822, 136708 (2021), arXiv:2105.05683[nucl-ex] doi: 10.1016/j.physletb.2021.136708
    [12] S. Acharya et al. (ALICE), Phys. Lett. B 790, 22 (2019), arXiv:1809.07899[nucl-ex] doi: 10.1016/j.physletb.2018.12.033
    [13] S. Acharya et al. (ALICE), Phys. Lett. B 833, 137335 (2022), arXiv:2111.06611[nucl-ex] doi: 10.1016/j.physletb.2022.137335
    [14] S. Acharya et al. (ALICE), Eur. Phys. J. C 83, 340 (2023), arXiv:2205.15176[nucl-ex] doi: 10.1140/epjc/s10052-023-11476-0
    [15] L. Fabbietti, V. Mantovani Sarti, and O. Vazquez Doce, Ann. Rev. Nucl. Part. Sci. 71, 377 (2021), arXiv:2012.09806[nucl-ex] doi: 10.1146/annurev-nucl-102419-034438
    [16] S. Acharya et al. (ALICE), Phys. Rev. D 106, 052010 (2022), arXiv:2201.05352[nucl-ex] doi: 10.1103/PhysRevD.106.052010
    [17] S. Acharya et al. (ALICE), (2022), arXiv: 2211.02491[physics.ins-det]
    [18] K. Morita, T. Furumoto, and A. Ohnishi, Phys. Rev. C 91, 024916 (2015), arXiv:1408.6682[nucl-th] doi: 10.1103/PhysRevC.91.024916
    [19] A. Ohnishi, K. Morita, K. Miyahara et al., Nucl. Phys. A 954, 294 (2016), arXiv:1603.05761[nucl-th] doi: 10.1016/j.nuclphysa.2016.05.010
    [20] V. M. Sarti, A. Feijoo, I. Vidaña et al., (2023), arXiv: 2309.08756[hep-ph]
    [21] K. Morita, A. Ohnishi, F. Etminan et al., Phys. Rev. C 94, 031901 (2016), arXiv:1605.06765[hep-ph] doi: 10.1103/PhysRevC.94.031901
    [22] T. Hatsuda, K. Morita, A. Ohnishi et al., Nucl. Phys. A 967, 856 (2017), arXiv:1704.05225[nucl-th] doi: 10.1016/j.nuclphysa.2017.04.041
    [23] D. L. Mihaylov, V. Mantovani Sarti, O. W. Arnold et al., Eur. Phys. J. C 78, 394 (2018), arXiv:1802.08481[hep-ph] doi: 10.1140/epjc/s10052-018-5859-0
    [24] J. Haidenbauer, Nucl. Phys. A 981, 1 (2019), arXiv:1808.05049[hep-ph] doi: 10.1016/j.nuclphysa.2018.10.090
    [25] K. Morita, S. Gongyo, T. Hatsuda et al., Phys. Rev. C 101, 015201 (2020), arXiv:1908.05414[nucl-th] doi: 10.1103/PhysRevC.101.015201
    [26] Y. Kamiya, T. Hyodo, K. Morita et al., Phys. Rev. Lett. 124, 132501 (2020), arXiv:1911.01041[nucl-th] doi: 10.1103/PhysRevLett.124.132501
    [27] Y. Kamiya, K. Sasaki, T. Fukui et al., Phys. Rev. C 105, 014915 (2022), arXiv:2108.09644[hep-ph] doi: 10.1103/PhysRevC.105.014915
    [28] Y. Kamiya, T. Hyodo, and A. Ohnishi, Eur. Phys. J. A 58, 131 (2022), arXiv:2203.13814[hep-ph] doi: 10.1140/epja/s10050-022-00782-y
    [29] Z.-W. Liu, J.-X. Lu, and L.-S. Geng, Phys. Rev. D 107, 074019 (2023), arXiv:2302.01046[hep-ph] doi: 10.1103/PhysRevD.107.074019
    [30] I. Vidana, A. Feijoo, M. Albaladejo et al., Phys. Lett. B 846, 138201 (2023), arXiv:2303.06079[hep-ph] doi: 10.1016/j.physletb.2023.138201
    [31] M. Albaladejo, J. Nieves, and E. Ruiz-Arriola, Phys. Rev. D 108, 014020 (2023), arXiv:2304.03107[hep-ph] doi: 10.1103/PhysRevD.108.014020
    [32] Z.-W. Liu, J.-X. Lu, M.-Z. Liu et al., Phys. Rev. D 108, L031503 (2023), arXiv:2305.19048[hep-ph] doi: 10.1103/PhysRevD.108.L031503
    [33] Z.-W. Liu, K.-W. Li, and L.-S. Geng, Chin. Phys. C 47, 024108 (2023), arXiv:2201.04997[hep-ph] doi: 10.1088/1674-1137/ac988a
    [34] N. Ikeno, G. Toledo, and E. Oset, Phys. Lett. B 847, 138281 (2023), arXiv:2305.16431[hep-ph] doi: 10.1016/j.physletb.2023.138281
    [35] M. Albaladejo, A. Feijoo, I. Vidaña et al., (2023), arXiv: 2307.09873[hep-ph]
    [36] A. Feijoo, L. R. Dai, L. M. Abreu et al., Phys. Rev. D 109, 016014 (2024), arXiv:2309.00444 [hep-ph] doi: 10.1103/PhysRevD.109.016014
    [37] R. Molina, C.-W. Xiao, W.-H. Liang et al., Phys. Rev. D 109, 054002 (2024), arXiv:2310.12593 [hep-ph] doi: 10.1103/PhysRevD.109.054002
    [38] R. Molina, Z.-W. Liu, L.-S. Geng et al., (2023), arXiv: 2312.11993[hep-ph]
    [39] S. Sakai, L. Roca, and E. Oset, Phys. Rev. D 96, 054023 (2017), arXiv:1704.02196[hep-ph] doi: 10.1103/PhysRevD.96.054023
    [40] R. Aaij et al. (LHCb), Nature Phys. 18, 751 (2022), arXiv:2109.01038[hep-ex] doi: 10.1038/s41567-022-01614-y
    [41] R. Aaij et al. (LHCb), Nature Commun. 13, 3351 (2022), arXiv:2109.01056[hep-ex] doi: 10.1038/s41467-022-30206-w
    [42] N. Li, Z.-F. Sun, X. Liu et al., Phys. Rev. D 88, 114008 (2013), arXiv:1211.5007[hep-ph] doi: 10.1103/PhysRevD.88.114008
    [43] M.-Z. Liu, T.-W. Wu, M. Pavon Valderrama et al., Phys. Rev. D 99, 094018 (2019), arXiv:1902.03044[hep-ph] doi: 10.1103/PhysRevD.99.094018
    [44] Z.-M. Ding, H.-Y. Jiang, and J. He, Eur. Phys. J. C 80, 1179 (2020), arXiv:2011.04980[hep-ph] doi: 10.1140/epjc/s10052-020-08754-6
    [45] A. Feijoo, W. H. Liang, and E. Oset, Phys. Rev. D 104, 114015 (2021), arXiv:2108.02730[hep-ph] doi: 10.1103/PhysRevD.104.114015
    [46] X.-K. Dong, F.-K. Guo, and B.-S. Zou, Commun. Theor. Phys. 73, 125201 (2021), arXiv:2108.02673[hep-ph] doi: 10.1088/1572-9494/ac27a2
    [47] K. Chen, R. Chen, L. Meng et al., Eur. Phys. J. C 82, 581 (2022), arXiv:2109.13057[hep-ph] doi: 10.1140/epjc/s10052-022-10540-5
    [48] M. Albaladejo, Phys. Lett. B 829, 137052 (2022), arXiv:2110.02944[hep-ph] doi: 10.1016/j.physletb.2022.137052
    [49] M.-L. Du, V. Baru, X.-K. Dong et al., Phys. Rev. D 105, 014024 (2022), arXiv:2110.13765[hep-ph] doi: 10.1103/PhysRevD.105.014024
    [50] S. S. Agaev, K. Azizi, and H. Sundu, JHEP 06, 057 (2022), arXiv:2201.02788[hep-ph] doi: 10.1007/JHEP06(2022)057
    [51] H.-X. Chen, W. Chen, X. Liu et al., Rept. Prog. Phys. 86, 026201 (2023), arXiv:2204.02649[hep-ph] doi: 10.1088/1361-6633/aca3b6
    [52] G.-J. Wang, Z. Yang, J.-J. Wu et al., (2023), arXiv: 2306.12406[hep-ph]
    [53] C. Chen, C. Meng, Z.-G. et al., Chin. Phys. C 83, 043102 (2024), arXiv:2307.12069 [hep-ph] doi: 10.1088/1674-1137/ad1a0b
    [54] L. R. Dai, L. M. Abreu, A. Feijoo et al., Eur. Phys. J. C 83, 983 (2023), arXiv:2304.01870[hep-ph] doi: 10.1140/epjc/s10052-023-12159-6
    [55] L. R. Dai, J. Song, and E. Oset, Phys. Lett. B 846, 138200 (2023), arXiv:2306.01607[hep-ph] doi: 10.1016/j.physletb.2023.138200
    [56] J. P. Ader, J. M. Richard, and P. Taxil, Phys. Rev. D 25, 2370 (1982) doi: 10.1103/PhysRevD.25.2370
    [57] S. Zouzou, B. Silvestre-Brac, C. Gignoux et al., Z. Phys. C 30, 457 (1986) doi: 10.1007/BF01557611
    [58] J. Carlson, L. Heller, and J. A. Tjon, Phys. Rev. D 37, 744 (1988) doi: 10.1103/PhysRevD.37.744
    [59] M. Bando, T. Kugo, S. Uehara et al., Phys. Rev. Lett. 54, 1215 (1985) doi: 10.1103/PhysRevLett.54.1215
    [60] S. Furui, R. Kobayashi, and M. Nakagawa, Nuovo Cim. A 108, 241 (1995) doi: 10.1007/BF02816745
    [61] U. G. Meissner, Phys. Rept. 161, 213 (1988) doi: 10.1016/0370-1573(88)90090-7
    [62] H. Nagahiro, L. Roca, A. Hosaka et al., Phys. Rev. D 79, 014015 (2009), arXiv:0809.0943[hep-ph] doi: 10.1103/PhysRevD.79.014015
    [63] D.-L. Yao, P. Fernandez-Soler, F.-K. Guo et al., Phys. Rev. D 101, 034014 (2020), arXiv:1906.00727[hep-ph] doi: 10.1103/PhysRevD.101.034014
    [64] D. Gamermann, J. Nieves, E. Oset et al., Phys. Rev. D 81, 014029 (2010), arXiv:0911.4407[hep-ph] doi: 10.1103/PhysRevD.81.014029
    [65] J. Song, L. R. Dai, and E. Oset, Eur. Phys. J. A 58, 133 (2022), arXiv:2201.04414[hep-ph] doi: 10.1140/epja/s10050-022-00753-3
    [66] T. Hyodo, D. Jido, and A. Hosaka, Phys. Rev. C 85, 015201 (2012), arXiv:1108.5524[nucl-th] doi: 10.1103/PhysRevC.85.015201
    [67] T. Sekihara, T. Hyodo, and D. Jido, PTEP 2015, 063D04 (2015), arXiv:1411.2308[hep-ph] doi: 10.1093/ptep/ptv081
    [68] F. Aceti, L. R. Dai, L. S. Geng et al., Eur. Phys. J. A 50, 57 (2014), arXiv:1301.2554[hep-ph] doi: 10.1140/epja/i2014-14057-2
    [69] T. Hyodo, Int. J. Mod. Phys. A 28, 1330045 (2013), arXiv:1310.1176[hep-ph] doi: 10.1142/S0217751X13300457
    [70] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, (1992)
    [71] B. Efron and R. Tibshirani, Statist. Sci. 57, 54 (1986) doi: 10.1214/ss/1177013815
    [72] M. Albaladejo, D Jido, J. Nieves et al., Eur. Phys. J. C 76, 300 (2016), arXiv:1604.01193[hep-ph] doi: 10.1140/epjc/s10052-016-4144-3
  • 加载中

Figures(3) / Tables(7)

Get Citation
Hai-Peng Li, Jing-Yu Yi, Chu-Wen Xiao, De-Liang Yao, Wei-Hong Liang and Eulogio Oset. Correlation function and the inverse problem in the BD interaction[J]. Chinese Physics C. doi: 10.1088/1674-1137/ad2dc2
Hai-Peng Li, Jing-Yu Yi, Chu-Wen Xiao, De-Liang Yao, Wei-Hong Liang and Eulogio Oset. Correlation function and the inverse problem in the BD interaction[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ad2dc2 shu
Milestone
Received: 2024-02-07
Article Metric

Article Views(1258)
PDF Downloads(9)
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:

Correlation function and the inverse problem in the BD interaction

  • 1. Department of Physics, Guangxi Normal University, Guilin 541004, China
  • 2. School of Physics and Electronics, Hunan University, Changsha 410082, China
  • 3. Guangxi Key Laboratory of Nuclear Physics and Technology, Guangxi Normal University, Guilin 541004, China
  • 4. School of Physics, Central South University, Changsha 410083, China
  • 5. Departamento de Física Teórica and IFIC, Centro Mixto Universidad de Valencia-CSIC Institutos de Investigación de Paterna, Aptdo. 22085, 46071 Valencia, Spain

Abstract: We study the correlation functions of the B0D+,B+D0 system, which develops a bound state of approximately 40MeV, using inputs consistent with the Tcc(3875) state. Then, we address the inverse problem starting from these correlation functions to determine the scattering observables related to the system, including the existence of the bound state and its molecular nature. The important output of the approach is the uncertainty with which these observables can be obtained, considering errors in the B0D+,B+D0 correlation functions typical of current values in correlation functions. We find that it is possible to obtain scattering lengths and effective ranges with relatively high precision and the existence of a bound state. Although the pole position is obtained with errors of the order of 50% of the binding energy, the molecular probability of the state is obtained with a very small error of the order of 6%. All these findings serve as motivation to perform such measurements in future runs of high energy hadron collisions.

    HTML

    I.   INTRODUCTION
    • Femtoscopic correlation functions are emerging as a tool for understanding the interactions of hadrons at small relative momenta. Considerable experimental research has already been conducted within the strange sector [114] (see also review paper [15]). Incursions in the D sector have been investigated [16], and plans to extend studies to this sector are expected in the near future [17]. In future runs of the LHC, the ALICE collaboration will likely also access the bottom sector.

      The theoretical community is also devoting considerable research to the subject [1833], and a model independent analysis of the correlation functions was very recently proposed [3438], where instead of contrasting theoretical models with experimental data, the inverse path was followed. The data were used to determine the scattering observables of coupled channels, explore the possibility of having several bound states, and eventually determine the nature of these bound states as molecular states or otherwise.

      In the present study, we address the BD interaction, where a bound state was predicted for the BD system with isospin I=0 in Ref. [39]. The observation of the Tcc state [40, 41] and subsequent theoretical research support this state as a bound state of the D0D+, D+D0 channels in I=0 [4253]. Two recent studies investigating Tcc from a different perspective also concluded that the Tcc state is indeed a molecular state [54, 55]. In the first study, the scattering length and effective range of the D0D+, D+D0 channels, as well as the shape of the D0D0π+ mass distributions, were theoretically studied [54], concluding that Tcc is a molecular state of the D0D+, D+D0 components, essentially constructing an I=0 state. In the second study, a different path was followed, assuming that Tcc could correspond to a non molecular, genuine, or preexisting state and be dressed by the meson-meson components where it is observed [55]. It was found that, in principle, it would be possible to have a Tcc state with negligible molecular probability, but at the heavy cost of a small D0D+ scattering length and large effective range, which are far from the observed experimental values.

      If Tcc is bound, using arguments of heavy quark symmetry, the ˉBD state should also be bound, even more than Tcc because there is a general rule that the heavier the quarks, the stronger the interaction and binding energies [5658]. This rule is also satisfied when the extension of the local hidden gauge approach [5962] is used as a source of interaction, exchanging vector mesons between the heavy meson components. In Ref. [39], different B()D(), BˉD() pairs were found to be bound using the same regulator for the loops as in Ref. [45], particularly the BD state in I=0, which was found to be bound by approximately 1530MeV. The existence of one bound state in the BD system is also supported by the phase moment obained in Ref. [63].

      The finding of the Tcc state has been extremely useful for placing constraints on the meson-meson interaction and its range, as reflected in the cutoff used to regularize the loops [64, 65]. The information and arguments on heavy quark symmetry used in Ref. [39] consolidated their results. With confidence in the predictions of Ref. [39] and their inputs consistent with the information obtained from Tcc, we present here one work which should stimulate experimental measurements to corroborate these findings. First, we reproduce the results of Ref. [39] and evaluate the correlation functions of D0B+ and D+B0. Next, assuming that these correlation functions correspond to actual data, we address the inverse problem of determining from them the value of the scattering observables, scattering length, and effective range for these two channels; the existence of a bound state below the threshold; the molecular probability of each of the D0B+ and D+B0 components; and most importantly, the precision with which we can determine these magnitudes, assuming reasonable error bars for the correlation function data. This information is important to obtain an idea of what to expect given the experimental constraints when such experiments are performed.

    II.   FORMALISM

      A.   Interaction of the D0B+, D+B0 channels

    • The details can be found in Ref. [39]; hence, we do not cover them here. Light vector mesons are exchanged between the B and D components, as shown in Fig. 1. We can see that heavy quarks are spectators in the exchange of light vector mesons. Furthermore, assuming mρ=mω, the interaction for the B+D0(1), B0D+(2) channels is given by

      Figure 1.  Vector exchange between the D0B+ and D+B0 components (in brackets, the four momentum of each particle).

      Vij=14f2Cij(p1+p2)(p3+p4);f=93 MeV,

      (1)

      with Cij the element of matrix C,

      Cij=(1111),

      (2)

      and projected over S-wave

      (p1+p2)(p3+p4)12[3s2(m2B+m2D)(m2Bm2D)2s].

      (3)

      The scattering matrix is then given by

      T=[1VG]1V,

      (4)

      where G=diag(G1,G2), with Gi being the BD meson loop function regularized with a cutoff

      Gi(s)=|q|<qmaxd3q(2π)3ω1(q)+ω2(q)2ω1(q)ω2(q)×1s[ω1(q)+ω2(q)]2+iε,

      (5)

      where ω1(q)=q2+m2i, ω2(q)=q2+M2i, and mi,Mi are the masses of the D and B mesons in channel i. With the isospin states (B+,B0), (D+,D0), the I=0 combination is given by

      |BD,I=0=12|B+D0+B0D+,

      (6)

      with which we can obtain

      BD,I=0|V|BD,I=0=12(V11+V22+2V12)=2V11,

      (7)

      which is the result obtained in Ref. [39], with Vij given by Eqs. (1), (2), and (3). The I=1 combination is given by

      |BD,I=1,I3=0=12|B+D0B0D+.

      (8)

      We reproduce the results of Ref. [39] using the cutoff regularization with an extra form factor (Eq. (27) from Ref. [39]) stemming from the S-wave projection of a vector meson exchange. However, practically identical results are obtained by ignoring this form factor and decreasing qmax. To match the formalism of the correlation functions of Ref. [30], we ignore this form factor and take qmax of the order of 420MeV, as used in Ref. [45], to obtain the binding of the related Tcc state.

    • B.   Scattering observables

    • We can evaluate the scattering length a, r0 for the B+D0 and B0D+ channels by recalling the relationship between our T matrix and that used in quantum mechanics [64]

      T=8πsfQM8πs11a+12r0k2ik,

      (9)

      with

      k=λ1/2(s,m21,m22)2s,

      (10)

      from which we easily find

      1ai=(8πsT1ii)|sth,i,

      (11)

      r0,i=[2sμis(8πsT1ii+iki)]sth,i,

      (12)

      where μi is the reduced mass in channel i, and sth,i and ki are the square of the threshold mass and the center of mass momenta of the mesons for channel i, respectively.

    • C.   Couplings and probabilities

    • We find that there is a pole, which is below the threshold of the two channels and hence corresponds to a bound state. The couplings are obtained from the T matrix in the vicinity of the pole,

      Tij=gigjssp,

      (13)

      where sp is the square of the mass at the pole. Thus,

      g21=limssp(ssp)T11,

      (14)

      g1gj=limssp(ssp)T1j,

      (15)

      which determine the relative sign of g2 with respect to g1. Once the couplings are evaluated, we calculate the molecular probabilities of the B+D0 and B0D+ channels, as in Refs. [64, 66, 67].

      Pi=g2iGis|s=sp.

      (16)

      Another magnitude of relevance is the wave function at the origin in coordinate space, given by [64]

      ψi(r=0)=giGi|sth,i.

      (17)
    • D.   Correlation functions

    • We follow the formalism of Ref. [30] and write the correlation functions for the two channels as

      CB+D0(pD0)=1+4πθ(qmaxpD0)drr2S12(r)×{|j0(pD0r)+TB+D0,B+D0(E)˜G(B+D0)(r;E)|2+|TB0D+,B+D0(E)˜G(B0D+)(r;E)|2j20(pD0r)},

      (18)

      CB0D+(pD+)=1+4πθ(qmaxpD+)drr2S12(r)×{|j0(pD+r)+TB0D+,B0D+(E)˜G(B0D+)(r;E)|2+|TB+D0,B0D+(E)˜G(B+D0)(r;E)|2j20(pD+r)},

      (19)

      where pi is the momentum of the particles in the rest frame of the pair,

      pi=λ1/2(s,m2i,M2i)2s,

      (20)

      S12(r) is the source function, parameterized as a Gaussian normalized to 1,

      S12(r)=1(4πR)3e(r2/4R2),

      (21)

      and the ˜G(i)(r;E) function is defined as

      ˜G(i)(r;E)=d3q(2π)3ω1(q)+ω2(q)2ω1(q)ω2(q)j0(qr)s[ω1(q)+ω2(q)]2+iε,

      (22)

      where j0(qr) is the spherical Bessel function, and E=s=m2i+p2i+M2i+p2i.

    III.   INVERSE PROBLEM: MODEL INDEPENDENT ANALYSIS OF THE CORRELATION FUNCTIONS
    • Here, we assume that the correlation functions have already been measured and attempt to extract the maximum information available using a general framework in which no model assumptions are made. To perform the test, we use the correlation function with the model described in the previous sections, assuming errors at the order of ±0.02, which are slightly larger than those obtained in current measurements of correlation functions.

      We begin by assuming that there is an interaction between the coupled channels, given by

      V=(V11V12V12V22),

      (23)

      where Vij are unknown potentials to be determined, and the T matrix is given by Eq. (4), using the G function of Eq. (5) with an unknown qmax. We make no assumption on the isospin structure of a possible bound state but assume that the interaction is isospin symmetric, which implies, according to Eqs. (6), (7), and (8), that

      BD,I=0|V|BD,I=1=12(V11V22)=0,

      (24)

      V22=V11.

      (25)

      In addition, to consider possible sources of interaction originating from channels neglected in our approach, we introduce several energy dependent terms, as discussed in Refs. [68, 69] and studied in [3436]:

      V11=V11+αm2V(ssth,1),

      (26)

      V12=V12+βm2V(ssth,1),

      (27)

      where the factor m2V (with mV=800MeV) is introduced to make α,β dimensionless. Then, we have 4 free parameters for the interaction, as well as qmax and R, which is a total of 6 parameters to fit the two correlation functions. We must be aware that there are strong correlations between these parameters because the input used to obtain the correlation functions corresponds to an interaction at I=0; hence, what matters is the combination V11+V12. This means that the values we obtain for the parameters in fits to the pseudodata are not meaningful, and only the values of the observables obtained from them are significant. To manage with these correlations, we use the bootstrap or resampling method [7072], generating random centroids of the data with a Gaussian distribution weight and performing a large number of fits to the data with the new centroids and same errors. After each fit, the values of the observables are evaluated, and the average and dispersion for each are calculated.

    IV.   RESULTS
    • Next, we use the model in Section II with a cutoff regularization of qmax=420MeV , as in the study of the Tcc state in Ref. [45]. We obtain a pole at s=7110.41 MeV and the couplings given in Table 1. Similarly, the probabilities obtained and the wave functions at the origin are given in Table 2, and the scattering length and effective range are shown in Table 3. As shown, the probability obtained for the sum of the two B+D0,B0D+ channels is of the order of 96%. The small deviation from unity is due to the energy dependence of the original potential of Eqs. (1) and (3). We also observe that while a1,r0,1 are real, a2,r0,2 are complex because the B+D0 channel is open at the threshold of the B0D+ channel. The couplings are very similar, and so are the wave functions at the origin, indicating an I=0 state, according to Eq. (6). The results for the correlation functions of the two channels are shown in Fig. 2, calculated with R=1fm.

      sp g1 g2
      (7110.41+0i) 31636.8 31631.0

      Table 1.  Pole position and couplings with qmax=420MeV. [in units of MeV]

      P1 P2 ψ1(r=0) ψ2(r=0)
      0.52 0.44 14.75 13.61

      Table 2.  Probability Pi and wave function at the origin ψi(r=0) for channel i.

      a1 a2 r1 r2
      0.71 0.500.16i 0.61 1.221.77i

      Table 3.  Scattering length ai and effective range ri for channel i. [in units of fm]

      Figure 2.  (color online) Correlation functions for the B+D0 and B0D+ channels.

      Next, we discuss the results obtained from the resampling method fits to the data. The data with the assumed errors are shown in Fig. 3. A warning should be given: as discussed in detail in Ref. [36], the sign of V12 is undefined in the present procedure; however, we rely on arguments of heavy quark flavor symmetry to choose solutions with V11 and V12 of the same sign.

      Figure 3.  (color online) Correlation functions for the B+D0 (left) and B0D+ (right) channels, with 26 points in each curve with an error of ±0.02. The centroids of the red data follow the theoretical curve. In blue, we plot the centroids obtained in one of the resampling runs, with a random Gaussian generation of the centroids of each point.

      In Table 4, we show the average values and dispersion of the parameters obtained. As we discussed above, they are not significative owing to the existing correlations. One indication of these correlations is the relatively large parameter errors. Nevertheless, the important aspect is the value of the observables. These can be found in Tables 5, 6, and 7.

      V11 V12 α β qmax/MeV R/fm
      1537.54±918.93 1512.32±913.55 31.89±10.00 165.16±133.31 407.9±60.8 1.01±0.04

      Table 4.  Values obtained for the parameters V11,V12, α,β, qmax , and R.

      sp g1 g2
      7107.84±17.79 34623.08±14300.15 34506.64±14304.57

      Table 5.  Average values and dispersion of the pole position and couplings. [in units of MeV]

      P1 P2 ψ1(r=0) ψ2(r=0)
      0.49±0.03 0.42±0.03 13.61±2.65 12.61±2.423

      Table 6.  Average value and dispersion of the probability Pi and wave function at the origin ψi(r=0) for channel i.

      a1 a2 r1 r2
      0.72±0.03 (0.51±0.02)(0.17±0.01)i 0.61±0.19 (1.41±0.28)(1.65±0.07)i

      Table 7.  Average value and dispersion of the scattering length ai and effective range ri for channel i. [in units of fm]

      In Table 5, we show the value of the energy at which the bound state is found, together with the values of the couplings. As shown, a bound state is found around 7108MeV, compatible with the bound state obtained with the original model within uncertainties. Interestingly, the error obtained is of the order of 18MeV. This is not a small error for a binding energy of 39MeV; however, this is what can be achieved with the assumed precision of the correlation data. More positively, using the data of the correlations at the BD threshold, we are still able to predict that there is a bound state with an approximately 40MeV binding.

      The couplings g1,g2 obtained are also compatible with the original ones, and the errors are also not small. Yet, the approximately equal values of the couplings suggest that we are dealing with an I=0 state.

      It is interesting to analyze our obtained probabilities of the states, which are shown in Table 6. We again obtain numbers for the probabilities of the two channels compatible with those obtained from the original model, but once again, it is the precision by which they can be obtained that is important. The uncertainties are very small, that is, of the order of 6%. This may be surprising in view of the formula used to obtain Pi (Eq. (16)), which is proportional to g2i, and gi has large errors according to Table 5. If g2i is bigger in a fit because the state is more bound, Gi/s also decreases in strength and the product becomes more stable. This is an interesting and fortunate result of our analysis, which allows us to conclude that the application of the inverse method from the femtoscopic correlation functions would allow us to determine the nature of the bound state obtained with a high accuracy.

      As shown in Table 7, it is also rewarding to find that we can determine the scattering lengths with good precision and the effective ranges with smaller precision but significant values.

      We should also stress that the inverse method allows us to obtain the size of the source function with a relative accuracy of approximately 4% (see Table 4).

    V.   CONCLUSIONS
    • We address the problem of evaluating the correlation functions for the BD system using inputs extracted from a successful study of the Tcc(3875) state. We take two channels, B0D+ and B+D0, the small mass differences between which induce visible differences in the correlation functions. The system also develops a bound state of approximately 40 MeV. Once this is achieved, we address the inverse problem of obtaining the observables associated with this system, starting from the correlation functions, assuming errors as in current measurements. Although we obtain results compatible with those obtained from the original model, an important new result is the uncertainty by which we can obtain the observables of the system. We obtain the size of the source with a precision of approximately 4%. We also determine that there is a bound state using the data of the correlation functions above the threshold of the channels, although with an uncertainty of approximately 50% of the binding energy. However, remarkably, we can determine the molecular nature of the obtained state with a good precision of approximately 6%. The scattering lengths of the two channels are obtained with good precision and significant values, and the effective ranges are also obtained but with smaller precision. All these results indicate that the measurement of these correlation functions in the future will allow us to obtain valuable information on the BD interaction and the bound states associated with this interaction. This study could serve as motivation to conduct such measurements in the future.

Reference (72)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return