Excited states of spin-(3/2) doubly-heavy baryons within the QCD sum rules method

Within the framework of the QCD sum rules method, the masses and pole residues of the excited states 2S and 1P of the doubly-heavy spin-(3/2) baryons Ξcc∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _{cc}^*$$\end{document}, Ξbb∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _{bb}^*$$\end{document}, Ξbc∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _{bc}^*$$\end{document}, Ωcc∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _{cc}^*$$\end{document}, Ωbb∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _{bb}^*$$\end{document}, and Ωbc∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _{bc}^*$$\end{document} are investigated. The excited state masses and pole residues are predicted with reasonable accuracy. We expect that soon, the current predictions can be tested against future experimental data.

The SELEX Collaboration, on the other hand, has observed the doubly-heavy charmed baryon state + cc with spin-(1/2) and mass (3519 ± 1) MeV in the pD + K − decay channel [16], and then later was confirmed by the same experiment in 2005 [17]. Whereas, the LHCb Collaboration have recently reported that the ground state mass of ++ cc is (3624.40 ± 0.72 ± 0.14) MeV [2]. These findings stimulated the search for a new class of long-lived baryons called doubly-heavy baryons. a e-mail: thabit.barakat@tedu.edu.tr (corresponding author) The doubly-heavy baryons are hadrons containing two heavy quarks of (c and/or b), and one light quark of (d or u or s). The existence of such doubly-heavy baryons is considered a natural consequence of the quark model of hadrons.
Meanwhile, the LHCb has already begun the searching for the doubly-heavy bottom charmed baryon cb and is still working to improve its signal sensitivity [18].
The existence of such experiments has prompted theorists and phenomenologists to investigate the structural features of the doubly-heavy flavor baryons. The charmed and bottom baryons, which are composed of one (two) heavy quark(s), are of particular interest for studying the dynamics of light quarks in the presence of heavy quarks and may provide an exceptional platform for testing quark model estimations and heavy quark symmetries.
Motivated by the possibility of observing these new doubly-heavy baryons shortly, the theoretical study of such baryons could shed light on the dynamics of the heavy quarks at the hadronic scale, and hence might be represented as a truthful investigation of the interaction between perturbative and non-perturbative QCD pictures. In the theory, the study of doubly-heavy baryons could provide crucial inputs for future experimental discoveries, as well as possibly clarifying and revealing secrets about the nature of other heavy baryons.
One of the most important physical properties of the doubly-heavy baryons is their masses and residues, which can be accessed experimentally. Based on this, many authors have created and established various phenomenology methods and models to estimate the mass spectrum of doublyheavy baryons, such as the non-relativistic and relativistic quark models [19,20], the three-body Faddeev method [21], the QCD MIT bag model [22], effective field theory with the non-relativistic QCD potential method [23], the heavy quark spin symmetry method [24], the variational approach method [25], the full QCD sum rules method [26], and the lattice QCD method [27,28]. In addition to that, many other theoretical works inspired by the above mentioned discoveries also have been carried out to determine the characteristic properties like the masses, lifetimes, and the strong coupling constants of such baryons within various models [29][30][31][32][33][34][35][36].
Thus far, the QCD sum rules method is widely recognized as one of the most essential nonperturbative analytical methods for estimating the ground state masses of the heavy hadrons, and the numerical estimation of this method is extremely accurate. The success of this method made it as one of the most popular theoretical methods for testing the various properties of doubly-heavy baryons. The masses of spin-(1/2) doubly-heavy baryons, for example, were computed using the QCD sum rules approach in [37][38][39][40][41], whereas the ground state masses of spin-(3/2) doubly-heavy baryons were computed in [38][39][40][41][42][43][44][45]. The masses and pole residues of the ground state of spin-(3/2) doubly-heavy baryons * Q Qq , * Q Q q , * Q Qq , and * Q Q q , on the other hand, have been the focus of previous research [44]. Therefore, there is a lot of work to be done to have a full spectrum of spin-(3/2) doublyheavy baryons.
Due to its numerical performance, we were tempted to extend the calculations of the QCD sum rules in [44]. However, this time to understand how the QCD sum rules method can be applied to derive the excited state masses and the excited state pole residues of spin-(3/2) doubly-heavy baryons. To do so, in the present work, we'll use the most common type of interpolation current for spin-(3/2) doublyheavy baryons.
With this in mind, the form of this paper is arranged as follows: Sect. 2 presents the formulation of the QCD sum rules for spin-(3/2) doubly-heavy baryons. The analytical formulations of the QCD sum rules are written in such a way that the reader can use them without proceeding into their derivation. Section 3, focuses on the numerical analysis of the QCD sum rules for excited state masses and excited state pole residues of spin-(3/2) doubly-heavy baryons; 2S and 1P of * cc , * bb , * bc , * cc , * bb , and * bc . Finally, in the last section of this paper, we presented a summary of our findings and some concluding remarks.

