Combined analysis on nature of X ( 3960 ) , χ c 0 ( 3930 ) , and X 0 ( 4140 )

In this work, a study of the D s ¯ D s interaction and its couplings to the channels D ¯ D and J /ψφ is performed in a quasipotential Bethe–Salpeter equation approach. The D + s D − s and D + D − invariant mass spectra in three-body B decays are investigated in order to understand the origin of X ( 3960 ) , χ c 0 ( 3930 ) , and X 0 ( 4140 ) structures reported at LHCb. With the help of effective Lagrangians, the potential kernel can be constructed with meson exchanges, from which the scattering amplitudes can be obtained. By inserting it into the three-body decay processes, the invariant mass spectra can be calculated with an additional Breit–Wigner resonance introduced. The D + s D − s invariant mass spectrum in the decay process B + → D + s D − s K + is well reproduced, and the X ( 3960 ) structure can be explained as a molecular state from the D s ¯ D s interaction. The state also exhibits as a peak in the D + D − invariant mass spectrum in the decay B + → D + D − K + , however, it is too narrow to reproduce the experimental structure around 3930 MeV. The dip around 4140 MeV in the D + s D − s invariant mass spectrum can be reproduced by the additional resonance with interference effect. However, a small bump instead of a dip is produced from the D s ¯ D s − J /ψφ coupled-channel effect, which suggests that the X 0 ( 4140 ) cannot be interpreted as the D s ¯ D s −


Introduction
Recently, the LHCb collaboration reported two near-threshold structures named X(3960) and X 0 (4140) in the D + s D − s invariant mass spectrum of the decay B + → D + s D − s K + .The mass, width and quantum numbers of these two structures were measured to be: M X(3960) = 3956 ± 5 ± 10 MeV, Γ X(3960) = 43 ± 13 ± 8 MeV, M X 0 (4140) = 4133 ± 6 ± 6 MeV, Γ X 0 (4140) = 67 ± 17 ± 7 MeV, and both with spin parity a Corresponding author: junhe@njnu.edu.cn and charge parity J PC = 0 ++ [1].The resonance peak of X(3960) is very close to the D + s D − s threshold, making it a good candidate of a D + s D − s hadronic molecule.As suggested by the LHCb collaboration, the X 0 (4140) might be caused by either a new resonance or by the D + s D − s − J/ψφ coupledchannel effect, but a determinative conclusion have not been drawn [1].
The X(3960) is certainly not the first state observed in the energy region of 3.9-4.2GeV.Over the past few decades, several experimental candidates of charmonium-like states have been observed, such as the X(3915), the χ c0 (3930) and the χ c2 (3930).The X(3915) was originally observed at Belle in its ωJ/ψ decay mode, and was produced in both B decay [2] and γγ collisions [3] with J PC determined as 0 ++ [4].In 2020, the LHCb collaboration reported the D + D − decay mode of the χ c0 (3930) and χ c2 (3930) using B decays and determined their J PC to be 0 ++ and 2 ++ , respectively [5].
In the current version of the Review of Particle Physics (PDG), the X(3915) decaying to ωJ/ψ and the χ c0 (3930) decaying to D + D − are assumed to be both the χ c0 (3915) with 0 ++ [6].The X(3960) seems to be the same particle as χ c0 (3930) if one considers that the mass and width of the X(3960) state are consistent with those of the χ c0 (3930) meson within 3σ measured by the LHCb Collaboration [1,5].However, such assumption lead to a partial width ratio Γ(X → D + D − )/Γ(X → D + s D − s ) = 0.29 ± 0.09 ± 0.10 ± 0.08.It contradicts the expectation that Γ(X → D + D − ) should be considerably larger than Γ(X → D + s D − s ) if X does not have any intrinsic s s content [1].Hence, the X(3960) and χ c0 (3930) are neither the same resonance, nor the same nonconventional charmonium-like state.It is also puzzled to denote X(3960) as a conventional charmonium-like state because the mass of χ c0 (3P) is around 4131-4292 MeV, which is far away from X(3960).
Many theoretical efforts have been made into understand the origin and structure of the X(3960).Since its mass is close to the D + s D − s threshold, the hadronic molecular state is a promising picture to explain the X(3960).In Refs.[7,8], the molecular state interpretation was proposed to assign the X(3960) as a D + s D − s state with J PC = 0 ++ in the QCD sum rules approach.In Ref. [9], the X(3960) can be well described by either a bound or a virtual state below the D + s D − s threshold.Similar conclusion can be found in Ref. [10], where calculation was performed in an effective Lagrangian approach to study the production rate of X(3960) in the B decays utilizing triangle diagrams and assuming the X(3960) as a bound/virtual state from the D + s D − s interaction.Coupled-channel analysis has also been proposed to reveal the nature of the X(3960) state, which are often combined with the analysis of the χ c0 (3930).Bayar et al. performed a coupled-channel calculation of the interaction D D − D s Ds in the chiral unitary approach and concluded that the X(3960) in the D + s D − s mass distribution and the X(3930) in the D + D − mass distribution are the same state [11].Chen et al. perfromed an analysis on decays χ c0 (3930) → D D, X(3960) → D s Ds , and X(3915) → J/ψω using both a K-matrix approach and a model of Flattelike parameterization.It was suggested that the X(3960) has probably the mixed nature of a cc confining state and D s Ds continuum [12].
There are also several theoretical and experimental works about the X 0 (4140) but its origin is still under debate.The LHCb collaboration interpreted it as a new resonance with the 0 ++ assignment or the D s Ds -J/ψφ coupled-channel effect [1].Badalian et al. considered X 0 (4140) and X(3960) as exotic four-quark states and concluded that the X 0 (4140) is formed in the transitions J/ψφ into D * + s D * − s and back in the framework of the extended recoupling model [13].Agaev et al. studied the process D s Ds − J/ψφ by QCD three-point sum rule analyses, and got the conclusion that the structure X 0 (4140) may be interpreted as a hadronic molecule D + s D − s , whereas the resonance X(3960) can be identified as a tetraquark [cs][ cs] [14].Since different interpretations were proposed in the literature, it is important to come up with new ideas that can help us reveal the origin of X 0 (4140).
In our previous studies, the D + s D − s molecular states were studied in a quasipotential Bethe-Salpeter equation (qBSE) approach [15].Our calculation favors the existence of hidden-heavy bound state D + s D − s with J P = 0 + at a value of cutoff Λ=1.6 GeV.Inspired by the new experimental observation of X(3960) and X 0 (4140), we will explore the D s Ds interaction and its couplings to channels D D and J/ψφ, and compare the theoretical results with the experimental data to discuss the origin of the X(3960), χ c0 (3930) and X 0 (4140).
This work is organized as follows.After the introduction, the formalism of three-body decay will be presented in Sec. 2. Lagrangians used to construct the potential kernels, the qBSE approach, and formula of the invariant mass spectrum will be also presented.The numerical results of the corresponding invariant mass spectrums will be given, and the origins of X(3960), χ c0 (3930), and X 0 (4140) will be discussed in Sec. 3. A summary of the whole work will be given in Sec. 4.
2 Formalism of three-body decay

Mechanism of three-body decay
In this work, the three-body decay mechanism as illustrated in Fig. 1  coupling will be also considered but only in the rescattering process (which is not shown in the diagram) in order to be consistent with experimental analysis [1].

T
Fig. 1 The diagram for the three-body decay and rescattering process of First, we need to treat the direct vertex B + → D + (s) D − (s) K + as shown as a blue full circle in Fig. 1.As in Ref. [16], the amplitudes can be constrained by Lorentz invariance to have a form of , where c 1 and c 2 are coupling constants for processes B + → D s + D s − K + and B + → D + D − K + , respectively.The values of c 1 and c 2 can be estimated by the decay width of each process.It can be obtained by multiplying their branching fractions [17] with the decay width of B meson obtained from its mean life Γ B = 1/τ = 6.105 × 10 11 s −1 = 4.017 × 10 −10 MeV [6].The values of c 1 and c 2 can be obtained as c 1 = 6.027 × 10 −5 and c 2 = 5.397 × 10 −5 .
Here the rescattering effect is not considered, which effect will be discussed later.

