We present two new methods for performing ab initio calculations of excited states for closed-shell systems within the in-medium similarity renormalization group (IMSRG) framework. Both are based on combining the IMSRG with simple many-body methods commonly used to target excited states, such as the Tamm-Dancoff approximation (TDA) and equations-of-motion (EOM) techniques. In the first approach, a two-step sequential IMSRG transformation is used to drive the Hamiltonian to a form where a simple TDA calculation (i.e., diagonalization in the space of $1\mathrm{p}1\mathrm{h}$ excitations) becomes exact for a subset of eigenvalues. In the second approach, EOM techniques are applied to the IMSRG ground-state-decoupled Hamiltonian to access excited states. We perform proof-of-principle calculations for parabolic quantum dots in two dimensions and the closed-shell nuclei $^{16}\mathrm{O}$ and $^{22}\mathrm{O}$. We find that the TDA-IMSRG approach gives better accuracy than the EOM-IMSRG when calculations converge, but it is otherwise lacking the versatility and numerical stability of the latter. Our calculated spectra are in reasonable agreement with analogous EOM-coupled-cluster calculations. This work paves the way for more interesting applications of the EOM-IMSRG approach to calculations of consistently evolved observables such as electromagnetic strength functions and nuclear matrix elements, and extensions to nuclei within one or two nucleons of a closed shell by generalizing the EOM ladder operator to include particle-number nonconserving terms.
We present two new methods for performing ab initio calculations of excited states for closed-shell systems within the in-medium similarity renormalization group (IMSRG) framework. Both are based on combining the IMSRG with simple many-body methods commonly used to target excited states, such as the Tamm-Dancoff approximation (TDA) and equations-of-motion (EOM) techniques. In the first approach, a two-step sequential IMSRG transformation is used to drive the Hamiltonian to a form where a simple TDA calculation (i.e., diagonalization in the space of 1p1h excitations) becomes exact for a subset of eigenvalues. In the second approach, EOM techniques are applied to the IMSRG ground-state-decoupled Hamiltonian to access excited states. We perform proof-of-principle calculations for parabolic quantum dots in two dimensions and the closed-shell nuclei 16O and 22O. We find that the TDA-IMSRG approach gives better accuracy than the EOM-IMSRG when calculations converge, but it is otherwise lacking the versatility and numerical stability of the latter. Our calculated spectra are in reasonable agreement with analogous EOM-coupled-cluster calculations. This work paves the way for more interesting applications of the EOM-IMSRG approach to calculations of consistently evolved observables such as electromagnetic strength functions and nuclear matrix elements, and extensions to nuclei within one or two nucleons of a closed shell by generalizing the EOM ladder operator to include particle-number nonconserving terms.
Abstract We present a pedagogical introduction to the in-medium similarity renormalization group (IMSRG) framework for ab initio calculations of nuclei. The IMSRG performs continuous unitary transformations of the nuclear many-body Hamiltonian in second-quantized form, which can be implemented with polynomial computational effort. Through suitably chosen generators, it is possible to extract eigenvalues of the Hamiltonian in a given nucleus, or drive the Hamiltonian matrix in configuration space to specific structures, e.g., band- or block-diagonal form. Exploiting this flexibility, we describe two complementary approaches for the description of closed- and open-shell nuclei: the first is the multireference IMSRG (MR-IMSRG), which is designed for the efficient calculation of nuclear ground-state properties. The second is the derivation of non-empirical valence-space interactions that can be used as input for nuclear shell model (i.e., configuration interaction (CI)) calculations. This IMSRG+shell model approach provides immediate access to excitation spectra, transitions, etc, but is limited in applicability by the factorial cost of the CI calculations. We review applications of the MR-IMSRG and IMSRG+shell model approaches to the calculation of ground-state properties for the oxygen, calcium, and nickel isotopic chains or the spectroscopy of nuclei in the lower sd shell, respectively, and present selected new results, e.g., for the ground- and excited state properties of neon isotopes.
We present a new ab initio method that uses similarity renormalization group (SRG) techniques to continuously diagonalize nuclear many-body Hamiltonians. In contrast with applications of the SRG to two- and three-nucleon interactions in free space, we perform the SRG evolution "in medium" directly in the A-body system of interest. The in-medium approach has the advantage that one can approximately evolve 3,…,A-body operators using only two-body machinery based on normal-ordering techniques. The method is nonperturbative and can be tailored to problems ranging from the diagonalization of closed-shell nuclei to the construction of effective valence-shell Hamiltonians and operators. We present first results for the energies of 4He, 16O, and 40Ca, which have accuracies comparable to coupled-cluster calculations.
Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.
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.
Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.
The goal of the present paper is twofold. First, a novel expansion many-body method applicable to superfluid open-shell nuclei, the so-called Bogoliubov in-medium similarity renormalization group (BIMSRG) theory, is formulated. This generalization of standard single-reference IMSRG theory for closed-shell systems parallels the recent extensions of coupled cluster, self-consistent Green's function or many-body perturbation theory. Within the realm of IMSRG theories, BIMSRG provides an interesting alternative to the already existing multi-reference IMSRG (MR-IMSRG) method applicable to open-shell nuclei. The algebraic equations for low-order approximations, i.e., BIMSRG(1) and BIMSRG(2), can be derived manually without much difficulty. However, such a methodology becomes already impractical and error prone for the derivation of the BIMSRG(3) equations, which are eventually needed to reach high accuracy. Based on a diagrammatic formulation of BIMSRG theory, the second objective of the present paper is thus to describe the third version (v3.0.0) of the ADG code that automatically (1) generates all valid BIMSRG(n) diagrams and (2) evaluates their algebraic expressions in a matter of seconds. This is achieved in such a way that equations can easily be retrieved for both the flow equation and the Magnus expansion formulations of BIMSRG. Expanding on this work, the first future objective is to numerically implement BIMSRG(2) (eventually BIMSRG(3)) equations and perform ab initio calculations of mid-mass open-shell nuclei.