The scalar exotic resonances X(3915), X(3960), X(4140)

The scalar resonances $X(3915), X(3960), X(4140)$ are considered as exotic four-quark states: $cq\bar c \bar q, cs\bar c \bar s, cs\bar c\bar s$, while the $X(3863)$ is proved to be the $c\bar c, 2\,^3P_0$ state. The masses and the widths of these resonances are calculated in the framework of the Extended Recoupling Model, where a four-quark system is formed inside the bag and has relatively small size ($\la 1.0$~fm). Then the resonance $X(3915)$ appears due to the transitions: $J/\psi\omega$ into $D^{*+}D^{*-}$ (or $D^{*0}\bar D^{*0})$ and back, while the $X(3960)$ is created due to the transitions $D_s^+D_s^-$ into $J/\psi\phi$ and back, and the $X(4140)$ is formed in the transitions $J/\psi\phi$ into $D_s^{*+}D_s^{*-}$ and back. The characteristic feature of the recoupling mechanism is that this type of resonances can be predominantly in the $S$-wave decay channels and has $J^P=0^+$. In two-channel case the resonance occurs to be just near the lower threshold, while due to coupling to third channel (like the $c\bar c$ channel) it is shifted up and lies by (20--30)~MeV above the lower threshold. The following masses and widths are calculated: $M(X(3915))=3920$~MeV, $\Gamma(X(3915))=20$~MeV; $M(X(3960))=3970$~MeV, $\Gamma(X(3960)=45(5)$~MeV, $M(X(4140))= 4120(20)$~MeV, $\Gamma(X(4140))=100$~MeV, which are in good agreement with experiment.


Introduction
In the region (3.9-4.2) GeV there are now three scalar resonances and the X(3915) was the first, observed by the Belle in the e + e − → J/ψωK process [1].Later this resonance was confirmed by the BaBar [2] and in several other experiments [3]), in particular, in two-photon collisions [4,5].For some years this resonance was assumed to be the conventional cc meson -χ co (2P ), although this interpretation has called out some doubts [6,7] (see discussion in the reviews [8,9]) and does not agree with predictions in different relativistic potential models (RPM) [10]- [13].The experimental masses of the X(3915) and χ c2 (2P ) were found to be almost equal, while in the RPMs a smaller mass, M(2 3 P 0 ) ∼ = 3870 ± 30 MeV, and much larger mass difference, δ 20 (2P ) = M(χ c2 (2P ) − M(χ c0 (2P ) ∼ = (70 − 100) MeV, were predicted.Notice that large mass difference δ 20 is kept even if the coupling of the χ c0 (2P ) to open channels is taken into account [14,15].Such theoretical expectations were supported by the Belle observation of the wide scalar X(3860) resonance [16], both in e + e − → J/ψD + D − and e + e − → J/ψD 0 D0 decays, which has the mass M = 3862 +26 −32 +40 −82 MeV and large width Γ ∼ = 200 MeV.The existence of the scalar X(3860) resonance is confirmed by the analysis of two-photon production, γγ → D D in [17].
Very recently the LHCb [18] has observed two more scalar resonances X(3960), X 0 (4140) in the D + s D − s mass spectrum in the B + → D + s D − s K + decays with the parameters: 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, both with J P C = 0 ++ .These new scalar resonances evidently look as exotic states and the X(3960) was interpreted as the molecular D + s D − s state within the QCD sum rules approach [19,20] and in a coupledchannel model [21]; in [22] it appears due to the triangle singularity, while in [23] the parameters of the X(3960), as a diquark-antidiquark state, were obtained in a good agreement with experiment, using the QCD sum rules approach.Notice that the masses of the X(3960) and X(4140) resonances lie by ∼ 20 MeV above the thresholds: D + s D − s and J/ψφ, respectively.In our paper we assume that the X(3915) and both the X(3960), X 0 (4140) belong to exotic four-quark states cqcq and cscs and to define their parameters we will use the Extended Recoupling Model (ERM), recently suggested in [24], which develops the Recouplimg Model, presented earlier [25].The ERM allows to calculate the mass and width of a scalar four-quark states, however, within suggested mechanism such resonances cannot exist in the systems with two identical mesons, like D + s D + s , D * + s , D * + s .This theoretical prediction is supported by the Belle experiment [26].In the ERM the system of two mesons, e.g.(J/ψ + φ), can transfer into another pair of the mesons (D + s , D − s ) by rearranging confining strings and back in the infinite chain of transformations, like J/ψφ → (D + s D− s ) → J/ψφ → ....Note that such sequences can also be treated, for example, in the standard OBE approximation with the meson exchanges, which, however, does not produce the singularities near the thresholds.In the coupled-channel models (CCM) [27,28] the interaction between hadrons, like D + s D − s , J/ψφ, is usually neglected, while in the ERM such interaction is taken into account, introducing the fourquark bag.It is important that all hadrons involved have rather small sizes, ∼ = (0.40 − 0.55) fm and only ω(1S) has a bit larger r.m.s.∼ 0.7 fm.We would like to underline the characteristic features of the ERM [24]: first, due to the string rearrangement of a four-quark system the singularity lies close to the lower threshold; second, this mechanism produces the resonance in the S-wave hadron-hadron system and therefore, the quantum numbers of these resonances J P C = 0 ++ , 1 ++ , 2 ++ ; third, a resonance does not appear, if hadrons are identical.
In the literature there are still a controversy, concerning the X(3915), and different interpretations were proposed.This resonance was considered in the tetraquark model within the Born-Oppenheimer approach in [29,30,31,32], due to the triangle singularity [22] and the threshold effects [33], as the molecular D s Ds bound state [34], or the lightest cscs state [35] and as the diquark-antidiquark state, using the QCD sum rule method [23,36].In contrast to a molecular structure of four-quark states in the ERM these systems are assumed to be compact systems, similar to the diquark-antidiquark states studied in [37].In such compact systems their wave functions at the origin are not small and therefore they can be produced in the γγ transitions.At this point one can assume a possible existence of at least two different but subsidiary mechanisms, producing resonances in the four-quark and multiquark systems: first, the resonances, which are formed inside a common multiquark bag and connected with external independent channels.As a result these resonances could be seen in all external channels.The theory of this type of approach was suggested long ago in [38].Within the diquark-qntidiquark model the compact Q 2 Q2 resonances were already predicted in 1988 [37].Second type of multiquark resonances refers to the channel-coupling resonances where the internal multiquark region is only needed to connect different external channels with sufficient probability and the considered here Extended Recoupling Model belongs to this second type.One can easily imagine the existence of mixed type models and mechanisms where two these dynamics interfere with each other.In what follows we shall consider only the ERM mechanism.
In our paper we will shortly discuss the higher scalars, X(4500), X(4700), observed by the LHCb [39], which admit different interpretations.
The structure of the paper is as follows.In next section we shortly remind the basic formulas in two-channel case and give the values of the parameters, needed to define the masses and widths of the recoupled four-quark resonances.In section 3 more general matrix representation of the ERM is presented.In section 4 we calculate the transition amplitudes and give the masses and widths of the scalar resonances, and compare them with experimental data.In section 4 the masses of high X(4500), X(4700) resonances, as the cc states, are discussed.Our conclusions are presented in section 5.

The two-channel approach in the Extended
Recoupling Model We study the experimental process where, among other products, two hadrons are produced and one pair of hadrons (the pair 1) can transfer into another pair of hadrons (the pair 2).In [24] the probability amplitude of this transition was denoted as V 12 (p 1 , p 2 ), with p 1 , p 2 -relative momenta of the hadrons, referring to the pair 1 and 2. If an infinite set of the transformations was supposed and the total production amplitude A 2 of the pair 2 was written as a product of the slowly varying function F (E) and the singular factor . This definition of the transition amplitude V 12 = V 21 differs of that in other approaches, where one or more the OBE diagrams with meson exchanges are taken.In the ERM [24] the process occurs through the intermediate stage of the Quark Compound Bag (QCB) [38,40], where all quarks and antiquarks of two hadrons are participating in the string recoupling and, possibly, the spin recoupling.Denoting the QCB wave functions as Φ(q i ) (i = 1, 2, 3, 4) and the two-hadron wave functions as Ψ i (h 1 , h 2 ), the amplitude V 12 can be written as, i.e. the amplitude V 12 = 1 1−N acquires the factorized form: ) with the factor N, written as ( Here z = z(E) can be called the transition probability, while I 1 (E), I 2 (E) are the following integrals (see [24]): where the hadron energies 2m ′ + m ′ , include corresponding thresholds E th i and the reduced masses µ i , namely, The result of the integration in I i (E) can be approximated by the form: with µ i , defined in (4), while ν i is expressed via the parameters of the hadron wave functions, which were calculated explicitly in [24].Here we would like to underline that the transition probability z(E) appears to be the only fitting parameter in the ERM.The whole series of the transitions from the pair 1 to 2 and back is summed up to the amplitude f 12 , where ν i are found from the four-quark wave functions, as in [37,40].The form of Eq. ( 6) takes place for the energies E > E 1 , E 2 , while for E < E 1 , E 2 , i.e. below thresholds, the amplitude . It is important that in the ERM the process proceeds with the zero relative angular momentum between two mesons, L = 0, otherwise the transition probability z 12 (E) is much smaller and a resonance may not appear.
Note also that if the recoupling mechanism is instantaneous, or the transition from one pair of the mesons to another proceeds instantaneously, then the transition amplitude V (12) does not factorize into V (1)V (2); such an assumption was used in the original Recoupling Model [25].However, in this approximation, e.g. for the T cc resonance agreement with experiment was not reached [25].On the contrary, in the ERM [24] the recoupling mechanism proceeds in two stages: at first stage the hadrons h 1 , h 2 collapse into common "compound bag" [38,40], where the four quarks are kept together by the confining interaction between all possible quark pairs.This compound bag has its own wave function Φ i (q 1 , q 2 , q 3 , q 4 ) and the probability amplitude of the h 1 , h 2 → Φ transition, which defines the factor V 1 (p 1 ) in Eq. ( 2).In a similar way the transition from the Bag state to the final hadrons h 3 , h 4 defines the factor V 2 (p 2 ) and we obtain the relation: and similar equation for v 2 (p 2 ), replacing h 1 , h 2 by h 3 , h 4 .From v i (p i ) the function I i (3) is defined and using (6), one obtains ν i .Now we give experimental data and corresponding the ERM parameters, referring to the four-quark systems, cqcq for X(3915) and cscs for the X(3960), X(4140).We give also the threshold energies E 1 , E 2 .
2) X(3960), J P = 0 + , Γ(exp .)= 43 (21) 3) X(4140), J P = 0 + , Γ(exp .)= 67 (24) Here q can be u, d quarks.To define the structure of the cross sections we start with the value of the recoupling probability z = 0.2 GeV 2 and the parameters from the item 1) to obtain the distribution |f 12 (E)| 2 ; the values of |f 12 (E)| 2 will be given in Section 4. In the amplitude f 12 (E) the resulting singularity can be found in the form of ( 6) and for equal threshold masses it produces a pole nearby thresholds; however, real distance between the thresholds is large, ∼ 100 MeV and the actual singularity structure can be more complicated.

