The nuclear many-body problem is among the most challenging many-body problems in physics. Our group develops and implements many-body methods to advance the understanding of nuclear properties from underlying effective interactions. These properties are important for understanding astrophysical processes, such as nucleosynthesis and supernovas, and for interpreting results expected from new radioactive beam facilities, such as FRIB.
High-performance computing has enabled significant progress in applying various ab initio methods to calculate properties of nuclei, but these methods are mostly limited to light nuclei. Self-consistent mean-field theories, such as density functional theory (DFT), is unique in providing a global theory of nuclei. However, extending DFT to include correlations beyond the mmean field remains challenging. The configuration-interaction (CI) shell model, a basic model of nuclear structure, provides an attractive framework as it accounts for both shell effects and correlations. However, the combinatorial growth of the dimension of the many-particle space with the number of valence orbitals and/or the number of valence nucleons hinders its application in heavy open-shell nuclei.
We have overcome in part this difficulty in applying the CI shell-model approach to mid-mass and heavy nuclei by using quantum Monte Carlo methods. In particular, we use the auxiliary-field Monte Carlo (AFMC) method, which is known in the context of the nuclear CI shell model as the shell model Monte Carlo (SMMC). For a recent review of AFMC in nuclei, see the chapter in Emergent Phenomena in Atomic Nuclei from Large-Scale Modeling (2017). SMMC has enabled CI shell model calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional diagonalization methods. While fermionic Monte Carlo methods are often limited by the so-called sign problem, the dominant components of nuclear interactions have a good sign in SMMC and often suffice for realistic calculations of statistical and collective properties. The smaller bad-sign components of the interactions can be treated by a method we introduced in PRL 1994.
The SMMC method is particularly useful for the microscopic calculation of statistical properties of nuclei, such as level densities. Level densities are required to estimate transition rates through Fermi’s golden rule and are important in the Hauser-Feshbach theory of compound-nucleus reactions, and therefore appear in numerous nuclear physics applications. However, their calculation in the presence of correlations is a difficult many-body problem. We have developed state-of-the-art SMMC methods for the microscopic calculation of level densities, including their dependence on good quantum numbers such as parity, spin and isospin (PRL 1997, PRL 1999, PRL 2000, PRL 2007, PRC Rapid Comm 2008, PRL 2008, PRL 2013, PRL 2014).
We have applied SMMC to nuclei as heavy as the lanthanides (PRL 2008, PRL 2013) and in ongoing work we are extending it to the actinides (arXiv:2509.26571), carrying out some of the largest SMMC calculations to date. Heavy nuclei exhibit various types of collectivity that are well described by phenomenological models, but a microscopic description has been lacking. Using SMMC, we have provided the first microscopic description of the crossover from vibrational to rotational collectivity (PRL 2013).
In a finite-size system such as the nucleus, it is crucial to use the canonical ensemble with a fixed number of particles. However, the projection on an odd number of particles leads to a new sign problem, which has prevented applications of SMMC to odd-mass and odd-odd nuclei. We developed a method that circumvents this odd-particle sign problem, enabling accurate calculations of the ground-state energy of odd-particle systems (PRL 2012). Using this method we calculated pairing gaps from odd-even mass differences and accurate level densities of odd-mass nuclei (PRC Rapid Comm. 2013). This method becomes time consuming in heavy nuclei and we have recently introduced a novel method, the partition function extrapolation method (PRC Letters 2024), that is suitable for extracting the ground-state energy of heavy odd-mass nuclei.
We have also used a configuration-space Monte Carlo (CSMC) method to solve exactly for the pairing Hamiltonian constructed from a density-dependent contact interaction in the framework of DFT (PRC 2011). This improved over the bare DFT+BCS results as measured by the rms deviation between the theoretical and experimental values of ~450 neutron pairing gaps.
Deformation is a key concept in our understanding of heavy nuclei. It arises naturally in the framework of the mean-field approximation, but it breaks the rotational invariance of the underlying nuclear Hamiltonian. We have introduced a novel method to study deformation in the rotationally invariant framework of the CI shell model (PRL 2014). We have used SMMC to calculate the axial quadrupole distribution in the laboratory frame, and showed that this distribution carries a model-independent signature of deformation. We have applied the method to isotope chains of lanthanides (PRC 2018). Using a Landau-like expansion of the logarithm of the quadrupole distribution in the so-called quadrupole invariants, we have determined its dependence on intrinsic deformation without invoking an intrinsic frame or a mean-field approximation. We can then calculate the dependence of statistical properties of nuclei on intrinsic deformation (PRC 2018b). This dependence is an important input to models of shape dynamics such as fission.
Most calculations of statistical properties of nuclei are based on mean-field approximations. However, their performance has not been tested against a theory that takes the full correlations into account. In PRC 2016 we benchmarked mean-field results against SMMC results. Finite-temperature mean-field theories are formulated in the grand-canonical ensemble, and it is necessary to project on fixed particle number. We introduced a novel projection formula in mean-field theories in which pairing correlations violate particle-number conservation (PRC 2017).
Using both mean-field and SMMC methods, we explained a long-standing problem in nuclear physics concerning the large enhancement (two orders of magnitude) of level densities in the crossover from closed-shell nuclei to mid-shell nuclei. Our SMMC level densities are in excellent agreement with recent experiments by the University of Oslo group for a chain of neodymium isotopes (PLB 2021).
Another central statistical property of nuclei is the gamma strength function, measuring the decay probability of the compound nucleus through the emission of gamma ray radiation. In SMMC, it is only possible to calculate the imaginary-time response function, whose inverse Laplace transform is the gamma strength function. However, this transform is numerically ill-defined. The standard method to carry out numerically the analytic continuation is the maximum-entropy method whose success depends crucially on a good choice of a prior strength function. We found out that the strength function in the static path plus random-phase approximation (SPA+RPA) (PRC Letters 2024) as well as in the SPA alone (PRC 2024, PRC 2025) constitute good prior strength functions for the maximum-entropy method.