Collective Excitations of Transactinide Nuclei in a Self-Consistent Mean Field Theory
We applied the ATDHFB approach for study of properties of collective quadrupole states in several transactinide nuclei: U, Pu, Pu, Cm, Cm, Cf and Cf. Calculated energies and B(E2) transition probabilities are in a reasonable agreement with experimental data. We present also results concerning superdeformed collective states in the second minimum of potential energy of the Pu nucleus.
Received (October 31, 2007)
Revised (revised date)
Fission properties of heavy and superheavy nuclei have been recently a subject of intensive studies within the frame of a self-consistent mean field theory. Most papers treated height of fission barriers but some results on half-lives of superheavy nuclei have also been published. The standard method of calculating the half-lives is the WKB(J) approach requiring knowledge of the potential and mass parameter(s), which are usually obtained from the Adiabatic Time Dependent Hartree-Fock-Bogolyubov (ATDHFB) or the Generating Coordinate Method (GCM) theory. These general theories can provide a frame for a description of other collective phenomena, e.g. quadrupole rotational-vibrational excitations. However such applications of ATDHFB or GCM in the region of transactinides are rather scarce. On the other hand collective levels and E2 transition probabilities of lighter transactinides (U,Pu,Cm) are well known experimentally. Moreover thanks to recent advances in experimental techniques the region of heavier nuclei accessible for measurements is growing rapidly. In this paper we present several results of application of the ATDHFB theory with the Skyrme interaction for describing quadrupole collective properties of U, Pu, Pu, Cm, Cm, Cf and Cf nuclei, for which there is an abundant experimental data. The main part of the paper is devoted to normally deformed states, i.e., built around the first minimum of the potential energy but Subsection 3.3 contains also results of calculation of superdeformed (connected with the second minimum of the potential) states in the Pu nucleus.
The main tool used in microscopic theory of full five-dimensional quadrupole collective dynamics is the ATDHFB theory, which leads to the generalized Bohr Hamiltonian when appropriate collective variables are chosen. The collective Hamiltonian contains potential energy and mass parameters (including moments of inertia) calculated solely from a microscopic input, i.e., nucleon-nucleon interaction. Details of the method we use and results obtained for lighter nuclei can be found in. Below we recall only some of main points.
Collective quadrupole variables , are defined by components of the second rank tensor of a nuclear mass distribution:
Hence in the axial case for a given nucleus is strictly proportional to the mass quadrupole moment. The self-consistent nuclear mean field is obtained from HF+BCS calculations with a double constraint:
Then the ATDHFB theory allow to calculate mass parameters and collective potential energy that enter the Bohr Hamiltonian. Its eigenvalues are interpreted as collective excitation energies and its eigenfunctions can be used, among others, to obtain E2 transition probabilities.
2.1 Details of calculations
We have chosen as a microscopic interaction the SkM* version of Skyrme forces, proposed already long time ago but still regarded as a good choice especially in the case of barrier heights and other fission properties. We have performed also calculations with the SIII Skyrme interaction which previously gave good results for lighter nuclei and a few tests with the more recent SLy4 forces. In the particle-particle channel we have taken the simplest form of the pairing interaction i.e. of seniority type (constant ). Strength of the pairing was fixed by comparing values of the 5-point formula pairing gap with minimal quasiparticle energies at points corresponding to a minimum of the potential in U and Pu nuclei. Final results for the strength are , MeV, MeV. As we checked in some test cases the results from a state dependent () pairing interaction are almost the same as from the constant force, provided the strength of the interaction is fixed using the same method.
The calculations were made for 220 points in the sextant (, ) of the deformation plane. These point form a regular grid with the distance of 0.05 and between them. The maximal value of (i.e. ) corresponds for the considered nuclei to a quadrupole moment around 100 b. For the lowest normally deformed collective states only the region is important but for studying the second minimum one must expand the range of considerably.
It must be stressed that we exclude the octupole deformation which is essential for fission processes, but for the considered nuclei it becomes important only for larger deformations than those studied by us. As we discussed previously in, see also, pairing vibrations and so called Thouless-Valatin corrections can have a noticeable influence on the mass parameters. In the present work we took a simplified approach introducing an average factor 1.3 by which we multiply values of the mass parameters obtained from the ATDHFB formulas.
For each nuclei we calculate seven functions (the potential energy and six mass parameters, including moments of inertia). In Fig. 1 we show only a small sample, namely the potential energy and the mass parameter for the Pu nucleus.
3.1 Energy levels
In Fig. 2 we show a calculated energy of the first level and of bandheads of and bands and we compare them with experimental data. Please note that the energies (theoretical and experimental) of the level are multiplied by 5 to make the figure more readable.
Keeping in mind the starting point of calculations which was the nucleon-nucleon effective interaction and the fact that we did not fit any parameter to analyzed data, the results shown in Fig. 2 are quite good, even despite too large energies of the theoretical and bandheads.
Below we present also analogous results for the SIII Skyrme interaction (Fig. 3). As it can be seen the and bandheads are now on average closer to experimental values, but their behavior with an increasing mass number does not follow the experimental one.
Detailed inspection of calculated potential energy and mass parameters shows that differences between the results of the SkM* and SIII variant of Skyrme interaction stem mainly from different behavior of the potential in the region of small deformations. It can be readily seen in Fig. 4, where we plot the potential energy for axial shapes of the Pu nucleus obtained using the SkM*, SIII and SLy4 forces. This figure shows also that dependence of on deformation for the SkM* and SLy4 forces is very similar.
Let us mention another remarkable point. In many papers the energy of the first level is estimated from the formula , where is the moment of inertia calculated in the minimum of the potential, in other words assuming a perfect rotor behavior. Within the frame of our approach we can compare the value obtained in such a way with the respective eigenvalue of the Bohr Hamiltonian . It appears that is grater typically by 16–20 keV than . This correction comes from vibrational degrees of freedom and is very small if compared with the energy of, loosely speaking, and phonons but is quite large in relation to .
3.2 E2 electromagnetic transitions
In this subsection we compare theoretical transition probabilities with the experimental ones, taken from. Fig. 5 shows transitions in the considered nuclei while Fig. 6 contains results concerning transitions within the ground state band in the Cm nucleus. Moreover in Table 1 we present some inter-band transitions in the Cf nucleus.
The conclusion is that theoretical results follow experiment very closely. We recall again that we do not use here any additional parameters (effective charge etc.). And again the SIII interaction performs slightly better than SkM*.
3.3 Superdeformed collective states
Theory presented in the previous sections gives also possibility to study superdeformed collective states i.e. built in the vicinity of the second minimum of the potential energy. It needs however a large extension of the basis used to solve the eigenproblem of the Bohr Hamiltonian. We presented such calculations in the approach with the phenomenological Nilsson potential in Ref. , see also papers by Libert et al. using Gogny forces. In the present paper we show results for the Pu nucleus, for which there are many experimental data, see also review in Ref. Moreover in this nucleus an outer barrier is sufficiently high (see also Fig. 4) so as we do not need to worry about e.g. coupling with the continuum.
Fig. 7 shows the probability distribution for the ground state and the lowest state that can be unambiguously identified as a superdeformed one. More precisely we plotted here the product , where is the metric tensor in the collective space.
Table 2 contains more detailed information on the calculated properties of the lowest normally and super deformed states in the Pu nucleus. We show here also the average values and .
Obtained qualitative agreement with experiment looks rather encouraging. Moreover, one must remember that we have not included in our calculation the so called rotational correction often discussed in papers on fission barriers. Its magnitude increases with a deformation and in consequence superdeformed states can be lowered by 0.8-0.9 MeV. However in our opinion a consistent introduction of this correction into our formalism (keeping in view also allowed triaxial shapes of a nucleus) needs more careful treatment.
The ATDHFB approach offers a considerable extension of the area of applicability of the self-consistent mean field theory, namely on low energy collective states, for which there is often rich spectroscopic data. Our study shows that starting with standard Skyrme interactions one can obtain reasonable description of energies and B(E2)’s also for very heavy nuclei. However some problems still remain open, e.g. questions concerning the Thouless-Valatin corrections, pairing vibrations, rotational and other so called zero point energy corrections.
Present work was partially supported by the Academy of Finland and Jyväskylä University within the FiDiPro programme.
- 1. M. Warda, J. L. Egido, L. M. Robledo, K. Pomorski, Phys. Rev. C 66, 014310 (2002).
- 2. L. Bonneau, P. Quentin, D. Samsoen, Eur. Phys. J. A 21, 391 (2004).
- 3. T. Burvenich, M. Bender, J. A. Maruhn, P.-G. Reinhard, Phys. Rev. C 69, 014307 (2004).
- 4. A. Staszczak, J. Dobaczewski, W. Nazarewicz, Int. J. Mod. Phys. E 15, 302 (2006).
- 5. J.-P. Delaroche, M. Girod, H. Goutte, J. Libert, Nucl. Phys. A 771, 103 (2006).
- 6. A. Baran, A. Staszczak, J. Dobaczewski, W. Nazarewicz, Int. J. Mod. Phys. E 16, 443 (2007).
- 7. L. Prochniak, P. Quentin, D. Samsoen, J. Libert, Nucl. Phys. A 730, 59 (2004).
- 8. L. Prochniak, Acta Phys. Pol. B 38, 1605 (2007).
- 9. J. Bartel, P. Quentin, M. Brack, C. Guet, H.-B. Hakansson, Nucl. Phys. A 386, 79 (1982).
- 10. L. Prochniak, K. Zajac, K. Pomorski, S. G. Rohozinski, J. Srebrny, Nucl. Phys. A 648, 181 (1999).
- 11. L. Prochniak, Int. J. Mod. Phys. E 15, 379 (2006).
- 12. J. Libert, M. Girod, J.-P. Delaroche, Phys. Rev. C 60, 054301 (1999).
- 13. Data extracted using the NNDC On-Line Data Service from the ENSDF database, file revised as of Sept. 10, 2007.
- 14. L. Prochniak, K. Zajac, K. Pomorski, S. G. Rohozinski, J. Srebrny, Acta Phys. Pol. B 33, 405 (2002).
- 15. D. Pansegrau, P. Reiter, D. Schwalm, H. Bauer, J. Eberth, D. Gassmann, D. Habs, T. Hartlein, F. Kock, H. G. Thomas, Phys. Lett. B 484, 1 (2000).
- 16. D. Gassmann, P. G. Thirolf, E. Mergel, D. Habs, M. Chromik, J. Domscheit, A. Gorgen, K. Hauschild, H. Hubel, M. Hunyadi, A. Krasznahorkay, A. Lopez-Martens, S. Neumann, A. Neusser, D. Pansegrau, P. Reiter, H. Scheit, G. Schonwasser, D. Schwalm, Phys. Lett. B 497, 181 (2001).
- 17. M. Hunyadi, D. Gassmann, A. Krasznahorkay, D. Habs, P. G. Thirolf, M. Csatlos, Y. Eisermann, T. Faestermann, G. Graw, J. Gulyas, R. Hertenberger, H. J. Maier, Z. Mate, A. Metz, M. J. Chromik, Phys. Lett. B 505, 27 (2001).