QCD sum rules for doubly-heavy spin-(3/2) baryon states
One of the most important methods for non-perturbative QCD is the QCD sum rules method. In this method, hadrons are represented by their interpolating quark currents, and the major object of this method is the so-called correlation function ( p 2 ), which can be expressed in terms of the interpolating quark currents. The correlation function ( p 2 ) is an analytic function of p 2 that can be defined at both negative (space-like) and positive (time-like) values of p 2 . When p 2 is shifted from large negative values towards positive values, the correlation function starts receiving contributions from the long-distance quark-gluon interactions. In this region, quarks start to form hadrons, and the phenomenological (or hadronic) side of the correlation function is investigated, and can be expressed in terms of hadronic degrees of freedom by inserting a complete set of intermediate hadronic states with the same quantum numbers as the current operators. In this representation, the correlation function leads to a masssum rule, which is in turn related to the physically observed hadrons.
On the other hand, when p 2 0 (space-like), the correlation function is studied in terms of the quarks and gluons and their interactions with the QCD vacuum. This is referred to as theoretical (or QCD) representation. In this representation, the Wilson operator product expansion (OPE) method is usually used to separate the short-and long-distance quark-gluon interactions, by expressing the correlation function in terms of the QCD degrees of freedom, such as quark condensate, gluon condensate, etc.
Finally, equating the two sides of these representations and applying the double Borel transformations to the momentum of the state to suppress the contribution of the higher states and continuum, the QCD sum rules for the physical observables such as masses and residues can be estimated by matching the coefficients of the same structure in both representations of the correlation function.
Consequently, we start the analysis of spin-(3/2) doublyheavy baryon states by taking into account the two-point correlation function μν ( p 2 ) between vacuum states, with no initial and final baryons: where T represents the time ordering operator, η μ is the interpolating current of the doubly-heavy baryon with spin-(3/2) that injects quarks into the QCD vacuum at point x, and p is the four-momentum of the doubly-heavy baryon. The conventional expression of the interpolating current, which may be coupled to both the positive parity as well as to negative parity spin-(3/2) doubly-heavy baryons is chosen to be of the following form: where q denotes the light quark, Q and Q are the two heavy quarks. T stands for transposition, C for charge conjugation, a, b, and c are for color indices.
The systematic procedure of the QCD sum rules method is started by saturating the above correlation function with full sets of the baryon states with equal quantum numbers as the interpolating current, and by summing on for the ground state and associated excited states, we obtain the physical part or phenomenological representation of the correlation function as: where the summation here only includes the required hadronic states produced by the quark current η μ , and... includes contributions from higher states as well as the continuum.
The matrix elements of the interpolating current evolved into Eq. (3), on the other hand, are parameterized in terms of the three pole residues for the three different baryonic states, 1S, 2S, and 1P, where the λ's are the pole residues of spin-(3/2) doubly-heavy baryon states, and u( p, r ) is the Rarita-Schwinger spinor satisfying the identity, where, m B 1S ≡ m g , m B 2S ≡ m, and m B 1P ≡m are to be the masses of spin-(3/2) doubly-heavy baryon states, as well as, λ B 1S ≡ λ g , λ B 2S ≡ λ, and λ B 1P ≡λ are to be the pole residues of spin-(3/2) doubly-heavy baryon states. Before going further, a word of caution is in order. The above-mentioned interpolating current η μ of Eq. (2) couples to both the spin-(1/2) and spin-(3/2) baryons with positive and negative parities. However, because we are only interested in studying spin-(3/2) doubly-heavy baryon states of both positive and negative parities in this work. Therefore, the unwanted spin-(1/2) contributions of both parities should be removed from this study. It is clear from Eq. (5) that the unwanted spin-(1/2) structures are proportional to either γ μ γ ν , p μ p ν , and ( p μ γ ν − p ν γ μ ).
Overall, Eq. (5) shows that there are only two independent Lorentz structures left, / pg μν and g μν , which are coupled only to spin-(3/2) doubly-heavy baryons with both positive and negative parities. On the whole, the hadronic side (HAD) therefore contains the (1S + 2S + 1P) spin-(3/2) doublyheavy baryonic states, and the excited states are considered as (2S) and (1P) states, In this regard, the operator product expansion approach (OPE) in the deep Euclidean region − p 2 m 2 Q on the other hand is used to determine the QCD side of the correlation function. The relevant correlation function in terms of the heavy (light) quark propagator may be obtained by putting the explicit expressions of the interpolating currents into the correlation function of Eq. (1) and then, by doing some contractions of all quark fields through the Wicks theorem, we obtain: whereS i j is the heavy (light) quark propagator in the coordinate space, and the explicit mathematical expressions of the heavy (light) quark propagators are, where K 1 and K 2 are the modified Bessel functions of the second kind, qq is the quark condensates, and m 2 0 is the conventional parametrization constant in the quark-gluon condensates.
At this point, we make the necessary substitutions of the heavy (light) quark propagators into Eq. (7), the twopoint correlation function can be then re-expressed in terms of some spectral densities ρ(s)'s with the two independent structures pg μν and g μν , respectively, where, The spectral densities ρ 1 (s) and ρ 2 (s) in Eq. (10) are obtained from the imaginary part of the correlator function; ρ i (s) =  2) is related to the structures / pg μν and g μν . Further, s 0 is the continuum threshold, which is not fully arbitrary and it is associated with the energy of the exited states, s min = (m Q + m Q ) 2 , and finally, the spectral densities ρ 1 (s) and ρ 2 (s) are written as (see also [44]): where, After calculating both the HAD and QCD sides of the correlation function, we matched the coefficients of the same structures from both sides of Eqs. (6) and (9), and finally applying the Borel transformation to the variable p 2 to eliminate the subtraction terms and suppress the contributions from higher states and the continuum. As a result, the desired sum rules for the states 1S + 2S + 1P are obtained: where M 2 is the Borel mass parameter, and QC D The derived sum rules in Eq. (13) contain six unknown parameters for the ground state and excited state baryons. However, Eq. (13) is broken down into two sum rule equations, each containing six unknowns: three masses and three pole residues. To solve this system of equations linearly, where, The whole set of six equations above makes the system's solution extremely difficult. Rather, to solve this system of equations, we picked up an appropriate value of s 0 at which the physical part of the correlation function is saturating only the ground state, and the other states, 2S and 1P are still remain in the continuum zoom. To put it differently, the initial guess of s 0 solved the system of equations and reproduced From there, we can see that the values of the ground state masses (m g ) and the corresponding pole residues (λ g ) of spin-(3/2) doublyheavy baryons * cc , * bb , * bc , * cc , * bb , and * bc have begun to appear and converge quickly within 5 GeV 2 ≤ M 2 ≤ 30 GeV 2 . The obtained results for (m g ) and (λ g ) with the corresponding s 0 values are also given in [44].
By doing so, we reduced the number of equations and the number of unknowns to four: two excited state masses m,m, and two excited state pole residues λ,λ, while we are still keeping the states 2S and 1P in the continuum zoom. To put it differently, we will solve the resulting four equations algebraically in terms (m g ) and (λ g ).
While solving the remaining four equations algebraically, we gradually increased the values of continuum zoom s 0 , and we are still keeping the Borel mass parameter M 2 running within the interval, 5 GeV 2 ≤ M 2 ≤ 30 GeV 2 . Bearing in mind, the results of the excited state masses and the excited state pole residues of the spin-(3/2) doubly-heavy baryons in Eq. (14) should also not depend on the Borel mass parameter M 2 . In fact, within the interval, 5 GeV 2 ≤ M 2 ≤ 30 GeV 2 we observed that there are different working regions where the results are supposedly acceptable and depend only moderately on the Borel mass parameter. For example, in Fig. 1,  Fig. 5 The dependency of the excited states 2S, 1P masses, and pole residues of the doubly-heavy of spin-(3/2) * cc baryon on the Borel mass M 2 at different values of the continuum threshold s 0 if we set the working region to be varied in the interval 18 GeV 2 ≤ M 2 ≤ 22 GeV 2 , clearly from this figure, we can still see that the excited state masses are rather stable with the variations of the Borel mass parameter.
However, we noticed that the excited state masses and the corresponding pole residues of the spin-(3/2) doubly-heavy baryons began to appear and are rapidly converge within the interval 25 GeV 2 ≤ M 2 ≤ 27 GeV 2 . The obtained results for the excited states 2S and 1P are presented in a series of Figs. 2, 3, 4, 5, 6 and 7 with respect to M 2 , along with the appropriate s 0 values, so that the reader might, if so motivated, repeat our calculations.

