Radiative decays of the p-wave charmed heavy baryons

The radiative decays of the p-wave charmed heavy baryons to the ground state baryon states are studied in the framework of the light cone QCD sum rules method. Firstly, the transition form factors that describe these transitions are estimated, and then using these form factors the corresponding decay widths are calculated. A comparison of our results on the decay widths with those predicted by the other approaches existing in literature is performed.


Introduction
Last fifteen years is characterized (marked) by the impressive progress in hadron physics experiments, especially on the charmed hadron physics. Many new particles have been discovered [1][2][3][4][5]. Part of the new meson states could not be described by the conventional qq picture. Initiated by these observations remarkable amount of theoretical works appear which lead to the conclusion that these states could be interpreted as molecules [6][7][8][9][10], tetraquarks [11,12] or hybrids [13]. Similar situation exists in charmed baryon sector. Some of the newly discovered charmed baryon states [14][15][16] as the J P = 1 2 − Λ c (2595) or the J P = 3 2 − Λ c (2625) can be interpreted as the meson-baryon molecule. At present time many negative parity baryons are discovered in the charmed sector, while only Σ b baryon observed in the beauty sector. the study of the modes of the newly discovered negative parity heavy baryons can be useful for establishing their nature, i.e., the conventional qqQ structure, or molecular picture, or more exotic structure. The radiative decays of these hadrons into their ground states may constitute significant part of the total width, if the hadronic states are suppressed by the phase-volume or if the coupling constants are small. The radiative decays can also be useful in the determination of the quantum numbers of the negative parity states, as well as for understanding their internal structures.
The present work is devoted to the study of the radiative decays of the negative parity heavy baryons to the ground state heavy ones in framework of the light cone QCD sum rules method (LCSR). In our analysis we have used the background field formalism, that was introduced in [17] to calculate the Σ → pγ decay. Note that the radiative decays between negative (positive) parity heavy baryons is studied in the same framework in [18] ( [19]). It should also be emphasized that the decays of the heavy hadrons are studied in framework of the LCSR in the leading order of the heavy quark effective theory [20].
2 Transition form factors of the B i → B f γ decays at Q 2 = 0 In the present section we will examine the transition form factors of the B i → B f γ decays at Q 2 = 0 (for the real photon). Before giving the details of the calculations we first present the definition of the transition matrix element between the initial (B i ) and final (B f ) baryons states induced by the electromagnetic current which can be written as, where F 1 and F 2 are the Dirac and the Pauli form factors, and q = p 1 − p f . Using this matrix element for the B 1 → B 2 γ decay, we can easily calculate the decay width whose expression is as, where is the magnitude of the photon three-momentum in the rest frame of the initial baryon, and α is the fine structure constant. It follows from Eq. (2) that the radiative decay width is described only by the form factor F 2 (Q 2 = 0). The process B i → B f γ can be described by the correlation function where η 1 and η 2 are the interpolating currents of the initial and final heavy baryons, j el µ = e qq γ µ q + e QQ γ µ Q is the electromagnetic current, e q and e Q are the electric charges for the light and heavy quarks, respectively.
According to the SU(3) classification the heavy baryons belong to the symmetric sextet and antitriplet representations with respect to the light quark exchanges. Interpolating currents corresponding to these representations are constructed in [20] which are given as, where a, b, c are the color indices, and β is the arbitrary auxiliary parameter. The light quark contents of the heavy baryons in 6 and3 representations are given in Table 1. The decay constant of the baryon interpolating current is defined by, Introducing an electromagnetic background field of the plane wave F µν = i ε (λ) µ q ν e iqx , the correlation function given in Eq. (3)can be rewritten as, where the subscript F means that all calculations are performed in the background field F µν . The background field method was introduced in [21] in order to calculate the decay width of the of the Σ → pγ process. Let us now consider the transitions involving heavy hadrons. Along the lines as is suggested by the usual sum rules method, we insert a complete set of heavy baryons between the interpolating currents η 1 and η 2 , and the electromagnetic current in the correlation function given in Eq. (3). In doing so, there appears a problem which is absent in the case of mesons, namely the chosen form of the interpolating current couples not only with the positive parity ground state baryons J P = 1 2 + , but also with a heavier baryon resonance with the negative parity J P = 1 2 − . Many of the heavy baryons, especially the ones involving c-quark have already been observed in the experiments. Hereafter, the negative parity baryons and their masses will be denoted as B * and m B * , respectively.
The mass difference between the negative parity and ground state positive parity baryons is about 300 MeV , and similar mass difference is expected for the baryons with b-quark. Consequently, we should take into account the contribution of the negative parity baryons in the correlation function. Keeping these remarks in mind, from the hadronic side, the correlation function containing double poles can be written as, where summation over α and β is performed over the positive and the negative parity baryons, and dots indicate the contributions of higher states and continuum. The matrix elements in Eq. (6) are defined as, where λ 2 λ * 2 are the residues of the baryons B 2 and B * 2 , and F 1 and F 2 are the form factors for the corresponding transitions, respectively.
The matrix elements As we already noted that the conservation of the electromagnetic current leads to the result that the radiative decay is described only by the form factor F 2 .
Using the relation for the summation of the Dirac spinors over the spin, and using the transversality condition q ·ε = 0 for the photon field, one can easily show that only the (p·ε) term is proportional to F 2 .
Separating the terms proportional to F 2−+ , we get the following result for the correlation function from the phenomenological part, where As has already been noted, in this work we study the radiative decays of the transitions of the negative parity heavy baryons to positive heavy baryons. The invariant function A −+ in Eq. (9) describes these radiative decays under consideration. Therefore the other invariant functions should be eliminated. It follows from Eq. (9) that we have four independent invariant functions, and hence we need four equations in order to solve for the invariant function A −+ . The four equations can be obtained from the coefficients of the structures (p·ε)/ p/ q, (p·ε)/ p, (p·ε)/ q, and (p·ε)I.
In order to obtain the sum rules for the form factors F 2−+ (0) the expression of the correlation function from the QCD side is needed. This correlation function can be obtained can be calculated using the operator product expansion (OPE). The OPE is performed in the deep Eucledian region for the variables p 2 , (p + q) 2 ≪ m 2 Q . In LCSR method the OPE is performed over the twists of the nonlocal operators. The expansion of the nonlocal operators up to twist-4 is performed in [17], where the four-particle contributions are neglected. Following this work we also neglect them. The matrix elements of the nonlocal operators between the vacuum and the photon state are parametrized in terms of the photon distribution amplitudes (DAs). The photon DAs are studied comprehensively in [22], and therefore we do not present them in this work.
Equating the coefficients of the structures (p · ε)/ p/ q, (p · ε)/ p, (p · ε)/ q, (p · ε)I from the OPE side to the double hadronic dispersion relations, we obtain the following four linearly independent equations from which we cal solve for the form factor F 2−+ , Solving the set of four equations in Eq. (11), we can easily find A −+ . Performing the double Borel transformation over the variables (p + q) 2 → M 2 1 , p 2 → M 2 2 , we obtain the LCSR for the form factor F 2−+ at the point Q 2 = 0 whose expression is given as, where Π B i denote the Borel-transformed invariant functions in the coefficients of the aforementioned Lorentz structures. The expressions of the invariant functions Π B i are quite lengthy and therefore we do not present them in this work. The function ρ h (s 1 , s 2 ) in the last term is the hadronic spectral density of all excited and continuum states with the quantum numbers of B * . The hadronic spectral density is estimated by using the quark-hadron duality ansatz. The residues of the corresponding positive (negative) parity baryons are calculated in [23], [24] ( [25]). In our calculations for residues we have used the results of these works.
For the processes under consideration we choose M 2 1 = M 2 2 = 2M 2 due to the fact that the masses of the initial and final baryons are quite close to each other. The continuum subtraction can be carried out with the help of the transformation, For the terms that are proportional to the zeroth or negative power of M 2 or negative powers of M 2 the continuum subtraction is not performed because these contributions are negligibly small (see for details [26]).

