S-wave single heavy baryons with spin-3/2 at finite temperature

The thermal behavior of the spectroscopic parameters of the S-wave single heavy baryons ΣQ∗,ΞQ∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Sigma _{Q}^{*}, \Xi _{Q}^{*}$$\end{document} and ΩQ∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Omega _{Q}^{*} $$\end{document} with spin-3/2 are investigated in QCD at finite temperature. We analyze the variations of the mass and residue of these baryons taking into consideration the contributions of QCD thermal condensates up to dimension eight in Wilson expansion. At finite temperature, due to the breakdown of the Lorentz invariance by the choice of reference frame and presence of an extra O(3) symmetry, some new four-dimensional operators come out in the form of the fermionic and gluonic parts of the energy momentum tensor that are taken into account in the calculations. Our analyses show that at lower temperatures, the parameters of baryons under consideration are not affected by the medium. These parameters, however, show rapid variations with respect to temperature at higher temperatures near to a pseudo-critical temperature, after which the baryons are melted. The results of the masses and residues at T→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T\rightarrow 0 $$\end{document} limit are compatible with the available experimental data and predictions of other theoretical studies.


Introduction
With the rising number of experimental data on charmed and bottom baryons, the interest in the investigation of heavy baryons has increased, considerably. Before giving the details of the experimental studies on heavy baryons, it would be useful to give some theoretical information. The Quark Model is one of the most successful tools to classify the mesons and baryons. The traditional single heavy baryons (Qqq) consist of one heavy (Q = b or c) and two light quarks (q = u, d or s). The mass of heavy quark is very large compared to the light quark masses and the light degrees of freedom form a diquark qq, which orbits the nearly static heavy Q quark. Therefore, infinitely heavy mass limit (m Q → ∞) a e-mail: kazem.azizi@ut.ac.ir (corresponding author) for the heavy quark is utilized to classify the single heavy baryons [1][2][3]. In this case, for the two light quarks, the total flavor-spin wave function has to be symmetric because their color wave function is antisymmetric. Hence there are two different representations for the S-wave heavy baryons (3 ⊗ 3 = 3 ⊕ 6): antisymmetric 3 or symmetric 6. The antitriplet (3) of baryons contain only spin-1/2 states while the sextet (6) of baryons contain both spin-1/2 and spin-3/2 states. In this study, we investigate the thermal properties of the single heavy bottom/charmed spin-3/2 sextet states: The members for charmed baryons are shown in Fig. 1.
Experimentally  [19]. Roncaglia et al. [20,21] estimated the heavy baryon masses with one/two heavy quark/quarks in the framework of Feynman-Hellman theorem. Jenkins [22] studied the masses of heavy baryons in the 1/m Q and 1/N c expansions. The 1/m corrections to heavy baryon masses were calculated by Dai et al. [23] in the framework of the HQET. QCD sum rule for heavy baryons at leading order in 1/m Q and at next to the leading order in α s were evaluated by Groote et al. [24]. Wang et al. [25] improved the analysis for the Q and Q baryon masses to order QC D /m Q from QCD sum rule. Mathur et al. [26] predicted the mass spectrum of charmed and bottom baryons from Lattice QCD. Wang and Huang [27] studied the mass, coupling constant, and Isgur-Wise function for ground-state heavy baryons within the framework of HQET by taking into account both the two and three-point correlation functions. Ebert et al. computed heavy baryon masses in the heavy-quark light-diquark approximation in the framework of constituent quark model [28]. Garcilazo et al. [29] solved exactly the three quark problem via Faddeev method in momentum space. Zang and Huang [30] calculated the charm and bottom baryon masses up to operator dimension six in operator product expansion (OPE) by the help of the QCD sum rule approach. The mass and residue of * c and * b with spin parity 3/2 + were studied by Wang via QCD sum rule [31]. A quark model was applied to the spectrum of baryons containing one heavy baryon by Roberts and Pervin [32]. Bottom baryon spectra were investigated using Faddeev method in momentum space by Valcarce et al. [33]. Liu et al. [34] performed a systematic study of the masses of bottom baryons up to 1/m Q in HQET. Groote et al. [35] computed the NLO perturbative corrections for the static properties of heavy baryons. The mass of Q and ( * ) Q baryons were calculated by Zhang and Huang [36] via QCD sum rule taking into account operators up to dimension six. Using the coupled channel formalism, Gerasyuta and Matskevich calculated the S-wave bottom baryons masses [37]. Karliner et al. [38] investigated the b-baryons in the quark model. In two-point and light cone QCD sum rule methods Aliev et al. studied the mass and magnetic moments of single heavy baryons with spin-3/2 [39]. Lewis and Shyn [40] predicted the bottom baryon masses based on a 2 + 1 flavor dynamical lattice QCD simulation. The spin-3/2 + heavy and doubly heavy baryon states [41] were investigated by subtracting the contributions from the corresponding negative parity by Wang. The mass spectra of heavy baryons were studied by the help of the motivated relativistic quark model by Ebert [42]. Kim et al. investigated the single heavy baryon mass based on the self-consistent Chiral quark soliton model [43]. Finally, Azizi and Er studied the in-medium properties of spin-3/2 heavy baryons in nuclear matter using QCD sum rule in a dense medium [44].
Theoretical investigations of spectroscopic parameters of the single heavy baryons at finite temperature will help us better understand and analyze the results of heavy-ion collision experiments and gain valuable information on the internal structures of these baryons, behavior of these baryons near to a pseudo-critical temperature, possible phase transition/ crossover [45,46] to/with quark gluon plasma (QGP) (adopted as a new phase of matter ) as well as the perturbative and nonperturbative dynamics of QCD. At extreme temperatures, two different possibilities can be considered: crossover and phase transition. Many Lattice calculations predict that crossover occurs at T pc ≈ 155 MeV [47,48]. For the QGP phase transition, we need greater temperature values and there is no unique temperature to the phase transition of QGP. At short distances, to describe the strong interaction QCD is a suitable theory. However, the calculations of hadronic parameters including nonperturbative effects (occur in low energy scale) usually need some nonperturbative phenomenological models. Many phenomenological models are available in the literature: QCD sum rule is one of the powerful ones among them. This method firstly suggested by Shifman, Vainshtein and Zakharov to investigate the vacuum properties of mesons [49,50] and then Ioffe [51] applied this method for baryons. The thermal version of the QCD sum rule was extended by Bochkarev and Shaposhnikov [52]. In addition to the vacuum expectation values of quark and gluon condensates, their thermal forms and some new operators appear in the thermal version.
In this study, we investigate the temperature effects on the spectroscopic parameters of the ground state sextet baryons including single heavy quark and with spin-3/2 via thermal QCD sum rule method. Taking into account the additional operators coming from OPE due to breaking of the Lorentz invariance by the choice of the thermal rest frame, condensates up to dimension eight are considered. The article is arranged in the following form. In Sect. 2, the in-medium sum rules for the mass and residues of the * Q , * Q and * Q single heavy baryons are obtained. In Sect. 3 the numerical analysis for the spectroscopic parameters under consideration is performed. The last section includes the summary and our concluding remarks.

