BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
H. Hergert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $β$ decay.
We use chiral interactions and several {\em ab initio} methods to compute the nuclear matrix elements (NMEs) for ground-state to ground-state double Gamow-Teller transitions in a range of isotopes, and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the {\em ab initio} calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double-beta decay.
Recently, ab initio techniques have been successfully connected to the traditional valence-space shell model. In doing so, they can either explicitly provide ab initio shell-model effective Hamiltonians or constrain the construction of empirical ones. In the present work, the possibility to follow a similar path for the nuclear energy density functional (EDF) method is analyzed. For this connection to be actualized, two theoretical techniques are instrumental: the recently proposed ab initio PGCM-PT many-body formalism and the MR-IMSRG pre-processing of the nuclear Hamiltonian. Based on both formal arguments and numerical results, possible new lines of research are briefly discussed, namely to compute ab initio EDF effective Hamiltonians at low computational cost, to constrain empirical ones or to produce them directly via an effective field theory that remains to be invented.
We report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-beta ($0νββ$) decays in light nuclei with mass number ranging from $A=6$ to $A=22$. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two-plus-three-nucleon interaction from chiral effective field theory is employed, and both isospin-conserving ($ΔT=0$) and isospin-changing ($ΔT=2$) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of $ΔT=0$ transitions are in good agreement among all calculations, at the level of 10%. For $ΔT=2$, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps towards improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross-checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the $0νββ$ decays in heavier candidate nuclei, where quasi-exact benchmarks are not available.
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.
Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. We use our scheme to compute the matrix element that governs the neutrinoless double beta decay of $^{48}$Ca to $^{48}$Ti, and find it to have the value $0.61$, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as $^{76}$Ge, $^{130}$Te, and $^{136}$Xe.
We compute nuclear matrix elements of neutrinoless double-beta decay mediated by light Majorana-neutrino exchange in the A = 6 system. The goal is to benchmark two many-body approaches, the No-Core Shell Model and the Multi-Reference In-Medium Similarity Renormalization Group. We use the SRG-evolved chiral N3LO-EM500 potential for the nuclear interaction, and make the approximation that isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-beta decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei.
We present ab initio predictions for ground and excited states of doubly open-shell fluorine and neon isotopes based on chiral two- and three-nucleon interactions. We use the in-medium similarity renormalization group, in both flow-equation and Magnus formulations, to derive mass-dependent sd valence-space Hamiltonians. The experimental ground-state energies are reproduced through neutron number N=14, beyond which a new targeted normal-ordering procedure improves agreement with data and large-scale multi-reference calculations. For spectroscopy, we focus on neutron-rich 23-26F and 24-26Ne isotopes near N=14,16 magic numbers. In all cases we find an agreement with experiment competitive with established phenomenology. Moreover, yrast states are well described in 20Ne and 24Mg, providing an ab initio description of deformation in the medium-mass region.
We present a pedagogical introduction to the In-Medium Similarity Renormalization Group (IM-SRG) framework for ab initio calculations of nuclei. The IM-SRG 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 IM-SRG (MR-IM-SRG), which is designed for the efficient calculation of nuclear ground-state properties. The second is the derivation of nonempirical valence-space interactions that can be used as input for nuclear Shell model (i.e., configuration interaction (CI)) calculations. This IM-SRG+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-IM-SRG and IM-SRG+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 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 {\it 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, {\it 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.
We study theoretically the electric dipole transitions of even Ni isotopes at low energies, using the self-consistent quasi-particle random-phase approximation (RPA) with the D1S Gogny interaction and a continuum-RPA model with the SLy4 Skyrme force. We analyze isoscalar states, isovector states, and the dipole polarizability. We define a reference value for the polarizability, to remove a trivial dependence on the mass number. We compare our results with data and other calculations, with a focus on collective states, shell effects, and threshold transitions. Our results support the presence of a strong isoscalar transition, with little or moderate E1 strength, as a universal feature of ordinary nuclei. In moderately neutron-rich Ni isotopes, namely 68Ni and neighboring isotopes, this transition is found bimodal due to couplings with surface neutrons. An adequate treatment of the continuum states appears essential for describing suprathreshold E1 strength, especially beyond 68Ni. Very exotic isotopes (N>50) are found highly polarizable, with practically all their E1 strength in the continuum. The dipole polarizability and the neutron skin thickness are influenced by shell structure in different ways, so they can appear anticorrelated. A comparison with existing results for lighter (Ca) and heavier (Sn) nuclei suggests that the so-called pygmy dipole strength is influenced strongly by shell effects and that, partly for that reason, its isospin structure depends on the mass region.
The concept of single-nucleon shells constitutes a basic pillar of our understanding of nuclear structure. Effective single-particle energies (ESPEs) introduced by French and Baranger represent the most appropriate tool to relate many-body observables to a single-nucleon shell structure. As briefly discussed in [T. Duguet, G. Hagen, Phys. Rev. C {\bf 85}, 034330 (2012)], the dependence of ESPEs on one-nucleon transfer probability matrices makes them purely theoretical quantities that "run" with the non-observable resolution scale $λ$ employed in the calculation. Given that ESPEs provide a way to interpret the many-body problem in terms of simpler theoretical ingredients, the goal is to specify the terms, i.e. the exact sense and conditions, in which this interpretation can be conducted meaningfully. State-of-the-art multi-reference in-medium similarity renormalization group and self-consistent Gorkov Green's function many-body calculations are employed to corroborate the formal analysis. This is done by comparing the behavior of several observables and of non-observable ESPEs (and spectroscopic factors) under (quasi) unitary similarity renormalization group transformations of the Hamiltonian parameterized by the resolution scale $λ$. The non-observable nature of the nuclear shell structure, i.e. the fact that it constitutes an intrinsically theoretical object with no counterpart in the empirical world, must be recognized and assimilated. Eventually, practitioners can refer to nuclear shells and spectroscopic factors in their analyses of nuclear phenomena if, and only if, they use consistent structure and reaction theoretical schemes based on a fixed resolution scale they have agreed on prior to performing their analysis and comparisons.
The structure of the $^{24}$F nucleus has been studied at GANIL using the $β$ decay of $^{24}$O and the in-beam $γ$-ray spectroscopy from the fragmentation of projectile nuclei. Combining these complementary experimental techniques, the level scheme of $^{24}$F has been constructed up to 3.6 Mev by means of particle-$γ$ and particle-$γγ$ coincidence relations. Experimental results are compared to shell-model calculations using the standard USDA and USDB interactions as well as ab-initio valence-space Hamiltonians calculated from the in-medium similarity renormalization group based on chiral two- and three-nucleon forces. Both methods reproduce the measured level spacings well, and this close agreement allows unidentified spins and parities to be consistently assigned.
We present a comprehensive review of the In-Medium Similarity Renormalization Group (IM-SRG), a novel ab inito method for nuclei. The IM-SRG employs a continuous unitary transformation of the many-body Hamiltonian to decouple the ground state from all excitations, thereby solving the many-body problem. Starting from a pedagogical introduction of the underlying concepts, the IM-SRG flow equations are developed for systems with and without explicit spherical symmetry. We study different IM-SRG generators that achieve the desired decoupling, and how they affect the details of the IM-SRG flow. Based on calculations of closed-shell nuclei, we assess possible truncations for closing the system of flow equations in practical applications, as well as choices of the reference state. We discuss the issue of center-of-mass factorization and demonstrate that the IM-SRG ground-state wave function exhibits an approximate decoupling of intrinsic and center-of-mass degrees of freedom, similar to Coupled Cluster (CC) wave functions. To put the IM-SRG in context with other many-body methods, in particular many-body perturbation theory and non-perturbative approaches like CC, a detailed perturbative analysis of the IM-SRG flow equations is carried out. We conclude with a discussion of ongoing developments, including IM-SRG calculations with three-nucleon forces, the multi-reference IM-SRG for open-shell nuclei, first non-perturbative derivations of shell- model interactions, and the consistent evolution of operators in the IM-SRG. We dedicate this review to the memory of Gerry Brown, one of the pioneers of many-body calculations of nuclei.
We use the newly developed Multi-Reference In-Medium Similarity Renormalization Group to study all even isotopes of the calcium and nickel isotopic chains, based on two- plus three-nucleon interactions derived from chiral effective field theory. We present results for ground-state and two-neutron separation energies and quantify their theoretical uncertainties. At shell closures, we find excellent agreement with Coupled Cluster results obtained with the same Hamiltonians. Our results highlight the importance of the chiral 3N interaction to obtain a correct reproduction of experimental energy trends, and their subtle impact on the location of the neutron drip lines in the Ca and Ni chains. At the same time, we uncover and discuss deficiencies of the input Hamiltonians which need to be addressed by the next generation of chiral interactions.
Starting from the matrix elements of the nucleon-nucleon interaction in momentum space we present a method to derive an operator representation with a minimal set of operators that is required to provide an optimal description of the partial waves with low angular momentum. As a first application we use this method to obtain an operator representation for the Argonne potential transformed by means of the unitary correlation operator method and discuss the necessity of including momentum dependent operators. The resulting operator representation leads to the same results as the original momentum space matrix elements when applied to the two-nucleon system and various light nuclei. For applications in fermionic and antisymmetrized molecular dynamics, where an operator representation of a soft but realistic effective interaction is indispensable, a simplified version using a reduced set of operators is given.
Two different experimental approaches were combined to study the electric dipole strength in the doubly-magic nucleus 48Ca below the neutron threshold. Real-photon scattering experiments using bremsstrahlung up to 9.9 MeV and nearly mono-energetic linearly polarized photons with energies between 6.6 and 9.51 MeV provided strength distribution and parities, and an (α,α'γ) experiment at E_α=136 MeV gave cross sections for an isoscalar probe. The unexpected difference observed in the dipole response is compared to calculations using the first-order random-phase approximation and points to an energy-dependent isospin character. A strong isoscalar state at 7.6 MeV was identified for the first time supporting a recent theoretical prediction.
We study the low-energy dipole (LED) strength distribution along the Sn isotopic chain in both the isoscalar (IS) and the isovector (IV, or E1) electric channels, to provide testable predictions and guidance for new experiments with stable targets and radioactive beams. We use the self-consistent Quasi-particle Random-Phase Approximation (QRPA) with finite-range interactions and mainly the Gogny D1S force. We analyze also the performance of a realistic two-body interaction supplemented by a phenomenological three-body contact term. We find that from N=50 and up to the N=82 shell closure (132Sn) the lowest-energy part of the IS-LED spectrum is dominated by a collective transition whose properties vary smoothly with neutron number and which cannot be interpreted as a neutron-skin oscillation. For the neutron-rich species this state contributes to the E1 strength below particle threshold, but much more E1 strength is carried by other, weak but numerous transitions around or above threshold. We find that strong structural changes in the spectrum take effect beyond N=82, namely increased LED strength and lower excitation energies. Our results with the Gogny interaction are compatible with existing data. On this basis we predict that a) the summed IS strength below particle threshold shall be of the same order of magnitude for N=50-82, b) the summed E1 strength up to approximately 12 MeV shall be similar for N=50-82 MeV, while c) the summed E1 strength below threshold shall be of the same order of magnitude for N ~ 64 - 82 and much weaker for the lighter, more-symmetric isotopes. We point out a general agreement of our results with other non-relativistic studies, the absence of a collective IS mode in some of those studies, and a possibly radical disagreement with relativistic models.
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 24O. Similarly, natural-parity spectra in 21,22,23,24O are in agreement with experiment, and we present predictions for excited states in 25,26O. The results exhibit a weak dependence on the harmonic-oscillator basis parameter and reproduce spectroscopy within the standard sd valence space.
We formulate the In-Medium Similarity Renormalization Group (IM-SRG) for open-shell nuclei using a multi-reference formalism based on a generalized Wick theorem introduced in quantum chemistry. The resulting multi-reference IM-SRG (MR-IM-SRG) is used to perform the first ab initio study of even oxygen isotopes with chiral NN and 3N Hamiltonians, 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.
We use the recently proposed In-Medium Similarity Renormalization Group (IM-SRG) to carry out a systematic study of closed-shell nuclei up to $\nuc{Ni}{56}$, based on chiral two- plus three-nucleon interactions. We analyze the capabilities of the IM-SRG by comparing our results for the ground-state energy to Coupled Cluster calculations, as well as to quasi-exact results from the Importance-Truncated No-Core Shell Model. Using chiral two- plus three-nucleon Hamiltonians whose resolution scales are lowered by free-space SRG evolution, we obtain good agreement with experimental binding energies in $\nuc{He}{4}$ and the closed-shell oxygen isotopes, while the calcium and nickel isotopes are somewhat overbound.
Recent theoretical work has not led to a consensus regarding the nature of the low-energy E1 strength in the 40,44,48Ca isotopes, for which high-resolution (gamma,gamma') data exist. Here we revisit this problem using the first-order quasi-particle random-phase approximation (QRPA) and different interactions. First we examine all even Ca isotopes with N=14-40. All isotopes are predicted to undergo dipole transitions at low energy, of large and comparable isoscalar strength but of varying E1 strength. Provided a moderate and uniform energetic shift is introduced to the results, QRPA with the Gogny D1S interaction is able to account for the (gamma,gamma') data, because, up to N=28, it yields a rather pure isoscalar oscillation. A neutron-skin oscillation is anticipated for N larger or equal to 30. This contradicts existing predictions that 44,48Ca develop a neutron-skin mode. Which theoretical result is correct cannot be resolved conclusively using the available data. We propose that alpha-scattering, possibly followed by an electroexcitation experiment, could resolve the situation and thereby help to improve the different models aspiring to describe reliably the low-energy dipole strength of nuclei.
We have developed a fully consistent framework for calculations in the Quasiparticle Random Phase Approximation (QRPA) with $NN$ interactions from the Similarity Renormalization Group (SRG) and other unitary transformations of realistic interactions. The consistency of our calculations, which use the same Hamiltonian to determine the Hartree-Fock-Bogoliubov (HFB) ground states and the residual interaction for QRPA, guarantees an excellent decoupling of spurious strength, without the need for empirical corrections. While work is under way to include SRG-evolved 3N interactions, we presently account for some 3N effects by means of a linearly density-dependent interaction, whose strength is adjusted to reproduce the charge radii of closed-shell nuclei across the whole nuclear chart. As a first application, we perform a survey of the monopole, dipole, and quadrupole response of the calcium isotopic chain and of the underlying single-particle spectra, focusing on how their properties depend on the SRG parameter $λ$. Unrealistic spin-orbit splittings suggest that spin-orbit terms from the 3N interaction are called for. Nevertheless, our general findings are comparable to results from phenomenological QRPA calculations using Skyrme or Gogny energy density functionals. Potentially interesting phenomena related to low-lying strength warrant more systematic investigations in the future.
Microscopic input to a universal nuclear energy density functional can be provided through the density matrix expansion (DME), which has recently been revived and improved. Several DME implementation strategies are tested for neutron drop systems in harmonic traps by comparing to Hartree-Fock (HF) and ab initio no-core full configuration (NCFC) calculations with a model interaction (Minnesota potential). The new DME with exact treatment of Hartree contributions is found to best reproduce HF results and supplementing the functional with fit Skyrme-like contact terms shows systematic improvement toward the full NCFC results.
We investigate the influence of phenomenological three-nucleon interactions on the systematics of ground-state energies and charge radii throughout the whole nuclear mass range from 4-He to 208-Pb. The three-nucleon interactions supplement unitarily transformed two-body interactions constructed within the Unitary Correlation Operator Method or the Similarity Renormalization Group approach. To be able to address heavy nuclei as well, we treat the many-body problem in Hartree-Fock plus many-body perturbation theory, which is sufficient to assess the systematics of energies and radii, and limit ourselves to regularized three-body contact interactions. We show that even with such a simplistic three-nucleon interaction a simultaneous reproduction of the experimental ground-state energies and charge radii can be achieved, which is not possible with unitarily transformed two-body interactions alone.