Numerical analysis
In order to perform the numerical analysis for the form factor F 2−+ (0) we first introduce the input parameters which we shall use in further calculations. The main ingredient in any of the LCSR approaches in calculating the form factors is the set of photon DAs of the particle under consideration (in our case photon DAs). The photon DAs are obtained in [22], and for completeness we present their expressions in Appendix A. The mass of the negative and positive parity baryons are taken from [28]. The mass of the virtual c quark appearing in the heavy quark propagator is set to its value calculated in the MS scheme which is given as:m c (m c ) = 1.28 ± 0.03 GeV [25]. In our analysis for the values of the quark condensates we use qq (1 GeV ) = −(246 +28 −19 MeV ) 3 [29,30], ss (1 GeV ) = 0.8 qq (1 GeV ), m 2 0 = (0.8 ± 0.2) GeV 2 [31]. The value of the magnetic susceptibility we have used is χ(1 GeV ) = −0.285 GeV −2 [32].
The sum rule for the form factors F 2−+ (0), in addition to the aforementioned input parameters, contain also three auxiliary parameters. These parameters are the Borel mass parameter M 2 , the continuum threshold s 0 , and the arbitrary parameter β appearing in the expressions of the interpolating currents. According to the QCD sum rules methodology we need to find the so-called working regions of these parameters, where F 2−+ (0) be insensitive to the variation of these parameters in their working regions. Using the standard criteria of the sum rules which requires that both power corrections and continuum contributions should sufficiently be suppressed. The result of the analysis of the sum rules for the negative parity heavy baryons which has been studied in [30], leads to the following intervals for the Borel parameter M 2 : The range of the values of the continuum threshold s 0 is determined by imposing the condition that the sum rules should reproduce the measured values of the form factors to within 10-15% accuracy, from which we obtain 11 GeV 2 ≤ s 0 ≤ 14.0 GeV 2 . Our final attempt is focused on the determination of the working region of the arbitrary parameter β. The working domain of the parameter β is obtained by demanding that the form factor F 2−+ (0) shows good stability with respect to the variation of β. We have gone through analysis for all possible transitions, and the form factors F 2−+ (0) seems to be practically insensitive to the variation in cos θ in the domain −1 ≤ cos θ ≤ −0.5 and observed that this working region of cos θ is practically common for all transitions. Our final results for the form factor F 2−+ (0) for all transition channels are given below, The errors in the results are estimated by changing various input parameters within their working regions, and by also taking into consideration the resulting separate errors of the form factors F 2−+ (0) quadratically.
Having obtained the values of the form factors F 2−+ (0), we can calculate the decay widths of the decays under consideration. The values of the decay widths are presented in Table 2. For completeness we also present the values of the decay widths for the same transitions obtained in framework of the other approaches such as the coupled channel dynamically generated model [33], relativistic quark model [34], light cone sum rules method [35], chiral perturbation theory [36], bound state picture [37], and constituent quark model [38].   Reviewing the values of the decay widths given in Table 1, we see that between our results and predictions coming from other approaches on many channels, there are considerable differences. Therefore, measurement of decay widths can serve as a good tool for "choosing the right model". From Table 1, we conclude that, our approach predict quite large values, especially for the decay channels and these channels ( can be observed in the new future in experiments conducting at accelerators. Our final remark on this section is that, the decay widths for corresponding beauty baryons can be obtained from these calculations by performing relevant replacements

Conclusion
In this work we study the radiative decays of the negative parity heavy baryons to the ground state positive parity heavy baryons in framework of the LCSR method. We obtain that the transitions Σ ++ * c → Σ ++ c γ, Ξ ′0 * c → Ξ ′0 c γ, Σ + * c → Λ + c γ, Ξ ′+ * c → Ξ + c γ, Ξ + * c → Ξ + c γ, and Λ + * c → Λ + c γ have sizable decay width values. Our results show that the decay widths calculated in this work are quite different from the ones predicted by the other approaches. The apparently sizable values of the decay widths we obtain indicate that these radiative decays can serve a very useful tool for gaining information about the properties of the negative parity heavy baryons.