Numerical analysis
Using the QCD sum rules Eqs. (13) and (14), the masses and pole residues of the excited states of spin-(3/2) doublyheavy baryons * cc , * bb , * bc , * cc , * bb , and * bc are calculated. Usually, in numerical calculations, the QCD sum rules approach involves a set of fixed input parameters.  [48]. Moreover, the QCD sum rules technique also includes two auxiliary parameters whose Fig. 6 The dependency of the excited states 2S, 1P masses, and pole residues of the doubly-heavy of spin-(3/2) * bc baryon on the Borel mass M 2 at different values of the continuum threshold s 0 working regions must be determined, specifically, the continuum threshold s 0 , and the Borel parameter M 2 .
The physical working regions for these parameters should, however, be managed in such a way that the involved physical quantities are roughly independent of these parameters. The working regions of s 0 and M 2 are obtained by the sum rules analysis using the conventional criteria of the QCD sum rules approach, means that, the perturbative contributions must be greater than the condensate contributions, and the pole contributions must be greater than the continuum contributions. As a result, the bounds of M 2 are obtained such that M 2 must be large enough to satisfy the convergence criterion of OPE. Therefore, the working regions of the continuum threshold s 0 and the Borel mass M 2 must be found with the condition that the masses and pole residues are stable with respect to variations of these parameters.
In Figs. 2, 3, 4, 5, 6 and 7, we showed our numerical findings for the masses and the pole residues of spin-(3/2) doublyheavy excited state baryons. The masses and the pole residues of the excited states 2S and 1P of spin-(3/2) doubly-heavy baryons are represented as a function of M 2 for different values of s 0 . The expected masses and the pole residues of the excited states 2S and 1P showed good stability with respect to the obtained physical working regions of M 2 , and s 0 .
As discussed in the previous section, according to our findings, these requirements are met when M 2 is changed from 25 GeV 2 ≤ M 2 ≤ 27 GeV 2 for all the considered spin-(3/2) doubly-heavy baryons, while for example in Fig. 1, s 0 is modified in the range: s 0 = (22 − 26) GeV 2 to saturate  Table 1 The masses m andm of the excited states (2S) and (1P) of spin-(3/2) doubly-heavy baryons * cc , * bc , and * bb (in units of GeV)  the masses and the pole residues of * cc (2S), and * cc (1P) baryon excited states.
Also, for example in Figs. 2 and 3, considering the heaviest baryons, namely the doubly-heavy baryons including the heavy quark (b), the values of the continuum zoom s 0 should be increased to a sufficiently value in order to saturate again the excited states 2S and 1P.
Finally, in Tables 1, 2, 3 and 4 we displayed our predictions along with the other theoretical approaches for comparison. Our findings are in line with the findings of other methodologies.