Calculations
In this section, QCD sum rules for the spectroscopic parameters of the spin-3/2 * Q , * Q and * Q baryons are obtained at finite temperature. To this end, we start with the following two-point thermal correlation function: where q is the four-momentum of the chosen baryon, is the ground state of the hot medium, T denotes the time-ordering operator and J μ (x) is the interpolating current of the single heavy baryon, B S H . As the standard procedures of the QCD sum rule, the correlation function given above can be calculated at different contexts. At large distances, it is evaluated in terms of the hadronic parameters such as the mass and residue of hadron. We call it the physical or hadronic representation of the correlator. The same correlator can be expressed in terms of the quark, gluon and mixed condensates by the help of the OPE at q 2 << 0 region. The computations in this way contain short distance effects. This representation, is generally called the OPE or QCD side of the correlation function. Finally, we match the two windows and compare the coefficients of the same Lorentz structures from both sides. To remove the unwanted contributions coming from the higher states and continuum, Borel transformation as well as continuum subtraction, supplied by the quark-hadron duality assumption at finite temperature, are performed. These procedures bring some auxiliary parameters, which we fix them before making any numerical estimations on the physical quantities.
To obtain the physical side of the correlator, a complete set of intermediate state with the same quantum numbers and quark content as the chosen current is inserted between the interpolating currents in correlation function. This is followed by the integral over four-x, which leads to where m B S H (T ) is the temperature-dependent mass of the ground state of B S H . The matrix element where u μ (q, s) is the Rarita-Schwinger spinor. The final form of the physical side can be obtained by inserting Eq.
(3) into Eq. (2) and summing over the spins of the B S H . The summation over Rarita-Schwinger spinors is performed using By using the above behest, we recast the physical side as Phys where . It should also be specified that the interpolating current J μ (x) couples to both the spin-1/2 and spin-3/2 states. In this study, we only consider the contribution of spin-3/2 heavy baryons and we need to comb out the pollution of spin-1/2 state. These unwanted contributions can be eliminated in two different ways: (1) For spin-3/2 state, it should be introduced a projection operator which destroys the spin-1/2 contributions, (2) By a specific ordering of the Dirac matrices and remove the terms corresponding to the spin-1/2 particles (for more details see for instance [53]). The contribution of the spin-1/2 states can be traced using where κ 1 and κ 2 are some constants. By applying the condition J μ γ μ = 0 (for more details see [54]), we get κ 1 in terms of κ 2 . Hence, As is seen from Eq. (7), the pollution coming from spin-1/2 resonances are commensurate to either q μ or γ μ . To remove these contributions, the Dirac matrices are ordered as γ μ qγ ν and terms proportional to q μ or q ν , also those beginning with γ μ or ending with γ ν are set to zero. Finally, the clean physical side of the correlator, in the Borel scheme, is obtained as Table 1 The light quark flavors for the single heavy baryons with spin-3/2 and the value of normalization constant where M 2 is the Borel parameter and dots denote the contributions of other structures as well as the higher states and continuum. The next step is to calculate the OPE side of the correlation function. In deep Euclidean region, the correlation function is evaluated in terms of the quark and gluon degrees of freedom by the help of Wilson expansion. To achieve this goal, the basic point it to choose a suitable interpolating current defining the particles under study. The interpolating current for spin-3/2 B S H in a compact form can be written as [55][56][57] where A is the normalization constant, abc is the antisymmetric Levi-Civita tensor, a, b, c are color indices, q 1(2) denotes the light quark (u, d or s), Q is the bottom (b) or charm (c) quark and C is the charge conjugation operator. The normalization constant A and the q 1(2) quark for the considered baryons are given in Table 1. By inserting the explicit form of the interpolating current into the correlator and contracting all heavy and light quark fields via Wick's theorem, we get the corelation function in the case of q 1 = q 2 in terms of the thermal light(heavy) quark propagators, S q(Q) , as Some extra contractions arise because of the identical particles in the case of q 1 = q 2 = q, and the correlator is obtained as where S i j To go further in the calculations, the thermal light quark propagator in coordinate space is selected as (see also [58,59]) which includes the thermal quark and gluon condensates ( qq T and g 2 s G 2 T ), gluon fields in thermal bath, mixed condensate (m 2 0 qq T = qg s σ Gq ) as well as new operators containing the energy momentum tensor, μν . For the heavy quark, the following propagator including the thermal gluon condensate and gluon fields in hot medium is used [60]: In Eqs. (12) and (13), m q(Q) denotes the light(heavy) quark mass.
The thermal quark condensates, qq T (for q = u, d) and ss T are parameterized in terms of the vacuum condensates, 0|qq|0 and 0|ss|0 . For these quantities, we use the following parametrizations in terms of temperature, which are based on the lattice QCD predictions [61]. Note that in this study the temperature dependence of these quantities are given up to a temperature T = 300 MeV. However, we parameterize them up to T pc ≈ 155 MeV, which is considered as the pseudocritical temperature for the crossover phase transition at zero chemical potential. We get, and where A 1 = −6.534 × 10 −4 and A 2 = −2.169 × 10 −5 .
As we previously mentioned, because of the choice of the thermal rest frame in Wilson expansion, the Lorentz invariance is broken. To restore that the four-velocity vector of the medium u μ = (1, 0, 0, 0) is introduced, which implies u 2 = 1 and q · u = q 0 . In the rest frame of heat bath, u μ f,g μν u ν = u f,g u = f,g 00 = f,g , as well. In thermal version, as also mentioned above, new operators representing the fermionic and gluonic parts of the energymomentum tensor arises in OPE. The fermionic part f μν appears explicitly in the light-quark propagator, while the gluonic part of the energy-momentum tensor g λσ appears in the expansion of the trace of two-gluon field strength tensor in heat bath [62]: The temperature dependent gluon condensate G 2 T is parameterized in terms of the vacuum gluon condensate 0|G 2 |0 [61] as: For the temperature-dependent strong coupling [65,66] we utilize where, M S T pc /1.14.
Alike to the physical part, the correlation function on the OPE side is expanded in terms of the Lorentz structures as where O P E 1(2) is the coefficient of the selected Lorentz structure. These functions can be expressed by the help of following dispersion integral: where s min = (m q 1 +m q 2 +m Q ) 2 , ρ O P E 1(2) (s, T ) is the spectral density obtained via the imaginary part of the perturbative correlation function ( pert in the following equation stands for the perturbative contributions) and nonpert 1 (2) represents the contributions coming from all the nonperturbative effects. In this step, our main aim is to calculate the spectral densities, corresponding to the perturbative effects in the present study, as well as the nonperturbative contributions to the QCD side. To this end, the explicit forms of the heavy and light quark propagators are inserted into Eqs. (10) and (11). The next step is to perform the standard but lengthy calculations: These calculations contain Fourier integrals appearing in different forms, Borel transformation as well as continuum subtraction. By matching the coefficients of the selected structures from both the physical and OPE sides of the correlation function, we find the desired sum rules: and where the functionsB O P E 1(2) denote the O P E 1 (2) in Borel scheme and are given aŝ with s 0 (T ) being the temperature-dependent continuum threshold. We will use the above sum rules to extract the values of the mass and residue of the baryons under consideation as well as their thermal behavior in next section.
As examples, we would like to present the explicit forms of the ρ O P E 1 (s, T ) andB nonpert 1 for the * b baryon. They are obtained as and where stands for the unit-step function, L(s, z) = s z(1 − z) − m 2 b z and β = z − 1.

