## Brazilian Journal of Physics

*Print version* ISSN 0103-9733

### Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007

#### http://dx.doi.org/10.1590/S0103-97332007000100020

**X (3872) in QCD sum rules**

**R. D. Matheus**

Instituto de Física, Universidade de São Paulo, C.P. 66318, 05315-970, São Paulo, SP, Brazil

**ABSTRACT**

QCD spectral sum rules is used to test the nature of the meson X(3872), assumed to be an exotic four-quark $(c\bar{c}q\bar{q})$ state with *J ^{PC}* = 1

^{++}. For definiteness, the current proposed recently by Maiani et al [1] is used, at leading order in a

_{s}, considering the contributions of higher dimension condensates. The value

*M*= (3.94±0.11) GeV is found which is compatible, within the errors, with the experimental candidate X(3872). The uncertainties of our estimates are mainly due to the one from the

_{X}*c*quark mass.

**Keywords:** QCD sum rules; Tetraquarks; Axial vector mesons

**I. INTRODUCTION**

In august 2003, there's been a report by BELLE[2] of a narrow resonance in *B*^{+ }® *X*(3872)*K*^{+ }® *J*/y p^{+}p^{-} *K*^{+} decays. It has been soon after confirmed by D0 [3], CDF II [4] and BABAR [5]. The collective data from these experiments give an averaged mass of *m _{X}* = 3871,9 ±0,5 MeV [6]. The width has been stablished by BELLE at a upper limit of 2,3 MeV at 90% confidence level and the most probable quantum numbers are:

*J*= 1

^{PC}^{++}[6].

BELLE has also found the resonance in the channel: *X*(3872) ® J/y p^{+}p^{-}p^{0} . The relative strength of this decay mode is given by[7]:

This strong isospin violation makes it difficult to understand the *X* as a charmoniun. Besides, the charmoniun spectrum state with correct quantum numbers, (3925), has both the mass and the decay width ( » 16 MeV) too high to be gracefully identified with the observed resonance[6].

The anomalous nature of the *X* has led to many speculations: tetraquark [1, 8], cusp [9], hybrid [10], or glueball [11]. Another explanation is that the *X*(3872) is a *D*^{*} bound state [12-16], as predicted before its discovery.

In this work the QCD spectral sum rules (the Borel/Laplace Sum Rules (SR) [17-19] and Finite Energy Sum Rules (FESR) [19-21]) will be used to study the two-point functions of the axial vector meson, *X*(3872), assumed to be a four-quark state. In previous calculations, the Sum Rule (SR) approach was used to study the light scalar mesons [22-25] and the (2317) meson [26, 27], considered as four-quark states and a good agreement with the experimental masses was obtained. However, the tests were not decisive as the usual quark-antiquark assignments also provide predictions consistent with data and more importantly with chiral symmetry expectations [19, 23, 28, 29]. In the four-quark scenario, scalar mesons can be considered as S-wave bound states of diquark-antidiquark pairs, where the diquark was taken to be a spin zero color anti-triplet. Ref. [1] will be followed here, and the *X*(3872) will be considered as a *J ^{PC}* = 1

^{++}state with the symmetric spin distribution: [cq]

_{S = 1}[]

_{S = 0}+[cq]

_{S = 0}[]

_{S = 1}. Therefore, the corresponding lowest-dimension interpolating operator for describing

*X*is given by:

_{q} where *a, b, c*, ... are color indices, *C* is the charge conjugation matrix and *q* denotes a *u* or *d* quark.

**II. THE QCD EXPRESSION OF THE TWO-POINT CORRELATOR**

The SR are constructed from the two-point correlation function

Since the axial vector current is not conserved, the two functions, P_{1} and P_{0}, appearing in Eq. (3) are independent and have respectively the quantum numbers of the spin 1 and 0 mesons.

