Recent Advances in the Microscopic Calculations of Level Den- sities by the Shell Model Monte Carlo Method

Y. Alhassid, M. Bonett-Matiz, S. Liu, A. Mukherjee, Hitoshi Nakada

Published Apr 1, 2014• 3 citations

Abstract

The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (i) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (ii) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes 59−64Ni and of a heavy deformed rare-earth nucleus 162Dy and found them to be in close agreement with various experimental data sets.

DOI: 10.1051/epjconf/20136900010/pdf

Related papers

Shell model Monte Carlo approach: the heavy nuclei

Y. Alhassid, L. Fang, Hitoshi Nakada

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.

Journal of Physics Conference SeriesJan 1, 20111 citation

Circumventing the odd particle-number sign problem in the shell model Monte Carlo

Y. Alhassid, P. Fanto, C. Özen

The shell model Monte Carlo (SMMC) method is a powerful method for calculating exactly (up to statistical errors) thermal observables and statistical properties of atomic nuclei. However, its application has been limited by a sign problem at low temperatures that arises from the projection onto odd particle number even for good-sign interactions. Here, we develop a technique - the partition function extrapolation method (PFEM) - to extract the ground-state energy of an odd-mass nucleus from the excitation partition function calculated at temperatures at which this sign problem is moderate. We validate the PFEM in heavy even-mass nuclei and systematically calculate ground-state energies for isotopic chains of heavy odd-mass nuclei. The PFEM can be extended to other finite-size quantum many-body systems.

arXiv:2308.14698Aug 28, 20230 citations

Extraction of Spectra in the Shell Model Monte Carlo Method Using Imaginary-Time Correlation Matrices

Y. Alhassid, M. Bonett-Matiz, C. N. Gilbreth +1 more

Nuclear energy levels are usually calculated using conventional diagonalization methods in the framework of the configuration-interaction (CI) shell model but these methods are prohibited in very large model spaces. The shell model Monte Carlo (SMMC) method is a powerful technique for calculating thermal and ground-state observables of nuclei in very large model spaces, but it is challenging to extract nuclear spectra in this approach. We present a novel method to extract within SMMC low-lying energy levels for given values of a set of good quantum numbers such as spin and parity. The method is based on imaginary-time correlation matrices (ITCMs) of one-body densities that satisfy asymptotically a generalized eigenvalue problem. We validate the method in a light nucleus that allows comparison with exact diagonalization results of the CI shell-model Hamiltonian. The method is broadly applicable to quantum many-body systems in other disciplines.

Physical Review LettersOct 28, 20241 citation

Signatures of phase transitions in nuclei at finite excitation energies

Y. Alhassid, C. Özen, H. Nakada

The mean-field approximation predicts pairing and shape phase transitions in nuclei as a function of temperature or excitation energy. However, in the finite nucleus the singularities of these phase transitions are smoothed out by quantal and thermal fluctuations. An interesting question is whether signatures of these transitions survive despite the large fluctuations. The shell model Monte Carlo (SMMC) approach enables us to calculate the statistical properties of nuclei beyond the mean-field approximation in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We have extended the SMMC method to heavy nuclei and used it to study the transition from vibrational (spherical) to rotational (deformed) nuclei in families of rare-earth isotopes. We have calculated collective enhancement factors of level densities as a function of excitation energy and found that the decay of the vibrational and rotational enhancements is well correlated with the pairing and shape phase transitions, respectively.

arXiv:1210.3643AIP conference proceedingsJan 1, 20120 citations

Signatures of phase transitions in nuclei at finite excitation energies

Y. Alhassid, C. Özen, Hitoshi Nakada

Abstract. The mean-field approximation predicts pairing and shape phase transitions in nuclei as a function of temperature or excitation energy. However, in the finite nucleus the singularities of these phase transitions are smoothed out by quantal and thermal fluctuations. An interesting question is whether signatures of these transitions survive despite the large fluctuations. The shell model Monte Carlo (SMMC) approach enables us to calculate the statistical properties of nuclei beyond the mean-field approximation in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We have extended the SMMC method to heavy nuclei and used it to study the transition from vibrational (spherical) to rotational (deformed) nuclei in families of rare-earth isotopes. We have calculated collective enhancement factors of level densities as a function of excitation energy and found that the decay of the vibrational and rotational enhancements is well correlated with the pairing and shape phase transitions, respectively.

Aug 17, 20160 citations

Extracting spectra in the shell model Monte Carlo method using imaginary-time correlation matrices

Y. Alhassid, M. Bonett-Matiz, C. N. Gilbreth +1 more

Conventional diagonalization methods to calculate nuclear energy levels in the framework of the configuration-interaction (CI) shell model approach are prohibited in very large model spaces. The shell model Monte Carlo (SMMC) is a powerful technique for calculating thermal and ground-state observables of nuclei in very large model spaces, but it is challenging to extract nuclear spectra in this approach. We present a novel method to extract low-lying energy levels for given values of a set of good quantum numbers such as spin and parity. The method is based on imaginary-time one-body density correlation matrices that satisfy asymptotically a generalized eigenvalue problem. We validate the method in a light nucleus that allows comparison with exact diagonalization results of the CI shell model Hamiltonian. The method is applicable to other finite-size quantum many-body systems that can be described within a CI shell model approach.

arXiv:2401.00613Dec 31, 20230 citations