Numerical results
In this section, we analyze the obtained sum rules for the masses and residues. They includes some input parameters such as the heavy and light quark masses, m 2 0 , quark and gluon condensates in vacuum and energy of the quasi-particle in medium, q 0 . Their numerical values are presented in Table 2.
In addition, we also need to have the gluonic and fermionic parts of the energy density. Based on the lattice QCD results on the thermal behavior of the energy-momentum tensor given in [63], their parametrizations, up to the pseudo-critical point under consideration in the present study, are obtained as which, we are going to use them in our numerical computations. The next problem is to obtain the parametrization of s 0 (T ) as a function of temperature. This function shall reduce to the vacuum threshold, s 0 , at zero temperature. We parameterize it as such that at T → 0 limit, f (T ) → 1. Hence, we should first determine s 0 based on the standard prescriptions of the method, afterwards we will extract the function f (T ) from the calculations. Besides the continuum threshold in vacuum the sum rules obtained in previous section include another auxiliary parameter, Borel parameter M 2 , which should also be fixed. We need to determine the working regions of s 0 and M 2 such that the physical quantities under consideration show mild dependence on these parameters. The continuum threshold s 0 is not totally free but it is related to the energy of the first excited state in the same channel. Thanks to the experiments that have provided many new results not only on the ground states but also on the excited states of some single heavy baryons, recently [70]. In view of PDG, we see that the excited states generally have energies about 300 MeV above the ground states masses. In choosing the working window for the s 0 , we also look after the pole dominance and OPE convergence in our sum rules. These considerations leads to the window: The upper and lower limits of the Borel parameter are fixed consider the criteria of the QCD sum rule method. To find the lower limit, we apply the criterion of the OPE convergence at the chosen window for the continuum threshold. To this end, we demand that the perturbative part exceeds the total nonperturbative contributions and the slogan of the higher the dimension of the nonperturbative operator the lower its contribution is satisfied. Our calculations show that the operators having eight dimensions, the higher dimension that we include into the analyses, constitute only one percent of the total contribution at lower value of M 2 , i.e. Figure 2 shows the perturbative and nonperturbative contributions to total OPE as well as the contributions of different nonperturbative operators with various mass dimensions, separately. This figure depicts a nice convergence of the OPE in our calculations. As it is clear, the perturbative contribution dominates over nonperturbative contributions and it is about 53% of the total at M 2 min = 6 GeV 2 . The main contribution in nonperturbative part belongs to the quark condensate, qq .
To obtain M 2 max , we utilize the condition of the pole dominance as As a result, we get the working region of the Borel parameter as M 2 ∈ [6, 10] GeV 2 . We plot, as an example, a 3D graphic of the mass of * b baryon as functions of M 2 and s 0 at T = 0 in Fig. 3. As is seen the mass shows good stability against the variations of the auxiliary parameters in the selected windows. Now, we proceed to find the function f (T ) and the temperature dependent mass m B S H (T ) and residue λ B S H (T ) of the single heavy spin-3/2 baryons. To this end, we use the two sum rules in Eqs. (23) and (24) Fig. 4 at average value of the vacuum continuum threshold. This figure shows that the spectroscopic parameters of the * b , * b and * b baryons are stable against the changes in temperature until a certain tem-perature but after that, they start to decrease with increasing the temperature. Our analyses show that the charmed baryons present similar behavior, as well. The points that the stability starts to break down for mass and residue are T ∼ = 0.14 GeV and T ∼ = 0.13 GeV, respectively. After these points the mass and residue starts to diminish. The mass and residue fall substantially near to the pseudo-critical temperature. The amount of decrements at T pc are 75% (66-71%) for the mass of bottom (charmed) and 71-80% (42-50%) for the residue of bottom (charmed) baryons, respectively compared to their vacuum values. These behavior of baryons can be interpreted as substantial melting of the heavy baryons near to the pseudo-critical temperature.
At the end of this section, we would like to present our results for the masses of the single heavy spin-3/2 baryons at T → 0 limit. This is done in Table 3. For comparison, we also present the existing theoretical predictions in the literature and experimental data in the same table. With a quick glance in this table, we see that our predictions, within the errors, are overall consistent with other theoretical predictions made using different methods and approaches. Our predictions are also well consistent with the existing experimental data for five members within the presented uncertainties. * b baryon is only missing member, which has not been discovered in the experiment. We hope that, our result together with other  theoretical predictions will help experimental group in the course of search for this particle.

