Magic Numbers of Ultraheavy Nuclei V

Home , Magic number (physics)

Physics of Atomic Nuclei, Vol. 68, No. 7, 2005, pp. 1133–1137. Translated from Yadernaya Fizika, Vol. 68, No. 7, 2005, pp. 1179–1184. Original Russian Text Copyright c 2005 by Denisov. NUCLEI Theory

Magic Numbers of Ultraheavy Nuclei V. Yu. Denisov Institute for Nuclear Research, National Academy of Sciences of Ukraine, pr. Nauki 47, Kiev, 03028 Ukraine Received March 30, 2004; in ﬁnal form, May 12, 2004

Abstract—For nuclei where the number of protons lies in the range 76 ≤ Z ≤ 400, proton and neutron shell corrections are calculated along the beta-stability line described by Green’s formula. The magic numbers of protons and neutrons are determined for ultraheavy nuclei. Alpha-decay half-lives and ﬁssion barriers are estimated for ultraheavy doubly magic nuclei. c 2005 Pleiades Publishing, Inc.

INTRODUCTION local minimum in the vicinity of a magic number [4, 9–14]. Calculating shell corrections for spherical nu- Magic numbers corresponding to the filling of nuclei involving various numbers of protons and neu- clear shells of single-particle levels have been known trons, one can therefore determine magic numbers since the middle of the past century, and the role from the positions of deep local minima in the proton that they play in nuclear physics is of crucial impor- and neutron shell corrections. It should be noted that tance [1–4]. Nuclei where the numbers of nucleons nuclei lying along the beta-stability line and having are magic are more stable and have a higher binding empirically known magic numbers of nucleons (Z = energy than their neighborsandaremoreabundant 8, 20, 28, 50, 82, N =8, 20, 28, 50, 82, 126) are in nature than them [1–4]. Many quantities, such as spherical [4, 11]. In the following, we will therefore the energies required for the separation of one and two also explore shell corrections in spherical nuclei. nucleons, the energies of alpha and beta transitions, pairing energies, and the excitation energies of low- lying vibrational states, undergo discontinuities upon CALCULATION OF SHELL CORRECTIONS passing a magic number [1–4]. Figure 1 shows the proton (δP ), neutron (δN ), The magic numbers Z =82and N = 126 are the and total (δP + δN ) shell corrections calculated for greatest magic numbers that have been empirically 76 ≤ Z ≤ 400 even–even spherical nuclei lying along confirmed to date for protons and neutrons, respec- the beta-stability line approximated by Green’s for- tively. It should be noted, however, that, at matter mula [15], from which it follows that a nucleus involv- densities close to nuclear-matter densities, there oc- ing Z protons and NGreen(Z)=(2/3)Z +(5/3) × curs a transition in neutron stars from nuclei, neutron (10 000 + 40Z + Z2)2 − 500/3 neutrons corresponds drops, neutrons, and protons to fusing and decaying to the beta-stability valley [15]. Green’s formula heavy nuclei [5, 6], which transform, as the density describes well the beta-stability line, which is asso- increases further, into a more complicated state of ciated with a specific relation between the numbers nuclear matter [5–7]. Therefore, very heavy nuclei can of protons and neutrons for nuclei known to date. be formed in neutron stars. The existence of neutron- Let us assume that that the relation between the rich nuclei where the number of neutrons is about 3 5 numbers of protons and neutrons in beta-stable 10 –10 and of supercharged nuclei where the num- nuclei that is described by Green’s formula is valid ber of protons is about 1600 is discussed in [6]. It for heavier nuclei inclusive. Figure 1 shows the shell would be interesting to find magic numbers in ultra- corrections calculated for nuclei involving an even heavy nuclei where the number of nucleons falls with- number Z of protons in the range from 76 to 400 ≤ ≤ in the range 300 A 1200. In neutron stars, the andanevennumberN of neutrons in the range relative production rate for ultraheavy nuclei involv- NGreen − 10 ≤ N ≤ NGreen +10,whereNGreen is the ing a magic number of nucleons would be enhanced even number closest to NGreen(Z).Thenumbersof because of their higher stability. The results presented neutrons and nucleons in nuclei were varied in the in [8] also provoke interest in studying magic numbers range 102 ≤ N ≤ 820 and in the range 178 ≤ A ≤ in the region of superheavy nuclei. 1218, respectively. It is well known that magic numbers correspond The energies of single-particle levels of nucleons to the filling of nucleon shells in spherical beta-stable were calculated for the nucleon mean field in the nuclei [1–4, 9–12]. The shell correction has a deep form of the Woods–Saxon potential with allowance