The matrix approach in the ERM
In previous Section we have presented the ERM equations in the case of two channels, which are convenient to define the mass of a resonance.However, they do not allow to study some details of the process, or to consider a larger number of channels, which can have a influence at the properties of a fourquark system.Therefore here we present a more general representation of the amplitude using the unitarity relation, when the standard form of the transition amplitudes f ij (E) (for L = 0) is or the unitarity relation can be realized through the M-matrix representation, where f , M , k are the matrices in the channel numbers [28].In some cases instead of the M it is more convenient to use the K matrix, M = − K−1 , where the matrix elements (m.e.) M ik (E) are the real analytic functions of E with the dynamical cuts.For two-channel system fM can be written as Here One can easily establish the relation between the equations ( 10)-( 12) and the amplitude f 12 (ERM) (6) in two-channel case, which is a partial case of these equations: and One can see that for z > 0 the values ν i = M ii are real analytic functions of E. In the ERM [24] ν i were positive constants (defined via the parameters of the compound bag model), while in general case Eqs. ( 12)-( 14) include other transition m.e.s f ik .Later in our analysis we will be interested only in the denominator D(E) (12) and the factors in ( 13), (14), which fully define the position of a resonance.
The value of z, in principle, can be calculated within the ERM, however, it can depend on many unknown parameters, and at the present stage we prefer to keep z as a single fitting parameter.It can be shown that z depends on the width of a resonance, but weakly depends on the resonance position.Now we consider three channels case to study more realistic case and choose the situation, when a resonance lies above the threshold 3.Here we do not need to specify the channel 3, which for example, may be a conventional cc state with J P C = 0 ++ .We introduce the 3 × 3 amplitude fM (E) with three thresholds E i (i = 1, 2, 3) and the momenta Here m 1i , m 2i are the masses of two hadrons in the channel i.In this case the form of Eq. ( 9) is kept, 16) For the energy E below the thresholds, 1 and 2, −ik 1 = |k 1 |, −ik 2 = |k 2 |, and the factor ∆M is a real function of E. For the threshold 3 below thresholds of 1 and 2 one can define the poles of the amplitude f3 , or the zeroes of D 3 (E), and rewrite the Eq. ( 15) as, where the transition probability z(E) One can see that z(E) acquires imaginary part, which can be of both signs.Therefore the influence of the third (or more) open channels, lying below the thresholds E 1 , E 2 in the 2 × 2 matrix f 12 (E), may be important in some cases.The channel 3 can be taken into account, introducing complex values of z(E), which can depend on the energy as in Eq. ( 18).

