In-medium similarity renormalization group for nuclei

A Similarity Renormalization Group Approach to Green's Function Methods

The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.

Ab initio predictions link the neutron skin of 208Pb to nuclear forces

Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 109 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape.

Ab Initio Structure Factors for Spin-Dependent Dark Matter Direct Detection

We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: ^{19}F, ^{23}Na, ^{27}Al, ^{29}Si, ^{73}Ge, ^{127}I, and ^{129,131}Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon currents at the one-body level, including effects from axial-vector two-body currents. We then apply the in-medium similarity renormalization group to construct effective valence-space Hamiltonians and consistently transformed operators of nuclear responses. Combining the recent advances of natural orbitals with three-nucleon forces expressed in large spaces, we obtain basis-space converged structure factors even in heavy nuclei. Generally results are consistent with previous calculations but large uncertainties in ^{127}I highlight the need for further study.

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg [J. Chem. Phys. 107, 432 (1997)]. This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple-ζ basis set predicts the spin splittings of [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ to be -35.7 and -17.1 kcal mol^-1, respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

Ab Initio Computation of Charge Densities for Sn and Xe Isotopes

We present the first ab initio calculations for open-shell nuclei past the tin isotopic line, focusing on Xe isotopes as well as doubly magic Sn isotopes. We show that, even for moderately hard interactions, it is possible to obtain meaningful predictions and that the NNLO_{sat} chiral interaction predicts radii and charge density distributions close to the experiment. We then make a new prediction for ^{100}Sn. This paves the way for ab initio studies of exotic charge density distributions at the limit of the present ab initio mass domain, where experimental data is becoming available. The present study closes the gap between the largest isotopes reachable by ab initio methods and the smallest exotic nuclei accessible to electron scattering experiments.

Unmixing Symmetries

The low-lying spectra of atomic nuclei display diverse behaviors, for example, rotational bands, which can be described phenomenologically by simple symmetry groups such as spatial SU(3). This leads to the idea of dynamical symmetry, where the Hamiltonian commutes with the Casimir operator(s) of a group, and is block diagonal in subspaces defined by the group's irreducible representations or irreps. Detailed microscopic calculations, however, show these symmetries are in fact often strongly mixed and the wave function fragmented across many irreps. More commonly, the fragmentation across members of a band are similar, which is called a quasidynamical symmetry. In this Letter I explicitly, albeit numerically, construct unitary transformations from a quasidynamical symmetry to a dynamical symmetry, adapting the similarity renormalization group (SRG) in order to transform away the symmetry-mixing parts of the Hamiltonian. The standard SRG produces unsatisfactory results, forcing the induced dynamical symmetry to be dominated by high-weight irreps irrespective of the original decomposition. Using spectral distribution theory to rederive and diagnose standard SRG, I introduce a new form of SRG. The new SRG transforms a quasidynamical symmetry to a dynamical symmetry, that is, unmixes the mixed symmetries, with intuitively more appealing results.

Physics of Nuclei: Key Role of an Emergent Symmetry

We show through first-principles nuclear structure calculations that the special nature of the strong nuclear force determines highly regular patterns heretofore unrecognized in nuclei that can be tied to an emergent approximate symmetry. This symmetry is ubiquitous and mathematically tracks with a symplectic symmetry group. This, in turn, has important implications for understanding the physics of nuclei: we find that nuclei are made of only a few equilibrium shapes, deformed or not, with associated vibrations and rotations. It also opens the path for ab initio large-scale modeling of open-shell intermediate-mass nuclei without the need for renormalized interactions and effective charges.

Multireference Theories of Electron Correlation Based on the Driven Similarity Renormalization Group

The driven similarity renormalization group (DSRG) provides an alternative way to address the intruder state problem in quantum chemistry. In this review, we discuss recent developments of multireference methods based on the DSRG. We provide a pedagogical introduction to the DSRG and its various extensions and discuss its formal properties in great detail. In addition, we report several illustrative applications of the DSRG to molecular systems.

78Ni revealed as a doubly magic stronghold against nuclear deformation

Nuclear magic numbers correspond to fully occupied energy shells of protons or neutrons inside atomic nuclei. Doubly magic nuclei, with magic numbers for both protons and neutrons, are spherical and extremely rare across the nuclear landscape. Although the sequence of magic numbers is well established for stable nuclei, experimental evidence has revealed modifications for nuclei with a large asymmetry between proton and neutron numbers. Here we provide a spectroscopic study of the doubly magic nucleus ⁷⁸Ni, which contains fourteen neutrons more than the heaviest stable nickel isotope. We provide direct evidence of its doubly magic nature, which is also predicted by ab initio calculations based on chiral effective-field theory interactions and the quasi-particle random-phase approximation. Our results also indicate the breakdown of the neutron magic number 50 and proton magic number 28 beyond this stronghold, caused by a competing deformed structure. State-of-the-art phenomenological shell-model calculations reproduce this shape coexistence, predicting a rapid transition from spherical to deformed ground states, with ⁷⁸Ni as the turning point.

Structure of the Lightest Tin Isotopes

We link the structure of nuclei around ^{100}Sn, the heaviest doubly magic nucleus with equal neutron and proton numbers (N=Z=50), to nucleon-nucleon (NN) and three-nucleon (NNN) forces constrained by data of few-nucleon systems. Our results indicate that ^{100}Sn is doubly magic, and we predict its quadrupole collectivity. We present precise computations of ^{101}Sn based on three-particle-two-hole excitations of ^{100}Sn, and we find that one interaction accurately reproduces the small splitting between the lowest J^{π}=7/2^{+} and 5/2^{+} states.

Ab Initio Description of Open-Shell Nuclei: Merging No-Core Shell Model and In-Medium Similarity Renormalization Group

We merge two successful ab initio nuclear-structure methods, the no-core shell model (NCSM) and the multireference in-medium similarity renormalization group (IM-SRG) to define a new many-body approach for the comprehensive description of ground and excited states of closed and open-shell nuclei. Building on the key advantages of the two methods-the decoupling of excitations at the many-body level in the IM-SRG and the access to arbitrary nuclei, eigenstates, and observables in the NCSM-their combination enables fully converged no-core calculations for an unprecedented range of nuclei and observables at moderate computational cost. We present applications in the carbon and oxygen isotopic chains, where conventional NCSM calculations are still feasible and provide an important benchmark. The efficiency and rapid convergence of the new approach make it ideally suited for ab initio studies of the complete spectroscopy of nuclei up into the medium-mass regime.

An integral-factorized implementation of the driven similarity renormalization group second-order multireference perturbation theory

We report an efficient implementation of a second-order multireference perturbation theory based on the driven similarity renormalization group (DSRG-MRPT2) [C. Li and F. A. Evangelista, J. Chem. Theory Comput. 11, 2097 (2015)]. Our implementation employs factorized two-electron integrals to avoid storage of large four-index intermediates. It also exploits the block structure of the reference density matrices to reduce the computational cost to that of second-order Møller-Plesset perturbation theory. Our new DSRG-MRPT2 implementation is benchmarked on ten naphthyne isomers using basis sets up to quintuple-ζ quality. We find that the singlet-triplet splittings (ΔST) of the naphthyne isomers strongly depend on the equilibrium structures. For a consistent set of geometries, the ΔST values predicted by the DSRG-MRPT2 are in good agreements with those computed by the reduced multireference coupled cluster theory with singles, doubles, and perturbative triples.

Nucleus-Dependent Valence-Space Approach to Nuclear Structure

We present a nucleus-dependent valence-space approach for calculating ground and excited states of nuclei, which generalizes the shell-model in-medium similarity renormalization group to an ensemble reference with fractionally filled orbitals. Because the ensemble is used only as a reference, and not to represent physical states, no symmetry restoration is required. This allows us to capture three-nucleon (3N) forces among valence nucleons with a valence-space Hamiltonian specifically targeted to each nucleus of interest. Predicted ground-state energies from carbon through nickel agree with results of other large-space ab initio methods, generally to the 1% level. In addition, we show that this new approach is required in order to obtain convergence for nuclei in the upper p and sd shells. Finally, we address the 1^{+}/3^{+} inversion problem in ^{22}Na and ^{46}V. This approach extends the reach of ab initio nuclear structure calculations to essentially all light- and medium-mass nuclei.

Towards numerically robust multireference theories: The driven similarity renormalization group truncated to one- and two-body operators

The first nonperturbative version of the multireference driven similarity renormalization group (MR-DSRG) theory [C. Li and F. A. Evangelista, J. Chem. Theory Comput. 11, 2097 (2015)] is introduced. The renormalization group structure of the MR-DSRG equations ensures numerical robustness and avoidance of the intruder-state problem, while the connected nature of the amplitude and energy equations guarantees size consistency and extensivity. We approximate the MR-DSRG equations by keeping only one- and two-body operators and using a linearized recursive commutator approximation of the Baker-Campbell-Hausdorff expansion [T. Yanai and G. K.-L. Chan, J. Chem. Phys. 124, 194106 (2006)]. The resulting linearized MR-DSRG scheme with one- and two-body operators [MR-LDSRG(2)] contains only 39 terms and scales as O(N(2)NP (2)NH (2)) where NH, NP, and N correspond to the number of hole, particle, and total orbitals, respectively. Benchmark MR-LDSRG(2) computations on the hydrogen fluoride and molecular nitrogen binding curves and the singlet-triplet splitting of p-benzyne yield results comparable in accuracy to those from multireference configuration interaction, Mukherjee multireference coupled cluster theory, and internally contracted multireference coupled cluster theory.

Radii and Binding Energies in Oxygen Isotopes: A Challenge for Nuclear Forces

We present a systematic study of both nuclear radii and binding energies in (even) oxygen isotopes from the valley of stability to the neutron drip line. Both charge and matter radii are compared to state-of-the-art ab initio calculations along with binding energy systematics. Experimental matter radii are obtained through a complete evaluation of the available elastic proton scattering data of oxygen isotopes. We show that, in spite of a good reproduction of binding energies, ab initio calculations with conventional nuclear interactions derived within chiral effective field theory fail to provide a realistic description of charge and matter radii. A novel version of two- and three-nucleon forces leads to considerable improvement of the simultaneous description of the three observables for stable isotopes but shows deficiencies for the most neutron-rich systems. Thus, crucial challenges related to the development of nuclear interactions remain.

Ab initio description of p-shell hypernuclei

We present the first ab initio calculations for p-shell single-Λ hypernuclei. For the solution of the many-baryon problem, we develop two variants of the no-core shell model with explicit Λ and Σ(+),Σ(0),Σ(-) hyperons including Λ-Σ conversion, optionally supplemented by a similarity renormalization group transformation to accelerate model-space convergence. In addition to state-of-the-art chiral two- and three-nucleon interactions, we use leading-order chiral hyperon-nucleon interactions and a recent meson-exchange hyperon-nucleon interaction. We validate the approach for s-shell hypernuclei and apply it to p-shell hypernuclei, in particular to (Λ)(7)Li, (Λ)(9)Be, and (Λ)(13)C. We show that the chiral hyperon-nucleon interactions provide ground-state and excitation energies that generally agree with experiment within the cutoff dependence. At the same time we demonstrate that hypernuclear spectroscopy provides tight constraints on the hyperon-nucleon interactions.

Auxiliary-field quantum Monte Carlo simulations of neutron matter in chiral effective field theory

We present variational Monte Carlo calculations of the neutron matter equation of state using chiral nuclear forces. The ground-state wave function of neutron matter, containing nonperturbative many-body correlations, is obtained from auxiliary-field quantum Monte Carlo simulations of up to about 340 neutrons interacting on a 10(3) discretized lattice. The evolution Hamiltonian is chosen to be attractive and spin independent in order to avoid the fermion sign problem and is constructed to best reproduce broad features of the chiral nuclear force. This is facilitated by choosing a lattice spacing of 1.5 fm, corresponding to a momentum-space cutoff of Λ=414  MeV/c, a resolution scale at which strongly repulsive features of nuclear two-body forces are suppressed. Differences between the evolution potential and the full chiral nuclear interaction (Entem and Machleidt Λ=414  MeV [L. Coraggio et al., Phys. Rev. C 87, 014322 (2013).

Nonperturbative shell-model interactions from the in-medium similarity renormalization group

We present the first ab initio construction of valence-space Hamiltonians for medium-mass nuclei based on chiral two- and three-nucleon interactions using the in-medium similarity renormalization group. When applied to the oxygen isotopes, we find experimental ground-state energies are well reproduced, including the flat trend beyond the drip line at (24)O. Similarly, natural-parity spectra in (21,22,23,24)O are in agreement with experiment, and we present predictions for excited states in (25,26)O. The results exhibit a weak dependence on the harmonic-oscillator basis parameter and reproduce spectroscopy within the standard sd valence space.

Coupled-cluster computations of atomic nuclei

In the past decade, coupled-cluster theory has seen a renaissance in nuclear physics, with computations of neutron-rich and medium-mass nuclei. The method is efficient for nuclei with product-state references, and it describes many aspects of weakly bound and unbound nuclei. This report reviews the technical and conceptual developments of this method in nuclear physics, and the results of coupled-cluster calculations for nucleonic matter, and for exotic isotopes of helium, oxygen, calcium, and some of their neighbors.

New applications of renormalization group methods in nuclear physics

We review recent developments in the use of renormalization group (RG) methods in low-energy nuclear physics. These advances include enhanced RG technology, particularly for three-nucleon forces, which greatly extends the reach and accuracy of microscopic calculations. We discuss new results for the nucleonic equation of state with applications to astrophysical systems such as neutron stars, new calculations of the structure and reactions of finite nuclei, and new explorations of correlations in nuclear systems.

Isotopic chains around oxygen from evolved chiral two- and three-nucleon interactions

We extend the formalism of self-consistent Green's function theory to include three-body interactions and apply it to isotopic chains around oxygen for the first time. The third-order algebraic diagrammatic construction equations for two-body Hamiltonians can be exploited upon defining system-dependent one- and two-body interactions coming from the three-body force, and, correspondingly, dropping interaction-reducible diagrams. The Koltun sum rule for the total binding energy acquires a correction due to the added three-body interaction. This formalism is then applied to study chiral two- and three-nucleon forces evolved to low momentum cutoffs. The binding energies of nitrogen, oxygen, and fluorine isotopes are reproduced with good accuracy and demonstrate the predictive power of this approach. Leading order three-nucleon forces consistently bring results close to the experiment for all neutron rich isotopes considered and reproduce the correct driplines for oxygen and nitrogen. The formalism introduced also allows us to calculate form factors for nucleon transfer on doubly magic systems.

Ab Initio calculations of even oxygen isotopes with chiral two-plus-three-nucleon interactions

We formulate the in-medium similarity renormalization group (IM-SRG) for open-shell nuclei using a multireference formalism based on a generalized Wick theorem introduced in quantum chemistry. The resulting multireference IM-SRG (MR-IM-SRG) is used to perform the first ab initio study of all even oxygen isotopes with chiral nucleon-nucleon and three-nucleon interactions, from the proton to the neutron drip lines. We obtain an excellent reproduction of experimental ground-state energies with quantified uncertainties, which is validated by results from the importance-truncated no-core shell model and the coupled cluster method. The agreement between conceptually different many-body approaches and experiment highlights the predictive power of current chiral two- and three-nucleon interactions, and establishes the MR-IM-SRG as a promising new tool for ab initio calculations of medium-mass nuclei far from shell closures.