Using the QCD sum rules method, we estimate the mass and residues of the first radial excitations of octet baryons. The contributions coming from the ground state baryons are eliminated by constructing the linear combinations of the sum rules corresponding to different Lorentz structures. Our predictions of the masses of the first radial excitations of octet baryons are in good agreement with the data.
1. Introduction
At present time, many radial excitations of the mesons and baryons which carry the same spin-parity quantum numbers as the ground states are observed in experiments [1]. The identification and investigation of the properties of the radial excitations on the background of the hadronic continuum are a quite difficult task in the experiments. These resonances are strongly coupled to the two, three hadrons, which leads to the large decay widths. Theoretical study of the radial excitations is also a challenging problem. One of the powerful theoretical methods for studying the properties of the hadrons is the sum rules method. The QCD sum rules method introduced in [2] was extremely useful for studying the properties of mesonic systems. This method was extended to ground state baryons in pioneering work [3]. It is well known that the choice of interpolating current which carries the same quantum numbers as appropriate baryon is not uniquely defined. The most general form of the interpolating currents for octet baryons is found in [4]. It is contemplated whether this method can be applied to the radial excitations of baryons. Note that the radial excitations of mesons within finite energy rules have been studied in [5, 6]. Recently, radial excitations of heavy-light mesons have been investigated in detail via the QCD sum rules in [7]. Recently, this method was applied for the determination of the masses of first radial excitation of mesons and nucleon by using the least square fitting method [8]. In this paper, we calculate the mass and residues of the first excited states of octet baryons within the QCD sum rules method. We modify the spectral representation from the hadronic part by taking into account two poles which correspond to the ground and first excited state baryons, respectively. Then, the method for eliminating the ground state contributions is employed for the determination of the mass and residues of the excited state baryons.
The paper is organized as follows. In Section 2, we derive the mass sum rules for octet baryons by including the first radial excitation baryons. Section 3 is devoted to the numerical analysis of obtained rules. Section 4 contains the summary of our results and conclusions.
2. Mass Sum Rules for Octet Baryons with Their Radial Excitation
For determination of the mass and residues of the mesons and baryons, usually two-point correlation function is used. Following the sum rules method strategy for extracting the mass and residues of radially excited octet baryons, we consider the following two-point correlation function:(1)Πq=i∫eipx0∣Tηxη¯0∣0d4x,where η(x) is the interpolating currents of the octet baryons. The most general forms of the interpolating currents for the octet baryons are [4](2)ηpx=2ϵabc∑l=12uaTCA1ldbA2luc,ηnx=ηpu⟷d,ηΣ+x=ηpd⟶s,ηΣ-x=ηnu⟶s,ηΞ0x=ηΣ+u⟷s,ηΞ-x=ηΣ-d⟷s,ηΣ0x=2ϵabc∑l=12uaTCA1lsbA2ldc+daTCA1lsbA2luc,where A11=I, A12=A21=γ5, and A22=β in which β is arbitrary parameter and a, b, and c are color indices.
The phenomenological part of the correlation function can be obtained by inserting a full set of baryons carrying the same quantum numbers as the interpolating current. Isolating the ground state and its first radial excitation from phenomenological side for the correlation functions we have(3)Πphys=λ2p+mp2-m2+λ12p+mp2-m12+⋯.
Here λ(λ1) and m(m1) are the residue and mass of the ground (first radial excitation) state baryon and ⋯ stands for the contributions of higher states and continuum. In derivation of (3), we used(4)0∣η∣Bp≡λup.
In order to suppress the higher states and continuum contributions the Borel transformation is applied. After performing Borel transformation from (3) we have(5)Bp2Πphys≡λp2p+me-m2/M2+λ12p+m1e-m12/M2+⋯.
The correlation function in terms of quark-gluon degrees of freedom (theoretical part) is calculated in the deep Euclidean domain (p2≪0) with the help of the operator product expansion (OPE), which contains the perturbative and nonperturbative (vacuum condensate) contributions. For structures p and I it can be written as(6)ΠiOPEp2=Πipertp2+Πinon-pertp2,where the invariant functions Π1(Π2) correspond to the coefficient of the structure p(I).
The expressions of various terms entering to the right side of (6) are calculated in numerous works (see, e.g., [9, 10]). Performing Borel transformation over p2 in theoretical part of the correlation function and equating the coefficients of the Lorentz structures p and I, one can obtain the sum rules for the mass.
The sum rules for the structures p and I can be written as(7)Π1M2=λ2e-m2/M2+λ′2e-m′/M2=C1E2M6+C2msq¯qE0M2+C3msm02q¯q+C4qq2+C5m02qq2M2+C6αSπG2M2,Π2M2=λ2me-m2/M2+λ′2m′e-m′2/M2C1′msM6E2x+C2′q¯qM4E1+C3′m02q¯qE0M2+C4′msqq2+C5′msαsπG2M2+C6′q¯qαsπG2.In (7), we consider the contributions of operators up to dimension 6. Obviously, Ci for different members of the octet baryons are different. Note that, in the derivation of (7), we set the light quark masses to zero; mu=md=0.
The continuum subtraction is done using the quark-hadron duality ansatz. In (7) the function E2s0/M2 is described by the higher states contributions and continuum contributions and determined as E2x=1-e-x∑xn/n!, where s0 is the effective continuum threshold. The coefficients Ci in (7) are given as follows (see [9, 10]).
For N(8)C1=1256π45+2β+5β2,C2=0,C3=0,C4=166-1+β2+-1+β2,C5=-12412-1+β2+-1+β2,C6=1256π25β2+2β+5.
For Σ(9)C1=1256π45+2β+5β2,C2=132π25+2β+5β2γ-12-1+β2,C3=196π2γ4+4β+4β2+21β2-1+316π2-1+β2γE-lnM2μ2,C4=166-1+β2γ+-1+β2,C5=-12412-1+β2γ+-1+β2,C6=1256π25β2+2β+5.
For Ξ(10)C1=1256π45+2β+5β2,C2=316π2-2-1+β2+1+β2γ,C3=196π215-γβ2-10γβ-15+f+316π2-1+β2γE-lnM2μ2,C4=f66-1+β2+-1+β2γ,C5=-124γ12-1+β2+-1+β2γ,C6=1256π25β2+2β+5.
For Λ(11)C1=1256π45+2β+5β2,C2=196π235+2β+5β2γ+41+4β-5β2,C3=-196π241+β+β2γ+-5β2+4β+1-116π21-β2γE-lnM2μ2,C4=-1181-β2γ1+5β+13+11β,C5=1721-β81+2βγ+25+23β,C6=1256π25β2+2β+5.Moreover, the coefficients Ci′ for the members of the octet baryons are given as follows.
For N(12)C1′=0,C2′=-116π27β2-2β-5,C3′=316π2β2-1,C4′=0,C5′=0,C6′=-128819β2+10β-29.
For Σ(13)C1′=164π4β-12,C2′=-116π26+γβ2-2γβ-6-f,C3′=316π2β2-1,C4′=165-3γβ2+2β+5+3γ,C5′=-1128π21-β2,C6′=-128824-5γβ2+10γβ-24+5γ.
For Ξ(14)C1′=332π4β2-1,C2′=-116π26γ+1β2-2β-6γ-1,C3′=3γ16π2β2-1,C4′=γ23-γβ2+2β+3+γ,C5′=364π2β2-1,C6′=-128824γ-5β2+10β-24γ+5.
For Λ(15)C1′=1192π411β2+2β-13,C2′=-148π210+11γβ2+-8+2γβ-2+13γ,C3′=-116π2-1-2γβ2+1+2γ,C4′=11815-5γβ2+6+4γβ+15+γ,C5′=1384π213β2-2β-11,C6′=-18644+53γβ2+40-10γβ-44+4γ,where γE is the Euler constant, γ=〈s¯s〉/〈u¯u〉, and μ is the renormalization scale parameter whose value is taken as μ=0.5GeV. Now using (5) and (7) one can determine the masses and residues of the octet baryons. We have four unknowns (two masses and two residues). Therefore we need four equations to find these unknowns.
The first two equations are obtained by equating the coefficients of the structures p and I in both representations of the correlation functions; that is,(16)λ2e-m2/M2+λ′2e-m′2/M2=Π1,λ2me-m2/M2+λ′2m′e-m′2/M2=Π2.
The remaining two equations can be obtained by taking derivatives with respect to -1/M2 from both sides of (16). Then we have(17)m2λ2e-m2/M2+λ′2m′2e-m2/M2=Π1′,λ2m3e-m2/M2+λ′2m′3e-m2/M2=Π2′.
From these equations we get(18)m′2=Π2′-mΠ1′Π2-mΠ1,λ′2=1m′2-m2em′2/M2Π1′-m2Π1.
To obtain the mass and residues of the radial excitations from the sum rules, we take the mass of the ground state as an input parameter. Note that the zero width approximation is assumed.
3. Numerical Analysis
The input parameters used in the above coefficients include the strange quark mass and quark condensates. In our numerical calculations, we use the following values of these parameters (see [3, 11, 12]):(19)m02=0.8±0.2GeV2,q¯q2GeV=-277-10+12MeV3,ms2GeV=95±5MeV,γ=s¯sq¯q=0.8±0.2,καsq¯q2=5.8±1.8×10-4GeV6.In numerical calculations we take κ=1.
The sum rules for the mass and residues of the radially excited baryons contain three auxiliary parameters in addition to these input parameters: the Borel mass M2, arbitrary parameter β (in expressions of the interpolating current), and continuum threshold s0. Usually the continuum threshold is related to the energy of the first excited state. In our calculations for s0 we have used s0=mground+ΔGeV where Δ varies between 0.3 and 0.8. These values of s0 include only the mass of the first radial excitation, while the higher excitations are included in the continuum states. The working interval of M2 is obtained by the following way. The lower bound of M2 is determined by demanding the convergence of the operator product expansion, while the upper bound is obtained by requiring that the continuum contribution remains subleading. Our calculations lead to the following working region of M2:(20)N:1.2GeV2≤M2≤2.2GeV2,Σ:1.4GeV2≤M2≤2.4GeV2,Ξ:1.8GeV2≤M2≤3GeV2,Λ:1.4GeV2≤M2≤2.4GeV2.
Having determined the working regions for Borel mass parameter M2, we can calculate the mass and residues of the first radial excitation of octet baryons. As an example, in Figures 1 and 2, we present the dependence of the mass of the radial excitation of N and Λ baryons on M2 at fixed values of β and s0. From these figures, we obtained that the masses of the radial excitations of N and Λ exhibit good stability with respect to the variation of M2 in the working region. We perform the same analysis for another fixed values of s0 and find that the results change about 3.5%. We also performed analysis for the other members of the octet baryons and obtained that the dependence of the mass of radial excitation of octet baryons on M2 is rather weak.
Dependence of mN′2 on the Borel mass parameter M2 for the fixed value of the continuum threshold s0=3.3GeV2 is depicted for the several fixed values of parameter β.
Dependence of mΛ′2 on the Borel mass parameter M2 for the fixed value of the continuum threshold s0=4.2GeV2 is depicted for the several fixed values of parameter β.
In order to determine the optimal working interval of the parameter β, we study the dependence of the mass and residues of octet baryons and their radial excitation on cosθ, where β=tanθ. Note that we use cosθ by the following reason. Exploring the whole region in β (-∞,∞) is equivalent to the very restricted domain (-1,1).
In Figures 3 and 4, we present the dependence of mN′2 and mΛ′2 on cosθ for two fixed values of M2 and at fixed value s0. From these figures, we observe that in the domain -0.6≤cosθ≤0.8 the mass of N′ and Λ′ baryons remains stable with respect to the variation of cosθ and we deduce the following values for their masses:(21)mN′2=2.1±0.2GeV2,mΛ′2=2.7±0.1GeV2.
Dependence of mN′2 on cosθ at the fixed value of the continuum threshold s0=3.3GeV2 is shown for various values of the Borel mass parameter M2.
Dependence of mΛ′2 on cosθ at the fixed value of the continuum threshold s0=4.2GeV2 is shown for various values of the Borel mass parameter M2.
Performing similar analysis for the mass of the Σ and Ξ baryons we obtained(22)mΣ′2=2.8±0.1GeV2,mΞ′2=3.4±0.2GeV2.
Finally, we can determine the residues of the radial excitation baryons. For this aim, we used (18). Using the working region for M2 and at given values of s0 we studied the dependency of the residue square on cosθ and we obtained the following values:(23)λN′2=2.0±0.5×10-3,λΣ′2=5.0±2.0×10-3,λΛ′2=5.0±2.0×10-3,λΞ′2=9.0±4.0×10-3.
Comparing our predictions on the mass of the radial excitation with the octet baryons we see that the results are in good agreement with the experimental data. Remember that experimentally the masses of the radial excitations of octet baryons are mN′=(1.43±0.02)GeV, mΛ′=(1.63±0.07)GeV, mΣ′=(1.66±0.03)GeV, and mΞ′=(1.950±0.015)GeV.
We observed that our predictions on mass for the radial excitations are in good agreement with the experimental data. We also calculate the residues for the radial excitations octet baryons.
Our final remark to this section is as follows. We perform our analysis without taking into account the radiative corrections to the two-point correlation function. The one-loop radiative corrections to the two-point correlation function for nucleon for Ioffe current (β=-1) is calculated in [13]. When we take into account these corrections, our results change by (5–8)%.
4. Conclusion
In the present work, we estimate the masses and residues of the first radial excited states of the octet baryons. The QCD sum rules are modified by including the contribution of the first radial excited states in addition to the ground state. Our predictions on the masses of the first radial excited states baryons are in good agreement with the experimental data. Hence, QCD sum rules work quite well not only for ground state baryons but also for the radial excitations. The obtained results for the residues can be checked by studying the electromagnetic or strong transitions of radial excitation to the ground state baryons.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
OliveK.Review of Particle Physics201438909000110.1088/1674-1137/38/9/090001ShifmanM. A.VainshteinA. I.ZakharovV. I.QCD and resonance physics. Theoretical foundations1979147538544710.1016/0550-3213(79)90022-12-s2.0-0001020047IoffeB. L.Calculation of baryon masses in quantum chromodynamics19811882317341Erratum to: Nuclear Physics B, vol. 191, pp. 591-592, 198110.1016/0550-3213(81)90259-5ChungY.DoschH. G.KremerM.SchallD.Baryon sum rules and chiral symmetry breaking19821971557510.1016/0550-3213(82)90154-72-s2.0-0002158303KrasnikovN.PivovarovA.The influence of the analytical continuation effects on the value of the QCD scale parameter Λ extracted from the data on charmonium and upsilon hadron decays19821162-316817010.1016/0370-2693(82)91001-2GorishnyS. G.KataevA. L.LarinS. A.Next-next-to-leading perturbative QCD corrections and light quark masses19841355-645746210.1016/0370-2693(84)90315-02-s2.0-0001534307GelhausenP.KhodjamirianA.PivovarovA. A.RosenthalD.Radial excitations of heavy-light mesons from QCD sum rules201474, article no. 297910.1140/epjc/s10052-014-2979-z2-s2.0-84908286499JiangJ.-F.ZhuS.-L.Radial excitations of mesons and nucleons from QCD sum rules20159207400210.1103/physrevd.92.074002LeeF. X.LiuX.Predictive ability of QCD sum rules for excited baryons200266110.1103/physrevd.66.014014AlievT.ÖzpineciA.SavciM.Octet baryon magnetic moments in light cone QCD sum rules200266101600210.1103/physrevd.66.016002NarisonS.Strange quark, tachyonic gluon masses and Vu s from hadronic tau decays20056261–410111010.1016/j.physletb.2005.08.085NarisonS.Power corrections to αsMT, Vu sand m-s20096731303610.1016/j.physletb.2009.01.062DrukarevE. G.RyskinM. G.SadovnikovaV. A.On the loop approximation in nucleon QCD sum rules201578784985310.1134/s1063778815070054