Flavor wave functions and Lagrangians involved
The potential kernel of rescattering process as shown in Fig. 1 should be constructed with effective Lagrangians and flavor wave functions to find the pole in the complex energy plane within qBSE approach and to calculate the invariant mass spectrum.In Ref. [18], it was explained explicitly how to construct a potential for states with definite isospins under SU(3) symmetry with the corresponding flavor wave functions.For the S-wave D s + D s − and D + D − states, the wave functions can be constructed as The meson exchanges were adopted to achieve the potential kernel for the interactions considered.According to the chiral and heavy quark symmetries, the Lagrangians for heavy mesons interacting with light mesons read [19,20,21,22,23] where being the mass of the initial or final heavy meson.The V denotes the vector matrix as The parameters involved here were determined in the literature as β = 0.9, g s = 0.76 and g V = 5.9 [24,25,26,27].In Refs.[28,29], it was suggested that the J/ψ exchange is important in the D D * interaction to reproduce the Z c (3900).In this work, the couplings of heavy-light charmed mesons to J/ψ are also considered and written with the help of heavy quark effective theory as [23,30], where the couplings are related to a single parameter g 2 as To investigate whether the structure found around 4140 MeV is caused by the J/ψφ − D s + D s − coupled-channel effect, which should be taken into consideration in rescattering process.Pseudoscalar Ds and D s mesons and vector D * s and D * s mesons are exchanged during the coupling process J/ψφ − D s + D s − .To describe this process, apart from Eq. (2) and Eq. ( 5), we also need the Lagrangians as [23,30] where √ m ψ , and the value of g 2 have been mentioned under Eq.( 5).

Potential kernel
The potential interaction can be constructed by the meson exchange as [15,31,32], Here the vertices Γ can be obtained by the standard Feynman rule from the Lagrangians given in the above.In the current work, the SU(3) symmetry is considered.Hence, with the help of the wave function in Eq. ( 1) and the flavor part in Lagrangians and the matrix in Eq. ( 4), a flavor factor can attracted for a state with a certain exchanged meson based on the Feynman diagram of the meson exchanges [15,18,29,32].Except the flavor factor I ρ = 3/2, I ω = 1/2 for the D D state, all other flavor factor has values of 1 as The value of the flavor factor I K * equals to 1 theoretically, but we take it equals to 0.3, which will be discussed later.The propagators are defined as usual as with q being the momentum of the exchanged meson.A form factor f (q 2 ) = Λ 2 e /(Λ 2 e − q 2 ) is introduced to compensate the off-shell effect of exchanged meson with a cutoff Λ e .The current form factor can avoid overestimation of the contribution of J/ψ exchange.