Summary and concluding remarks
In this study, we have performed two-point thermal QCD sum rule analyses for * Q , * Q and * Q single heavy baryons which are the members of the spin-3/2 sextet family. In the OPE, operators up to dimension eight were taken into account which lead to a good OPE convergence as well as pole dominance. We included the thermal effects by two ways: We replaced the vacuum condensates by their thermal versions and considered the extra operators, appearing in the forms of the fermionic and gluonic parts of the energy momentum tensor due to the restoration of the Lorentz invariance. We fixed the auxiliary parameters entering the calculations by the standard prescriptions of the method. By simultaneous solving of the two sum rules obtained together with an extra equation derived from one of the sum rules, we found three unknowns: Thermal continuum threshold, temperature-dependent mass and temperature-dependent residue. We discussed the thermal behavior of the mass and residue for all the bottom and charmed baryon members having the spin-3/2. We observed that the spectroscopic parameters remain unchanged up to a certain temperature, after which they start to diminish considerably near to the pseudo-critical temperature. The decrements order in the mass and residue of the considered baryons near to the pseudo-critical temperature are obtained as (66-75)% and (42-80)%, respectively, representing substantial melting of the heavy baryons near to the pseudo-critical temperature. In the literature, there are no other studies on the thermal behavior of single heavy baryons to make a compression with our predictions. However, there are some studies on the temperature dependence of the masses of light baryons, In Refs. [58,71] the authors investigated the light octet and decuplet baryons using the thermal QCD sum rule, but considering a pseudo-critical temperature of T pc = 197 MeV. They obtained that the shifts in the masses of the considered baryons are overall about 80%. The pole mass of the octet and decuplet baryons were also evaluated in Ref. [72] via the chiral perturbation theory. The authors observed that a 20% mass shift occurs around the temperature T 150 MeV, where the freeze-out in the relativistic heavy-ion collision is expected to be formed. Using the many-body techniques at finite temperature, all baryonic states of the octet and decuplet flavors were examined in Ref. [73]. They obtained that the baryon masses decrease with the temperature and there are strong dependencies on the melting (or deconfinement) temperature depending on the flavor content of the baryons. In the framework of the thermal QCD sum rule, the masses of the decuplet baryons were also investigated in Ref. [74]. According to this study, the masses of the decuplet baryons show very little temperature dependence below T = 0.11 GeV and the melting or hadron-quark phase transition occurs at a temperature T ≥ 0.11 GeV. Our results indicate that this point is T = 0.14 GeV for heavy baryons, after which the masses start to decrease with the increasing of the temperature and the dependence of the masses on temperature near to the critical temperature is very strong. These information on the behavior of the masses of different baryons may help experimental groups in the analyses of the results of the in-medium and heavy ion collision experiments, despite the statistical hadronization model claims that any thermal modification of masses is negligibly small at pseudo-critical temperature and the in-medium mass shifts at T pc would be excluded.
We extracted the values of the masses for both the bottom and charmed baryons at T → 0 limit and compared with the predictions of other phenomenological models and experimental data. The obtained results are well consistent with existing experimental data. Our result on the mass of * b baryon as the only undiscovered member together with other predictions may help the experimental group to hunt this particle and measure its parameters.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All the numerical and mathematical data have been included in the paper and we have no other data regarding this paper.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .