-
In 2018, the BELLE Collaboration made an exciting announcement regarding the discovery of the
Ω(2012) hyperon. This discovery was based on theΩ∗−→Ξ0K− andΩ∗−→Ξ−K0s decay channels, with a measured mass ofm=2012.4±0.7(stat)±0.6(sys)MeV and decay width ofΓtot=6.4+2.5−2.0(stat)±1.6(sys)MeV [1]. However, knowing only the mass of the state is not sufficient to determine the quantum numbers of a state. For instance, within the QCD sum rule method, the mass of theΩ(2012) baryon is estimated, assuming it to be either the 1P or 2S excitation state [2]. Both assumptions yield the same mass value, although the estimated residues differ. Thus, additional physical quantities, such as the decay width, are necessary to identify the quantum numbers of newly discovered particles.In a previous study [3], the
Ω(2012)→Ξ0K− transition was investigated, and its corresponding decay width was estimated by considering two possible scenarios forΩ(2012) : either a1P or2S state. A comparison of the total decay widths obtained in this work led to the conclusion thatΩ(2012) is itself aJP=32− state. Moreover, predictions from various theoretical models also converge on the likely quantum numbersJP=32− for the observed state [4−15].In this study, considering
Ω(2012) as theJP=32− state, the strong couplings of the SU(3) partners of this state are investigated within the framework of light cone sum rules (LCSRs) using the distribution amplitudes (DAs) of the octet baryon. Note that this problem was also investigated in [16] using the flavor SU(3) symmetry approach.The structure of this paper is as follows. Section II introduces the LCSRs for the strong couplings of the transition
32−→12+ + pseudoscalar mesons. Section III provides a numerical analysis of the LCSRs, focusing on the relevant strong couplings. Within this section, we also present the computed values of the decay widths based on the obtained coupling constants. Additionally, we compare our results with those obtained using the flavor SU(3) symmetry method. Finally, our conclusions are summarized in Section IV. -
To calculate the strong couplings of
SU(3) partners, denoted as32− states in the following discussions, we introduce the vacuum-to-octet baryon correlation function:Πμν(p,q)=i∫d4xeiqx⟨0|T{ημ(0)Jν(x)}|O(p)⟩,
(1) where
ημ represents the interpolating current of the decuplet baryons,Jν=ˉq1γνγ5q2 is the interpolating current of the pseudoscalar mesons, and|O(p)⟩ represents the octet baryon state.The interpolating current of the decuplet baryons can be written as
ημ=εabcA{(qaT1Cγμqb2)qc3+(qaT2Cγμqb3)qc1+(qaT3Cγμqb1)qc2},
(2) where
a,b,c are the color indices, C is the charge conjugation operator, and A is the normalization factor. The quark content of the decuplet baryons and the normalization factor A are presented in Table 1.A q1 q2 q3 Δ+ √1/3 u u d Σ+(3/2) √1/3 u u s Σ0(3/2) √2/3 u d s Σ− (3/2)√1/3 d d s Ξ0 (3/2)√1/3 s s u Ξ−(3/2) √1/3 s s d Ω− 1 s s s Table 1. Quark content of the decuplet baryons and the normalization factor A.
To derive the LCSRs for the strong coupling constants, the correlation function is computed in two ways: in terms of hadrons and in terms of quark-gluon fields within the deep Euclidean domain. By applying the quark-hadron duality ansatz, the relevant sum rules can be derived.
Strong coupling constants appear in the double dispersion relation for the correlation function given in Eq. (1). Hence, to calculate these constants, the double dispersion relation for the correlation function must be calculated. The double dispersion relation is obtained via analytical continuation of the imaginary part of the corresponding invariant amplitudes with respect to the variables
p′2 andq2 in the spin-3/2 and pseudoscalar meson channels, respectively.Before delving into the details of the calculations, it is important to highlight the following aspect: the interpolating current for the decuplet baryons interacts not only with the ground positive parity states
JP=32+ , but also with the negative parity statesJP=32− and even with the statesJP=12− .To eliminate the contributions from unwanted states,
JP=32+ andJP=12− , a technique involving the linear contributions of different Lorentz structures is employed (for more details about this approach, refer to [17]).Following the standard procedure, we insert the total set of baryons with
JP=32 into the correlation function along with the corresponding pseudoscalar mesons. Then, we obtainΠμν(p,q)=∑i=±⟨0|ημ|32i(p′)⟩m2i−p′2⟨32i(p′)P(q)|O(p)⟩m2P−q2×⟨0|Jν(x)|P(q)⟩,
(3) where summation is over positive and negative states, and
mP is the mass of the corresponding pseudoscalar mesonP with momentum q. The matrix elements in the above equation are defined as⟨0|ημ|32+(p′)⟩=λ+uμ(p′),⟨0|ημ|32−(p′)⟩=λ−γ5uμ(p′),⟨32+(p′)P(q)|O(p)⟩=g+ˉuα(p′)u(p)qα,⟨32−(p′)P(q)|O(p)⟩=g−ˉuα(p′)γ5u(p)qα,⟨0|Jν|P(q)⟩=ifPqν,
(4) where
λ± are the residues of the related32± baryons,g± represents the coupling constants of theJP=32± baryons with the octet baryons and pseudoscalar mesons,fP is the decay constant of the pseudoscalar meson and q denotes its 4-momentum, anduμ(p′) andu(p) are the Rarita-Schwinger and Dirac spinors, respectively. Performing summation over the spins of the Rarita-Schwinger spinors using the formula∑s′uμ(p′,s′)ˉuα(p′,s′)=−(⧸p′+m)[gμα−13γμγν−2p′μp′α3m2+p′μγα−p′αγμ3m],
(5) and using Eqs. (3) and (4), we can obtain an expression for the correlation function from the hadronic part. It should be reminded that the interpolating current interacts not only with spin
32 states, but also with spin12 states.Using the condition
γμημ=0 , it can easily be shown that⟨0|ημ|12(p′)⟩∼[αγμ−βp′μ]u(p′).
(6) It follows from this equation that any structure containing
γμ orp′μ is "contaminated" by the contributions of spin12 -states. Hence, to remove the contributions of spin12 -states, such structures are all discarded.Another problem is all Dirac structures not being independent of each other. To overcome this issue, Dirac structures must be arranged in a specific order. In this study, we choose the ordering
γμ⧸p′⧸qγν .Keeping this in mind, and using Eqs. (3), (4), and (5), we obtain the correlation function from the phenomenological part as follows:
Πμν=λ+g+(−⧸q+m++mO)qμqνfP(m2+−p′2)(m2P−q2)u(p)+λ−g−(⧸q+m−−mO)qμqνfP(m2−−p′2)(m2P−q2)u(p)+...,
(7) where
mO is the mass of the relevant octet baryon, andm+(m−) is the mass of the spin-32 positive (negative) parity baryon. Here,... indicates the contributions of the excited states and continuum.As a final step, we must eliminate the contributions of
JP=32+ states. For this purpose, the linear combinations of the invariant functions corresponding to different Lorentz structures are considered.We now turn our attention to the calculation of the correlation function using operator product expansion (OPE) in the deep Euclidean region for the variables
p′2=(p−q)2 andq2≪0 . To calculate OPE, the explicit forms of the interpolating current are placed in the correlator, and possible contractions are performed between quark fields using Wick's theorem. As an example, for the correlation function ofΣ0(3/2)→NK , we obtainΠμν=√23∫d4xeiqxϵabc(Cγμ)αβ(γνγ5)ρσ×{⟨0|uaα(0)ubσ(x)dcβ(0)|N⟩Sγρ(−x)+⟨0|ucγ(0)ubσ(x)dcα(0)|N⟩Sβρ(−x)+⟨0|uaβ(0)ubσ(x)dcγ(−x)|N⟩Sαρ(−x)},
(8) where
S(−x) is the strange quark propagator. From this expression, it follows that the OPE results are obtained via convolution of the quark propagator to the sum of the nucleon DAs, obtained from theϵabc⟨0|uaα(0)ubβ(x)dcγ(0)|N⟩ matrix element. A diagrammatic description of Eq. (8) is given in Fig. 1. Once we use the explicit expressions of the quark propagators and the definition of the DAs of octet baryons, the following master integral appears in the coefficients of different Lorentz structures:Figure 1. Diagrammatic representation of the correlation function. The wavy lines denote the external currents, solid lines correspond to the quark fields, and shaded regions correspond to the DAs of the nucleon.
In,k=∫duuk[m2−(pu−q)2]n;n=1,2,3.
For the calculation of the double spectral densities, it is sufficient to find the double spectral representations of the master integrals. The details of the spectral density calculations for the
n=1 case are presented in Appendix A. The cases ofn=2 andn=3 are calculated in a similar manner.The invariant amplitudes are related to the spectral densities via the double dispersion relation as follows:
Π[(p−q)2,q2]=∫ds1∫ds2ρ(s1,s2)[s1−(p−q)2](s2−q2)+⋯
(9) The spectral density can be obtained from
Π[(p−q)2,q2] by applying two subsequent double Borel transformations (for more details of the calculation, see Appendix A.)Matching the OPE results with the double dispersion relations for the relevant Lorentz structures of the hadrons, applying the quark-hadron duality ansatz, and performing double Borel transformation with respect to the variables
−(p−q)2 and−q2 , we obtain the LCSRs for the relevant coupling constants whose explicit form can be written asg−=em2−/M21em2P/M22fPλ−(m++m−)1π2∫s00ds1×∫min(s′0,t2(s1))t1(s1)ds2e−s1/M21e−s2/M22Ims1Ims2×{Π1(m+−mO)+Π2},
(10) where
Π1 andΠ2 are the invariant functions of the Lorentz structures⧸qqμqν andqμqν , respectively, andt1,2=s1+m2O∓2mO√s1−m2,
where m is the corresponding mass of light quarks. Here,
s′0 is the continuum threshold in the pseudoscalar meson channel. The continuum thresholds′0 is chosen as a mass square of the first radial excitation of the corresponding pseudoscalar meson. Finally, note that in themO→0 limit, the applied method must be modified (for more details, see [18,19]). -
This section is devoted to the numerical analysis of the coupling constants derived in the previous section within the LCSRs. The main nonperturbative input of the considered LCSRs is the DAs of the octet baryons, namely, N, Σ, and Ξ. The explicit expressions of the relevant DAs are obtained in [20−23]. The DAs contain the normalization constants f,
λ1 , andλ2 , which are determined from the analysis of mass sum rules as well as lattice QCD [24, 25]. The normalization constant of the leading twist f, (for N, Σ, and Ξ baryons) is defined via the matrix element of the local current (all quark fields are at the same point).ϵabc⟨0|(qa1(0)C⧸nqb2(0))γ5⧸nqc3(0)|O(p)⟩=f(pn)⧸nu(p).
(11) Moreover, the DAs of higher twist contributions involve two additional normalization constants,
λ1 andλ2 , which are defined as the matrix elements of local three quark twist-four operators,ϵabc⟨0|(qa1(0)Cγμqb2(0))γ5γμqc3(0)|O(p)⟩=λ1mOu(p),ϵabc⟨0|(qa1(0)Cσμνqb2(0))γ5σμνqc3(0)|O(p)⟩=λ2mOu(p),
(12) where n is the light-like vector, and
u(p) is the Dirac bispinor.The normalization constants f,
λ1 , andλ2 for the Λ baryon can be obtained from Eqs. (11) and (12) via the following replacements:C⧸n→Cγ5⧸nγ5⧸n→⧸nfor f
(13) Cγμ→Cγ5γμγ5γμ→γμfor λ1
(14) Cσμν→Cγ5γ5σμν→1for λ2
(15) In our analysis, we use the parameter values obtained from lattice QCD that are presented in Table 2 for completeness. The masses of the
SU(3) partners ofΩ(2012) are obtained in [16] and presented below.f λ1 λ2 N 3.54 −44.9 93.4 Σ 5.31 −46.1 85.2 Ξ 6.11 −49.8 99.5 Λ 4.87 −42.2 98.9 Table 2. Numerical values of the parameters
f,λ1 , andλ2 , given in units of10−3GeV2 .m−={1700±90MeVforΔ,1805±100MeVforΣ(3/2),1910±110MeVforΞ(3/2),2012.4±0.9MeVforΩ[26].
These mass values are used in our numerical analysis. For the masses of the ground state baryons, we adapt values from the PDG [26]. In addition, the value of the quark condensate is taken as
⟨ˉqq⟩=−(246+28−19MeV)3 [17] and⟨ˉss⟩=0.8⟨ˉqq⟩ [27].The residues of the negative parity
JP=32− baryons are related to the residues of the radial excitations of the decuplet baryons as follows:λ−=λrad√m−−m+m−+m+.
The residues of the radial excitations of the decuplet baryons are calculated in [2]. Using these results, we can easily determine the residues of the
JP=32− baryons.The working regions of the Borel mass parameters and continuum thresholds,
s0 ands′0 , used in the numerical analysis are presented in Table 3. Determination of the working regions of the Borel parameters is based on the criteria that both power corrections and continuum contributions should be suppressed. Moreover, the continuum thresholds0 is obtained under the condition that the mass of the considered states reproduces the experimental values with an accuracy of approximately 10%.Borel mass parameters Continuum threshold Continuum threshold M21/GeV2 M22/GeV2 s0/GeV2 s′0/GeV2 Δ→Nπ 3÷4 0.25÷0.35 5.0±0.2 1.7 Σ(3/2)→NK 3÷4 0.25÷0.35 5.5±0.2 2.0 Σ(3/2)→Λπ 3÷4 0.42÷0.44 5.5±0.2 1.7 Σ(3/2)→Σπ 3÷4 0.42÷0.44 5.5±0.2 1.7 Ξ(3/2)→ΛK 3÷4 0.45÷0.47 6.0±0.2 2.0 Ξ(3/2)→ΣK 3÷4 0.60÷0.65 6.0±0.2 2.0 Ξ(3/2)→Ξπ 3÷4 0.50÷0.60 6.0±0.2 1.7 Ω→ΞK 3÷4 0.55±0.65 6.5±0.2 2.0 Table 3. Working regions of the Borel mass parameters and continuum threshold
s0 .Having the values of all input parameters at hand, we can perform the numerical analysis of the relevant coupling constants. As an example, in Fig. 2, we present the dependency of the coupling constant on
M21 at fixed values of the continuum thresholdss0 ,s′0 , andM22 for theΩ→ΞK transition, because this transition has already been discovered. From this figure, we observe that there is good stability of the coupling constant whenM21 varies in its working region (Table 3). The obtained coupling constants are presented in Table 4. The errors in the results for the coupling constants can be attributed to the uncertainties in the input parameters as well as errors to the Borel mass parametersM21 andM22 and continuum thresholds0 ands′0 .Figure 2. (color online) Dependency of the coupling constant of the
Ω(2012)→Ξ−K+ transition on the Borel mass parameterM21 at several fixed values of the Borel parameterM22 and the continuum thresholds0=6.5GeV2 .Decay channels g−/GeV−1 Γ/MeV (This study)Γ/MeV [16]Δ→Nπ 12±3 71.6×(1.0±0.5) 39−58 Σ→NK 6±2 11.1×(1.0±0.6) 7−12 Σ→Λπ 9±3 23.7×(1.0±0.6) 11−18 Σ→Σπ 5±1 4.5×(1.0±0.4) 4−7 Ξ→ΛK 10±2 15.5×(1.0±0.4) 5−10 Ξ→ΣK 6±2 2.7×(1.0±0.6) 2−5 Ξ→Ξπ 7±2 6.9×(1.0±0.5) 5−9 Ω→ΞK 12±3.5 7.4×(1.0±0.6) - Table 4. Decay widths of the
JP=32− baryons.Having determined the coupling constants, we can calculate the decay widths of the corresponding transitions. Using the matrix elements for the considered
32−→12++pseudoscalar meson transitions, the decay width can be written asΓ=g2−24πm2−[(m−−mO)2−m2P]|→p|3,
(16) where
|→p|=12m−√m4−+m4O+m4P−2m2−m2O−2m2−m2P−2m2Om2P,
is the momentum of octet baryon, and
mO andmP are the mass of the octet baryon and pseudoscalar meson, respectively. Using the values of the coupling constants obtained within this study, we estimate the decay widths of the relevant transitions summarized in Table 4. For comparison, we also present the results of the decay widths obtained from the flavorSU(3) analysis [16]. We would like to make the following remark at this point. From the expression of the decay width, it is evidently sensitive to the mass splitting among theSU(3) partners of theΩ(2012) and ground state baryons. Thus, for a fair comparison, we use the same mass values as in [16].Finally, we compare our results with the values obtained within the framework of the flavor
SU(3) method [16]. In this analysis, the coupling constant forΩ→ΞK is taken as the input parameter, and all the remaining couplings are expressed in terms of this coupling usingSU(3) symmetry relations. Using the experimental value of the decay widthΩ→ΞK , we can determine the coupling constant of this transition via Eq. (16), and hence all the other coupling constants can be determined. When we compare our results on the coupling constants and decay widths of the considered decays with those obtained within the flavorSU(3) analysis, we find that they are compatible within the uncertainties of the the model predictions. Small deviations in the results can be attributed to theSU(3) violation effects and uncertainties of the input parameters of the theory. Furthermore, our prediction of the decay width forΩ→ΞK is compatible with those observed by the BELLE Collaboration within the uncertainties [1]. Moreover, note that the coupling constant, and hence the decay width ofΩ→ΞK , within LCSRs method was calculated using the DAs of pseudoscalar mesons in [3]. However, in this study, we recalculate these quantities within the same framework using the DAs of the Ξ baryon. In this method, the calculations of the theoretical part of the sum rules can be achieved using only one quark propagator; however, in [3], two quark propagators were required, making the calculations difficult because each quark propagator contains many terms. Another advantage of the present method lies in dealing with the contributions of baryons with different parities, especially when mass splittings are small. In this method, no pollution arise due to negative parity baryons. However, with the methods used in [3], the problem of the separation of the contributions of positive baryons remains unsolved. Another difference between the two methods is that in this study, we consider both Borel mass parametersM21 andM22 but in [3],M21=M22 was considered. The uncertainties of the parameters entering the DAs of baryons are larger than those of meson DAs. Once the errors are minimized in the determination of these parameters, more precise results can be obtained. When we compare our results on the coupling constant forΩ→ΞK , we find that our result is consistent with that in [3] within the uncertainties. -
In conclusion, we employ the LCSR method to compute the strong coupling constants and decay widths for the
SU(3) partners of theΩ(2012) baryon in32−→12++pseudoscalar meson transitions. The "contamination" caused by theJP=32+ baryons are eliminated by considering the linear combinations of the sum rules obtained from different Lorentz structures. By comparing our decay width results with the findings of [16], we ascertain the compatibility of our decay width predictions with the outcomes of the flavorSU(3) symmetry analysis. The small discrepancy between the predictions of the two methods may be attributed to theSU(3) violation effects. Moreover, our estimated decay width for theΩ→ΞK transition is also compatible with the measurement of the BELLE Collaboration within the uncertainty limits. In addition, our result on the coupling constant forΩ→ΞK calculated using the DAs of Ξ is consistent with the prediction in [3], where the DAs of pseudoscalar mesons are used.Our results on the branching ratios can provide useful hints about the nature of the
SU(3) partners of theΩ(2012) baryon. -
We are grateful to Y. M. Wang for discussions on numerical analysis. We also thank A. Ozpineci for useful remarks.
-
Here, we provide the detailed derivation of the spectral density (see also [28]).
After applying the double Borel transformation over the variables
−p′2 and−q2 to Eq. (9), we obtainΠB1(M21,M22)=∫ds1∫ds2e−s1/M21−s2/M22ρ(s1,s2).
(A1) Before implementing the second double Borel transformation, we introduce the new variables
σ1=1M2i . The second Borel transformation can be performed over the new Borel parameterτi using the relationBτe−sσ=δ(1τ−s).
(A2) As a result, we have
Bτ1Bτ2ΠB1(M21,M22)=ρ(1τ1,1τ2).
(A3) Hence, the double spectral density can be obtained as follows:
ρ(s1,s2)=B1s1(σ1)B1s2(σ2)ΠB(1σ1,1σ2).
Let us now focus on the double spectral density for the
n=1 case. Using−(pu−q)2=−u(p−q)2−ˉuq2+uˉum2O,
where
ˉu=1−u ,I1,k can be written asI1,k=∫duuk[m2−u(p−q)2−ˉuq2+ˉuum2O]=∫duukD,
where m is the corresponding quark mass. Using the Schwinger representation for the denominator and performing the first double Borel transformation over the variables
−(p−q)2 and−q2 , we obtainI1,k=σk2(σ1+σ2)k+1exp[−m2Oσ1σ2σ1+σ2−m2(σ1+σ2)],=σk2(σ1+σ2)k+1exp[m2Oσ21+σ222(σ1+σ2)−(m2+m2O2)(σ1+σ2)],
where
σi=1M2i . To perform the second double Borel transformation, we use the relation√σ1+σ22π∫+∞−∞dxiexp[−σ1+σ22x2i−σimOxi]=exp[m2Oσ2i2(σ1+σ2)]. Then, we obtain
IB1,k=12π∫+∞−∞dx1∫+∞−∞dx2σk2(σ1+σ2)kexp[−σ1(m2+(mO+x1)2+x222)−σ2(m2+(mO+x2)2+x212)]=12π1Γ(k)∫+∞−∞dx1∫+∞−∞dx2∫∞0dttk−1σk2exp[−σ1(m2+(mO+x1)2+x222+t)−σ2(m2+(mO+x2)2+x212+t)]=12π1Γ(k)∫+∞−∞dx1∫+∞−∞dx2∫∞0dttk−1exp[−σ1(m2+(mO+x1)2+x222+t)]×(−∂∂t)kexp[−σ2(m2+(mO+x2)2+x212+t)].
After performing the second Borel transformation, we obtain the the spectral density corresponding to
I1,k :ρ1,k(s1,s2)=12π1Γ(k)(−∂∂s2)k∫+∞−∞dx1∫+∞−∞dx2∫∞0dttk−1δ[s1−(m2+(mO+x1)2+x222+t)]×δ[s2−(m2+(mO+x2)2+x212+t)]=12π1Γ(k)(−∂∂s2)k∫+∞−∞dx1∫+∞−∞dx2∫∞0dttk−1δ[s1−(m2+(mO+x1)2+x222+t)]×δ[s2−(m2+(mO+x2)2+x212+t)].
Using two Dirac delta functions, we can easily perform integrals over t and
x2 :ρ1,k(s1,s2)=12πΓ(k)mO(−∂∂s2)k∫+∞−∞dx1[s1−(m2+(mO+x1)2+x222)]k−1Θ[s1−(m2+(mO+x1)2+x222)],
where
x2=s2−s1mO+x1,
and
Θ(x) is the Heaviside step function, which restricts the integral overx1 between the limitsy±(s1,s2) , wherey±(s1,s2)=−m2O+s1−s2±√Δ2mO,
and
Δ=−m4O−(s1−s2)2+2m2O(−2m2+s1+s2).
As a result of the above summarized calculations, the spectral density can take the following form:
ρ1,k(s1,s2)=12π1Γ(k)1mO(−∂∂s2)k∫y+y−dx[(y+−x)(x−y−)]kΘ(Δ).
To evaluate the x integral, we introduce a new variable through the relation
x=(y+−y−)y+y−,
so that the spectral density can be written as
ρ1,k(s1,s2)=12πΓ(k)Γ(2k)1m2kO(−∂∂s2)k[Δk−12Θ(Δ)].
(A4) The double spectral densities for
I2,k andI3,k can be calculated using the following relations:I2,k=(−∂∂m2)I1,k,and,
I3,k=12(−∂∂m2)2I1,k,
(see also [17] for the calculation of the spectral densities
I2,k andI3,k ).
Analysis of strong decays of SU(3) partners of Ω(2012) baryon
- Received Date: 2024-03-26
- Available Online: 2024-08-15
Abstract: We estimate the coupling constants and decay widths of the