We extend the shell model Monte Carlo approach to heavy deformed nuclei using a new proton-neutron formalism. The low excitation energies of such nuclei necessitate low-temperature calculations, for which a stabilization method is implemented in the canonical ensemble. We apply the method to study a well-deformed rare-earth nucleus, 162Dy. The single-particle model space includes the 50-82 shell plus 1f_{7/2} orbital for protons and the 82-126 shell plus 0h_{11/2}, 1g_{9/2} orbitals for neutrons. We show that the spherical shell model reproduces well the rotational character of 162Dy within this model space. We also calculate the level density of 162Dy and find it to be in excellent agreement with the experimental level density, which we extract from several experiments.
The Shell Model Monte Carlo (SMMC) approach has been applied to calculate level densities and partition functions to temperatures up to ∼ 1.5 − 2 MeV, with the maximal temperature limited by the size of the configuration space. Here we develop an extension of the theory that can be used to higher temperatures, taking into account the large configuration space that is needed. We first examine the configuration space limitation using an independent-particle model that includes both bound states and the continuum. The larger configuration space is then combined with the SMMC under the assumption that the effects on the partition function are factorizable. The method is demonstrated for nuclei in the iron region, extending the calculated partition functions and level densities up to T ∼ 4 MeV. We find that the back-shifted Bethe formula has a much larger range of validity than was suspected from previous theory. The present theory also shows more clearly the effects of the pairing phase transition on the heat capacity.
We apply the configuration-interaction method to calculate the spectra of two-component Fermi systems in a harmonic trap, studying the convergence of energies at the unitary interaction limit. We find that for a fixed regularization of the two-body interaction the convergence is exponential or better in the truncation parameter of the many-body space. However, the conventional regularization is found to have poor convergence in the regularization parameter, with an error that scales as a low negative power of this parameter. We propose a new regularization of the two-body interaction that produces exponential convergence for systems of three and four particles. We estimate the ground-state energy of the four-particle system to be (5.045 +/- 0.003) variant Planck's constant over 2 pi omega.
We introduce a simple model to calculate the nuclear moment of inertia at finite temperature. This moment of inertia describes the spin distribution of nuclear levels in the framework of the spincutoff model. Our model is based on a deformed single-particle Hamiltonian with pairing interaction and takes into account fluctuations in the pairing gap. We derive a formula for the moment of inertia at finite temperature that generalizes the Belyaev formula for zero temperature. We show that a number-parity projection explains the strong odd-even effects observed in shell model Monte Carlo studies of the nuclear moment of inertia in the iron region.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.