1063-7788/05/6807-1133$26.00 c 2005 Pleiades Publishing, Inc. 1134 DENISOV

δ δ P + N, MeV 12

–6

–12

–18 200 400 600 800 1000 1200 A δ N, MeV

–4 228 184 –8 126 644 772 308 524 –12 406 100 300 500 700 N δ P, MeV

114 210 –4

164 354 82 274 –8 100 200 300 400 Z

Fig. 1. Proton (δP ), neutron (δN ) and total (δP + δN ) shell corrections for even–even spherical nuclei. for spin–orbit and Coulomb interactions [2–4, 11, levels in light, medium-mass, heavy, and superheavy 16]. We employed a “universal” set of parameters of spherical and deformed nuclei. Also, it was success- the Woods–Saxon potential [16]. This set makes it fully used to calculate various properties of nuclei [14, possible to describe well the spectra of single-particle 16, 17]. The residual pairing interaction of nucleons

PHYSICS OF ATOMIC NUCLEI Vol. 68 No. 7 2005 MAGIC NUMBERS OF ULTRAHEAVY NUCLEI 1135 was taken into account in the Lipkin–Nogami ap- estimate the alpha-decay half-lives and fission bar- 298 proximation [13, 18] , the coupling constant being riers for the super- and ultraheavy nuclei 114184, 472 616 798 set to rmic =3.30 fm [13]. In order to calculate shell 164308, 210406,and 274524. corrections, we employed a basis formed by oscillator The alpha-decay half-lives for these nuclei will be wave functions for the lowest 35 shells of an axially found with the aid of the phenomenological Viola– deformed harmonic oscillator and took into account Seaborg formula [19], which relates the alpha-decay the lowest 2330 single-particle levels. The degree of half-life to the energy of alpha particles and the charge the correcting polynomial was chosen to be six. This of the primary nucleus. For the constants of the phe- choice is conventional in calculating shell corrections nomenological Viola–Seaborg formula, the authors for medium-mass, heavy, and superheavy nuclei [11]. of [20] found values that made it possible to repro- The energies of single-particle levels of nucleons were duce faithfully the experimental half-lives of 58 nuclei calculated with the aid of the WSBETA code [16], heavier than 208Pb . Knowing the total shell cor- which was refined in order to take into account greater 126 rections calculated here and the macroscopic binding numbers of shells and levels. energies of nuclei as calculated by means of the mass From the data in Fig. 1, it follows that the pro- formula from [13], we determine the energies of alpha ton shell corrections have deep local minima at Z = particles (in MeV) emitted by the aforementioned 82, 114, and 164 and that the neutron shell cor- super- and ultraheavy nuclei. We have rections have deep local minima at N = 126, 184, 298 472 and 228. It should be noted that Z =82, 114, and Q( 114184) ≈ 9.4,Q( 164308) ≈ 13.1, 164 and N = 126, 184, and 228 are or are assumed Q(616210 ) ≈ 20.9,Q(798274 ) ≈ 35.0. to be magic numbers [4]. For example, the empir- 406 524 ically known magic numbers Z =82 and N = 126 After that, we determine the half-lives with respect correspond to the doubly magic spherical nucleus to the alpha decay of these nuclei with the aid of the 208 Pb126. The values that we found for the shell cor- modified Viola–Seaborg formula [20]. The results are 208 rections in the doubly magic nuclei Pb126 and (in s) 298 114184 are in good agreement with their counter- T (298114 ) ≈ 1.1 × 102, parts calculated in [13, 14]. Thus, our calculations 1/2 184 472 4 reproduce known results and make it possible to per- T1/2( 164308) ≈ 2.3 × 10 , form an extrapolation to the region of heavier nuclei. 616 −6 T1/2( 210406) ≈ 4.2 × 10 , The proton shell corrections have three deep local 798 ≈ × −21 minima in the region 164

PHYSICS OF ATOMIC NUCLEI Vol. 68 No. 7 2005 1136 DENISOV

E, MeV 16

β 0 0.2 0.4 0.6 2

298 Fig. 2. Deformation energy for the super- and ultraheavy doubly magic nuclei (solid curve) 114184, (dashed curve) 472 616 798 164308 , (dotted curve) 210406 , and (dash-dotted curve) 274524 .

472 616 From the data in Fig. 2, it follows that the fission 164308 and 210406 could therefore readily be 298 barrier in the 114184 nucleushasatwo-humped measured if they were formed. shape. The identical fission-barrier shape was ob- The nucleus involving the magic number Z = 164 298 tained in [14] for the 114184 nucleus. The inner- of protons can be formed, for example, in the fusion and outer-barrier heights determined from the da- of two lead nuclei, whereas the nucleus involving the ta in Fig. 2 are larger than their counterparts cal- magic number N = 308 of neutrons can be generated 252 culated in [14], this being due to the fact that the in the fusion of two Cf154 nuclei. However, nuclei multipole surface deformations β2, β3, β4, ..., β8 arising in these reactions are rather far off the beta- were taken into account in [14]. However, the inner- stability line. Ultraheavy doubly magic nuclei can be fission-barrier heights determined from Fig. 2 and produced in a collision of two heavy neutron-rich in [14] are rather close, because, at small values of β2, nuclei that is accompanied by the absorption of many the effect of deformations of higher multipole orders neutrons. on the barrier height is weak. This makes it possible to estimate the shape and the height of fission barriers in ultraheavy doubly magic nuclei. DISCUSSION OF THE RESULTS From Fig. 2, one can see that the ultraheavy dou- AND CONCLUSION 472 616 798 bly magic nuclei 164308, 210406,and 274524 Recently, there appeared the article of Zhang have rather high but narrow one-humped fission bar- et al. [22], who studied magic numbers for 100 ≤ Z ≤ riers. With increasing nuclear charge, the fission bar- 140 nuclei within the relativistic continuum Hartree– rier becomes narrower and occurs at smaller values of Bogolyubov approximation, employing various ver- β2. An increase in the barrier height with increasing sions of microscopic forces. Various nuclear shapes mass number in these doubly magic nuclei is due to were also taken into account in that article. It is an increase in the amplitude of the shell correction well known that the inclusion of nuclear deformations (see Fig. 1). leads to the emergence of additional local minima in 298 The fission half-life of the 114184 nucleus is the dependence of shell corrections on the number about 1010 times longer than its alpha-decay half- of nucleons [9–11, 14]. These minima, which are life [14]. The lifetimes of the 114 ≤ Z ≤ 120 neigh- associated with the filling of shells in deformed nuclei, boring nuclei are also determined by their alpha- correspond to “quasimagic” numbers. It should be decay periods [14, 17]. Similarly, the lifetimes of the noted that deformed nuclei having filled shells—that 472 616 798 164308, 210406,and 274524 nuclei are related is, quasimagic numbers of nucleons—also possess to their alpha-decay half-lives, since the fission bar- enhanced stability and other properties inherent in riers in these nuclei are rather high. As was indi- spherical beta-stable nuclei involving magic numbers cated above, the half-lives of the doubly magic nuclei of nucleons. However, the amplitude of magicity

