-
The
X(6900) peak observed by the LHCb Collaboration in the di-J/ψ invariant mass spectrum [1, 2], and later in theJψψ(3686) invariant mass spectrum [3], has stimulated many discussions of theoretical aspects (see for example Ref. [4] for an incomplete list of references). Moreover,X(6900) is close to the threshold ofJ/ψψ(3770) ,J/ψψ2(3823) ,J/ψ ψ3(3842) , andχc0χc1 , whereasX(7200) is close to the threshold ofJ/ψψ(4160) andχc0 χc1(3872) ). Inspired by this, Ref. [5] studied the properties ofX(6900) andX(7200) by assuming theX(6900) coupling toJ/ψJ/ψ ,J/ψψ(3770) ,J/ψψ2(3823) ,J/ψψ3(3842) , andχc0χc1 channels and theX(7200) coupling toJ/ψJ/ψ ,J/ψψ(4160) , andχc0χc1(3872) channels. For the S-waveJ/ψJ/ψ coupling, the pole counting rule (PCR) [6], which has been applied to the studies of "XYZ " physics in Refs. [7−10], was employed to analyze the nature of the two structures. It was found that the di-J/ψ data alone are not sufficient to judge the intrinsic properties of the two states. It was also pointed out thatX(6900) is unlikely a molecule ofJ/ψψ(3686) [5], a conclusion drawn before the discovery of Ref. [3]. More recently, Refs. [4, 11, 12] investigatedX(6900) using a combined analysis of di-J/ψ andJ/ψψ(3686) data and concluded thatX(6900) cannot be aJ/ψψ(3686) molecule.Nevertheless, as already stressed in Ref. [4], even though
X(6900) is very unlikely a molecule ofJ/ψψ(3686) , this does not mean that it has to be an "elementary state" (i.e., a compactˉcˉccc tetraquark state). It was pointed out that it is possible thatX(6900) be a molecular state composed of other particles, such asJ/ψψ(3770) , which form thresholds closer toX(6900) if the channel coupling is sufficiently large1 .This note will discuss a possible mechanism for the enhancement of the
J/ψψ(3770) channel coupling. TheDˉD (orD∗ˉD∗ ) component insideψ(3770) may play an important role, so far ignored in the literature, in explaining theX(6900) resonant peak through the anomalous threshold emerged from the triangle diagram generated by the D (D∗ ) loop, as depicted in Fig. 1.Noticing that
ψ(3770) orψ″ couples dominantly toDˉD , we start from the Feynman diagram as depicted in Fig. 1 by assuming that it contributes toJ/ψψ″ elastic scatterings near theJ/ψψ″ threshold2 . Assuming an interaction Lagrangian3 L=−ig(D0∂μˉD0−ˉD0∂μD0)ψ″μ−ig(D+∂μD−−D−∂μD+)ψ″μ+g′D0ˉD0J/ψμJ/ψμ+g′D+D−J/ψμJ/ψμ ,
(1) after performing the momentum integration, the amplitude as depicted by Fig. 1 is
iM=(−16g2g′)(ϵ2⋅ϵ4)ϵμ1ϵν3∫dDk×kμkν((k−p1)2−m2D)((k−p3)2−m2D)(k2−m2D)≡(−16g2g′)(ϵ2⋅ϵ4)ϵμ1ϵν3Aμν
(2) where
Aμν=−i16π2∫10dx∫1−x0dy{pμ3pν1xyΔ−gμν4Γ(ϵ)1Δϵ}≡−i16π2{pμ3pν1B+∫10dx∫1−x0dy(−gμν4Γ(ϵ)1Δϵ)} ,
(3) and
B(t)=∫10dxx2M2(2(M2(1−2x)+tx)×ArcTan(M2−txΛ(t,x))−ArcTan(M2(2x−1)−txΛ(t,x))Λ(t,x)+lnm2+t(x−1)xm2+M2(x−1)x) .
(4) Here,
Λ(t,x)=√4m2M2−M4+4M2tx2−2M2tx−t2x2 ,Δ=M2(x2+y2)+(2M2−t)xy−M2(x+y)+m2 , andΓ(ϵ)1Δϵ=1ϵ−lnΔ−γ+ln4π+O(ϵ) . On the right hand side of Eq. (3), only theB term will be considered since the rest will be absorbed by the contact interactions to be introduced latter. M is the mass ofψ(3770) , and m is the mass ofD(ˉD) . Parameter g is the coupling strength of theψ″DˉD three point vertex, andg′ is the coupling strength of theJ/ψJ/ψDˉD four point vertex; {t=(p2−p4)2 }. Parameter g can be determined by theψ″→D0ˉD0 decay process ,iMψ″DD=igϵ(ψ″)⋅[p(D0)−p(ˉD0)] ,Γ=18π|iMψ″DD|2q(DD)M2ψ″=16πg2q(DD)3M2ψ″ ,
(5) where
q(DD) is the norm of the three-dimensional momentum ofD0 orˉD0 in the final state. The PDG valueΓψ″→D0ˉD0∼27.2×52%×10−3 [16] determinesg∼12 4 . Parameterg′ is unknown and is left as a free parameter. The amplitude, Eq. (4), contains a rather complicated singularity structure, especially the well-known anomalous threshold, which was discovered by Mandelstam who used it to explain the looseness of the deuteron wave function [18]. The anomalous threshold is located atsA=4m2−(M2−2m2)2m2 .
(6) Considering the mass of the
ψ″ andD0 mesons, one obtainssA=−1.28 GeV2 (for theD+ loop, –0.98 GeV2 ). Numerically, functionB is plotted in Fig. 2(a), where one clearly sees the anomalous threshold beside the normal one att=4m2 . Note that ifM2 is smaller than2m2 , the anomalous branch point is located below the physical threshold, but on the second sheet. It touches the physical threshold4m2 and turns up to the physical sheet if the value ofM2 increases to2m2 . With a further increase inM2 , the anomalous threshold moves towards the left on the real axis, passes the origin whenM2=4m2 , and finally reaches the physical value, i.e., –1.28 GeV2 . The situation is depicted in Fig. 2(b). Note thatsA here is negative, contrary to what occurs with the deuteron, because the latter is a bound state with a normalizable wave function, whereasψ″ is an unstable resonance.Figure 2. (color online) Left: triangle diagram contribution (the y axis label is arbitrary). Right: the trajectory of the anomalous threshold with respect to the variation of
M2 .To proceed, one needs to make the partial wave projection of
M and obtainTJμ1μ2μ3μ4(s)=132π12q2(s)∫0−4q2(s)dtdJμμ′×(1+t2q2(s))Mμ1μ2μ3μ4(t) ,
(7) where the channel momentum square reads
q2(s)=[s−(M−MJ)2][s−(M+MJ)2]/4s ,MJ is the mass ofJ/ψ ,μi denotes the corresponding helicity configuration, andμ=μ1−μ2 ,μ′=μ3−μ4 . The key observation is that the integral interval in Eq. (7) will coversA if√s>6.96 GeV (for theD+ loop,√s>6.94 GeV). In other words, the partial wave amplitude will be enhanced in the vicinity of theX(6900) peak by the anomalous threshold enhancement in the t channel, as can be observed in Fig. 3.Based on the above observation, it is suggested that the
X(6900) peak may at least partly be explained by the anomalous threshold generated by the triangle diagram depicted in Fig. 1. Furthermore, to obtain theL=0 (s wave) amplitudes, we need the relation between the s wave and helicity amplitudes (see Refs. [11, 12] for further discussions):T0L=0(s)=13[2T0++++(s)+2T0++−(s)−2T0++00(s)−2T000++(s)+T00000(s)],
(8) T2L=0(s)=115[T2++++(s)+T2++−(s)]+√615[T2+++−(s)+T2+−++(s)+T2++−+(s)+T2−+++(s)]+√315[T2+++0(s)+T2+0++(s)+T2++0+(s)+T20+++(s)+T2++−0(s)+T2−0++(s)+T2++0−(s)+T20−++(s)]+15[T2+00+(s)+T20++0(s)+T2+00−(s)+T20−+0(s)+T2+0+0(s)+T20+0+(s)+T20+0−(s)+T2+0−0(s)]+√25[T2+−+0(s)+T2+0+−(s)+T2+−0+(s)+T20++−(s)+T2−++0(s)+T2+0−+(s)+T2−+0+(s)+T20+−+(s)]+215[T2++00(s)+T200++(s)+T20000(s)]+2√615[T2+−00(s)+T200+−(s)]+25[T2+−+−(s)+T2+−+(s)]+2√315[T2+000(s)+T200+0(s)+T20+00(s)+T2000+(s)] ,
(9) In practice, it is found that the anomalous enhancement gives a more prominent effect to the
J=2 amplitude than theJ=0 amplitude. Furthermore, to estimate the triangle diagram contribution, a combined fit withJ/ψJ/ψ andJ/ψψ(3686) data is made. A coupled–channel K-matrix unitarization scheme is employed includingJ/ψJ/ψ ,J/ψψ(3686) , andJ/ψψ(3770) . The tree level amplitudes are also taken into account from the following contact interaction Lagrangian [19]:Lc=c1VμVαVμVα+c2VμVαVμV′α+c3VμV′αVμV′α+c4VμV′μVαV′α+c5VμVαVμV′′α+c6VμV′′αVμV′′α+c7VμV′′μVαV′′α+c8VμV′αVμV′′α+c9VμV′μVαV′′α,
(10) and these tree level amplitudes are as follows:
iMJ/ψJ/ψ→J/ψJ/ψ=i8c1(ϵ1μϵ2αϵ†μ3ϵ†α4+ϵ1μϵ2αϵ†α3ϵ†μ4+ϵ1μϵμ2ϵ†3αϵ†α4),
(11) iMJ/ψJ/ψ→J/ψψ(2S)=i2c2(ϵ1μϵ2αϵ†μ3ϵ′†α4+ϵ1μϵμ2ϵ†3αϵ′†α4+ϵ1αϵ2μϵ†μ3ϵ′†α4),
(12) iMJ/ψψ(2S)→J/ψψ(2S)=i4c3(ϵ1μϵ′2αϵ†μ3ϵ′†α4)+i2c4(ϵ1μϵ′μ2ϵ†3αϵ′†α4+ϵ1μϵ′α2ϵ†3αϵ′†μ4),
(13) iMJ/ψJ/ψ→J/ψψ(3770)=i2c5(ϵ1μϵ2αϵ†μ3ϵ′′†α4+ϵ1μϵμ2ϵ†3αϵ′′4†α+ϵ1αϵ3μϵ†μ3ϵ′′4),
(14) iMJ/ψψ(3770)→J/ψψ(3770)=i4c6(ϵ1μϵ′′2αϵ†μ3ϵ′′+α4)+i2c7(ϵ1μϵ′′μ2ϵ†3αϵ′′†α4+ϵ1μϵ′′α2ϵ†3αϵ′′†μ4),
(15) iMJ/ψψ(2S)→J/ψψ(3770)=ic9(ϵ1μϵ′μ2ϵ†3αϵ′′†α4+ϵ1αϵ′μ2ϵ†3μϵ′′†α4)+i2c8(ϵ1μϵ′2αϵ†μ3ϵ′′α4).
(16) After the same partial wave projection process of Eqs. (7) – (9), the coupled–channel partial wave amplitudes at the tree level,
MJ,ijL=0(J=0,2 andi,j=1,2,3 ), are determined. By taking into account K-matrix unitarization and final state interaction, the unitarized partial wave amplitude isFJi(s)=3∑k=1αk(s)TJ,kiL,U(s) ,
(17) where
TJL,U(s)=[1−iρKJL]−1 .
(18) αk(s) is the real polynomial function in general and is set to be constant here, and we setα1(s)2=1 . Particularly, as for the casei=j=3 andJ=2 , the triangle diagram needs to be taken into consideration. That is to sayKJ=2,i=3j=3L=0=MJ=2,i=3j=3L=0,tree+TJ=2triangle , in which the triangle diagram contribution comes from Fig. 15 . For other cases,KJ,ijL=0=MJ,ijL=0,tree . Further, to fit the experimental data, one hasdEventsid√s=Nipi(s)|Fi|2,
(19) where
pi(s) refers to the abs of three momenta for the corresponding channel. According to partial wave convention, for theJ=0 andJ=2 case, they have a total scale factor [11, 12]:|Fi|2=|FJ=0i|2+5|FJ=2i|2 .
(20) The fit is overdone since there are many parameters. One solution is shown in Fig. 4, and the fit parameters are listed in Table 1 for illustration. In this fit, we set
c1,c2,⋯,c4 to be negligible simply because they are not directly related to theJ/ψ(3770) channel and the fit can be performed reasonably well without them. The error band in Fig. 4 is rather large; this is due to the two normalization factors,N1 andN2 , which contain rather large error bars.Parameter χ2/d.o.f g′ N1 N2 α2 α3 Fit 0.93 −18.3±1.6 23±10 0.34±0.23 −4.94±0.02 3.97±0.05 c5 c6 c7 c8 c9 11.34±0.05 50.79±0.05 −35.01±0.19 −64.71±0.07 1.529±0.002 Table 1. Fit parameters of Fig. 4. The
ci parameters are defined in Eq. (10). The errors are statistical only.During the fit, many solutions exist. Nevertheless, it was found that the triangle diagram contributions are all small. This behavior is unclear; however, one possible reason could be that the peak position through the anomalous threshold contribution, as shown in Fig. 3, is approximately 60 – 80 MeV above the
X(6900) peak; hence, the fit becomes difficult. One possible way to resolve the problem is to adopt another parameterization in which the background contributions are more flexible to be tuned. Thus, the interference between the background and the anomalous enhancement can lead to the shift of the peak position by a few tens of MeV. Another possible mechanism for the suppression of the triangle diagram is that theψ(3770)DˉD vertex is in the p-wave form; hence, it may provide another suppression factor due to (non-relativistic) power counting. [20]. We defer this investigation for future studies.
HTML
-
We would like to thank De-Liang Yao and Ling-Yun Dai for the very helpful discussions.