The fundamental assumption of the sum rules approach is the principle of duality. Specifically, we assume that there is an interval over which the correlation function may be equivalently described at both the quark and the hadron levels. Therefore, on one hand, we calculate the correlation function at the quark level in terms of quark and gluon fields. On the other hand, the correlation function is calculated at the hadronic level introducing hadron characteristics such as masses and coupling constants. At the quark level, the complex structure of the QCD vacuum leads us to employ the Wilson's operator product expansion (OPE). The calculation of the phenomenological side proceeds by inserting intermediate states for the meson *X*. Parametrizing the coupling of the axial vector meson 1^{++}, *X*, to the current, *j*_{µ}, in Eq. (2) in terms of the meson decay constant *f _{X}* as:

the phenomenological side of Eq. (3) can be written as

where the Lorentz structure projects out the 1^{++} state. The dots denote higher axial-vector resonance contributions that will be parametrized, as usual, through the introduction of a continuum threshold parameter *s*_{0}.

The OPE side will be evaluated at leading order in a_{s} and the contributions of condensates will be considered up to dimension five.

The correlation function, P_{1}, in the OPE side can be written as a dispersion relation:

where the spectral density is given by the imaginary part of the correlation function: pr(*s*) = Im[(*s*)]. After making an inverse-Laplace (or Borel) transform of both sides, and transferring the continuum contribution to the OPE side, the sum rule for the axial vector meson *X* up to dimension-five condensates can be written as:

where r(*s*) is the sum of the r^{Dim}(*s*) below:

where the integration limits are given by a_{min} = (1-)/2, a_{max} = (1+)/2 and (b_{min} = a)/( *s*a-). The dominant contributions from the dimension-five condensates have also been included :

**III. LSR PREDICTIONS OF M_{X}**

In order to extract the mass *M _{X}* without worrying about the value of the decay constant

*f*, one must take the derivative of Eq. (7) with respect to 1/

_{X}*M*

^{2}, divide the result by Eq. (7) and obtain:

In the numerical analysis of the sum rules, the values used for the quark masses and condensates are (see e.g. [19, 30-32]): *m _{c}* = (1.23 ±0.05) GeV,

*m*= 2.3 MeV ,

_{u}*m*= 6.4 MeV,

_{d}*m*= (

_{q}*m*+

_{u}*m*)/2 = 4.3 MeV, á

_{d}*q*ñ = -(0.23±0.03)

^{3}GeV

^{3}, ágs.Gqñ = áqñ with = 0.8 GeV

^{2}and ág

^{2}G

^{2}ñ = 0.88 GeV

^{4}. The sum rules were evaluated in the range 1.6

__<__

*M*

^{2}

__<__2.8 for three values of s

_{0}: = 4.1 GeV, = 4.3 GeV and = 4.5 GeV.

Comparing the relative contribution of each term in Eqs. (8) and (9), to the right hand side of Eq. (7) there is a quite good OPE convergence for *M*^{2} > 1.9 GeV^{2}. This analysis allows us to determine the lower limit constraint for *M*^{2} in the sum rules window [33].

The upper limit constraint for *M*^{2} is determined by imposing that the QCD continuum contribution should be smaller than the pole contribution. The maximum value of *M*^{2} for which this constraint is satisfied depends on the value of *s*_{0}. The comparison between pole and continuum contributions gives *M*^{2} < 3.0 for = 4.5 GeV, *M*^{2} < 2.6 for = 4.3 GeV and *M*^{2} < 2.3 GeV^{2} for = 4.1 GeV.

Figure 1 shows the *X* meson mass obtained from Eq. (10), in the relevant sum rules window, with the upper and lower validity limits indicated. From Fig. 1 we see that the results are reasonably stable as a function of *M*^{2}. In the numerical analysis the range of *M*^{2} values from 2.1 GeV^{2} until the one allowed by the sum rule window criteria shall then be considered for each value of *s*_{0}.

Varying the QCD parameters inside the limits shown before and taking into account the range of values for *M*^{2} and *s*_{0} we get:

The error is due to the combined effect of *M*^{2}, *s*_{0}, áqñ, and *m _{c}* (

*m*having the greatest effect).

_{c}

**IV. FESR PREDICTION FOR M_{X}**

