-
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 [1–14] (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 [18–33], and a model independent analysis of the correlation functions was very recently proposed [34–38], 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 theBD system with isospinI=0 in Ref. [39]. The observation of theTcc state [40, 41] and subsequent theoretical research support this state as a bound state of theD0D∗+ ,D+D∗0 channels inI=0 [42–53]. Two recent studies investigatingTcc from a different perspective also concluded that theTcc state is indeed a molecular state [54, 55]. In the first study, the scattering length and effective range of theD0D∗+ ,D+D∗0 channels, as well as the shape of theD0D0π+ mass distributions, were theoretically studied [54], concluding thatTcc is a molecular state of theD0D∗+ ,D+D∗0 components, essentially constructing anI=0 state. In the second study, a different path was followed, assuming thatTcc 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 aTcc state with negligible molecular probability, but at the heavy cost of a smallD0D∗+ 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 thanTcc because there is a general rule that the heavier the quarks, the stronger the interaction and binding energies [56–58]. This rule is also satisfied when the extension of the local hidden gauge approach [59–62] is used as a source of interaction, exchanging vector mesons between the heavy meson components. In Ref. [39], differentB(∗)D(∗) ,B∗ˉD(∗) pairs were found to be bound using the same regulator for the loops as in Ref. [45], particularly theBD state inI=0 , which was found to be bound by approximately15−30MeV . The existence of one bound state in theBD 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 fromTcc , 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 ofD0B+ andD+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 theD0B+ andD+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 theB+D0(1) ,B0D+(2) channels is given byFigure 1. Vector exchange between the
D0B+ andD+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[3s−2(m2B+m2D)−(m2B−m2D)2s].
(3) The scattering matrix is then given by
T=[1−VG]−1V,
(4) where
G=diag(G1,G2) , withGi being theBD meson loop function regularized with a cutoffGi(√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 , andmi,Mi are the masses of the D and B mesons in channel i. With the isospin states(B+,B0) ,(D+,−D0) , theI=0 combination is given by|BD,I=0⟩=−1√2|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). TheI=1 combination is given by|BD,I=1,I3=0⟩=−1√2|B+D0−B0D+⟩.
(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 takeqmax of the order of420MeV , as used in Ref. [45], to obtain the binding of the relatedTcc state. -
We can evaluate the scattering length a,
r0 for theB+D0 andB0D+ channels by recalling the relationship between our T matrix and that used in quantum mechanics [64]T=−8π√sfQM≃−8π√s1−1a+12r0k2−ik,
(9) with
k=λ1/2(s,m21,m22)2√s,
(10) from which we easily find
−1ai=(−8π√sT−1ii)|sth,i,
(11) r0,i=[2√sμi∂∂s(−8π√sT−1ii+iki)]sth,i,
(12) where
μi is the reduced mass in channel i, andsth,i andki 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=gigjs−sp,
(13) where
sp is the square of the mass at the pole. Thus,g21=lims→sp(s−sp)T11,
(14) g1gj=lims→sp(s−sp)T1j,
(15) which determine the relative sign of
g2 with respect tog1 . Once the couplings are evaluated, we calculate the molecular probabilities of theB+D0 andB0D+ channels, as in Refs. [64, 66, 67].Pi=−g2i∂Gi∂s|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πθ(qmax−pD0)∫drr2S12(r)×{|j0(pD0r)+TB+D0,B+D0(E)˜G(B+D0)(r;E)|2+|TB0D+,B+D0(E)˜G(B0D+)(r;E)|2−j20(pD0r)},
(18) CB0D+(pD+)=1+4πθ(qmax−pD+)∫drr2S12(r)×{|j0(pD+r)+TB0D+,B0D+(E)˜G(B0D+)(r;E)|2+|TB+D0,B0D+(E)˜G(B+D0)(r;E)|2−j20(pD+r)},
(19) where
pi is the momentum of the particles in the rest frame of the pair,pi=λ1/2(s,m2i,M2i)2√s,
(20) S12(r) is the source function, parameterized as a Gaussian normalized to1 ,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, andE=√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 unknownqmax . 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(V11−V22)=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 [34–36]:
V11=V′11+αm2V(s−sth,1),
(26) V12=V′12+βm2V(s−sth,1),
(27) where the factor
m2V (withmV=800MeV ) is introduced to makeα,β dimensionless. Then, we have4 free parameters for the interaction, as well asqmax and R, which is a total of6 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 atI=0 ; hence, what matters is the combinationV11+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 [70–72], 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 theTcc 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 twoB+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 whilea1,r0,1 are real,a2,r0,2 are complex because theB+D0 channel is open at the threshold of theB0D+ channel. The couplings are very similar, and so are the wave functions at the origin, indicating anI=0 state, according to Eq. (6). The results for the correlation functions of the two channels are shown in Fig. 2, calculated withR=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.50−0.16i −0.61 1.22−1.77i Table 3. Scattering length
ai and effective rangeri for channel i. [in units of fm]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 withV11 andV12 of the same sign.Figure 3. (color online) Correlation functions for the
B+D0 (left) andB0D+ (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.
V′11 V′12 α β 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
V′11,V′12 ,α,β ,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 rangeri 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 of18MeV . This is not a small error for a binding energy of39MeV ; however, this is what can be achieved with the assumed precision of the correlation data. More positively, using the data of the correlations at theBD threshold, we are still able to predict that there is a bound state with an approximately40MeV 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 anI=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 tog2i , andgi has large errors according to Table 5. Ifg2i 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 theTcc(3875) state. We take two channels,B0D+ andB+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 theBD interaction and the bound states associated with this interaction. This study could serve as motivation to conduct such measurements in the future.