Conclusion
In the context of the QCD sum rules method, we numerically calculated the excited state masses and the pole residues 2S and 1P of spin-(3/2) doubly-heavy * cc , * bb , * bc , * cc , * bb , and * bc baryons. The originality of the present work is that the contributions of 1S, 2S, and 1P states in Eqs. (13) and (14) are considered simultaneously.
On the whole, there should be a total of six unknowns: three masses and three pole residues. However, we only have two equations for two independent Lorentz structures pg μν and g μν . So, we get the remaining four equations by doing a differentiation for all sides of the sum rule equations with respect to , and then we solve these equations simultaneously.
We have used a slightly different method to solve this system of equations. We set a suitable value of s 0 and then we saturated the physical part of the correlation function with ground state contributions only. The spectroscopic parameters were then used as input parameters for the remaining set of equations. As a result, we have worked out a system of four equations as an alternative to six equations. Then, we have calculated the masses and the pole residues of the excited states 2S and 1P.
From figures, we have noticed that the masses of the 1P states is slightly exceeding the masses of the 2S states in the doubly-heavy baryons that are containing one light (u or d) quarks. Whereas, the masses of the 1P states is slightly less than the masses of the 2S states in the doubly-heavy baryons that are containing the strange quark. On the other hand, from figures, if one look at the values of the pole residues of the 2S states and that of the 1P states in the same region of s 0 values, one can notice that the values of the pole residues are of the same order to each other, and this is due to the reason that we have mentioned in Sect. 2. That is, the interpolating current in this work have been chosen in such a way that this current can couple to doubly-heavy baryons of both positive and negative parities.
By comparing the QCD sum rules results in this work with the results of the hyper-central constituent quark models [49,50], we have noticed that the 2S and 1P mass states of the doubly-heavy baryons are exceeding the mass states of the doubly-heavy baryons obtained within the potential models. Finally, we hope that future LHC experimental results can be compared to the current estimates.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data underlying the results are available as part of the article and no additional source data are required.] 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 copy-right holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .