application/xmlStructure of 33Mg sheds new light on the [formula omitted] island of inversionR. KanungoC. NociforoA. ProchazkaY. UtsunoT. AumannD. BoutinD. Cortina-GilB. DavidsM. DiakakiF. FarinonH. GeisselR. GernhäuserJ. GerlR. JanikB. JonsonB. KindlerR. KnöbelR. KrückenM. LantzH. LenskeY. LitvinovK. MahataP. MaierbeckA. MusumarraT. NilssonT. OtsukaC. PerroC. ScheidenbergerB. SitarP. StrmenB. SunI. SzarkaI. TanihataH. WeickM. WinklerNucleon removal or knockout reactionsShell structureMagic numberRadioactive beamsEikonal modelPhysics Letters B 685 (2010) 253-257. doi:10.1016/j.physletb.2010.02.008journalPhysics Letters BCopyright © 2010 Elsevier B.V. All rights reserved.Elsevier B.V.0370-26936854-58 March 20102010-03-08253-25725325710.1016/j.physletb.2010.02.008http://dx.doi.org/10.1016/j.physletb.2010.02.008doi:10.1016/j.physletb.2010.02.008http://vtw.elsevier.com/data/voc/oa/OpenAccessStatus#Full2014-01-01T00:14:32ZSCOAP3 - Sponsoring Consortium for Open Access Publishing in Particle Physicshttp://vtw.elsevier.com/data/voc/oa/SponsorType#FundingBodyhttp://creativecommons.org/licenses/by/3.0/JournalsS300.3PLB26523S0370-2693(10)00169-310.1016/j.physletb.2010.02.008Elsevier B.V.ExperimentsFig. 1Schematic setup of the experiment for one-neutron removal from 33Mg using a carbon reaction target at the central focal plane of the FRS.Fig. 2The longitudinal momentum distribution data for 33Mg→32Mg with carbon reaction target (open squares) and without the carbon reaction target (open circles).Fig. 3The longitudinal momentum distribution data for 33Mg→32Mg (filled circles). The curves are eikonal model calculations (a) with the core 32Mg in its ground state coupled to a neutron in the 1f7/2 orbital (dotted line), 1d3/2 orbital (solid line) and 2p3/2 orbital (dashed line), respectively, (b) dashed line (solid line) with the MCSM (MCSMmod) configurations (Table 1). The MCSMmod has the 2p3/2 orbital lowered by 1 MeV.Fig. 4The longitudinal momentum distribution data for 33Mg→32Mg (filled circles). The solid line is the eikonal model calculation from the configuration mixing (see text) with a 3/2− ground state for 33Mg. The individual significant configurations obtained from best fit χ2 search are labeled as (1) Mg21+32+ϕn(2p3/2), (2) Mg1−32+ϕn(2s1/2), (3) Mg21+32+ϕn(1f7/2), (4) Mg1−32+ϕn(1d3/2) and (5) Mggs32+ϕn(2p3/2).Fig. 5The longitudinal momentum distribution data for 33Mg→32Mg (filled circles). The solid line is the eikonal model calculation from the configuration mixing (see text) with a 3/2+ ground state for 33Mg. The individual significant configurations obtained from best fit χ2 search are labeled as (1) Mg3−32+ϕn(2p3/2), (2) Mg22+32+ϕn(2s1/2), (3) Mg3−32+ϕn(1f7/2).Table 1Shell model spectroscopic factors for 33Mg→32Mg using the SDPF-M interaction (MCSM) and modified with lower 2p3/2 orbital (MCSMmod). Jc is the spin of the 32Mg core state and Energy is its excitation energy. Configurations with S⩾0.1 are shown.JcEnergy (MeV)Nljσ−nsp (mb)S(MCSM)S(MCSMmod)S(3/2−)fitS(3/2+)fit0+0.002p3/240.120.120.6−0.5+0.30+0.001d3/222.390.0−0.0+0.521+0.892p3/234.130.140.230.5−0.3+0.721+0.891f7/220.870.320.240.5−0.5+0.241+2.321f7/217.931.130.883−2.862p3/224.602.2−0.5+0.23−2.861f7/217.231.1−0.5+0.122+3.012p3/225.820.190.2022+3.012s1/224.810.1−0.1+0.011−3.812s1/220.940.220.210.1−0.1+0.011−3.811d3/214.040.150.141.5−1.3+1.321−3.761d3/213.300.610.5722−3.942s1/220.910.290.273−3.991d3/213.610.440.393+4.002p3/221.960.1242+4.391f7/215.940.230.18Structure of 33Mg sheds new light on the N=20 island of inversionR.Kanungoa⁎ritu@triumf.caC.NociforobA.ProchazkabcY.UtsunodT.AumannbD.BoutincD.Cortina-GileB.DavidsfM.DiakakigF.FarinonbcH.GeisselbR.GernhäuserhJ.GerlbR.JanikiB.JonsonjB.KindlerbR.KnöbelbcR.KrückenhM.LantzjH.LenskecY.LitvinovbkK.MahatabP.MaierbeckhA.MusumarralmT.NilssonjT.OtsukanC.PerroaC.ScheidenbergerbB.SitariP.StrmeniB.SunbI.SzarkaiI.TanihataoH.WeickbM.WinklerbaAstronomy and Physics Department, Saint Mary's University, Halifax, NS B3H 3C3, CanadabGSI Helmholtzzentrum für Schwerionenforschung, D-64291 Darmstadt, GermanycPhysics Department, Justus-Liebig University, 35392 Giessen, GermanydJapan Atomic Energy Agency, Tokai, Ibaraki 319-1195, JapaneUniversidad de Santiago de Compostela, E-15706 Santiago de Compostela, SpainfTRIUMF, Vancouver, BC V6T 2A3, CanadagNational Technical University, Athens, GreecehPhysik Department E12, Technische Universität München, D-85748 Garching, GermanyiFaculty of Mathematics and Physics, Comenius University, 84215 Bratislava, SlovakiajFundamental Physics, Chalmers University of Technology, SE 412-916 Göteborg, SwedenkMax-Planck-Institut für Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, GermanylUniversita' di Catania, 95153 Catania, ItalymINFN-Laboratori Nazionali del Sud, 95123 Catania, ItalynCenter for Nuclear Study, University of Tokyo, Wako Shi, Saitama 351-0198, JapanoRCNP, Osaka University, Mihogaoka, Ibaraki, Osaka 567 0047, Japan⁎Corresponding author.Editor: V. MetagAbstractThe first reaction spectroscopy on the ground state structure of 33Mg through the measurement of the longitudinal momentum distribution from the one-neutron removal reaction using a C target at 898A MeV is reported. The experiment was performed at the FRS, GSI. The distribution has a relatively narrow width (150±3 MeV/c (FWHM)) and the one-neutron removal cross-section is 74±4 mb. An increased contribution from the 2p3/2 orbital is required to explain the observation showing its lowering compared to existing model predictions. This provides new information regarding the configuration of 33Mg and the island of inversion.KeywordsNucleon removal or knockout reactionsShell structureMagic numberRadioactive beamsEikonal modelExploring the changes in shell structure of neutron-rich nuclei has been one of the major foci of nuclear physics in recent times. The so-called island of inversion [1] around neutron-rich N=20 nuclei continues to be a region of great interest. The increased separation energy of the neutron-rich Na isotopes was the first signature for the breakdown of the conventional N=20 shell gap [2]. The observations of a low-lying first 2+ excited state in 32Mg [3–5] indicated the presence of deformation. The lowering of the first 2+ excited state is found to continue progressively for 34Mg [6] and 36Mg [7]. Large BE(2) values were also observed for lighter N=20 isotones 31Na [8] and 30Ne [9]. The low-lying (2+) excited state of 32Ne [10] was a signature of extension of the island of inversion to N=22 for the Ne isotopes.The cause for the quenching of the N=20 shell gap has been theoretically studied in the shell model framework [11–14] where fp shell contributions are taken into account. It has been discussed that in this region due to the reduced strength of the T=0 attractive monopole interaction, the 1d3/2 and 1f7/2 orbitals come closer to each other leading to deformation through np–nh neutron excitations [12,13]. The deformation increases the binding of the neutron-rich N=20 isotones.For a complete understanding on the evolution of the island of inversion it is crucial to understand how the nuclear orbitals across the broken N=20 gap are arranged. We report here the first investigation of the structure of 33Mg through the longitudinal momentum distribution from one-neutron removal from 33Mg. The shape of the distribution is much narrower than predictions from the shell model with the SDPF-M interaction, suggesting a larger occupancy of neutrons in the 2p3/2 orbital. This signals a lowering of the 2p3/2 orbital. The increased p-wave occupancy has not been discussed in the recent works that debate the ground state spin of 33Mg [15–17]. The present finding therefore opens a new dimension in our attempt to understand the island of inversion. A recent calculation using a Wood–Saxon potential for neutron-rich Mg isotopes discusses the near-degeneracy of the 1f7/2 and 2p3/2 orbitals in a spherical potential to be a basic element in producing deformation [18].In-beam and beta-delayed gamma spectroscopy experiments have tried to assess the boundary of the island of inversion. Among the excited states observed for 25–29Ne, those for 28Ne have been interpreted to show intruder signatures, while 25–27Ne were in agreement with normal sd-shell predictions [19]. On the other hand proton inelastic scattering has identified low-lying intruder states in 27Ne [20] that are in agreement with a study of the 26Ne(d,p) reaction [21]. The level scheme of 29,30Na constructed on the basis of the beta decay of 29,30Ne reflected the dominance of intruder configurations in the low-lying states [22,23]. These experiments therefore, provide signatures of intruder configurations in the excited states around N=20.It is worthwhile to point out here that the island of inversion refers to intruder orbitals in the ground states of nuclei [1,11,13]. Intruder configurations were reported for the ground state of 28Ne [26]. The 3/2− intruder ‘core’ state was found to have a spectroscopic factor (SF) of 0.32(4), in fair agreement with SDPF-M shell model predictions. The momentum distribution for this intruder component was explained by a combination of l=0 and 1 contributions. A measurement of the neutron knockout from 32Mg found that the ground state consists dominantly of 2p–2h N=20 cross-shell configurations [27]. The SF's for the intruder 2p3/2 and the 1f7/2 orbitals were found to be 0.59(11) and 1.19(36) respectively. The abnormal spin of 1/2+ observed for 31Mg showed a deformed ground state with an intruder configuration [28].On the neutron-rich side of the N=20 gap, the beta decay of 33Mg suggested an abnormal spin of 3/2+ due to a 1p–1h configuration resulting from an inversion of the d3/2 and f7/2 orbitals [15]. The intermediate energy Coulomb excitation of 33Mg populated a state at 485 keV that was attributed to a rotational excitation and suggested a 5/2+ spin for this state [29]. Such an excitation would support 33Mg having similar deformation as observed for other nuclei in the island of inversion. The measurement of a negative magnetic moment for the 33Mg ground state has however suggested a spin of 3/2− with a 2p–2h configuration [16]. In this configuration, it was noted that there are an odd number of neutrons in the fp shell. Very recently the large branching ratio observed for the beta decay of 33Mg to the 5/2+ ground state of 33Al suggests the 33Mg ground state to be 3/2+[17]. It was proposed to have an admixture of 1p–1h with two neutrons in the f7/2 orbital and 3p–3h with four neutrons in the f7/2 orbital configurations for 33Mg. In the proposed configurations (Fig. 4 [17]) the odd neutron occupies the d3/2 orbital and no occupancy in the 2p3/2 orbital is discussed. In this work the ground state configuration of 33Mg is investigated by measuring the momentum distribution after removal of one neutron.The experiment was performed at the fragment separator, FRS, at GSI, Darmstadt, Germany. The schematic view of the experiment is shown in Fig. 1 where the first three foci F1–F3 of the FRS are dispersive while the last F4 is achromatic. The 33Mg secondary beam was produced via the fragmentation of a 1A GeV48Ca primary beam on a 6.347 g/cm2 thick Be target. The nuclei produced were then separated and identified by transporting them to the second dispersive focus of the FRS. Detectors placed at F1 and before the target at F2 served to identify the secondary beam on target. The plastic scintillators at F1 and F2 measure the time-of-flight (TOF) between these locations. The position on the reaction target was determined by tracking using the time projection chambers (TPC) at F2. The momentum of the incident particles can be determined from the position and the central magnetic rigidity of the first half of the FRS. In addition, it can be derived from the TOF measurement. The energy-loss (ΔE) measured using the multi-sampling ionization chamber (MUSIC) provided the charge (Z) of the nuclei. The combined information from Bρ–ΔE-TOF allowed a complete event-by-event identification of the 33Mg secondary beam with less than 1% contamination.A 4.05 g/cm2 carbon reaction target was placed at the F2 dispersive mid-plane. The 32Mg fragments from the one-neutron removal reaction were transported to F4 using the second half of the FRS. Here they were identified using the same principle as above with the TPC, MUSIC and scintillator detectors placed at F4.To determine the momentum distribution due to one-neutron removal at the target, the FRS was operated in dispersion-matched mode. In this mode, the position distribution at F4 is a measure of the momentum distribution due to the reaction, independent of the initial momentum spread of the 33Mg secondary beam. The energy loss in and the shape of the reaction target at F2 leads to some deviation from this ideal achromatic image condition as observed through a correlation of the horizontal position at the reaction target at F2 and the horizontal position at the F4 focus. This correlation could then be accounted for to improve the achromatic condition by making the position at F4 independent of that at F2.Reactions occurring in the residual matter (detectors, vacuum windows) were separately measured in a setting where the carbon reaction target was removed. The magnetic rigidity of the second half of FRS was set to transport 32Mg as the central fragment at F4. To subtract these background contributions, the number of 32Mg fragments detected without the carbon target were normalized by the ratio of the incident beam (33Mg) particles measured with and without the target. Fig. 2 shows the normalized background distribution (open circles). This distribution was then subtracted from the distribution with the carbon target (Fig. 2, open squares). The background subtracted momentum distribution is shown in Figs. 3–5. The momentum acceptance of the FRS is shown by the vertical dashed lines in Fig. 2.The momentum resolution was determined in a separate setting, where the second half of the FRS was set to transport the unreacted 33Mg beam after the reaction target to F4. The observed momentum distribution fitted by a Gaussian function gives a measure of the resolution as 29.0±0.2 MeV/c (FWHM). The width of the background subtracted distribution for 33Mg→32Mg, shown by the black data points in Fig. 3, was 150±3 MeV/c (FWHM). Accounting for the momentum resolution, the intrinsic width (assuming gaussian shape) is 147±2 MeV/c (FWHM). The one-neutron removal cross-section is 74±4 mb. The error bars shown in Fig. 3 include both statistical and systematic errors. The systematic errors arise from the 1.2% uncertainty in the measured target thickness and the estimated 5% uncertainty in transmission.The data are interpreted in the framework of eikonal model calculations [30] based on a Mg32+n description of 33Mg. The distributions calculated in this model have been folded with the experimental momentum resolution. The wavefunctions for the neutrons were obtained from a Wood–Saxon binding potential, the depth of which was adjusted to reproduce the effective neutron separation energy, including the excitation energy of the 32Mg core for each configuration.Following the conventional shell model with a shell closure at N=20, the last neutron in 33Mg should occupy the 1f7/2 orbital. The dotted line in Fig. 3(a) shows the distribution for the Mggs32+n(1f7/2) configuration with the spectroscopic factor S=1. It can be seen that this configuration does not describe the data either in shape or in magnitude thereby confirming that neutron(s) are distributed differently between the sd-shell and the fp-shell. This rules out a 7/2− spin assignment in a 0p–0h description of the ground state of 33Mg.The beta decay studies suggest the ground state of 33Mg to be 3/2+ with an odd number of neutron(s) in the 1d3/2 orbital either in a 1p–1h or a 3p–3h configuration. The solid line in Fig. 3(a) shows the calculated momentum distribution for Mggs32+n(1d3/2) with a spectroscopic factor of 1.0. This distribution fails as well in reproducing the shape and magnitude of the data.The relatively narrow width of the momentum distribution suggests neutron removal from orbitals with low angular momentum. We consider neutron removal from the 2p3/2 orbital, that may be lowered both due to deformation as well as the decreasing neutron separation energy for neutron-rich nuclei. The momentum distribution for the configuration with the 32Mg core in its ground state and one neutron in the 2p3/2 orbital resulting in a ground state spin of 3/2−, is shown by the dashed line in Fig. 3(a) with a spectroscopic factor of unity. Its width of 108 MeV/c (FWHM) is narrower than the data, suggesting that the core 32Mg should also be in an excited state. This would increase the effective separation energy, thereby making the momentum distribution wider.The Monte Carlo shell model (MCSM) calculation with the SDPF-M interaction, that has been considered to be very successful in describing the island of inversion, predicts the ground state of 33Mg to be 3/2−. A highly mixed configuration is predicted with the configurations and spectroscopic factors listed in Table 1. Fig. 3(b). shows that the calculated momentum distribution (dashed line) with the MCSM configurations (Table 1) does not explain the data. The experiment clearly shows the necessity for an enhanced low angular momentum component in the ground state of 33Mg. The modified MCSM calculations with the 2p3/2 orbital lowered by 1 MeV are shown in Table 1 (MCSMmod) and in Fig. 3(b) by the solid line. This leads to a better reproduction of the momentum distribution. The magnetic moment values predicted by the MCSM is −0.750 while that with the modified MCSM is −0.763 showing that it is minimally affected by the lowering of the 2p3/2 orbital. Both these values are consistent with the experiment [16]. The present data thus provides a sensitive way to understand the gap between the 1f7/2 and the 2p3/2 orbitals which is otherwise not known experimentally yet. The influence of the 2p3/2 orbital is also discussed to be important towards understanding the magnetic moment of 34Al [31].We performed a search for a best fit to the data from χ2 minimization with the spectroscopic factors as free parameters for the following configurations in the 3/2− ground state of 33Mg. The ground state is described by α1⋅Mg21+32ϕn(2p3/2)+α2⋅Mg1−32ϕn(2s1/2)+α3⋅Mg21+32ϕn(1f7/2)+α4⋅Mg1−32ϕn(1d3/2)+α5⋅Mggs32ϕn(2p3/2). Si=αi2 is the spectroscopic factor. The maximum allowed values of the spectroscopic factor are restricted to the respective sum rule values for each case.The 2s1/2 orbital in the 32Mg core is occupied in the ground state as well as in the 2+ excited states. The average occupation number predicted by MCSM is 1.899 and 1.891 for the first two 2+ states. This would mean that in the configuration for 33Mg there is not much possibility of placing the odd-neutron in the 2s1/2 orbital outside the 32Mg core. In the configuration mixing considered in the fitting we have therefore restricted the allowed occupancy for the valence neutron of 33Mg to be in the 2s1/2 orbital to be 0.1.Fig. 4 shows the best fit (solid) lines to the data with the individual components shown as S1 (dashed line), S2 (dashed–dotted line), S3 (dotted line) and S4 (dashed–double-dotted line). The resultant spectroscopic factors were found to be S1=0.5−0.3+0.7, S2=0.1−0.1+0.0, S3=0.5−0.5+0.2, S4=1.5−1.3+1.3 and S5=0.6−0.5+0.3. Inclusion of higher core excited states in the fitting did not lead to any strength in them.Considering the ground state spin of 33Mg to be 3/2+ we looked into a configuration mixing where the ground state can be described by α1⋅Mg3−32ϕn(2p3/2)+α2⋅Mg22+32ϕn(2s1/2)+α3⋅Mg3−32ϕn(1f7/2)+α4⋅Mggs32ϕn(1d3/2). The spectroscopic factors for the different configurations Si=(αi)2 were obtained from a χ2 minimization to the data. The values were found to be S1=2.2−0.5+0.2, S2=0.1−0.1+0.0, S3=1.1−0.5+0.1, S4=0.0−0.0+0.5. The resultant momentum distributions are shown in Fig. 5.Therefore, the important observation in this work is that the narrow observed momentum distribution is a clear signature of significant occupation of the 2p3/2 orbital. It is much larger than the prediction from the MCSM, which is unable to describe the data. The p-wave spectroscopic strength is larger than that observed for 32Mg from one-neutron removal [27]. In 32Mg the occupancy of the neutron in the 1f7/2 orbital (1.19(36)[27]) is larger than the 2p3/2 orbital (0.59(11)[27]). In 33Mg greater occupancy in 2p3/2 is observed (1.1−0.58+0.76 for Jgsπ=3/2− and 2.2−0.5+0.2 for Jgsπ=3/2+) compared to 1f7/2 (0.5−0.5+0.2 for Jgsπ=3/2− and 1.1−0.50.1 for Jgsπ=3/2+). This enhanced probability of neutrons residing in the 2p3/2 orbital shows the lowering of the 2p3/2 orbital thereby confirming that the island of inversion is well pronounced beyond N=20. The modified MCSM calculation with the 2p3/2 orbital lowered by 1 MeV leads to a better reproduction of the data lending support to this view.Although these nuclei are deformed, in this case we do not employ Nilsson orbitals to interpret the data due to several reasons. It is already known from the study of 11,12Be [24,25] where deformation plays a key role in the inversion of the 2s1/2 and 1p1/2 orbitals that the momentum distributions are well explained using spherical single-particle wavefunctions. It was also applied successfully to the study of 28Ne [26] and 32Mg [27]. This reflects the fact that these distributions clearly exhibit the orbital angular momentum of the single-particle orbital. Moreover, the deformation in these neutron-rich nuclei has been understood in the context of particle–hole excitations within the spherical shell model. A Nilsson model that includes orbital re-arrangement in spherical configurations for neutron-rich nuclei is yet to be developed. Furthermore, its extension to reaction theory also needs to be developed. Moreover the excited states of 32Mg do not exhibit the spectrum of a rigid rotor.In summary, the longitudinal momentum distribution of one-neutron removal from 33Mg was measured at 898A MeV at the FRS of GSI. The one-neutron removal cross-section was measured to be 74±4 mb. The narrow structure of the momentum distribution clearly supports significant occupancy in low angular momentum orbitals. Within the eikonal model, the data could be explained with a large spectroscopic strength for the 2p3/2 orbital. The data could not be described by the spectroscopic strengths predicted by the MCSM with the SDPF-M interaction. The modified MCSM calculation with the 2p3/2 orbital lowered by 1 MeV shows an improved agreement with the data. The magnetic moment is found to be less sensitive to this lowering and is consistent with experimental observation both for MCSM and modified MCSM. This shows the lowering of the 2p3/2 orbital due to deformation and weaker binding in 33Mg. Thus the island of inversion becomes more pronounced across the N=20 gap. The present data has some sensitivity to the relative placing of the 1f7/2 and 2p3/2 orbitals which is otherwise not known experimentally yet.AcknowledgementsThe authors are thankful for the support of the GSI accelerator staff and the FRS technical staff for an efficient running of the experiment. The support from NSERC for this work is gratefully acknowledged. R. Kanungo gratefully acknowledges the kind support of the Alexander von Humboldt foundation and the kind hospitality of GSI. This work was supported by the BMBF under contract 06MT238, by the DFG cluster of excellence Origin and Structure of the Universe. T. Nilsson is supported by a grant from the Knut and Alice Wallenberg Foundation.References[1]B.H.WildenthalW.ChungPhys. Rev. C2219802260[2]C.ThibaultPhys. Rev. C121975644[3]T.MotobayashiPhys. Lett. B34619959[4]B.V.PritychenkoPhys. Lett. B4611999322[5]H.ScheitNucl. Phys. A746200496[6]K.YonedaPhys. Lett. B4992001233[7]A.GadePhys. Rev. Lett.992007072502[8]B.V.PritychenkoPhys. Rev. C632000011305R[9]Y.YanagisawaPhys. Lett. B566200384[10]P.DoornenbalPhys. Rev. Lett.1032009032501[11]E.K.WarburtonPhys. Rev. C4119901147[12]Y.UtsunoPhys. Rev. C601990054315[13]Y.UtsunoPhys. Rev. C642001011301R[14]T.OtsukaPhys. Rev. Lett.872001082502[15]S.NumellaPhys. Rev. C642001054313[16]D.T.YordanovPhys. Rev. Lett.992007212501[17]V.TripathiPhys. Rev. Lett.1012008142504[18]I.HamamotoPhys. Rev. C762007054319[19]M.BelleguicPhys. Rev. C722005054316[20]Zs.DombradiPhys. Rev. Lett.962001182501[21]A.ObertelliPhys. Lett. B332006633[22]V.TripathiPhys. Rev. Lett.942005162501[23]V.TripathiPhys. Rev. C762007021301R[24]T.AumannPhys. Rev. Lett.84200035[25]A.NavinPhys. Rev. Lett.852000266[26]J.R.TerryPhys. Lett. B640200686[27]J.R.TerryPhys. Rev. C772008014316[28]G.NeyensPhys. Rev. Lett.942005022501[29]B.V.PritychenkoPhys. Rev. C652002061304[30]Y.OgawaNucl. Phys. A5711994784[31]P.HimpePhys. Lett. B6582008203