Abstract

In this work, the mass and pole residue of resonance is studied by using QCD sum rules approach at finite temperature. Resonance is described by a diquark-antidiquark tetraquark current, and contributions to operator product expansion are calculated by including QCD condensates up to dimension six. Temperature dependencies of the mass and the pole residue are investigated. It is seen that near a critical temperature , the values of and decrease to 87% and to 44% of their values at vacuum.

1. Introduction

Heavy quarkonia systems provide a unique laboratory to search the interplay between perturbative and nonperturbative effects of QCD. They are nonrelativistic systems in which low energy QCD can be investigated via their energy levels, widths, and transition amplitudes [1]. Among these heavy quarkonia states, vector charmonium and bottomonium sectors are experimentally studied very well, since they can be detected directly in annihilations. In the past decade, observation of a large number of bottomonium-like states in several experiments increased the interest in these structures [26]. However, these observed states could not be conveniently explained by the simple picture of mesons. The presumption of hadrons containing quarks more than the standard quark content ( or ) is introduced by a perceptible model for diquarks plus antidiquarks, which was developed by Jaffe in 1976 [7]. Later Maiani, Polosa, and their collaborators proposed that the X, Y, Z mesons are tetraquark systems, in which the diquark-antidiquark pairs are bound together by the QCD color forces [8]. In this color configuration, diquarks can play a fundamental role in hadron spectroscopy. Thus, probing the multiquark matter has been an intensely intriguing research topic in the past twenty years and it may provide significant clues to understand the nonperturbative behavior of QCD.

In 2007, Belle reported the first evidence of and first observation for decays near the peak of the state at [2]. Assigning these signals to , the partial widths of decays and were measured unusually larger (more than two orders of magnitude) than formerly measured decay widths of states. Following these unusually large partial width measurement, Belle measured the cross sections of and and reported that the resonance observed via these decays does not agree with conventional line shape. These observations led to the proposal of existence of new exotic hidden-beauty state analogous to broad resonance in the charmonium sector, which is a Breit-Wigner shaped resonance with mass () MeV/, and width () MeV/, and is called [5]. In literature, there are several approaches to investigate the structure of exotic resonance. In [9], is considered as a bound state with a highly large binding energy. In [10, 11], is interpreted as a tetraquark and its mass is estimated by using QCD sum rules at vacuum.

Moreover, it is likely that at very high temperatures within the first microseconds following the Big Bang, quarks and gluons existed freely in a homogenous medium called the quark-gluon plasma (QGP). One of the first quark-gluon plasma signals proposed in the literature is suppression of particles [12]. In 2011, CMS collaboration reported that charmonium states and melt or were suppressed due to interacting with the hot nuclear matter created in heavy-ion interactions [13, 14]. Following these observations in the charmonium sector, CMS collaboration also reported suppression of bottomonium states, and relative to the ground state [15, 16]. The dissociation temperatures for the states are expected to be related to their binding energies and are predicted to be and for the , and mesons, respectively, where is the critical temperature for deconfinement [1618].

In addition to these studies, one of the first investigations on thermal properties of exotic mesons was done in [19], in which the authors investigated in medium properties of under the hypothesis that it is a state or state, but making no assumption on its structure. They estimated that the mass of the molecular state decreases with increasing temperature; the mass of charmonium or tetraquark state is almost stable.

Inspired by these findings and motivated by the aforementioned discussions, we focus on the resonance and its thermal behavior. This paper is organized as follows. In Section 2, theoretical framework of Thermal QCD sum rules (TQCDSR) and its application to are presented, and obtained analytical expressions of the mass and pole residue of are given up to dimension six operators. Numerical analysis is performed and results are obtained in Section 3. Concluding remarks are discussed in Section 4. The explicit forms of the spectral densities are written in Appendix.

2. Finite Temperature Sum Rules for Tetraquark Assignment

QCD sum rules (QCDSR) approach is based on Wilson’s operator product expansion (OPE) which was adapted by Shifman, Vainshtein, and Zakharov [24] and applied with remarkable success to estimate a large variety of properties of all low-lying hadronic states [2527]. Later, this model is extended to its thermal version that is firstly proposed by Bochkarev and Shaposnikov and led to many successful applications in QCD at [2835]. Very recently, in [36], the authors claimed that any QCDSR study on tetraquark states should contain contributions to OPE, which are unknown. It is a very important and strong argument for studying tetraquark states within QCDSR; however, those terms and their thermal behaviors are not known, and their calculation is beyond the scope of this work. Thus we follow the traditional sum rules to investigate the thermal behavior of hadronic parameters of , which were used very successfully in predicting properties of tetraquark states as well. In this work, we proceed with traditional sum rules analysis by using a tetraquark current, following several successful applications to exotic hadrons [3742].

In this section, the mass and pole residue of the exotic resonance are studied by interpreting it as a bound tetraquark via TQCDSR technique which starts with the two point correlation functionwhere represents the hot medium state, is the interpolating current of the state, and denotes the time ordered product. The thermal average of any operator in thermal equilibrium is given aswhere is the QCD Hamiltonian, and is inverse of the temperature, and is the temperature of the heat bath. Chosen current must contain all the information of the related meson, like quantum numbers, quark contents and so on. In the diquark-antidiquark picture, tetraquark current interpreting can be chosen as [10]where , is the charge conjugation matrix and are color indices. Shorthand notations and are also employed in (3).

In TQCDSR, the correlation function given in (1) is calculated twice, as in two different regions corresponding two perspectives, namely the physical side (or phenomenological side) and the QCD side (or OPE side). By equating these two approaches, the sum rules for the hadronic properties of the exotic state under investigation are achieved. To derive mass and pole residue via TQCDSR, the correlation function is calculated in terms of hadronic degrees of freedoms in the physical side. A complete set of intermediate physical states possessing the same quantum number as the interpolating current are inserted into (1), and integral over is handled. After these manipulations, the correlation function is obtained ashere is the temperature-dependent mass of meson. Temperature-dependent pole residue is defined in terms of matrix element aswhere is the polarization vector of the satisfyingAfter employing polarization relations, the correlation function is written in terms of Lorentz structures in the formwhere dots denote the contributions coming from the continuum and higher states. To obtain the sum rules, coefficient of any Lorentz structure can be used. In this work, coefficients of are chosen to construct the sum rules and the standard Borel transformation with respect to is applied to suppress the unwanted contributions. The final form of the physical side is obtained ashere is the Borel mass parameter. In the QCD side, is calculated in terms of quark-gluon degrees of freedom and can be separated into two parts over the Lorentz structures as where and are invariant functions connected with the scalar and vector currents, respectively. In the rest framework of (), can be expressed as a dispersion integral, where corresponding spectral density is described as The spectral density can be separated in terms of operator dimensions asIn order to obtain the expressions of these spectral density terms, the current expression given in (3) is inserted into the correlation function given in (1) and then the heavy and light quark fields are contracted, and the correlation function is written in terms of quark propagators as where are the full quark propagators and is used. The quark propagators are given in terms of the quark and gluon condensates [25]. At finite temperatures, additional operators arise due to the breaking of Lorentz invariance by the choice of thermal rest frame. Thus, the residual invariance brings additional operators to the quark propagator at finite temperature. The expected behavior of the thermal averages of these new operators is opposite of those of the Lorentz invariant old ones [43]. The heavy-quark propagator in coordinate space can be expressed asand the thermal light quark propagator is chosen aswhere is the external gluon field, with Gell-Mann matrices, symbolizes color indices, implies the strange quark mass, is the four-velocity of the heat bath, is the temperature-dependent light quark condensate, and is the fermionic part of the energy-momentum tensor. Furthermore, the gluon condensate related to the gluonic part of the energy-momentum tensor is expressed via relation [43]:where and coefficients areIn order to remove contributions originating from higher states, the standard Borel transformation with respect to is applied in the QCD side as well. By equating the coefficients of the selected structure in both physical and QCD sides, and by employing the quark hadron duality ansatz up to a temperature-dependent continuum threshold , the final sum rules for are derived asTo find the mass via TQCDSR, one should expel the hadronic coupling constant from the sum rules. It is commonly done by dividing the derivative of the sum rule given in (18) with respect to to itself. Following these steps, the temperature-dependent mass is obtained aswhere the thermal continuum threshold is related to continuum threshold at vacuum via relation[44, 45]. For compactness, the explicit forms of spectral densities are presented in Appendix.

3. Phenomenological Analysis

In this section, the phenomenological analysis of sum rules obtained in (18) and (19) is presented. First, the input parameters and temperature dependance of relevant condensates are given. Following, the working regions of the obtained sum rules at vacuum are analyzed. The behavior of QCD sum rules at is used to test the reliability of our analysis.

3.1. Input Parameters

During the calculations, input parameters given in Table 1 are used. In addition to these input parameters, temperature-dependent quark and gluon condensates, and the energy density expressions are necessary. The thermal quark condensate is chosen aswhere is the light quark condensate at vacuum and which is credible up to a critical temperature . The expression given in (21) is obtained in [46, 47] from the Lattice QCD results given in [48, 49]. The temperature-dependent gluon condensate is parameterized via [46, 50]where is the gluon condensate in vacuum state and . Additionally, for the gluonic and fermionic parts of the energy density, the following parametrization is used [46]which is extracted from the Lattice QCD data in [51].

3.2. Analysis of Sum Rules at

In order to get reliable results, obtained sum rules should be tested at vacuum, and the working regions of the parameters and should be determined. Within the working regions of and , convergence of OPE and dominance of pole contributions should be assured. In addition, the obtained physical results should be independent of small variations of these parameters. Convergence of the OPE is tested by the following criterion. The contribution of the highest order operator in the OPE should be very small compared to the total contribution. In Figure 1, the ratio of the sum of the terms up to the specified dimension to the total contribution is plotted to test the OPE convergence. It is seen that all higher order terms contribute less than the perturbative part for GeV. On the other hand, dominance of the pole contribution is tested as follows. The contribution coming from the pole of the ground state should be greater than the contribution of the continuum. In this work, the aforementioned ratio iswhen as can be seen in Figure 1. After checking these criteria, the working regions of the parameters and are determined as which is also consistent with GeV norm [10]. Within these working regions, the variations of the mass of with respect to and are plotted in Figure 2. It is seen that the mass is stable with respect to variations of and .

4. Results and Discussions

Following the analysis presented in previous section, the mass of the ground state estimated by the tetraquark current given in (3) at is presented together with other results from literature in Table 2. It is seen that our results agree with other theoretical estimates and also with the experimental data on [10, 11, 20]. Thus, the broad resonance can be described by the tetraquark current given in (3), and our analysis can be extended to finite temperatures.

To analyze the thermal properties of resonance, the temperature dependencies of the mass and the pole residue of are plotted in Figure 3. It is seen that the mass and the pole residue of stay monotonous until . However, after this point, they begin to decrease promptly with increasing temperature. At the vicinity of the critical (or so called deconfinement) temperature, the mass reaches nearly of its vacuum value. On the other hand, the pole residue decreased to of its value at vacuum as shown in Figure 3.

Our predictions presented in Figure 3 are in good agreement with other QCD sum rules analysis on thermal behaviors of conventional or exotic hadrons [3538]. However in [19], authors predicted a decrease of in the mass of molecular state, and no change in the mass of charmonium state, even beyond Hagedorn temperature MeV. Since the decay properties also depend on temperature, and while the mass and pole residue diminish, the decay width might increase with increasing temperature [52] and decay widths at finite temperature should also be investigated. However the current status of resonance is very complicated, since it is very close to state. Thus studying its decays requires establishment of a good model in the hidden-beauty sector.

Finally, we would like to highlight the following remarks:(i)We observed that the mass (the pole residue) of exotic state starts to decrease near .(ii)Both quantities tend to diminish with increasing temperature up to critical temperature .(iii)Even though the sum rules at estimates the mass of consistent with experimental data, more theoretical efforts are required to discriminate and .(iv)In order to get more reliable results on tetraquarks from QCD sum rules, O() contributions suggested in [36] should be investigated.

In summary, we revisited the hidden-beauty exotic state and studied its properties at vacuum and finite temperatures. To describe the hot medium effects to the hadronic parameters of the resonance , TQCDSR method is used considering contributions of condensates up to dimension six. Our results for are in reasonable agreement with the available experimental data and other QCD studies in the literature. Numerical findings show that can be well described by a scalar-vector tetraquark current. In the literature, remarkable drop in the values of the mass and the pole residue in hot medium was regarded as the signal of the QGP, which is called the fifth state of matter, phase transition. We hope that precise spectroscopic measurements in the exotic bottomonium sector can be done at Super-B factories, and this might provide conclusive answers on the nature and thermal behaviors of the exotic states.

Appendix

Thermal Spectral Density for State

In this appendix, the explicit forms of the spectral densities obtained in this work are presented. The expressions for and are shown below as integrals over the Feynman parameters and , where is the step function, where explicit expressions of the functions and are The below definitions are used for simplicity:

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

J. Y. Süngü, A. Türkan, and E. Veli Veliev thank Kocaeli University for the partial financial support through the grant BAP 2018/082. H. Da acknowledges support through the Scientific and Technological Research Council of Turkey (TUBITAK) BIDEP-2219 grant.