PHYSICS OF ATOMIC NUCLEI Vol. 68 No. 7 2005 MAGIC NUMBERS OF ULTRAHEAVY NUCLEI 1137 effects in deformed nuclei is less than that in spher- fission barriers in them are rather large; therefore, ical nuclei. In the region 100 ≤ Z ≤ 140,Zhang searches for such nuclei are of great interest. et al. [22] found a greater number of magic numbers than the present author and the authors of [4, 12]. Some of the magic numbers found in [22] are REFERENCES quasimagic. However, the magic numbers Z = 114, ≤ ≤ 1. M. Goeppert-Mayer and J. Jensen, Elementary The- N = 184,andN = 228 for 100 Z 140 nuclei ory of Nuclear Shell Structure (Wiley, New York, were found here and in [4, 22] as well. 1955; Inostrannaya Literatura, Moscow, 1958). The shell corrections for 40 ≤ Z ≤ 200, 40 ≤ N ≤ 2. V. G. Solov’ev, Theory of Atomic Nuclei (Energoiz- 420 nuclei were calculated in [12] in the Hartree– dat, Moscow, 1981; Institute of Physics, Bristol, Eng- Fock approximation with various versions of Skyrme land, 1992). forces and on the basis of the relativistic mean-field 3. P. Ring and P. Schuck, The Nuclear Many-Body model. Within various versions of the calculations, Problem (Springer-Verlag, New York, 1980). magic numbers in the interval Z ≈ 114–126, an in- 4. S. G. Nilsson and I. Ragnarsson, Shapes and Shells terval around Z ≈ 164, the interval N ≈ 172–184, in Nuclear Structure (Cambridge Univ. Press, Cam- and an interval around N ≈ 308 were found in that bridge, 1995). study. It is worth noting that, in the figures presented 5. S. L. Shapiro and S. A. Teukolsky, Black Holes, in [12], deep minima in the shell corrections are seen White Dwarfs, and Neutron Stars: the Physics in the regions close to N ≈ 228 and N ≈ 406;unfor- of Compact Objects (Wiley, New York, 1983; Mir, tunately these minima are not discussed there. Moscow, 1985). 6. A. B. Migdal, Fermions and Bosons in Strong Because of the Coulomb repulsion of protons, a Fields (Nauka, Moscow, 1978) [in Russian]; Rev. region depleted in nucleons can arise at the center Mod. Phys. 50, 107 (1978). of super- and ultraheavy nuclei [12, 23]. This effect 7. A. B. Migdal, E. E. Saperstein, M. A. Troitsky, and was taken into account in [12]. However, the doubly D. N. Voskresensky, Phys. Rep. 192, 179 (1990). 472 magic nucleus 164308 found here also proved to 8. S. V. Adamenko and A. S. Adamenko, nucl- be doubly magic in the calculations performed in [12] ex/0307011. by using some parametrizations of Skyrme forces. 9. V. M. Strutinsky, Nucl. Phys A 95, 420 (1967). The Coulomb repulsion of protons at the center of 10. V. M. Strutinsky, Nucl. Phys. A 122, 1 (1968). a nucleus could affect the values of the magic num- 11. V. Brack, J. Damgaard, A. S. Jensen, et al.,Rev. bers for heavier nuclei, but a detailed investigation of Mod. Phys. 44, 320 (1972). this effect would require substantially changing the 12. M. Bender, M. Nazarewicz, and P.-G. Reinhard, parametrization of the nucleon-mean-field potential. Phys. Lett. B 515, 42 (2001). 13. P.Moller, J. R. Nix, W.D. Myers, and W.J. Swiatecki, It should be noted that a parametrization that would At. Data Nucl. Data Tables 59, 185 (1995). take into account a decrease in the density at the 14. R. Smolanczuk, Phys. Rev. C 56, 812 (1997). center of a nucleus has not yet been investigated. 15. A. E. S. Green, Nuclear Physics (McGraw-Hill, New Therefore, the effect of the Coulomb repulsion of pro- York, 1955). tons at the center of super- and ultraheavy nuclei was 16. S. Cwiok, H. Dudek, W. Nazarewicz, and J.Skalsi, not considered in the present study. Comput. Phys. Commun. 46, 379 (1987). In [12, 22], it was indicated that the values of 17. I. Muntian, Z. Patyk, and A. S. Sobczewski, Acta magic numbers in super- and ultraheavy nuclear re- Phys.Pol.B32, 691 (2001); 34, 2141 (2003). gions depend both on the choice of model and on the 18. H. C. Pradhan, Y. Nogami, and J. Law, Nucl. Phys. A choice of parameters of forces in microscopic calcu- 201, 357 (1973). lations. Therefore, it would be of interest to perform 19. V. E. Viola and G. T. Seaborg, J. Inorg. Nucl. Chem. investigations similar to those in [12, 22] for heavier 28, 741 (1966). 20. R. Smolanczuk, J. Skalski, and A. Sobiczewski, in nuclei and to compare the results of such investiga- Proceedigns of the International Workshop 24 on tion with the results obtained here. Gross Properties of Nuclei and Nuclear Excita- By studying shell corrections, we have determined tions, Ed. By H. Feldmeier et al. (GSI, Darmstadt, the proton magic numbers Z = 114, 164, 210, 274, 1996), p. 35. and 354 and the neutron magic numbers N = 184, 21. O. Moller and J. R. Nix, At. Data Nucl. Data Tables 228, 308, 406, 524, 644, and 772 in super- and ul- 26, 165 (1981). traheavy nuclear regions. Ultraheavy nuclei involving 22. W. Zhang, J. Meng, S. Q. Zhang, et al.,nucl- magic numbers of protons and neutrons are expected ex/0403021. to be more stable and to possess higher binding en- 23. O. Decharge,´ J.-F. Berger, M. Girod, and Dietrich, ergies than neighboring nuclei. The alpha-decay pe- Nucl. Phys. A 716, 55 (2003). riods of some doubly magic ultraheavy nuclei and Translated by A. Isaakyan

PHYSICS OF ATOMIC NUCLEI Vol. 68 No. 7 2005