Shell Model Monte Carlo Approach Research Articles

Level densities of heavy nuclei in the shell model Monte Carlo approach

Nuclear level densities are necessary input to the Hauser-Feshbach theory of compound nuclear reactions. However, the microscopic calculation of level densities in the presence of correlations is a challenging many-body problem. The configurationinteraction shell model provides a suitable framework for the inclusion of correlations and shell effects, but the large dimensionality of the many-particle model space has limited its application in heavy nuclei. The shell model Monte Carlo method enables calculations in spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods and has proven to be a powerful tool in the microscopic calculation of level densities. We discuss recent applications of the method in heavy nuclei.

The shell model Monte Carlo approach to level densities: Recent developments and perspectives

We review recent advances in the shell model Monte Carlo approach for the microscopic calculation of statistical and collective properties of nuclei. We discuss applications to the calculation of i) level densities in nickel isotopes, implementing a recent method to circumvent the odd-particle sign problem; ii) state densities in heavy nuclei; iii) spin distributions of nuclear levels; and iv) finite-temperature quadrupole distributions.

Direct microscopic calculation of nuclear level densities in the shell model Monte Carlo approach

Direct microscopic calculation of nuclear level densities in the shell model Monte Carlo approach

Nuclear state densities of odd-mass heavy nuclei in the shell model Monte Carlo approach

While the shell model Monte Carlo approach has been successful in the microscopic calculation of nuclear state densities, it has been difficult to calculate accurately state densities of odd-even heavy nuclei. This is because the projection on an odd number of particles in the shell model Monte Carlo method leads to a sign problem at low temperatures, making it impractical to extract the ground-state energy in direct Monte Carlo calculations. We show that the ground-state energy can be extracted to a good precision by using level counting data at low excitation energies and the neutron resonance data at the neutron threshold energy. This allows us to extend recent applications of the shell-model Monte Carlo method in even-even rare-earth nuclei to the odd-even isotopic chains of $^{149-155}$Sm and $^{143-149}$Nd. We calculate the state densities of the odd-even samarium and neodymium isotopes and find close agreement with the state densities extracted from experimental data.

Collective enhancement of nuclear state densities by the shell model Monte Carlo approach

The shell model Monte Carlo (SMMC) approach allows for the microscopic calculation of statistical and collective properties of heavy nuclei using the framework of the configuration-interaction shell model in very large model spaces. We present recent applications of the SMMC method to the calculation of state densities and their collective enhancement factors in rare-earth nuclei.

Recent Advances in the Application of the Shell Model Monte Carlo Approach to Nuclei

The shell model Monte Carlo (SMMC) method is a powerful technique for calculating the statistical and collective properties of nuclei in the presence of correlations in model spaces that are many orders of magnitude larger than those that can be treated by conventional diagonalization methods. We review recent advances in the development and application of SMMC to mid-mass and heavy nuclei.

Collectivity in Heavy Nuclei in the Shell Model Monte Carlo Approach

The microscopic description of collectivity in heavy nuclei in the framework of the configuration-interaction shell model has been a major challenge. The size of the model space required for the description of heavy nuclei prohibits the use of conventional diagonalization methods. We have overcome this difficulty by using the shell model Monte Carlo (SMMC) method, which can treat model spaces that are many orders of magnitude larger than those that can be treated by conventional methods. We identify a thermal observable that can distinguish between vibrational and rotational collectivity and use it to describe the crossover from vibrational to rotational collectivity in families of even-even rare-earth isotopes. We calculate the state densities in these nuclei and find them to be in close agreement with experimental data. We also calculate the collective enhancement factors of the corresponding level densities and find that their decay with excitation energy is correlated with the pairing and shape phase transitions.

Level densities of nickel isotopes: Microscopic theory versus experiment

We apply a spin-projection method to calculate microscopically the level densities of a family of nickel isotopes $^{59-64}$Ni using the shell model Monte Carlo approach in the complete $pfg_{9/2}$ shell. Accurate ground-state energies of the odd-mass nickel isotopes, required for the determination of excitation energies, are determined using the Green's function method recently introduced to circumvent the odd particle-number sign problem. Our results are in excellent agreement with recent measurements based on proton evaporation spectra and with level counting data at low excitation energies. We also compare our results with neutron resonance data, assuming equilibration of parity and a spin-cutoff model for the spin distribution at the neutron binding energy, and find good agreement with the exception of $^{63}$Ni.

Crossover from vibrational to rotational collectivity in heavy nuclei in the shell-model Monte Carlo approach.

Heavy nuclei exhibit a crossover from vibrational to rotational collectivity as the number of neutrons or protons increases from shell closure towards midshell, but the microscopic description of this crossover has been a major challenge. We apply the shell model Monte Carlo approach to families of even-even samarium and neodymium isotopes and identify a microscopic signature of the crossover from vibrational to rotational collectivity in the low-temperature behavior of ⟨J(2)⟩(T), where J is the total spin and T is the temperature. This signature agrees well with its values extracted from experimental data. We also calculate the state densities of these nuclei and find them to be in very good agreement with experimental data. Finally, we define a collective enhancement factor from the ratio of the total state density to the intrinsic state density as calculated in the finite-temperature Hartree-Fock-Bogoliubov approximation. The decay of this enhancement factor with excitation energy is found to correlate with the pairing and shape phase transitions in these nuclei.

Recent developments in the shell model Monte Carlo approach to nuclei

