skip to main content


Title: Nonempirical Interactions for the Nuclear Shell Model: An Update
The nuclear shell model has perhaps been the most important conceptual and computational paradigm for the understanding of the structure of atomic nuclei. While the shell model has been used predominantly in a phenomenological context, there have been efforts stretching back more than half a century to derive shell model parameters based on a realistic interaction between nucleons. More recently, several ab initio many-body methods—in particular, many-body perturbation theory, the no-core shell model, the in-medium similarity renormalization group, and coupled-cluster theory—have developed the capability to provide effective shell model Hamiltonians. We provide an update on the status of these methods and investigate the connections between them and their potential strengths and weaknesses, with a particular focus on the in-medium similarity renormalization group approach. Three-body forces are demonstrated to be important for understanding the modifications needed in phenomenological treatments. We then review some applications of these methods to comparisons with recent experimental measurements, and conclude with some remaining challenges in ab initio shell model theory.  more » « less
Award ID(s):
1713901
NSF-PAR ID:
10179961
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Annual Review of Nuclear and Particle Science
Volume:
69
Issue:
1
ISSN:
0163-8998
Page Range / eLocation ID:
307 to 362
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We evaluate the allowed $\beta^-$-decay properties of nuclei with $Z = 8$–$15$ systematically under the framework of the nuclear shell model using the valence space Hamiltonians derived from modern ab initio methods, such as in-medium similarity renormalization group and coupled-cluster theory. For comparison we also show results obtained with fitted interaction derived from chiral effective field theory and phenomenological universal $sd$-shell Hamiltonian version B interaction. We have performed calculations for O $\rightarrow$ F, F $\rightarrow$ Ne, Ne $\rightarrow$ Na, Na $\rightarrow$ Mg, Mg $\rightarrow$ Al, Al $\rightarrow$ Si, Si $\rightarrow$ P, and P $\rightarrow$ S transitions. Theoretical results for $B(GT)$, $\log ft$ values, and half-lives are discussed and compared with the available experimental data. 
    more » « less
  2. null (Ed.)
    Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitly treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born−Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron’s molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron’s chemical and physical behavior. 
    more » « less
  3. The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron’s vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron’s central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object. 
    more » « less
  4. The influence of functional traits on species survivorship has been evaluated in various contexts in both modern and ancient ecosystems, but an important direction for research is to integrate datasets that include both extinct and extant taxa. This approach can provide a more reliable understanding of the effects of functional traits on macroecological and macroevolutionary dynamics. Knowledge of the links between individual traits and survivorship is crucial for developing accurate extinction risk predictive models. Here we test the impact of numerous functional traits on the survival and extinction of species through time, using bivalve and gastropod species from the rich fossil record of the western Atlantic over the last ~3 million years, along with the associated extant biota. We also compare the impact of these organismic traits on survival relative to a group level trait: geographic distribution. Analyses use a dataset of 12 functional traits including life habit, feeding behavior and basal metabolic rate (BMR), for 115 species from 36 families. Traits were observed and measured from specimens in the collections of the Paleontological Research Institution, Florida Museum of Natural History, and University of Kansas, as well as surveys of the literature and online databases such as the Neogene Marine Biota of Tropical America (NMITA). Results derived from Principal Coordinates Analysis (PCoA) show there is a clear distinction between extinct and extant species, overall, when comparing them based on life habit, maximum body size, shell composition and BMR. Most traits showed little direct relation with survival, except BMR and associated maximum body size, supporting the Metabolic Theory of Ecology. Since many functional traits do not explain survival, their function may be mis- or over-interpreted, and traits posited to represent important organismic adaptations may not play a prominent role in long-term species survival, especially during the major climate changes over the last ~ 3 million years. Some traits do show significant interactions, and these were more fully explored using additional multivariate analyses. The relative importance of geographic range size suggests group-level characters may be the primary determinant of extinction patterns over macroevolutionary time scales. 
    more » « less
  5. Abstract Simulating quantum systems in a finite volume is a powerful theoretical tool to extract information about them. Real-world properties of the system are encoded in how its discrete energy levels change with the size of the volume. This approach is relevant not only for nuclear physics, where lattice methods for few- and many-nucleon states complement phenomenological shell-model descriptions and ab initio calculations of atomic nuclei based on harmonic oscillator expansions, but also for other fields such as simulations of cold atomic systems. This contribution presents recent progress concerning finite-volume simulations of few-body systems. In particular, it discusses details regarding the efficient numerical implementation of separable interactions and it presents eigenvector continuation as a method for performing robust and efficient volume extrapolations. 
    more » « less