arxiv: v1 [nucl-th] 28 Oct 2017

Transcription

1 Two-Proton Radioactivity with p halo in light mass nuclei A=8 G. Saxena a, M. Kumawat b, M. Kaushik c, S. K. Jain b, Mamta Aggarwal d, a Department of Physics, Government Women Engineering College, Ajmer-, India b Department of Physics, School of Basic Sciences, Manipal Univ., Jaipur-7, India c Department of Physics, Shankara Institute of Technology, Kukas, Jaipur-8, India d Department of Physics, University of Mumbai, Kalina Campus, Mumbai-98, India arxiv:7.v [nucl-th] 8 Oct 7 Abstract Two-proton radioactivity with p halo is reported theoretically in light mass nuclei A = 8-. We predict 9 Mg, Si, 6 S, Ar and Ca as promising candidates of ground state p-radioactivity with S p < and S p >. Observation of extended tail of spatial charge density distribution, larger charge radius and study of proton single particle states, Fermi energy and the wave functions indicate p halo like structure which supports direct p emission. The Coulomb and centrifugal barriers in experimentally identified p unbound Si show a quasi-bound state that ensures enough life time for such experimental probes. Our predictions are in good accord with experimental and other theoretical data available so far. Keywords: Relativistic mean-field theory; Nilson Strutinsky approach; Two-proton radioactivity; Oneand Two-proton separation energy; Halo nuclei, Proton drip-lines. Two-proton radioactivity, a simultaneous decay of two protons from exotic nuclei near or beyond proton drip line [,,, ], is a recently discovered [,, 6], a relatively lesser known exotic nuclear decay mode that has attracted a lot of attention in the recent times. For a p-radioactive nucleus S p > and S p < which means that p should not be an open decay channel [7] so that the two valence protons are emitted simultaneously where the role of pairing interaction becomes significant. This valence p weakly bound system in the decay channel might form a halo [8, 9,, ] like structure in the low density nuclear matter which is spatially far fromthe tightly bound coresystem andmaylead to enhanced probability of direct p decay. The pairing studies are expected to give new insights into the role of p halo of the valence proton system in the direct p emission. However p correlations [, ] in p decay and its mechanisms are still not clear and require further theoretical and experimental investigations []. So far p-radioactivity has been observed in the decay of ground state Fe [, ], Zn [, 6], 8 Ni [7, 8], 67 Kr [9], 9 Mg [], Ar [], Ar corresponding authors: Mamta Aggarwal, mamta.a@gmail.com [, ] Mg [,, 6], and excited states Mg and 9 Ag (for p and p decay) [7, 8] in experimental and theoretical [9, ] works have provided useful insights into this exotic decay mode of proton rich nuclei. Recent experimental observations of (i) β-delayed two-proton emission in the decay of Si [] (ii) p correlated emissions in 8,9 S but not in 7,8 P which indicated [8, 9, ] that the diproton correlations may result from proton configuration in the initial state like p halo rather than the deformed orbit (iii) and the speculation [] that the p halo plays an important role in the diproton emission, lack a precise theoretical description that has invoked the present study. The objective in this letter is to search for p halo and p unbound nuclear systems in even Z= to in a theoretical framework and present for the first time an analysis to probe the interdependence between the weakly bound p halolikestructureand the direct pemission. p-radioactivity which is essentially the energetically allowed simultaneous two-proton emission is investigated by employing two simple and effective well established theoretical formalisms (i) Relativistic mean-field plus state dependent BCS (RMF+BCS) approach [,,,, 6, 7, 8] Preprint submitted to Physics Letters B June, 8

