Abstract
In this article, we choose the type tetraquark current to study the hadronic coupling constants in the strong decays , , , , , , , , with the QCD sum rules based on solid quarkhadron quality. The predicted width is in excellent agreement with the experimental data from the Belle collaboration, which supports assigning the to be the type tetraquark state with . In calculations, we observe that the hadronic coupling constants , which is consistent with the observation of the in the mass spectrum, and favors the molecule assignment. It is important to search for the process to diagnose the nature of the , as the decay is greatly suppressed.
Strong decays of the as an vector tetraquark state in solid quarkhadron duality
[2mm] ZhiGang Wang ^{1}^{1}1Email: zgwang@aliyun.com.
Department of Physics, North China Electric Power University, Baoding 071003, P. R. China
PACS number: 12.39.Mk, 12.38.Lg
Key words: Tetraquark state, QCD sum rules
1 Introduction
In 2007, the Belle collaboration observed the and in the invariant mass distribution with statistical significances and respectively in the precess between threshold and using of data collected with the Belle detector at KEKB [1]. In 2008, the Belle collaboration observed the in the invariant mass distribution with a significance of in the exclusive process with an integrated luminosity of at the KEKB [2].
In 2014, the Belle collaboration measured the cross section from to with the full data sample of the Belle experiment using the ISR technique, and determined the parameters of the and resonances and superseded previous Belle determination [3]. The masses and widths are shown explicitly in Table 1. Furthermore, the Belle collaboration studied the invariant mass distribution and observed that there are two clusters of events around the masses of the and corresponding to the and , respectively.
In the scenario of conventional twoquark states, the structures of the and in the ideal mixing limit can be symbolically written as,
(1) 
While in the scenario of tetraquark states, the structures of the and in the ideal mixing limit can be symbolically written as [4, 5, 6],
(2) 
In Ref.[7], we take the nonet scalar mesons below as the twoquarktetraquark mixed states and study their masses and pole residues with the QCD sum rules in details. We determine the mixing angles, which indicate that the dominant components are the twoquark components. The maybe have constituent. The decay takes place, if the and are the same particle, the decay is OkuboZweigIizuka suppressed, there should be some rescattering mechanism to account for the decay.
The threshold of the is from the Particle Data Group [8], which is just above the mass from the Belle collaboration [3]. The can be assigned to be a molecular state [9, 10, 11] or a hadrocharmonium [12]. Other assignments, such as a 2P tetraquark state [13], a state [14], a state [15], a ground state Pwave tetraquark state [16, 17, 18, 19, 20, 21, 22] are also possible.
In Table 2, we list out the predictions of the masses based on the QCD sum rules [10, 11, 16, 17, 18, 19, 20, 21, 22], where the , , and denote the scalar (), pseudoscalar (), axialvector () and vector () diquark states. From the Table, we can see that it is not difficult to reproduce the experimental value of the mass of the with the QCD sum rules. However, the quantitative predications depend on the quark structures, the input parameters at the QCD side, the pole contributions of the ground states, and the truncations of the operator product expansion.
In the QCD sum rules for the hiddencharm (or hiddenbottom) tetraquark states and molecular states, the integrals
(3) 
are sensitive to the energy scales , where the are the QCD spectral densities, the are the Borel parameters, the are the continuum thresholds parameters, the predicted masses depend heavily on the energy scales . In Refs.[19, 23, 24], we suggest an energy scale formula with the effective quark mass to determine the ideal energy scales of the QCD spectral densities. The formula enhances the pole contributions remarkably, we obtain the pole contributions as large as , the largest pole contributions up to now. Compared to the old values obtained in Ref.[19], the new values based on detailed analysis with the updated parameters are preferred [20]. The energy scale formula also works well in the QCD sum rules for the hiddencharm pentaquark states [25].
For the correlation functions of the hiddencharm (or hiddenbottom) tetraquark currents, there are two heavy quark propagators and two light quark propagators, if each heavy quark line emits a gluon and each light quark line contributes a quark pair, we obtain a operator , which is of dimension , we should take into account the vacuum condensates at least up to dimension in the operator product expansion.
In Refs.[19, 20, 21, 26], we study the mass spectrum of the vector tetraquark states in a comprehensive way by carrying out the operator product expansion up to the vacuum condensates of dimension , and use the energy scale formula or modified energy scale formula to determine the ideal energy scales of the QCD spectral densities in a consistent way. In the scenario of tetraquark states, we observe that the preferred quark configurations are the and . In this article, we choose the quark configuration .
In Ref.[27], we assign the to be the diquarkantidiquark type axialvector tetraquark state, study the hadronic coupling constants , , with the QCD sum rules by taking into account both the connected and disconnected Feynman diagrams in the operator product expansion. We pay special attentions to matching the hadron side of the correlation functions with the QCD side of the correlation functions to obtain solid duality. The routine works well in studying the decays [28].
In this article, we assign the to be the vector tetraquark state, and study the strong decays , , , , , , , , with the QCD sum rules based on the solid quarkhadron duality, and reexamine the assignment of the .
Year 
Mass (MeV)  Width (MeV)  Experiment  

2007 
Belle [1]  
Belle [1]  
2008 
Belle [2]  
2014 
Belle [3]  
Belle [3] 
Structures  OPE (No)  mass(GeV)  References  

[10]  
[11]  

[16]  
[17]  

[18]  
[20]  

[19]  
[20]  

[19]  
[20]  

[21]  

[22] 
The article is arranged as follows: we illustrate how to calculate the hadronic coupling constants in the twobody strong decays of the tetraquark states with the QCD sum rules in section 2, in section 3, we obtain the QCD sum rules for the hadronic coupling constants , , , , , , , ; section 4 is reserved for our conclusion.
2 The hadronic coupling constants in the twobody strong decays of the tetraquark states
In this section, we illustrate how to calculate the hadronic coupling constants in the twobody strong decays of the tetraquark states with the QCD sum rules. We write down the threepoint correlation functions ,
(4) 
where the currents interpolate the tetraquark states , the and interpolate the conventional mesons and respectively,
(5) 
the , and are the pole residues or the decay constants.
At the phenomenological side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators , , into the threepoint correlation functions and isolate the ground state contributions to obtain the result [29, 30],
(6)  
where , the are the hadronic coupling constants defined by
(7) 
the four functions , , and have complex dependence on the transitions between the ground states and the higher resonances or the continuum states.
We rewrite the correlation functions at the hadron side as
(8)  
through dispersion relation, where the are the hadronic spectral densities,
(9) 
where the and are the thresholds, the , , are the continuum thresholds.
Now we carry out the operator product expansion at the QCD side, and write the correlation functions as
(10) 
through dispersion relation, where the are the QCD spectral densities,
(11) 
However, the QCD spectral densities do not exist,
(12)  
because
(13) 
Thereafter we will write the QCD spectral densities as for simplicity.
We math the hadron side of the correlation functions with the QCD side of the correlation functions, and carry out the integral over firstly to obtain the solid duality [27],
the denotes the thresholds . Now we write down the quarkhadron duality explicitly,
No approximation is needed, we do not need the continuum threshold parameter in the channel. The channel and channel are quite different, we can not set the continuum threshold parameters in the channel as , i.e. we can not set in the present case, where the denote the , , , , because the contaminations from the excited states , , , are out of control.
We can introduce the parameters , , and to parameterize the net effects,
(16) 
In numerical calculations, we take the relevant functions and as free parameters, and choose suitable values to eliminate the contaminations from the higher resonances and continuum states to obtain the stable QCD sum rules with the variations of the Borel parameters.
If the are charmonium or bottomnium states, we set and perform the double Borel transform with respect to the variables and , respectively to obtain the QCD sum rules,
(17) 
where the and are the Borel parameters. If the are opencharm or openbottom mesons, we set and perform the double Borel transform with respect to the variables and , respectively to obtain the QCD sum rules,
(18) 
where .
3 The width of the as vector tetraquark state
Now we write down the threepoint correlation functions for the strong decays , , , , , , , , , respectively, and apply the method presented in previous section to obtain the QCD sum rules for the hadronic coupling constants , , , , , , , .
For the twobody strong decays , , the correlation function is
(19) 
where
(20) 
For the twobody strong decay , the correlation function is
(21) 
where
(22) 
For the twobody strong decay , the correlation function is
(23) 
where
(24) 
For the twobody strong decay , the correlation function is
(25) 
where
(26) 
For the twobody strong decay , the correlation function is
(27) 
where
(28) 
For the twobody strong decay , the correlation function is
(29) 
For the twobody strong decay , the correlation function is
(30) 
At the phenomenological side, we insert a complete set of intermediate hadronic states with the same quantum numbers as the current operators into the threepoint correlation functions and isolate the ground state contributions to obtain the hadron representation [29, 30],
For the decays , ,