# Tetraquark interpretation of Belle data and BaBar data

###### Abstract

We summarize the main features of the spectroscopy, production and decays of the tetraquarks in the sector, concentrating on the lowest state called . The tetraquark framework is used to analyze the BaBar data on the cross section ( energy scan) between and 11.20 GeV and the Belle data on the processes near the peak of the resonance. The BaBar energy scan is consistent with an additional state at a mass of 10.90 GeV and a width of about 28 MeV, in broad agreement with the state GeV seen by Belle in the exclusive final states. We argue that the decay widths and the dipion invariant mass distributions measured by Belle are naturally explained by the tetraquark interpretation of .

Tetraquark interpretation of Belle data and BaBar data

Ahmed Ali^{†}^{†}thanks: Speaker.

Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, D-22607 Hamburg

E-mail: ahmed.ali@desy.de

\abstract@cs

## 1 Introduction

Experiments at the factories and Tevatron have in the past several years revived the interest in the spectroscopy of the Quarkonium-like exotic states. Labeled tentatively as , and , due to a lack of consensus on their interpretation, they have masses above the open charm threshold, with the being the lightest and the heaviest state observed so far [1]. There is also evidence for an exotic bound state having the quantum numbers , first observed in the initial state radiation (ISR) process , where is the scalar state. In the sector, Belle [2] has observed enhanced production for the processes in the center-of-mass energy between GeV and GeV, which does not agree with the conventional line shape [3]. The enigmatic features of the Belle data are the anomalously large decay widths for the mentioned final states and the dipion invariant mass distributions, which are strikingly different from the conventional QCD expectations for such dipionic transitions. A fit of the Belle data, using a Breit-Wigner resonance, yields a mass of MeV and a width of MeV [2]. This particle is given the tentative name . In [4, 5], is interpreted as a tetraquark state, which is a linear superposition of the flavour eigenstates and . The mass eigenstates (for the lighter) and (for the heavier) of the two are almost degenerate, with their small mass difference arising from isospin-breaking [6]. A dynamical model for the decay mechanisms of and the final state distributions measured by Belle was developed in [4] and refined in [5], yielding good fits of the Belle data. One anticipates that is also visible in the energy scan of the cross section, which was undertaken by the BaBar collaboration between GeV and 11.20 GeV [7]. A fit of the BaBar data on -scan is consistent with a structure around and yields a better than the fits without the tetraquark states. More data are required to resolve this and related structures in the line shape. This contribution summarizes the work done in [4, 5, 6] interpreting the Belle [2] and BaBar [7] data in terms of the tetraquark states.

## 2 Spectrum of bottom diquark-antidiquark states

The mass spectrum of tetraquarks with , , and can be described in terms of the constituent diquark masses, , spin-spin interactions inside the single diquark, spin-spin interaction between quark and antiquark belonging to two diquarks, spin-orbit, and purely orbital term [8], i.e., with a Hamiltonian

(2.0) |

where:

(2.0) |

Here is the coupling of the spin-spin interaction between the quarks inside the diquarks, are the spin-spin couplings ranging outside the diquark shells, is the spin-orbit coupling of diquark and characterizes the contribution of the total angular momentum of the diquark-antidiquark system to its mass.

The parameters involved in the above Hamiltonian (2.0) can be obtained
from the known meson and baryon masses by resorting to the constituent quark
model [9]:
,
where the sum runs over the hadron constituents. The coefficient depends on the flavour of the constituents , and on the particular
colour state of the pair. Using the entries in the PDG [3] for hadron masses along
with the assumption that the spin-spin interactions are independent of
whether the quarks belong to a meson or a diquark, the results for the
masses corresponding to the tetraquarks () were
calculated in [6]. The lowest eight
tetraquark states (, which
are all orbital excitations with , have the following spin and
orbital angular momentum eigenvalues:
,
,
,
and

.
Identifying the lowest lying state
with the measured by Belle, and using the estimates for the
other parameters entering in Eq. (2.0), fixes the diquark mass
GeV. The uncertainties on the masses of the other
six states () are higher, as they depend in addition on
the mass-splittings between the good and bad diquarks,
, estimated as
MeV [10, 11]. The central values of their masses
are:
MeV, MeV, and MeV.
Assuming isospin symmetry, the states and are degenerate for
each . Including isospin-symmetry breaking lifts this degeneracy with the mass difference
between the lighter and the heavier of the two states estimated as MeV, where is a mixing angle and the
mass eigenstates are defined as:
and [6].
The resulting mass differences are small. However, depending on , the
electromagnetic couplings of the and may turn out to be
significantly different from each other, and hence also their contributions to .

## 3 Decay Widths of and other tetraquarks

As the masses of all the eight tetraquark states lie above the thresholds for the decays , they decay readily into these final states. For the state (having a mass of 11257 MeV), also the decay is energetically allowed. In [6], the decay widths have been estimated (up to a tetraquark hadronic size parameter ) in terms of the corresponding partial decay widths , which can be calculated with the help of the entries in the PDG [3]. Specifically, the following relations are assumed

(3.0) |

and likewise for the and decays. Noting that the decays , as well as the decays are Zweig-forbidden, one expects, concentrating on the lowest mass state, . Using the PDG value [3] MeV, we get MeV for the total decay widths. Equating this decay width to the measured value of the total decay width MeV by Belle [2], one gets . This suggests that the tetraquarks have a hadronic size of the same order as that of the . The hadronic widths of the other tetraquarks are estimated as [6]: MeV, MeV and MeV.

To calculate the production cross sections, we have derived the corresponding Van Royen-Weisskopf formula for the leptonic decay widths of the tetraquark states made up of point-like diquarks [5]:

(3.0) |

where is the fine-structure constant, , are the diquark charges in units of the proton electric charge, and GeV [12] is the square of the derivative of the radial wave function for taken at the origin. Hence, the leptonic widths of the tetraquark states are estimated as [5]

(3.0) |

which are substantially smaller than the leptonic width of the [3]. The electronic widths of the mass eigenstates and depend, in addition, on the mixing angle .

## 4 Analysis of the BaBar data on -scan

BaBar has reported the cross section measured in a dedicated energy scan in the range GeV and GeV taken in steps of 5 MeV [7]. Their measurements are shown in Fig. 1 (left-hand frame) together with the result of the BaBar fit which contains the following ingredients: A flat component representing the -continuum states not interfering with resonant decays, called , added incoherently to a second flat component, called , interfering with two relativistic Breit-Wigner resonances, having the amplitudes , and strong phases, and , respectively. Thus,

(4.0) | |||||

with . The results summarized in their Table II for the masses and widths of the and differ substantially from the corresponding PDG values [3], in particular, for the widths, which are found to be MeV for the , as against the PDG value of MeV, and MeV for the , as compared to MeV in PDG. As the systematic errors from the various thresholds are not taken into account, this mismatch needs further study. The fit shown in Fig. 1 (left-hand frame) is not particularly impressive having a of approximately 2.

The BaBar -data is refitted in [6] by modifying the model in Eq. (4.0) by taking into account two additional resonances, corresponding to the masses and widths of and . Thus, formula (4.0) is extended by two more terms

(4.0) |

which interfere with the resonant amplitude and the two resonant amplitudes for and shown in Eq. (4.0). Using the same non-resonant amplitude and as in the BaBar analysis [7]. the resulting fit is shown in Fig. 1 (right-hand frame). Values of the best-fit parameters yield the masses of the and and their respective full widths which are almost identical to the values obtained by BaBar [7]. However, quite strikingly, a third resonance is seen in the -line-shape at a mass of GeV, tantalisingly close to the -mass in the Belle measurement of the cross section for , and a width of about MeV. In the region around 11.15 GeV, where the states are expected, our fits of the BaBar -scan do not show a resonant structure due to the larger decay widths of the states . The resulting with the 4 Breit-Wigners shown in Fig. 1 (right frame) is better than that of the BaBar fit [7].

The quantity is given by the ratio of the two amplitudes and , which also fixes the mixing angle . From the fit shown in the right-hand frame in Fig. 1, one obtains

(4.0) |

yielding

(4.0) |

for the mixing angle and the mass difference between the eigenstates, respectively. For the mass eigenstates and , the electronic widths and are given by [5] keV. With the above determination of and , we get

(4.0) |

## 5 Analysis of the Belle data on

With the for both and , the dipionic final state is allowed to have the quantum numbers and . There are only three low-lying resonances in the PDG which can contribute as intermediate states, namely, the two states, and , which we take as the lowest tetraquark states, and the meson state . All three states contribute for the final state . However, kinematics allows only the in the final state . In addition, a non-resonant contribution with a significant -wave fraction is needed by the data on these final states. This model accounts well the shape of the measured distributions, as shown in Fig. 2 for (left-hand frames) and for (right-hand frames). As the decays are Zweig-allowed, one expects larger decay widths for these transitions, typically of MeV [5], than the decay widths for the conventional dipionic transitions, such as , which are of order 1 keV [3]. Further tests of the tetraquark hypothesis involving the processes are presented in [5].

## References

- [1] S. L. Olsen, Nucl. Phys. A 827, 53C (2009) [arXiv:0901.2371 [hep-ex]]; A. Zupanc [for the Belle Collaboration], arXiv:0910.3404 [hep-ex].
- [2] K. F. Chen et al. [Belle Collaboration], Phys. Rev. Lett. 100, 112001 (2008); I. Adachi et al. [Belle Collaboration], arXiv:0808.2445 [hep-ex].
- [3] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008).
- [4] A. Ali, C. Hambrock and M. J. Aslam, Phys. Rev. Lett. 104, 162001 (2010).
- [5] A. Ali, C. Hambrock and S. Mishima, DESY 10-076 [arXiv:1011.4856 [hep-ph]].
- [6] A. Ali, C. Hambrock, I. Ahmed and M. J. Aslam, Phys. Lett. B 684, 28 (2010).
- [7] B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. 102, 012001 (2009).
- [8] N. V. Drenska, R. Faccini and A. D. Polosa, Phys. Lett. B 669, 160 (2008); N. V. Drenska, R. Faccini and A. D. Polosa, Phys. Rev. D 79, 077502 (2009).
- [9] A. De Rujula, H. Georgi and S. L. Glashow, Phys. Rev. D 12, 147 (1975).
- [10] R. L. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003).
- [11] R. L. Jaffe, Phys. Rept. 409, 1 (2005) [Nucl. Phys. Proc. Suppl. 142, 343 (2005)].
- [12] E. J. Eichten and C. Quigg, Phys. Rev. D 52, 1726 (1995).