2 and (ii) Macroscopic-microscopic approach (Mac- Mic) with Nilson Strutinsky (NS) prescription [, 9]. We calculate one- and two-proton separation energies of proton rich Mg, Si, S, Ar and Ca isotopes near proton drip-line and identify 9 Mg, Si, 6 S, Ar and Ca as potential candidates of p-radioactivity in agreement with the available experimental data [,, ] for 9 Mg, Si and Ar. The structural evidences for the existence of halo are searched in the above identified p emitters within the RMF+BCS approach by investigating spatial density distribution, larger charge radius, wave functions and single particle energies that reveal a halo like structure of weakly bound protons which is the highlight of this work. Another important feature here is to look for a quasi bound state using RMF+BCS to provide finite life time to enable experimental study in p unbound Si which has been recently identified a candidate of p-radioactivity experimentally [8]. Although various theoretical approaches [6,,,,, ] have been used to study p unbound nuclei, the RMF+BCS approach together with a realistic mean-field has proved a very useful tool especially for drip line nuclei as shown in our earlier works [7, 8, ] where it has been used to study magicity, weakly bound structure and halos. (Description of detailed theoretical formalisms that are adequately described in our earlier works, are avoided here. Kindly see references (RMF+BCS [7, 8] and Mac-Mic [, 9]) for details). Our calculations show consensus with experimental and other theoretical works and provide much needed insights on this subject. Two-and one- proton separationenergy(s p and S p ) computed using RMF with TMA [], NL* [6], NLSH [7] and PK [8] parameters and Mac- Mic approach [] are shown in Fig. (a) and (b) respectively where we find 9 Mg, Si, 6 S, Ar and Ca to be p unbound with S p < and S p >. Our calculated values of S p and S p and the prediction of the first two-proton unbound nuclei agree very well with the experimental data [9], results from various parameters of RMF and other theories like nonrelativistic approach viz. Skyrme- Hartree-Fock method with the HFB- functional (described in Ref. []) and FRDM []. However, wefindsmallanomalyfor Caand 6 Sforwhichall the theories and experiment show S p < whereas RMF shows Ca to be p bound with S p > and one parameter (NL*) of RMF shows 6 S to be p bound with S p >. This anomaly requires further S P [MeV] S P [MeV] Expt. Mg RMF(TMA) Si S Ar Ca RMF(NLSH) RMF(PK) Mic-Mac HFB- FRDM (a) Mg Si Expt. RMF(TMA) S Ar RMF(NLSH) Ca RMF(PK) Mic-Mac HFB- FRDM (b) Mass Number A Figure : Variation of(a) two- and (b) one- proton separation energy for Mg, Si, S, Ar and Ca isotopes along with available experimental data (with error bars) [9], HFB [] and FRDM [] data. investigation. To probe the structure of p unbound orbitals, we evaluate nuclear deformations using deformed RMF approach [, 6,,, ] and Mac-Mic approach [, 9] where Si and Ca are found to be spherical(β=.) in agreement with HFB [] and FRDM [] and 9 Mg, 6 S and Ar show large deformation which we study along with their neighbouring isotopes to check if the deformation shows any major structural changes in p unbound nuclei while moving from normally bound to unbound nuclei. Our estimated β values for 9 Mg (.,.9,.8), 6 8 S (.6,.6,.6) and Ar (.,.,.) show a gradual variation without any significant structural differences between deformed p emitters and their neighbouring isotopes. This suggests that the deformation appears to have no role in influencing p emission as speculated in recent works [8, 9, ], where, mainly a -body simultaneous decay in phase space and p sequential emission for odd-z 7 P while a diproton-type decay for even-z 8 S indicated that the p halo-like structure rather than large deformation is responsible for p correlations which may play an important role in the mechanism of p emission. However there may not be any taboo to deformations influencing p emission in other nuclei which needs investigation in various mass regions theoretically as well as experimentally to reach any conclusion on this. The p emission can happen either by sequential decay, direct diproton decay or -body demo-

