Sensitivity of single-particle energies to the spin-orbit coupling and to nuclear core structure in p-shell and sd-shell hypernuclei
We introduce a mean field model based on realistic 2-body baryon interactions and calculate spectra of a set of -shell and -shell hypernuclei - C, O, Ne, Si and Ca. The hypernuclear spectra are compared with the results of a relativistic mean field (RMF) model and available experimental data. The sensitivity of single-particle energies to the nuclear core structure is explored. Special attention is paid to the effect of spin-orbit interaction on the energy splitting of the single particle levels and . In particular, we analyze the contribution of the symmetric (SLS) and the anti-symmetric (ALS) spin-orbit terms to the energy splitting. We give qualitative predictions for the calculated hypernuclei.
keywords:hypernuclei, spin orbit splitting, interacion, mean field model
One of the major goals in hypernuclear physics is to obtain information on baryon-baryon interactions in a unified way. However, due to the limitation of hyperon () nucleon () scattering data, the potential models proposed so far, such as the Nijmegen models, have a large degree of ambiguity. Therefore, quantitative analyses of light hypernuclei, where the features of interactions appear rather clearly in observed level structures are indispensable. For this purpose, accurate measurements of -ray spectra Tamura00 (); Ajimura01 (); Akikawa02 (); Ukai04 (); Ukai06 () and high resolution reaction experiment Hotchi01 () have been performed systematically. These experiments are source of invaluable information about spin-dependent components of interaction. Useful constraints on the interaction components have been provided by shell model Millener85 () and few-body cluster Hiyama2009 (); Hiyama2012 () calculations. Among the spin-dependent components, spin-orbit terms are important, since they are intimately related to modeling the short-range part of the interaction. For example, it is well known that the antisymmetric spin-orbit (ALS) forces come out qualitatively different in one-boson-exchange (OBE) models model-D (); model-F (); soft-core (); Rijken99 () and in quark models Morimatsu84 (); QCM1 (). It was pointed out that the ALS force based on quark model QCM1 () is as strong as to cancel the symmetric spin-orbit (SLS) force, while the ALS force of Nijmegen potentials model-D (); model-F (); soft-core (); Rijken99 () is of a smaller strength. To extract information on these spin-orbit forces, specific -ray experiments Ajimura01 (); Akikawa02 () and reaction experiments Hotchi01 () have been performed. In the -ray experiments, spin-orbit dominated energy splittings of keV for the and doublet of levels in Be Akikawa02 () and keV for the doublet in C Ajimura01 () were measured. Considerably larger energy splittings of , and MeV were reported from the analysis of the -orbit, -orbit and -orbit peaks, respectively, observed in the Y spectrum of the reaction Hotchi01 (). However, these splittings are most likely caused by core-excited configurations and have little to do with the spin-orbit interaction Motoba2008 ().
The strength of the spin-orbit forces was studied within few-body Hiyama2000 () and shell model calculations Millener2012 () of Be and C. Reasonable strengths were introduced to reproduce the experimental data. It is worth mentioning that small energy splittings for Be and C were predicted by Hiyama et al. Hiyama2000 () using ALS forces based on a quark model. On the other hand, it was pointed out in the shell-model calculation by Millener Millener2001 () that a tensor contribution was as important as spin-orbit for these energy splittings.
More information about spin-orbit interaction in other hypernuclei is certainly needed. For this purpose, it is planned to explore the level structure of some medium heavy hypernuclei at JLab. Our aim in the present paper is to study the spin-orbit doublet states and in hypernuclei built on nuclear cores. It is expected that the spin-spin and tensor terms of the interaction are not effective as a consequence of the -cluster dominating structure, and we can thus safely assess spin-orbit forces. However, if clusters in these hypernuclei are broken, the spin-spin and the tensor interactions might contribute to the and energy splitting as well. In this paper, we investigate spin-spin and spin-orbit contributions to this energy splitting and give qualitative predictions of the splittings in C, O, Ne, Si and Ca. For this aim we apply the ESC08c potential, which has been recently proposed by Nijmegen group, using a mean field approach based on realistic baryon interactions HYPproc (). In the next section, we briefly describe our method. Selected results are presented and discussed in the third section, and a summary is given in the final section.
We calculate spectra of single hypernuclei using a Hamiltonian of the form
where is the kinetic term, () denotes interaction among nucleons ( and nucleons) and is a center of mass term. In our calculations we use a realistic interaction and a realistic interaction corrected by the density dependent term which simulates force. If we introduce the creation (annihilation) operators () for nucleons and () for , we can express the Hamiltonian (1) in the formalism of second quantization
with the kinetic matrix elements
antisymmetrized interaction matrix elements
and interaction matrix elements
all expressed in the harmonic oscillator basis. The harmonic oscillator basis depends on one parameter which defines the oscillator lengths and for the wave functions of nucleons and , respectively, due to the relations
In basis which is large enough the physical results do not depend on . In our calculations we use and MeV.
The mean field is constructed in our model as follows. First we solve the nuclear part of Hamiltonian (2) within the Hartree-Fock approximation. As a result we obtain the wave function as well as the nuclear density of the nuclear core of a hypernucleus. Then we calculate the single-particle energies by diagonalizing the matrix where
assuming that the hyperon interacts with the mean field of the core nucleus. The hypernuclear wave function at the level of mean-field approximation is defined as .
In our model, we used the realistic interaction NNLO n2lo (). It is a chiral next-to-next-to leading order potential with parameters optimized to minimize the effect of three-body interactions (their effect, however, still remains relatively important), which makes this force useful for many-body calculations (for more details about the optimization procedure see n2lo ()). Matrix elements of this interaction were generated by the CENS code Hjort (). However, the effect of three-body forces is still not negligible. If we perform the calculations purely with the two-body interactions we do not obtain correct distribution of the nuclear density. The nuclear density distribution is more compressed which leads to much smaller nuclear rms radii than are the experimental values. In general, this has influence on the single particle energies of , particularly on the splitting between the and states as can be deduced from the Bertlmann-Martin inequalities Bertl (). For this reason we add a corrective density dependent (DD) term of the form
where is the coupling constant and is the spin exchange operator. This density dependent interaction term was first introduced in Hergert (). It was shown Waroq () that it gives the same contribution to the Hartree-Fock energy as the contact three-body interaction
The term (8) is necessary for reasonable description of correct nuclear single particle spectra within the mean field calculations with realistic interactions Ves () and improves significantly the description of the nuclear density distributions and radii Ves2 (). In this paper we fix the values of for each hypernucleus independently to get the realistic values of nuclear radii as well as nuclear densities with respect to the available experimental data Radii () and the calculations within the RMF model Jarka (). Our future goal is to implement directly the chiral interaction instead of the density dependent term (8). In this case we would not need to fit any independent parameter .
For more details including the values of the parameters , , , see Isaka (). We represent the interaction (10) in the form of the interaction elements of Eq. (5). It should be noted that there are no tensor terms in the ESC08c version Isaka () used in this work.
The interaction depends explicitly on the Fermi momentum . We can either consider as a free parameter of our model and fit its value to the observed hypernuclear spectra or we can fix the value of within the Thomas-Fermi approximation through the relation
where is the density of the nuclear core and is the density of in the hypernucleus. Note that we have to perform the hypernuclear calculation to obtain and determine the value of Fermi momentum . In other words the value of has to be evaluated self-consistently because the equation (12) depends on the result of a calculation which itself depends on .
The symmetric (SLS) and antisymmetric (ALS) spin-orbit terms in the potential are included within the Scheerbaum approximation Scheer (). Due to this approximation we include the effect of the SLS and ALS terms directly at the mean field level. We add the following contribution into the matrix (7):
The value of in (15) is set to 0.7 fm Scheer (). The form of the function in (15) is identical for the SLS and ALS terms, they only differ by the values of the input parameters. It should be noted that the - coupling term in ESC08c is renormalized into the - part of G-matrix interaction, giving rise to an important part of the density dependence Isaka ().
In the RMF approach, the strong interactions among point-like hadrons are mediated by effective mesonic degrees of freedom. The formalism is based on the Lagrangian density of the form
Here, is the standard nuclear Lagrangian Jarka () and we used the NL-SH parametrization in this work NLSH (). The is the anomalous (tensor) coupling term. This term is crucial in order to get negligible spin-orbit splitting for larger values of the couplings required by constituent quark model jen (); Mares ().
For the coupling constants and we used the naive quark model values and was tuned so as to reproduce the binding energy of in the state in O Mares ().
In this section, we present selected results of our calculations of the hypernuclei C, O, Ne, Si and Ca. In the hypernuclei with the doubly magic nuclear core - O and Ca - the calculations within Hartree-Fock method are straightforward. However, in the case of C, Ne and Si, the ground states of the corresponding core nuclei have more complex structure Soyeur (); Miller () and it is necessary to take into account configuration mixing and perform calculations within a deformed basis in order to describe their structure properly. Nevertheless, even calculations within the spherical HO basis could provide interesting information about these hypernuclei if we consider various configurations of the corresponding nuclear cores. In this work, we performed calculations for the following configurations: and in C; , , and in Ne; , and in Si. We treated these configurations always symmetrically for both protons and neutrons. It is to be noted that more configurations could be realized in the ground states of the above nuclei. We selected just some of them in order to illustrate the effect of the wave function of the nuclear core on the energy splitting of the single particle levels and in the considered hypernuclei.
For each particular nucleus we first fixed the parameter to obtain reasonable density distribution and rms radius of the nuclear core, comparable with the available experimental data and RMF calculations within the NL-SH parametrization NLSH (). The values of the charge rms radii are summarized in Table 1.
We used the values MeVfm for C, MeVfm for O, MeVfm for Ne, MeVfm for Si, and MeVfm for Ca. In case of hypernuclei with the open-shell core we did not repeat tunning of the parameter for each configuration separately but we fixed it for one case (the configuration for C, for Ne and for Si). The corresponding radii for the remaining configurations differ from the values shown in Table 1 but the differences are much smaller than the differences between the values in the columns A and B for each given nucleus. We chose to fit the radii and radial distributions for the above configurations because they are the lowest configurations in energy due to the empirical ordering of the and levels.
In Fig. 1, we present the radial nuclear density distributions calculated within the mean field model based on realistic 2-body baryon interactions and the RMF model NL-SH NLSH (). The figure illustrates the importance of the DD term in the interaction (8) on selected nuclei – C (in the configuration), O, Ne (in the configuration), Si ( in the configuration), and Ca. Calculations performed without the DD term () yield unrealistically large central densities and, as a consequence, the corresponding rms radii are too small. After including the DD term the density distributions become more diffused and get closer to the RMF predictions which are in agreement with empirical density distributions Jarka (); neon () (compare also charge rms radii in Table 1).
The values of the Fermi momentum used in the present calculations were determined using the ADA approximation (Eqs. (11) and (12)). We obtained fm for C (in the configuration ), fm for O, fm for Ne (in the configuration ), fm for Si (in the configuration ), and fm for Ca.
Before focusing on the main objective of the present work, the energy splittings between and , we discuss on the single particle energies in the considered hypernuclei. Fig. 2 shows our hypernuclear spectra, calculated for the configurations for which the parameters and were tuned, and the spectra calculated within the RMF model NL-SH for the same configurations. We see that our results are in good agreement with the experimental data shown for comparison. However, it should be pointed out that the data in the figure are for O, Si, and Ca since the data for hypernuclei with the same and closed nuclear cores are not available at present (unlike the C case).
The most pronounced difference between the HF calculations based on realistic interactions and the RMF model NL-SH occurs in the spectrum of Si. Here the RMF model predicts significantly more binding for - the lowest level has nearly the same energy as the level in Ca. It will be demonstrated in Fig. 4 that the discrepancy between the calculated Si spectra is smaller for other nuclear core configurations and could be attributed to considerably different nuclear spin-orbit splittings in the two considered models. We expect that proper calculations allowing for nuclear core deformation and configuration mixing of the nuclear core wave function could decrease this discrepancy in predicted spectra in Si. It is to be noted that the energies calculated in both considered models could get closer to each other if we fine tuned e.g. the parameter for each particular configuration of the nuclear core in the mean field model based on realistic interactions and/or the RMF coupling separately for Si. However, being aware of the limitations of our current hypernuclear calculations we do not intend to do so in the present study.
Since our single particle energies reproduce the data reasonably well, let us discuss on the main subject of our study, the energy splitting between the single particle levels and . In Table 2, we present the calculated values of for the following options: the spin-orbit (SLS+ALS) forces completely switched off (A), only the SLS term included (B), both the SLS and ALS terms included (C). The results obtained within the RMF model are presented for comparison. The configurations for which the parameters and were tuned are given in bold face.
It is to be stressed that the energy splittings calculated within the mean field based on realistic baryon interactions and RMF models are of different origin. In the former case, the splitting is caused by the 2-body interaction with its various spin dependent terms. On the other hand, the splitting of the spin orbit partners in the RMF model results from a delicate balance between strong scalar and vector mean fields in the Dirac equation jen (). This explains qualitatively different predictions for within the two approaches. Clearly the calculations based on realistic baryon interactions give considerably larger variations of the values in the studied hypernuclei than the RMF approach.
The ESC08c potential gives a non-zero, negative energy splitting of the and levels in several configurations in C, Ne and Si even if the SLS and ALS forces are switched off (see column A in Table 2). The spin-spin term in the ESC08c interaction thus contributes significantly to the splitting in these configurations of hypernuclear cores. In the case of O and Ca (hypernuclei with the doubly magic nuclear core) and also in C for the configuration , in Ne for the configuration and in Si for the configuration the and levels are close to be degenerate. We checked that these states are roughly degenerate also in the other two hypernuclei with doubly magic core - Zr and Pb. When the SLS term is included (B), the splitting becomes positive in all considered hypernuclei except Si in the core configuration. The ALS term acts in an opposite way to the SLS term (thus weakens the effect of SLS) and its magnitude is of the SLS magnitude. The combined (SLS + ALS) term (column C) remains strong enough to cause positive splitting in most of the cases except Si in the core configurations and , and C in the core configuration . It is also to be pointed out that in the cases when the states and are close to be degenerate in the column A, the energy splitting comes almost entirely from the SLS and ALS spin-orbit interaction terms.
Finally, there is a missing RMF value for the Ne configuration in Table 2 – this configuration could not be calculated since the proton level was found to be unbound in the applied RMF parametrization.
The competition between the SLS and ALS forces can be illustrated with the help of the function defined as follows:
evaluated for (). The radial integral of determines the contribution of the SLS (ALS) force to the energy splitting between the and levels (see Eq. 13).
In Fig. 3, we compare for with for in O (left panel), Ne in the configuration (middle panel) and Si in the configuration (right panel). The difference between the areas delimited by SLS and -ALS curves above and below zero determines the spin-orbit (SLS + ALS) contribution to the energy splitting . This difference is relatively large in O while in the case of Ne, the ”negative” and ”positive” contributions compensate more each other. Even larger compensation effect is seen for the Si calculated within the configuration . Consequently, the (SLS + ALS) splitting in Si calculated in this configuration is about twice smaller than in O (compare Table 2).
In Fig. 4, we show the single particle energies in Si with three different configurations of the nuclear core, calculated within the mean-field model with realistic interactions (black lines) and the RMF model NL-SH (red lines). The levels are denoted by dotted lines. The experimental values for Si Pile (); Haseg () are shown for comparison (there are no data for Si). The figure demonstrates how the single particle spectrum is affected by the wave function of the nuclear core. We can see that the levels and switch their ordering in the case of mean-field based on realistic baryon forces. This is in contrast to the RMF model for which the energy splitting of both levels remains roughly constant. We found possible explanation for this effect by analyzing the nucleon single particle energies. While the RMF model is able to reproduce well the empirical spin-orbit splitting of the and states for protons and neutrons, the mean field model based on realistic baryon forces appears to be quite sensitive to different configurations considered for the ground state. We do not obtain realistic splitting of the nucleonic and states in the configurations and – in the former case we even get wrong ordering of these levels. Only for the configuration we obtain satisfactory agreement of the nucleonic spin-orbit splitting with the empirical values (and also the RMF values). Consequently, we also get standard ordering of the and levels (see Table 2).
In this work, we performed calculations of selected - and - shell hypernuclei, namely C, O, Ne, Si and Ca within the mean field model based on realistic 2-body baryon interactions and compared the results with the predictions of the RMF model NL-SH. We introduced the density dependent 2-body interaction term which mimics the effect of the 3-body force in order to get realistic charge radii and density distributions of the nuclear cores of the studied hypernuclei. This appeared important also in the calculations of hypernuclear spectra since the ESC08c interaction depends explicitly on the Fermi momentum which was determined using the averaged density approximation. Reasonable description of the density distributions in the studied (hyper)nuclear systems is thus crucial.
The main objective of the present calculations is to study the influence of SLS and ALS spin orbit terms on the energy splitting of the levels and . The splittings in O and Ca, calculated within the mean field model based on realistic baryon interactions and the RMF model NL-SH are very close to each other. In the case of C, Ne and Si it is desirable to perform calculations within deformed basis and take into account configuration mixing of the nuclear core wave function. Nevertheless, our calculations in the spherical HO basis, which considered several configurations of the nuclear core of these hypernuclei yielded valuable insight into the issue of the spin dependence of the interaction and the spin-orbit splitting in these open-shell hypernuclei. We found that the energy splittings of the levels and calculated using realistic and interactions depend strongly on the chosen configuration of the nuclear core, unlike the RMF approach. For the configurations which give the energy splitting close to zero when the spin-orbit interaction is switched off, the splitting is almost entirely due to the ALS and SLS terms and is in rough agreement with the RMF values. By comparing the results for the ESC08c model shown in columns B and C in Table 2 we conclude that the magnitude of the ALS term which acts in an opposite way to the SLS term is about 2/3 of the SLS magnitude.
Our results demonstrate that it is highly desirable to explore further the energy splitting of the and levels in - and -shell hypernuclei, both experimentally and theoretically, in order to extract important information about the spin-dependence of the interaction, as well as the inner structure of the hypernuclei under study.
There are several issues left for further improvements of the present calculations. First, we intend to develop the code which will allow to perform calculations within an axially symmetric single particle basis and allow for configuration mixing in the nuclear core wave function. In this case we would be able to calculate open-shell hypernuclei (such as Ne or Si) more precisely. Second, we intend to implement directly the 3-body forces instead of the 2-body density dependent term in our Hamiltonian. Third, it is desirable to incorporate the tensor terms and explore their contribution to the energy splitting.
Another extension is to include the core polarization effects. We intend to develop a scheme which couples single-particle states with one-phonon or possibly multi-phonon excitations of the core nucleus within an Equation of Motion Phonon Model (EMPM) EMPM () treating nuclear excitations in multiphonon basis. In this method the Tamm-Dancoff phonon operators are used to build Hilbert space spanned by one-, two- and three-phonon configurations. In this way we will get rather complex description of hypernuclei which includes not only core polarization effects but also beyond mean field correlations.
This work was partly supported by the GACR grant P203/15/04301S. Highly appreciated was the access to computing and storage facilities provided by the MetaCentrum under the program LM2010005 and the CERIT-SC under the program Center CERIT Scientific Cloud, part of the Operational Program Research and Development for Innovations, Reg. no. CZ.1.05/3.2.00/08.0144. P. V. thanks RIKEN for the kind hospitality during his stay. This work was partly supported by RIKEN iTHES Project.
- (1) H. Tamura et al., Phys. Rev. Lett. 84 (2000) 5963.
- (2) S. Ajimura et al., Phys. Rev. Lett. 86 (2001) 4255.
- (3) H. Akikawa et al., Phys. Rev. Lett. 88 (2002) 082501; H. Tamura et al., Nucl. Phys. A 754 (2005) 58c.
- (4) M. Ukai et al., Phys. Rev. Lett. 93 (2004) 232501.
- (5) M. Ukai et al., Phys. Rev. C 73 (2006) 012501(R).
- (6) H. Hotchi et al., Phys. Rev. C 64 (2001) 044302.
- (7) D. J. Millener, A. Gal, C. B. Dover, and R. H. Dalitz, Phys. Rev. C31 (1985) 499.
- (8) E. Hiyama and T. Yamada, Prog. Part. Nucl. Phys. 63 (2009) 339.
- (9) E. Hiyama, Few-Body. Syst. 53 (2012) 189.
- (10) M. M. Nagels, T. A. Rijken, and J. J . de Swart, Phys. Rev. D 12 (1975) 744.
- (11) M. M. Nagels, T. A. Rijken, and J. J. de Swart, Phys. Rev. D 20 (1979) 1633.
- (12) M. M. Nagels, T. A. Rijken and J. J. de Swart, Phys. Rev. D 17 (1978) 768; P. M. M. Maessen, T. A. Rijken and J. J. de Swart, Phys. Rev. C 40 (1989) 226.
- (13) Th. A. Rijken, V. G. J. Stoks and Y. Yamamoto, Phys. Rev. C 59 (1999) 21; V. G. J. Stoks and Th. A. Rijken, Phys. Rev. C 59 (1999) 3009.
- (14) O. Morimatsu, S. Ohta, K. Shimizu, and K. Yazaki, Nucl. Phys. A 420 (1984) 573.
- (15) Y. Fujiwara, C. Nakamoto and Y. Suzuki, Prog. Theor. Phys. 94 (1995) 65; ibid. 94 (1995) 215; ibid. 94 (1995) 353; Phys. Rev. Lett. 76 (1996) 2242; Phys. Rev. C 54 (1996) 2180.
- (16) T. Motoba, D. E. Lanskoy, D. J. Millener and Y. Yamamoto, Nucl. Phys. A 804, (2008) 99.
- (17) E. Hiyama, M. Kamimura, T. Motoba, T. Yamada and Y. Yamamoto, Phys. Rev. Lett. 85 (2000) 270.
- (18) D. J. Millener, Nucl. Phys. A 881 (2012) 298.
- (19) D. J. Millener, Nucl. Phys. A 691 (2001) 93c.
- (20) P. Veselý, J. Mareš, E. Hiyama, to appear in proc. HYP2015, Sendai, Sept. 7 -12, 2015.
- (21) A. Ekström, G. Baardsen, C. Forssén, G. Hagen, M. Hjorth-Jensen, G. R. Jansen, R. Machleidt, W. Nazarewicz, T. Papenbrock, J. Sarich, S. M. Wild, Phys. Rev. Lett. 110 (2013) 192502.
T. Engeland, M. Hjorth-Jensen, G. R. Jansen,
CENS, a Computational Environment for Nuclear Structure, in preparation,
- (23) R. J. Lombard, S. Marcos, J. Mareš, Phys. Rev. C 50 (1994) 2900.
- (24) H. Hergert, P. Papakonstantinou, R. Roth, Phys. Rev. C 83 (2011) 064317.
- (25) M. Waroquier, K. Heyde, H. Vincx, Phys. Rev. C 13 (1976) 1664.
- (26) D. Bianco, F. Knapp, N. Lo Iudice, P. Veselý, F. Andreozzi, G. De Gregorio, A. Porrino, J. Phys. G: Nucl. Part. Phys. 41 (2014) 025109.
- (27) F. Knapp, N. Lo Iudice, P. Veselý, F. Andreozzi, G. De Gregorio, A. Porrino, Phys. Rev. C 92 (2015) 054315.
- (28) I. Angeli, Atomic Data and Nuclear Data Tables 87 (2004) 185.
- (29) B. D. Serot, J. D. Walecka, Adv. Nucl. Phys. 16 (1986) 1.
- (30) T. Rijken, M. Nagels, Y. Yamamoto, Prog. Theor. Phys. Suppl. 185 (2010) 14.
- (31) M. Isaka, K. Fukukawa, M. Kimura, E. Hiyama, H. Sagawa, Y. Yamamoto, Phys. Rev. C 89 (2014) 024310.
- (32) Y. Yamamoto, T. Motoba, T. Rijken, Prog. Theor. Phys. Suppl. 185, (2010) 72.
- (33) R. R. Scheerbaum, Nucl. Phys. A 257 (1976) 77.
- (34) M. M. Sharma, M. A. Nagarajan, P. Ring, Phys. Lett. B 312 (1993) 377.
- (35) B. K. Jennings, Phys. Lett. B 246 (1990) 325.
- (36) J. Mareš, B. K. Jennings, Phys. Rev. C 49 (1994) 2472.
- (37) M. Soyeur, A. Zuker, Journal de Physique Colloques 32 (1971) C6.
- (38) H. G. Miller, H. P. Schroder, Z. Phys. A - Atoms and Nuclei 304 (1982) 273
- (39) Y. K. Gambir, P. Ring, A. Thimet, Ann. Phys. 198 (1990) 132.
- (40) M. Agnello et al., Phys. Lett. B 698 (2011) 219.
- (41) M. May, et. al., Phys. Rev. Lett. 78 (1997) 4343.
- (42) P. H. Pile, et al., Phys. Rev. Lett. 66 (1991) 2585.
- (43) T. Hasegawa, et al., Phys. Rev. C 53 (1996) 1210.
- (44) R. E. Chrien, et al., Nucl. Phys. A 478 (1988) 705c.
- (45) D. Bianco, F. Knapp, N. Lo Iudice, F. Andreozzi, A. Porrino, Phys. Rev. C 85 (2012) 014313.