As an alternative, one can use the FESR, which can be obtained from Eq. (7) by taking the limit 1/*M*^{2} ® 0 and equating the same power in 1/*M*^{2} in the two sides of the sum rules to get *n* equations:

Finally, dividing two subsequent equations (with *n* and *n*+1) one can obtain the mass *M _{X}* for any chosen value of

*n*(which, formally, is expected to be the same for any

*n*):

In contrast to the previous method, the FESR have the advantage of giving correlations between the mass and the continuum threshold *s*_{0}, which can be used to avoid inconsistencies in the determination of these parameters. Ideally, one looks at a minimum in the function *M _{X}*(

*s*

_{0}), which would provide a good criteria for fixing both

*s*

_{0}and

*M*. The results for different values of

_{X}*n*are very similar, therefore, in Fig. 2, only the results for

*n*= 0 and

*n*= 1 are shown. One can see in Fig. 2 that there is no stability in

*s*

_{0}, which presumably indicates the important role of the QCD continuum in the analysis.

The FESR results agree with LSR ones only for a small value of *s*_{0}: » 4.2 GeV, but since no stability was found we consider the LSR to be more reliable in this case.

**V. CONCLUSIONS**

We have presented a QCD spectral sum rules analysis of the two-point function of the *X*(3872) meson considered as a four quark state. We find that the sum rules result in Eq. (11) is compatible with experimental data. An improvement of this result needs an accurate determination of *m _{c}*.

Once the mass of the *X*(3872) is understood, it remains to explain why it is so narrow. There are presumably many multiquark states, but most of them are very broad and cannot be singled out from the continuum. In a recent study [34], based on the same interpolating field as the one used here, it was shown that, in order to explain the small width of the *X*(3872), one has to choose a particular set of diagrams contributing to its decay. However, it will be desirable if this investigation can be checked from alternative approaches, like e.g. lattice calculations. If confirmed, this method can be straightforwardly repeated to a variety of currents for understanding the width and the internal structure of the *X*(3872).

**Acknowledgments**

The author acknowledges M. Nielsen, S. Narison and J.-M. Richard for fruitful discussions. This work has been supported by FAPESP.

[1] L. Maiani, F. Piccinini, A.D. Polosa, and V. Riquer, Phys. Rev. D **71**, 014028 (2005). [ Links ]

[2] S.-L. Choi et al. [Belle Collaboration], Phys. Rev. Lett. **91**, 262001 (2003). [ Links ]

[3] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. **93**, 162002 (2004). [ Links ]

[4] D. Acosta et al. [CDF II Collaboration], Phys. Rev. Lett. **93**, (2004) 072001. [ Links ]

[5] B. Aubert et al. [ BaBar Collaboration], Phys. Rev. D **71**, 071103 (2005). [ Links ]

[6] Eric S. Swanson, Phys. Rept. **429**, 243 (2006). [ Links ]

[7] K. Abe et al., ''Evidence for X(3872) ® gamma J/psi and the sub-threshold decay X(3872) ® omega J/psi'', [hep-ex/0505037]. [ Links ]

[8] H. Høgaasen, J.-M. Richard, and P. Sorba, Phys. Rev. D **73**, 054013 (2006). [ Links ]

[9] D. Bugg, Phys. Lett. B **598**, 8 (2004). [ Links ]

[10] B.-A. Li, Phys. Lett. B **605**, 306 (2005). [ Links ]

[11] K. K. Seth, Phys. Lett. B **612**, 1 (2005). [ Links ]

[12] N. Törnqvist, [hep-ph/0308277]. [ Links ]

[13] F. E. Close and P. Page, Phys. Lett. B **578**, 119 (2004). [ Links ]

[14] C.Y. Wong, Phys. Rev. C **69**, 055202 (2004). [ Links ]

[15] S. Pakvasa and M. Suzuki, Phys. Lett. B **579**, 67 (2004). [ Links ]

[16] E. S. Swanson, Phys. Lett. B **588**, 189 (2004); [ Links ]Phys. Lett. B **598**, 197 (2004); [ Links ]Phys. Rept. **429**, 243 (2006). [ Links ]

[17] M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B **147**, 385 (1979). [ Links ]

[18] L.J. Reinders, H. Rubinstein, and S. Yazaki, Phys. Rept. **127**, 1 (1985). [ Links ]

[19] S. Narison, *QCD as a theory of hadrons, Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol.* **17**, (2002) 1-778 [hep-h/0205006]; * [ Links ]QCD spectral sum rules, World Sci. Lect. Notes Phys.* **26**, 1-527 (1989). [ Links ]

[20] R. A. Bertlmann, C. A. Dominguez, G. Launer, and E. de Rafaël, Nucl. Phys. B **250**, 61 (1985); [ Links ]R.A. Bertlmann, G. Launer, and E. de Rafael, Z. Phys. C **39**, 231 (1988). [ Links ]

[21] K. Chetyrkin, N.V. Krasnikov, and A.N. Tavkhelidze, Phys. Lett. B **76**, 83 (1978); [ Links ]N.V. Krasnikov, A.A. Pivovarov, and A.N. Tavkhelidze, Z. Phys. C **19**, 301 (1983). [ Links ]

[22] J. Latorre and P. Pascual, J. Phys. G **11**, L231 (1985). [ Links ]

[23] S. Narison, Phys. Lett. B **175**, 88 (1986). [ Links ]

[24] T. V. Brito, F. S. Navarra, M. Nielsen, and M. E. Bracco, Phys. Lett. B **608**, 69 (2005). [ Links ]

[25] H-J. Lee and N.I. Kochelev, Phys. Lett. B **642**, 358 (2006). [ Links ]

[26] M. E. Bracco, A. Lozea, R. D. Matheus, F. S. Navarra, and M. Nielsen, Phys. Lett. B **624**, 217 (2005). [ Links ]

[27] H.-Y. Cheng, W.S. Hou, Phys. Lett. B **566**, 193 (2003); [ Links ]A. Datta, P. J. O'Donnell, Phys. Lett. B **572**, 164 (2003). [ Links ]P. Colangelo, F. De Fazio, R. Ferrandes, Mod. Phys. Lett. A **19**, 2083 (2004); [ Links ]A. Hayashigaki, K. Terasaki, hep-ph/0411285; [ Links ]D. Becirevic, S. Fajfer, and S. Prelovsek, Phys. Lett. B **599**, 55 (2004); [ Links ]Z.-G. Wang, Shao-Long Wan, hep-ph/0602080; [ Links ]J.P. Pfannmoeller , hep-ph/0608213 [ Links ]

[28] S. Narison , hep-ph/0512256; [ Links ]S. Narison, Phys. Lett. B **605**, 319 (2005); [ Links ]S. Narison, Phys. Lett. B **210**, 238 (1988). [ Links ]

[29] W.A. Bardeen *et al.*, Phys. Rev. D **68**, 054024 (2003); [ Links ]S. Godfrey, Phys. Lett. B **568**, 254 (2003); [ Links ]M. Harada *et al.*, Phys. Rev. D **70**, 074002 (2004); [ Links ]Y.B. Dai *et al.* Phys. Rev. D **68**, 114011 (2003). [ Links ]

[30] S. Narison, Nucl. Phys. Proc. Suppl. **86**, 242 (2000); [ Links ]S. Narison, hep-ph/0202200; [ Links ]S. Narison, Phys. Rev. D **74**, 034013 (2006); [ Links ]S. Narison, Phys. Lett. B **341**, 73 (1994); [ Links ]H.G. Dosch and S. Narison, Phys. Lett. B **417**, 173 (1998); [ Links ]S. Narison, Phys. Lett. B **216**, 191 (1989). [ Links ]

[31] S. Narison, Phys. Lett. B **466**, 345 (1999). [ Links ]

[32] S. Narison, Phys. Lett. B **361**, 121 (1995); [ Links ]S. Narison, Phys. Lett. B **387**, 162 (1996). [ Links ]S. Narison, Phys. Lett. B **624**, 223 (2005). [ Links ]

[33] R. D. Matheus, S. Narison, M. Nielsen, and J. M. Richard, Phys. Rev. D **75**, 014005 (2007). [ Links ]

[34] F. S. Navarra, M. Nielsen, Phys. Lett. B **639**, 272 (2006). [ Links ]

Received on 29 September, 2006