3 Proton Density [fm - ] Proton S.P. Energy [MeV].. E- E- (a) (b) Radius [fm] d s / d / / Si Si 6 Si 8 Si Si Si Si Proton Density Si Si Si 6 Si 8 Si Si Si Si Proton Density HFB- Si d / Z= Weakened Gap s / for Si d / Proton Single Particle States (c) Neutron Number N Figure : Variation of charge density of Si vs. radius for (a) (b) HFB [] (c) Variation of proton single particle states (d /, s /, d / ) for Si using. Fermi energy is depicted by solid squares. cratic decay. However in view of the suggestion [] that p halo would lead to larger spectroscopic factor for a direct p decay than for a sequential decay and the speculations [8, 9,, ] that p halo plays an important role in the diproton emission, we look for the structural evidences for the existence and significance of p halo in p emitters. We study the charge density distribution and variation of proton single particle energies of sd shell (d /, s /, d / )within the sphericalframework of RMF with NL* parameter for the experimentally identified p emitter Si []. From charge density distribution curves plotted in (Fig. (a)) it is evident that Si shows an extended tail in the charge density distribution which is much more extended than its neighboring isotope Si. This indicates weakly bound halo like structure of last two protons in the outer most shell in Si that has more probability of decay by p emission than relatively normally bound Si shows dependence of p halo on p emission. The charge densities of Si calculated by Skyrme-Hartree-Fockmethod with the HFB- functional are also shown in Fig. (b) where loosely bound structure of protons in Si match in a good deal with our calculations of RMF and is a strong evidence in support of decay via two-proton emission. Fig. (c) shows proton single particle energies of sd shell where we note a large gap between d / and s / that produces a shell closure at Z= which is an area of recent interest as Si isotopes Si [] and 8 Si [6] have been identified as doubly magic nuclei in agreement withourcalculationsusingrmffor Si(seeinFig. (c)). The gap between d / and s / in Si is found to be quite large around 6 MeV which reduces to.6 MeV for Si while moving from N= to N=8 towards the drip-line. This gap decreases from7mevto.9mevand7.mevto.mevusing TMA and NLSH parameters respectively. Due to the reduced gap in Si, protons are partially occupied in higher state s / with non zero occupancy of. whereas the other isotopes 8 Si have zero occupancy of proton in s / state. The Fermi energy (shown in Fig. (c)) in Si seems to beveryclosetothisoccupieds / statewhichindicates somewhat loosely bound valence protons near Fermi level in Si comprising a p halo like structure which is more likely to decay by p emission as compared to other isotopes of Si. This analysis provides a strong evidence towards the significant role of p halo in p decay. The charge density Charge Density [fm - ] Mg 9 Mg Si -6 Si S S -6 Ar Ar Figure : Variation of charge density for few isotopes of Mg, Si, S and Ar near drip-line calculated by deformed RMF+BCS approach []. for deformed two-proton emitters 9 Mg, 6 S, Ar along with Si calculated by RMF+BCS approach including deformation [] is shown in Fig.. It is again remarkable to find the well spread density tail for these proton emitters in comparison to their neighboring isotopes indicating halo like structure. Due to extended density distribution of loosely bound protons, the charge radii is expected to grow larger. In order to check this, we compute charge radii of drip-line isotopes of Mg, Si, S, Ar and Ca shown in Fig. along with available experimental data [7]. We observe much larger charge radii of