The masses and widths of the scalar resonances
We start with the X(3915) resonance and consider the following recoupling process: J/ψω → D * D * .At first we look at two-channel situation and choose the recoupling parameter z 2 = 0.18 GeV 2 .For the X(3915) structure -cqcq the parameters µ i , ν i , E i are given in the item 1) of section 2. Then inserting all parameters to the Eq. ( 13), one obtains the distribution |f 12 (E)| 2 (f 2 ≡ f 12 ).Its values for different E are given in Table 1, which show that the maximum takes place at E = 3880 MeV, just near the lower threshold, and Γ 2 = Γ(2 − channels) ∼ = 15 MeV.In experiment for this resonance, observed by the Belle group in the process e + e − → e + e − J/ψω [1], the larger mass M(exp .)= (3918.4± 1.9) MeV and Γ(exp .)= (20 ± 5) MeV [3] were obtained.
In the case of 3-channels, when e.g. the coupling to the cc channel is taken into account, the factor z 3 (E) acquires an imaginary part.In this case we calculate the amplitude f 3 (E), taking z 3 = (0.18−i0.20)GeV 2 ; the values of |f 3 (E)| 2 are given in Tab. 1.  1 one can see that in the 3-channel case the peak is shifted up by ∼ 35 MeV and corresponds the mass E R ∼ = 3.915 GeV and the width Γ 3 ∼ = 20 MeV, which are in good agreement with the experimental mass and Γ(exp.)= 20 (5) MeV [3].
The scalar resonance X(3960) with J P C = 0 ++ was recently observed by the LHCb in the B + → J/ψφK + [18] and within the ERM it can be explained due to the infinite chain of the transitions: J/ψφ → D + s D − s and back.In two-channel approximation the X(3960) parameters (ν i , µ i , E i , (i = 1, 2) are given in the item 2) (Section 2), which are used to define the amplitude (13).First, we choose z 2 = 0.30 GeV 2 and calculate the transition amplitudes |f 12 (E)| 2 ; their values are given in the Table 2.In the two-channel approximation the numbers from Table 2 show the peak at E = 3940 MeV, near D + s D − s threshold, and Γ(2 − ch.) ∼ = 15 MeV.In the 3-channel case the mass of the X(3960) resonance is shifted up to the position M(3 − ch.) = 3970 MeV and the width increases to the value Γ(th.) ∼ = 45 (5) MeV; these values are in agreement with the experimental numbers: M(X(3960)) = 3956(15) MeV, Γ(X(3960)) = (43 ± 21) MeV [18].
In [18] the LHCb has reported about another, the X(4140) resonance, with (12) MeV is close to the J/ψφ threshold.We consider this resonance as the cscs system and first calculate the squared amplitudes |f 12 (E)| 2 in twochannel case, taking the parameters µ i , ν i , E i from the item 3) of Section 2. In this 2-channel case: J/ψφ and D * + s D * − s the transition probability z 2 = 0.35 is taken and the calculated values of |f 12 | 2 are given in Table 3.
In three-channel case the channel D + s D − s is added as the third one, then the values |f 3 | 2 are calculated for z 3 = 0.20 − i0.20 and given in Table 3.
From Table 3 one can see the peak at E R = (4.09± 0.01) GeV, Γ(th.) = 60 MeV in two-channel approximation and the peak at E R = (4.12±0.02)GeV  Our numbers in Tables 1-3 show that in two-channel case the resonance always lies just near the lower threshold, however, if the coupling to the third channel is taken into account, then it is shifted up and its position occurs to be close to the experimental number.The masses and widths of the exotic resonances, X(3915), X(3960), X(4140), defined in the ERM, are given in the Table 4 together with experimental data.5) 43 (21) [18] X(4140) 4120( 20) 4133 (12) 100 67 (24) [18] From Table 4 one can see that in the ERM the predicted masses and the widths of the scalar four-quark resonances are in good agreement with experiment, if besides two channels, which creates the resonance, the coupling of the resonance to third channel is taken into account.
Comparing our results with those in literature, one can notice that our conclusions on the four-quark structure of the X(3915), X(3960, X(4140)) also agree with the analysis in the paper [33], based on the coupled channel model of the cc and meson-meson systems.Notice that the general structure of the channel-coupling matrix elements in both approaches is similar.

The scalar X(4500), X(4700) resonances
High scalar resonances X(4500), X(4700), or χ c0 (4500), χ c0 (4700), [39], were studied in many papers and for them two interpretations were suggested.First, the X(4500) and X(4700) are considered as the cc states -4 3 P 0 and 5 3 P 0 and their masses were calculated in relativistic quark models, where coupling to open channels was taken into account [14,15,41].In [41] the influence of open channels is studied using the so-called screened potential [11], while in [13] the spectrum was calculated using the relativistic string Hamiltonian [42] with the flattened confining potential [43]; this flattening effect arises due to creation of virtual q q pairs.Notice that the flattened confining potential appears to be universal for all types of the mesons and it produces the hadronic shifts down ∼ (100 − 130) MeV for the 4P, 5P charmonium states and gives the masses of the 4 3 P 0 , 5 3 P 0 states in a reasonable agreement with experiment [13].On the contrary, in [44], within the 3 P 0 model, much smaller shifts due to the coupled-channel effects, < ∼ 30 MeV , were obtained for the 4 3 P 0 , 5 3 P 0 states, while in [41] these states acquire too large mass shifts for the chosen screened potential.
Model-independent analysis of the cc spectrum can also be done by means of the Regge trajectories, if they are defined not for the meson mass M(nL) but for the excitation energy: E(nL) = M(nL) − 2 mQ [45], where mQ is the current heavy quark mass [13]: (19) This Regge trajectory gives M(4 3 P 0 ) = 4.474 GeV and M(5 3 P 0) = 4.719 GeV, in good agreement with the LHCb data [39] (see Table 5).
In Table 5 the masses M(2 3 P 0 ) = 3863 MeV, M(4 3 P 0 ) = 4473 MeV and M(5 3 P 0 ) = 4719 MeV, show very good agreement with those of χ c0 (3862) [16], X(4500) and X(4700) [39].At present other high excitations with J P = 1 + , 2 + (n = 4,5) are not yet found and their observation would be very important to understand the fine-structure effects of high charmonium, in particular, the fine-structure splitting have to decrease for a screened GE potential.
Notice that the resonance X(4700) lies very close to the ψ(2S)φ threshold and this fact indicates a possible connection between the cc and the cscs states.The four-quark interpretation of the X(4500), X(4700) was discussed 4941 abs in different models [19], [46]- [49], where in the mass region (4.4-4.8)GeV the radial or orbital excitations of a diquark-antidiquark systems can exist.

Conclusions
In our paper the scalar resonances X(3915), X(3960), X(4140) are assumed to be the four-quark states, produced due to recoupling mechanism, when one pair of mesons can transform into another pair of mesons infinitely many times.These resonances do not exist in the cc spectrum.As the four-quark states they have several specific features: 1.The resonance appears only in the S-wave decay channel.
2. Within the ERM it lies rather close to the lower threshold.
3. The scalar four-quark resonance can be created in two channel case due to transitions between channels, but it can also be coupled to another channel 3, e.g. the cc channel.
4. These resonances have no large sizes, being the compact systems, and this fact may be important for their observation.In the case of the X(3915) this statement is confirmed by the Belle analysis of the Q 2 distribution of the X(3915) → J/ψω decays in [50].
The masses and widths of the X(3915), X(3960), X(4140), presented in Table 4, are obtained in a good agreement with experiment.The authors are grateful to N. P. Igumnova for collaboration.

Table 4 :
The ERM predictions for the masses and widths (in MeV) of exotic resonances with J P C = 0 ++

Table 5 :
The Regge trajectory predictions for the masses of the charmonium n 3 P 0 states (in MeV)