The shell model Monte Carlo (SMMC) approach provides a powerful method for the microscopic calculation of statistical and collective nuclear properties in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods. We discuss recent applications of the method to describe the emergence of collectivity in the framework of the configuration-interaction shell model and the crossover from vibrational to rotational collectivity in families of rare-earth nuclei. We have calculated state densities of these rare-earth nuclei and find their collective enhancement factors to be correlated with the pairing and shape phase transitions. We also discuss an accurate method to calculate the ground-state energy of odd-even and odd-odd nuclei, circumventing the sign problem that originates in the projection on an odd number of particles. We have applied this method to calculate pairing gaps in families of isotopes in the iron region.

Shell model Monte Carlo approach: the heavy nuclei

The auxiliary-field Monte Carlo (AFMC) method, also known as the shell model Monte Carlo (SMMC), enables us to calculate microscopically statistical and collective properties of nuclei in the presence of strong correlations. These calculations can be carried out in shell model spaces that are many orders of magnitude larger than spaces that can be treated with conventional diagonalization methods. A recent major development has been the extension of the AFMC approach to heavy nuclei. Such applications to heavy nuclei have been a major challenge. On the conceptual level, a crucial question is whether a truncated spherical shell model Hamiltonian can describe the proper collectivity observed in heavy nuclei and, in particular, the rotational character of strongly deformed nuclei. On the technical level, the low excitation energies make it necessary to perform calculations down to much lower temperatures. At such low temperatures the propagator becomes ill-conditioned and requires the introduction of stabilization methods.

Pairing Reentrance Phenomenon in Heated Rotating Nuclei in the Shell-Model Monte Carlo Approach

Rotational motion of heated 72Ge is studied within the microscopic shell-model Monte Carlo approach. We investigate the angular momentum alignment and nuclear pairing correlations associated with J^{π} Cooper pairs as a function of the rotational frequency and temperature. The reentrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests itself through the anomalous behavior of specific heat and level density.

Isospin-projected nuclear level densities by the shell model Monte Carlo method

We have developed an efficient isospin projection method in the shell model Monte Carlo approach for isospin-conserving Hamiltonians. For isoscalar observables this method has the advantage of being exact sample by sample. It allows us to take into account the proper isospin dependence of the nuclear interaction, thus avoiding a sign problem that such an interaction introduces in unprojected calculations. We apply the method to calculate the isospin dependence of level densities in the complete $\mathit{pf}+{g}_{9/2}$ shell. We find that isospin-dependent corrections to the total level density are particularly important for $N~Z$ nuclei.

Heavy Deformed Nuclei in the Shell Model Monte Carlo Method

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.

Spin Projection in the Shell Model Monte Carlo Method and the Spin Distribution of Nuclear Level Densities

We introduce spin projection methods in the shell model Monte Carlo approach and apply them to calculate the spin distribution of level densities for iron-region nuclei using the complete (pf + g9/2) shell. We compare the calculated distributions with the spin-cutoff model and extract an energy-dependent moment of inertia. For even-even nuclei and at low excitation energies, we observe a significant suppression of the moment of inertia and odd-even staggering in the spin dependence of level densities.

Nuclear level statistics: Extending shell model theory to higher temperatures

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.

Microscopic nuclear level densities by the shell model Monte Carlo method

We have developed ab initio methods to calculate nuclear level densities within the shell model Monte Carlo approach. We have applied these methods to nuclei in the mass region 50 ≲ A ≲ 70, and found remarkably good agreement with experimental data. Using projection methods in the Monte Carlo approach, we have also studied the parity-and isospin-dependence of level densities.

Signature of a Pairing Transition in the Heat Capacity of Finite Nuclei

The heat capacity of iron isotopes is calculated within the interacting shell model using the complete $({pf+0g}_{9/2})$ shell. We identify a signature of the pairing transition in the heat capacity that is correlated with the suppression of the number of spin-zero neutron pairs as the temperature increases. Our results are obtained by a novel method that significantly reduces the statistical errors in the heat capacity calculated by the shell model Monte Carlo approach. The Monte Carlo results are compared with finite-temperature Fermi gas and BCS calculations.

Unblocking of the Gamow-Teller strength in stellar electron capture on neutron-rich germanium isotopes

We propose a new model to calculate stellar electron capture rates for neutron-rich nuclei. These nuclei are encountered in the core-collapse of a massive star. Using the Shell Model Monte Carlo approach, we first calculate the finite temperature occupation numbers in the parent nucleus. We then use these occupation numbers as a starting point for calculations using the random phase approximation. Using the RPA approach, we calculate electron capture rates including both allowed and forbidden transitions. Such a hybrid model is particularly useful for nuclei with proton numbers Z<40 and neutron numbers N>40, where allowed Gamow-Teller transitions are only possible due to configuration mixing by the residual interaction and by thermal unblocking of $pf$-shell single-particle states. Using the even germanium isotopes Ge-68 to Ge-76 as examples, we demonstrate that the configuration mixing is strong enough to unblock the Gamow-Teller transitions at all temperatures relevant to core-collapse supernovae.

Rotational and pairing properties of 74Rb

Abstract We use a cranked shell model Monte Carlo (SMMC) approach to study rotational properties and pair correlationsin the odd-odd N = Z nucleus 74Rb. The calculation is performed in the complete 1p-0ƒ 5 2 -0g 9 2 model space with a residual interaction derived from the Paris potential. The calculated ground state is dominated by isovector J = 0 proton-neutron (pn) pairing. With increasing frequency the J = 0 pn correlations decrease to a constant value at around 〈J〉 = 3 ± 1.5 h while the isoscalar pn correlations (mainly J = 9) increase. Relatedly, the isospin decreases with frequency from its ground state value T = 1 as isoscalar correlations set in. This finding is in agreement with experiement, where at higher rotational frequency a T = 0 band becomes energetically favored over the T = 1 ground state band.