4 9 Mg, Si, 6 S, Ar (seen very evidently in Fig. ) which is another strong evidence in the support of halo like structure in p emitters and thier possible influence on p decay which has not been shown earlier. Here we used deformed RMF approach with NL*, TMA, NLSH and PK parameters that are in reasonable match with each other. Furthermore the difference between protonradius of Si and the daughter nucleus Mg is found very large which suggests that the removal of two loosely bound protons from Si after p emission results in a more confined configuration of Mg which is more stable with extra binding B/A = 6.7 MeV as compared to B/A = 6.8 MeV of Si. This shows that Si with spherical configuration is a loosely bound drip-line nucleus which decays via p emission and gains stability. Much larger charge radii in deformed p emitters 9 Mg, 6 S and Ar indicates loosely bound halo like structure in well deformed nuclei and is in conformity with the speculation [8, 9,, ] that p halo is responsible for p correlations rather than the large deformation. However, the result for Ca is found albeit in difference for RMF theory where the charge radii is not much enlarged and also S p > although the experimental data [9] and other theories show S p < and support energetically favoured p emission in Ca which needs further investigation. In order to explain this, we evaluate and plot nuclear potential, centrifugal barrier calculated for d state (l=), Coulomb barrier and total potential which is sum of all three in Fig.. The nuclear potential along with Coulomb and centrifugal barriers has sufficient barrier height for the trapping of wavestoformaquasi-boundstatewhichenable Si to acquire rather long mean life time even it is unbound against two-proton emission and eventually Si may decay by the process wherein two-protons tunnel through the barrier leading to observation of two-protons radioactivity (shown by thick solid line infig. (a)). Duetothebarriereffect, suchametastable state remains mainly confined to the region of the potential well and the wave function exhibits characteristics similar to that of a bound state. Indeed wavefunctions of proton valence shells d / and s / of Si are found similar to the wavefunctions of other bound isotopes (as seen in Fig. (b)) and remain in the potential region providing finite life time. For comparison we also show wavefunction of s / state of Si by solid line which significantly lies outside the potential region as seen in Fig. (b) that indicates an unbound nucleus. Interestingly,in O,themirrornucleusof Si, Coulomb barrier is absent for N= therefore neutron s / state is found with -. MeV energy that leads to a perfect bound nucleus. Charge Radii [fm] Si 9 Mg Mg Charge Radii using Deformed RMF 6 S Ar Ca S Isotopes Ca Isotopes Ar Isotopes Mg Isotopes Expt. NL* TMA NLSH PK Neutron Number N Figure : Variation of charge radii for Mg, Si, S, Ar and Ca isotopes near proton drip line. Available experimental data (with error bars) is also shown [7]. Another remarkable feature of this work is to show the quasi bound state in two-proton emitter Sistudiedrecentlyexperimentally[]whichindicates that it must have finite life time large enough to study such p unbound nucleus experimentally. Energy [MeV] Wave Function 8 Centrifugal Barrier Si Coulomb Barrier Nuclear Potential - Total Potential Centrifugal Barrier Potential+Coulomb Barrier Coulomb Barrier (a) Radius R [fm].6-8 Si. (s / ) Si. d /. Si -. s Si Proton Single Particle States / 6 Si -. 8 Si (b) 6 8 Radius R [fm] Figure : (a) The RMF potential energy, centrifugal barrier energy for proton resonant states d / and Coulomb barrier as a function of radius for Si. Total potential (sum of potential energy, Centrifugal and Coulomb barrier), and sum of potential and Coulomb barrier are also shown. (b) Wavefunction of proton single particle d / and s / states are shown for 8 Si. Wavefunction of proton s / state of Si is also shown.