Quasipotential Beth-Salpeter approach
In the above, we construct the potential of the interactions considered in the current work.The rescattering amplitude can be obtained with the qBSE [18,29,33,34,35].After the partial-wave decomposition, the qBSE can be reduced to a 1-dimensional equation with a spin-parity J P as [29], where T J P and VJ P are partial-wave rescattering amplitude and potential.The n, m, k denotes the independent helicities λ (2,3) of two constituents for the initial, final and intermediate mesons, respectively.If two constituents have helicities λ 1 = λ 2 = 0, such as the case of the D s Ds interaction in the current work, a factor f = 1/ √ 2 should be introduced for the partial wave amplitude T J P and potential VJ P , otherwise, f = 1.The explicit explanations about the independent amplitudes and the factor can be find in Ref. [29].G 0 (p ′′ ) is reduced propagator with the spectator approximation in the center-of-mass frame with P = (M, 0) as As required by the spectator approximation, the heavier particle (remarked as h) is put on shell, which has for the lighter particle (remarked as l) is then W − E h (p ′′ ) with W being the total energy of the system of 2 and 3 particles.Here and hereafter we define the value of the momentum p = | p|.
The potential kernel VJ P m,n /( can be obtained from the potential in Eq. ( 7) as where η = PP 2 P 3 (−1) J−J 2 −J 3 with P and J being parity and spin for system and constituent 2 or 3.
The initial and final relative momenta are chosen as p ′ = (0, 0, p ′ ) and p = (p sin θ, 0, p cos θ).The d J λλ ′ (θ) is the Wigner d-matrix.An exponential regularization was also introduced as a form factor into the reduced propagator as r with the m l and Λ r being the mass of light constituent and a cutoff, respectively [29].The cutoff Λ r and the cutoff Λ e for exchanged meson provide similar effect on the results, hence, we take them as a unified parameter Λ.
The rescattering amplitude T can be obtained by discretizing the momenta p ′ , p, and p ′′ in the integral equation (9) by the Gauss quadrature with a weight w(p i ).After such treatment, the integral equation can be transformed to a matrix equation [29], The propagator G is a diagonal matrix as with on-shell momentum p

Invariant mass spectrum
With the preparation above, the invariant mass spectrum of processes considered can be given with direct three-body decay amplitude A J P m,λ B = A B + →D (s) + D (s) − K + (here we remark it as ÂJ P m,λ B = f m A J P m,λ B with pseudoscalar B meson having a helicity λ B = 0) and the rescattering amplitude T J P as Ref. [33], with where M is the mass of the B meson.N 2 J = (2J + 1)/4π = 1/4π with J = 0 here.p cm 3 and p ′cm 3 are the center of mass momentum in the 23 and 2 ′ 3 ′ system, which can be expressed as p cm Eq. ( 15) can be abbreviated as matrix form M = A + TGA by the same discretizing in Eq. ( 12), where T satisfies T = V + VGT .The rescattering amplitude T can be solved as T = (1 − VG) −1 V. Since we focus on the pole of the rescattering amplitude, we need to find the position where |1 − VG| = 0 with z = E R + iΓ/2 corresponding to the total energy and width at the real axis.

Numerical results and discussions
With the above preparation, the invariant mass spectra of B decays considered in the current work will be calculated and compared with the experimental data.In this section, the D s Ds interaction and its couplings to channels D D and J/ψφ will be considered.First, we consider the simplest model, that is, only the D s Ds interaction is included in the rescattering to provide a basic picture of molecular state from the interaction.In such model, the intermediate and final products are all D s Ds .And in the current work, we only consider the S-wave state with spin parity J P = 0 + .The corresponding cutoff Λ is adjusted to 1.8 GeV, which is a little bigger than the one in our previous work [15], which makes the result fit better with the new experimental data.Besides the obvious peak near the D s Ds threshold, a structure can be seen around 4140 MeV.Obviously, such structure can not be produced from the D s Ds interaction as shown in Eq. ( 14).Here, as suggested by LHCb collaboration [1], we introduce a Breit-Wigner resonance, which is used to fit the dip around 4140 MeV as where M 0 and Γ 0 are the mass and width of the structure obtained from the experimental results.The strength factor a and phase angle θ are free parameters.The theoretical results for events can be obtained by multiplying the theoretical decay distribution on efficiency, which equals to N sig /(N B + B B + →D s + D s − K + ) with N sig , N B + and B being the numbers of signal events, the total number of B candidates, and the branching fraction, respectively.Unfortunately, though the N sig and B were provided by LHCb collaboration in Ref. [1,17], the N B + was not reported, so we have to introduce a parameter multiplying on the theoretical decay distribution to fit the experimental data.After such treatment, the comparison between the theoretical and experimental results can be carried out.We would like to remind that such treatment make the absolute values of the coupling constants c The black (full), red (dashed), and blue (dotted) curves are for the total, rescattering as M in Eq. ( 15), and Breit-Wigner resonance in Eq. ( 16).The red points with error bars are the data from the LHCb experiment cited from Ref. [1]. is found at 3932.8 GeV in the first Riemann sheet on the real axis near D s Ds threshold, which suggests a very weakly bound state from the D s Ds interaction and close to the mass of X(3960) measured by the LHCb Collaboration.To confirm their relation, the invariant mass spectrum in energy region of 3800 − 4800 MeV is calculated and compared with the experiment.In the lower panel of Fig. 2, a sharp peak arises near the D s Ds threshold, which is from the D s Ds rescattering.With increase of the energies, a dip can be found around 4140 MeV, which is due to the Breit-Wigner resonance introduced in Eq. ( 16).We evaluate the χ 2 between the experimental data and our theoretical results in two energy regions of 3800 − 4170 MeV and of 3800 − 4800 MeV.The former one is predominantly relevant to the X(3960) and X 0 (4140) structures, and the latter one is relevant to all energy region.The χ 2 in energy region of 3800 − 4170 MeV is equal to 13.45 for 13 data points, while the χ 2 in energy region of 3800 − 4800 MeV is equal to 53.06 for 43 data points, which indicates that our theoretical data are in a good agreement with experimental data, and favors the assumption of the X(3960) as a D s Ds molecular state.

Discussion about relation between X(3960) and χ c0 (3930)
In this subsection, we perform a study of the B meson decay with rescattering of D s Ds -D D in order to discuss the relation between X(3960) and χ c0 (3930).The intermediate and final products can be D s Ds and D D with spin parity J P = 0 + .The value of flavor factor I K * is adjusted to 0.3 instead of 1.It is interesting to found that if we take I K * = 1 as its theoretical value, the D D state will remain but the D s Ds state will disappear, only if decreasing the value of I K * will these two states appear, which is analogous to Ref. [11] and may cause by SU(3) symmetry breaking for the K * meson.A Breit-Wigner resonance is still added into the amplitude of D s Ds -D D rescattering to match the dip found around 4140 MeV in Refs.[1,17].The corresponding cutoff Λ will be also chosen as 1.8 GeV for both D s Ds and D D channels.The D + s D − s and D + D − invariant mass spectra are presentedn in Fig. 3 and compared with the experiment.
In the upper panel of Fig. 3, the poles leave the real axis and become two conjugate poles in the complex energy plane after including the coupled-channel effect.These two conjugate poles are located at 3933.6±0.25iMeV, which means a molecular state with a mass of 3933.6 MeV and a width of 0.5 MeV close to the corresponding D s Ds threshold.The D + s D − s invariant mass spectrum in the middle panel of Fig. 3 looks similar to the single-channel results in energy region of 3800 − 4170 MeV in Fig. 2. The χ 2 between the experimental data and our theoretical data in energy region of 3800 − 4170 MeV is equal to 14.88 for 13 data points, which is close to the above single-channel results.However, the χ 2 in energy region of 3800 − 4800 MeV is equal to 76.13 for 43 data points, which is a little larger than the Events Re(z)(GeV) Fig. 3 The pole (upper panel) and theD + s D − s (middle panel) and D + D − (lower panel) invariant mass spectra with the D D-D s Ds coupledchannel interaction.The black (full), red (dashed), and blue (dotted) curves are for the total, rescattering as M in Eq. ( 15), and Breit-Wigner resonance in Eq. ( 16).The red points with error bars are the data of X(3960) (middle panel) and X(3930) (lower panel) from the LHCb experiment [24].one obtained in subsection.3.1.Since in the higher energy region, the real mechanism should be more complex, the current result is still accepted if we focus on the structure near the D s Ds threshold.
In the lower panel of Fig 3, the D + D − invariant mass spectrum are presented and compared with the data of X(3930) from LHCb experiment [5].As shown in Fig 3, the peak near 3930 MeV is extremely narrow, seems too sharp compared with the experimental structure.If a wider peak in the D + D − invariant mass spectrum is reached, the peak in the D + s D − s invariant mass spectrum will disappear.Our calculation suggests that though peaks can be produced simultaneously near the D s Ds threshold in both D + s D − s and D + D − invariant mass spectrum, the explicit shape of the experimental structure can not be well-fitted simultaneously.In the experimental article [5], such structure is suggested to be formed by two states χ c0 (3930) and χ c2 (3930), and the χ c0 (3930) has a smaller width.Hence, it is still possible to assign the χ c0 (3930) as a D s Ds molecular state but with a very small width.If so, the χ c2 (3930) may play a more important role to form the structure around 3930 MeV in the D + D − invariant mass spectrum.

Events
Re(z)(GeV) Fig. 4 The D + s D − s invariant mass spectrum with the D s Ds -J/ψφ coupled-channel interaction.The yellow, red, black, blue and green curves are for the mass distribution with flavor factor (FF) equal to 2, 1.5, 1, 0.5 and 0.1, respectively.The red points with error bars are the data of X(3960) from the LHCb experiment [1].
The black curve is the mass distribution obtained by theoretical values of flavor factor of the exchanged mesons, that is, I D s =I Ds =I D * s =I D * s =1.The peak near 3960 MeV is analogous to Fig. 2 and Fig. 3, which also favors the assumption of the X(3960) as a D s Ds molecular state, whereas a small bump is obtained instead of a dip at around 4140 MeV.This bump-like structure in the D + s D − s invariant mass spectrum is caused by the D s Ds -J/ψφ coupling.However, we cannot reproduce a dip from this hump-like structure by adjusting the parameters.In the current model, the coupled-channel effect with coupling to the J/ψφ channel is constrained by the heavy quark symmetry as shown in the Lagrangians.In the experimental treatment by the LHCb collaboration [1], the free parameters were introduced to adjust the contributions of the coupled-channel effect.Hence, our results do not conflict with experimental analysis, but are obtained with more theoretical constraint.
Besides, adjusting other parameters, such as the value of flavor factor, doesn't work either.To prove this point, we take values of flavor factor (FF) of the exchanged mesons I D s , I Ds , I D * s and I D * s in this process all equal to 0.1, 0.5, 1, 1.5 and 2, respectively.From Fig. 4 we can see although the value of the peak for the bump-like structure around 4140 MeV gradually decrease with the decreasing of the value of FF, a dip structure still cannot be reproduced in our model.Based on these calculations, one can say that, given the fact that the coupled channel D s Ds -φJ/ψ cannot reproduce a dip around 4140 MeV in our model, the X 0 (4140) could not be easily interpreted as a D s Ds -φJ/ψ coupling.

Summary
Inspired by the newly observed X(3960), we study the B meson decay with D s Ds rescattering and couplings to channels D D and J/ψφ in a qBSE approach with spin parity 0 + .With the help of effective Lagrangian, the potential kernel can be constructed by meson exchanges and the invariant mass spectra are calculated and compared with the LHCb experiment.
With the D s Ds rescattering, a pole is found at 3932.8 MeV on the real axis near D s Ds threshold, which shows that a very weakly bound state can be produced from the D s Ds interaction.A peak structure near the D s Ds threshold is found in the D + s D − s invariant mass spectrum and have a good agreement with the experimental results, favoring the assumption of the X(3960) as a D s Ds molecular state.After adding D D as intermediate and final channel in the process, the theoretical D + s D − s invariant mass spectrum is still in agreement with the relevant experimental data.
After adding D D channel, an extremely sharp peak can be obtained near 3930 MeV in D + D − invariant mass spectrum.The χ c0 (3930) can be assigned as the same D s Ds state in the D + s D − s invariant mass spectrum, but with a very small width.The experimental structure around 3930 MeV in the D + D − invariant mass spectrum can not be reproduced in the current model only with this molecular state with 0 ++ from the D s Ds interaction.The experimental analysis of such structure suggests that it may be formed by states with spins 0 and 2, χ c0 (3930) and χ c2 (3930).The current results indict that the χ c2 (3930) may play a more important role to form the structure near 3930 MeV in the D + D − invariant mass spectrum [5].
With an additional Breit-Wigner resonance, the dip can be well reproduced around 4140 MeV in the D + s D − s invariant mass spectrum.The possible role of the D s Ds -J/ψφ coupled-channel effect on the appearance of the dip X(4140) is also discussed in the current model.However, a small bump instead of a dip is found around 4140 MeV.Moreover, the dip can not be reproduced by adjusting the parameters, which does not support the assignment of X 0 (4140) as the D s Ds -J/ψφ coupled-channel effect.It suggests that the X(4140) has an independent origin different from the peak X(3960) which can be well reproduced from the D s Ds interaction in the current work.

3. 1
D + s D − s invariant mass spectrum with D s Ds interaction

Fig. 2
Fig.2The pole (upper panel) and the D + s D − s invariant mass spectrum for decay B + → D + s D − s K + (lower panel) with the D s Ds interaction.The black (full), red (dashed), and blue (dotted) curves are for the total, rescattering as M in Eq. (15), and Breit-Wigner resonance in Eq. (16).The red points with error bars are the data from the LHCb experiment cited from Ref.[1].

3. 3
Discussion about origin of X(4140)As suggested by the LHCb collaboration, the additional structure found around 4140 MeV in the D + s D − s invariant mass spectrum might be caused either by a new resonance with the 0 ++ assignment or by a D s Ds -J/ψφ coupledchannel effect[1].In the above, we adopt the former assumption by introducing a Breit-Wigner resonance.For better understanding of the origin of this structure, the D s Ds -J/ψφ coupled-channel effects with spin parity J P = 0 + is introduced to replace the Breit-Wigner resonance.The corresponding cutoff Λ will be chosen as 1.8 GeV as above calculation.We also study effect of flavor factors of exchanged meson I D s , I Ds , I D * s and I D * s by changing their values.The result is shown in Fig.4.