5 To conclude, two-proton radioactivity, a relatively lesser known exotic nuclear decay mode along with p halo are investigated in light mass region around A = 8- by employing (i) RMF+BCS approach and (ii) Macroscopic-microscopic approach with Nilson Strutinsky prescription. We predict 9 Mg, Si, 6 S, Ar and Ca as p-radioactive in their ground state with S p < and S p >. Si and Ca are found to be spherical whereas 9 Mg, 6 S and Ar show large deformations. Structural evidencesfortheexistenceofphalointhepemitters are searched within the RMF+BCS approach. We observe long tail in extended spatial density distribution and larger charge radius in 9 Mg, Si, 6 S, Ar (except Ca) indicate p halo in p emitters. Study of single particle states reveals weakly bound halo like structure of protons due to weakening of shell gap and occupancy in proton s / state shownforexperimentallyidentifiednucleus Sifurther emphasizes the existence of p halo and its role in supporting p emission. The nuclear potential along with Coulomb and centrifugal barriers in Si appear to form a quasi-bound state to ensure a long delay in the decay of Si that may result its existence for finite life time and enable experimental probe. p-radioactivity and the link between p halo and diproton emission need more efforts on the experimental and theoretical fronts. The authors thank Prof. H. L. Yadav, BHU, India, and Prof. L. S. Geng, Beihang University, China for valuable guidance. Authors GS and MA thank financial support by DST, India under YSS//9 and WOS-A respectively. References [] M. Pfutzner et al., Eur. Phys. Jour. A () 79. [] B. Blank and M. J. G. Borge, Prog. Part. Nucl. Phys. 6 (8). [] A. A. Sonzogni, Nucl. Data Sheets 9 (). [] P. J. Woods and C. N. Davids, Annu. Rev. Nucl. Part. Sci. 7 (997). [] B. Blank and M. Ploszajczak, Rep. Prog. Phys. 7 (8) 6, [6] M. Pfutzner, et al., Rev. Mod. Phys. 8 () 67. [7] V. I. Goldansky, Nucl. Phys. 9 (96) 8. [8] X. X. Xu and C. J. Lin et al., Phys. Lett. B 77 () 6. [9] C. J. Lin et al., Phys. Rev. C 8 (9). [] X. X. Xu et al., Phys. Rev. C 8 () 7. [] C. J. Lin et al. Nuclear Physics Review (6) 6. [] T. Maruyama, T. Oishi, K. Hagino, and H. Sagawa, Phys. Rev. C 86 (). [] T. Oishi, M. Kortelainen, and A. Pastore, arxiv:66. (6). [] J. Giovinazzo et. al., Phys. Rev. Lett. 89 (). [] B. Blank et. al., Phys. Rev. Lett. 9 (). [6] P. Ascher et. al., Phys. Rev. Lett. 7 (). [7] C. Dossat et. al., Phys. Rev. C 7 (). [8] M. Pomorski et. al., Phys. Rev. C 8 () 6. [9] T. Goigoux et al., Phys. Rev. Lett. 7 (6) 6. [] I. Mukha et al., Phys. Rev. Lett. 99 (7) 8. [] I. Mukha et al., Phys. Rev. Lett. (). [] G. T. Koldste et al., Phys. Lett. B 77 () 8. [] A. A. Lis et al., Phys. Rev. C 9 () 69. [] J. P. Wallace et al., Phys. Lett. B 7 () 9. [] M. V. Lund et al., Eur. Phys. J. A (6). [6] L. J. Sun et al., Phys. Rev. C 9 (7). [7] Y. G. Ma et al., Phys. Lett. B 7 () 6. [8] I. Mukha et. al., Nature 9 (6) 98. Phys. Rev. Lett. 9 (). [9] Mamta Aggarwal, Phys. Scr. T (6) 78. [] Mamta Aggarwal, Phys. Lett. B 69 () 89. [] X. X. Xu et al., arxiv:6.89v (6). [] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 6 (986). [] J. Boguta and A. R. Bodmer, Nucl. Phys. A 9 (977). [] Y. Sugahara and H. Toki, Nucl. Phys. A 79 (99) 7. [] P. Ring, Prog. Part. Nucl. Phys. 7 (996) 9. [6] H. L. Yadav, M. Kaushik, and H. Toki, Int. Jour. Mod. Phys. E () 67. [7] G. Saxena, D. Singh, M. Kaushik, and S. Somorendro Singh, Canadian Journal of Physics 9 (). [8] G. Saxena, D. Singh, M. Kaushik, H. L. Yadav, and H. Toki, Int. Jour. Mod. Phys. E (). [9] Mamta Aggarwal, Phys. Rev. C 89 (). [] L. V. Grigorenko et al., Phys. Lett. B 677 (9). [] Q. Zhao, J. M.Dong, J. L. Song, and W. H. Long, Phys. Rev. C 9 () 6. [] Y. Lim, X. Xia, and Y. Kim, Phys. Rev. C 9 (6). [] D. S. Delion, R. J. Liotta, and R. Wyss, Phys. Rev. C 87 () 8. [] G. Saxena and M. Kaushik, Chin. J. Phys. (7) 9. [] L. S. Geng, H. Toki, A. Ozawa, and J. Meng, Nucl. Phys. A 7 () 8. [6] G. A. Lalazissis, S. Karatzikos, R. Fossion, D. Pena Arteaga, A. V. Afanasjev, and P. Ring, Phys. Lett. B 67 (9) 6. [7] G. A. Lalazissis, D. Vretenar, and P. Ring, Phys. Rev. C 6 (). [8] W. Long, J. Meng, N. Van Giai, and Shan-Gui Zhou, Phys. Rev. C 69 () 9. [9] M. Wang et al., Chin. Phys. C 6 () 6. [] S. Goriely et al., Phys. Rev. Lett. (9), [] P. Moller, A. J. Sierk, T. Ichikawa, and H. Sagawa, At. Data Nucl. Data Tables, 9 (6). [] D. Singh, G. Saxena, M. Kaushik, H. L. Yadav, and H. Toki, Int. Jour. Mod. Phys. E () 76. [] L. S. Geng, H. Toki, S. Sugimoto, and J. Meng, Prog. Theor. Phys. () 9. [] P. Ring, Y. K. Gambhir, and G. A. Lalazissis, Comput. Phys. Commun. (997) 77 and reference therein. [] A. Mutschler et al., Nature Physics (7). [6] J. J. Li, J. Margueron, W. H. Long, and N. Van Giai, Phys. Lett. B 7 (6) 97. [7] I. Angeli and K. P. Marinovab, At. Data Nucl. Data Tables 99 () 69.