Optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order

Improved Tensor Current Limit from ^{8}B β Decay Including New Recoil-Order Calculations

A precision measurement of the β^{+} decay of ^{8}B was performed using the Beta-decay Paul Trap to determine the β-ν angular correlation coefficient a_{βν}. The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, a_{βν}=-0.3345±0.0019_{stat}±0.0021_{syst}. This value agrees with the standard model value of -1/3 and improves uncertainties in ^{8}B by nearly a factor of 2. By combining results from ^{8}B and ^{8}Li, a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision β decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3σ. The present results are of comparable sensitivity and do not support this finding.

Description of the Proton-Decaying 0_{2}^{+} Resonance of the α Particle

The recent precise experimental determination of the monopole transition form factor from the ground state of ^{4}He to its 0_{2}^{+} resonance via electron scattering has reinvigorated discussions about the nature of this first excited state of the α particle. The 0_{2}^{+} state has been traditionally interpreted in the literature as the isoscalar monopole resonance (breathing mode) or, alternatively, as a particle-hole shell-model excitation. To better understand the nature of this state, which lies only ∼410  keV above the proton emission threshold, we employ the coupled-channel representation of the no-core Gamow shell model. By considering the [^{3}H+p], [^{3}He+n], and [^{2}H+^{2}H] reaction channels, we explain the excitation energy and monopole form factor of the 0_{2}^{+} state. We argue that the continuum coupling strongly impacts the nature of this state, which carries characteristics of the proton decay threshold.

NuHamil : A numerical code to generate nuclear two- and three-body matrix elements from chiral effective field theory

The applicability of nuclear ab initio calculations has rapidly extended over the past decades. However, starting research projects is still challenging due to the required numerical expertise in the generation of underlying nuclear interaction matrix elements and many-body calculations. To ease the first issue, in this paper we introduce the numerical code NuHamil to generate the nucleon-nucleon (NN) and three-nucleon (3N) matrix elements expressed in a spherical harmonic-oscillator basis, inputs of many-body calculations. The ground-state energies for the selected doubly closed shell nuclei are calculated with the no-core shell-model (NCSM) and in-medium similarity renormalization group (IMSRG). The code is written in modern Fortran, and OpenMP+MPI hybrid parallelization is available for the 3N matrix-element calculations.

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

Proton Distribution Radii of ^{16-24}O: Signatures of New Shell Closures and Neutron Skin

The root mean square radii of the proton density distribution in ^{16-24}O derived from measurements of charge changing cross sections with a carbon target at ∼900A  MeV together with the matter radii portray thick neutron skin for ^{22-24}O despite ^{22,24}O being doubly magic. Imprints of the shell closures at N=14 and 16 are reflected in local minima of their proton radii that provide evidence for the tensor interaction causing them. The radii agree with ab initio calculations employing the chiral NNLO_{sat} interaction, though skin thickness predictions are challenged. Shell model predictions agree well with the data.

Impact of Clustering on the ^{8}Li β Decay and Recoil Form Factors

We place unprecedented constraints on recoil corrections in the β decay of ^{8}Li, by identifying a strong correlation between them and the ^{8}Li ground state quadrupole moment in large-scale ab initio calculations. The results are essential for improving the sensitivity of high-precision experiments that probe the weak interaction theory and test physics beyond the standard model. In addition, our calculations predict a 2^{+} state of the α+α system that is energetically accessible to β decay but has not been observed in the experimental ^{8}Be energy spectrum, and has an important effect on the recoil corrections and β decay for the A=8 systems. This state and an associated 0^{+} state are notoriously difficult to model due to their cluster structure and collective correlations, but become feasible for calculations in the ab initio symmetry-adapted no-core shell-model framework.

Recent Progress in Gamow Shell Model Calculations of Drip Line Nuclei

The Gamow shell model (GSM) is a powerful method for the description of the exotic properties of drip line nuclei. Internucleon correlations are included via a configuration interaction framework. Continuum coupling is directly included at basis level by using the Berggren basis, in which, bound, resonance, and continuum single-particle states are treated on an equal footing in the complex momentum plane. Two different types of Gamow shell models have been developed: its first embodiment is that of the GSM defined with phenomenological nuclear interactions, whereas the GSM using realistic nuclear interactions, called the realistic Gamow shell model, was introduced later. The present review focuses on the recent applications of the GSM to drip line nuclei.

Quasifree Neutron Knockout Reaction Reveals a Small s-Orbital Component in the Borromean Nucleus ^{17}B

A kinematically complete quasifree (p,pn) experiment in inverse kinematics was performed to study the structure of the Borromean nucleus ^{17}B, which had long been considered to have a neutron halo. By analyzing the momentum distributions and exclusive cross sections, we obtained the spectroscopic factors for 1s_{1/2} and 0d_{5/2} orbitals, and a surprisingly small percentage of 9(2)% was determined for 1s_{1/2}. Our finding of such a small 1s_{1/2} component and the halo features reported in prior experiments can be explained by the deformed relativistic Hartree-Bogoliubov theory in continuum, revealing a definite but not dominant neutron halo in ^{17}B. The present work gives the smallest s- or p-orbital component among known nuclei exhibiting halo features and implies that the dominant occupation of s or p orbitals is not a prerequisite for the occurrence of a neutron halo.

Ab Initio Calculations of Low-Energy Nuclear Scattering Using Confining Potential Traps

A recently modified method to enable low-energy nuclear scattering results to be extracted from the discrete energy levels of the target-projectile clusters confined by harmonic potential traps is tested. We report encouraging results for neutron-α and neutron-^{24}O elastic scattering from analyzing the trapped levels computed using two different ab initio nuclear structure methods. The n-α results have also been checked against a direct ab initio reaction calculation. The n-^{24}O results demonstrate the approach's applicability for a large range of systems provided their spectra in traps can be computed by ab initio methods. A key ingredient is a rigorous understanding of the errors in the calculated energy levels caused by inevitable Hilbert-space truncations in the ab initio methods.

Elastic Antiproton-Nucleus Scattering from Chiral Forces

Elastic scattering of antiprotons off ^{4}He, ^{12}C, and ^{16,18}O is described for the first time with a consistent microscopic approach based on the calculation of an optical potential (OP) describing the antiproton-target interaction. The OP is derived using the recent antiproton-nucleon (p[over ¯]N) chiral interaction to calculate the p[over ¯]N t matrix, while the target densities are computed with the ab initio no-core shell model using chiral interactions as well. Our results are in good agreement with the existing experimental data and the results computed at different chiral orders of the p[over ¯]N interaction display a well-defined convergence pattern.

Extending the Southern Shore of the Island of Inversion to ^{28}F

Detailed spectroscopy of the neutron-unbound nucleus ^{28}F has been performed for the first time following proton/neutron removal from ^{29}Ne/^{29}F beams at energies around 230  MeV/nucleon. The invariant-mass spectra were reconstructed for both the ^{27}F^{(*)}+n and ^{26}F^{(*)}+2n coincidences and revealed a series of well-defined resonances. A near-threshold state was observed in both reactions and is identified as the ^{28}F ground state, with S_{n}(^{28}F)=-199(6)  keV, while analysis of the 2n decay channel allowed a considerably improved S_{n}(^{27}F)=1620(60)  keV to be deduced. Comparison with shell-model predictions and eikonal-model reaction calculations have allowed spin-parity assignments to be proposed for some of the lower-lying levels of ^{28}F. Importantly, in the case of the ground state, the reconstructed ^{27}F+n momentum distribution following neutron removal from ^{29}F indicates that it arises mainly from the 1p_{3/2} neutron intruder configuration. This demonstrates that the island of inversion around N=20 includes ^{28}F, and most probably ^{29}F, and suggests that ^{28}O is not doubly magic.

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.

A 21st Century View of Nuclear Structure

Exploiting exact and special symmetries to unmask simplicity within complexity, which remains the “holy grail” of nuclear physics, will be considered within its historical context and as evolving through 21st century ab initio methods, including emerging results linked to the internal structure of nucleons. Some exemplar results for very light to medium mass nuclei will be presented, and what these may portend for heavier systems, including species beyond known lines of stability, will be proffered.

Evidence for prevalent Z = 6 magic number in neutron-rich carbon isotopes

The nuclear shell structure, which originates in the nearly independent motion of nucleons in an average potential, provides an important guide for our understanding of nuclear structure and the underlying nuclear forces. Its most remarkable fingerprint is the existence of the so-called magic numbers of protons and neutrons associated with extra stability. Although the introduction of a phenomenological spin-orbit (SO) coupling force in 1949 helped in explaining the magic numbers, its origins are still open questions. Here, we present experimental evidence for the smallest SO-originated magic number (subshell closure) at the proton number six in 13-20C obtained from systematic analysis of point-proton distribution radii, electromagnetic transition rates and atomic masses of light nuclei. Performing ab initio calculations on 14,15C, we show that the observed proton distribution radii and subshell closure can be explained by the state-of-the-art nuclear theory with chiral nucleon-nucleon and three-nucleon forces, which are rooted in the quantum chromodynamics.

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.

Light-Nuclei Spectra from Chiral Dynamics

In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our knowledge of few- and many-body systems. In this Letter, we present Green's function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4-12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.

The equation of motion phonon method and its application in the neutron rich oxygen region

An equation of motion phonon method, developed for even-even nuclear systems and extended to odd nuclei, is applied to 22O and to its odd neighbors 23O and 23F. A calculation using the chiral potential NNLOopt is carried out in a space encompassing up to two phonons. The computed dipole cross section in 22O and the spectra of 22O and 23O are in a satisfactory agreement with the experimental data. However, the calculation describes poorly the spectrum of 23F. This discrepancy originates from the strong coupling between the odd proton and the 22O phonons of neutron nature. This coupling pushes down in energy several states enhancing the level density at low energy. We suggest that a viable route for the solution of this problem could be the inclusion of the three-body interaction using the new chiral potential NNLOsat.

Can Tetraneutron be a Narrow Resonance?

The search for a resonant four-neutron system has been revived thanks to the recent experimental hints reported in [1]. The existence of such a system would deeply impact our understanding of nuclear matter and requires a critical investigation. In this work, we study the existence of a four-neutron resonance in the quasistationary formalism using ab initio techniques with various two-body chiral interactions. We employ no-core Gamow shell model and density matrix renormalization group method, both supplemented by the use of natural orbitals and a new identification technique for broad resonances. We demonstrate that while the energy of the four-neutron system may be compatible with the experimental value, its width must be larger than the reported upper limit, supporting the interpretation of the experimental observation as a reaction process too short to form a nucleus.

Proton Distribution Radii of ^{12-19}C Illuminate Features of Neutron Halos

Proton radii of ^{12-19}C densities derived from first accurate charge changing cross section measurements at 900A  MeV with a carbon target are reported. A thick neutron surface evolves from ∼0.5  fm in ^{15}C to ∼1  fm in ^{19}C. The halo radius in ^{19}C is found to be 6.4±0.7  fm as large as ^{11}Li. Ab initio calculations based on chiral nucleon-nucleon and three-nucleon forces reproduce the radii well.

Evidence of soft dipole resonance in ^{11}li with isoscalar character

The first conclusive evidence of a dipole resonance in ^{11}Li having isoscalar character observed from inelastic scattering with a novel solid deuteron target is reported. The experiment was performed at the newly commissioned IRIS facility at TRIUMF. The results show a resonance peak at an excitation energy of 1.03±0.03 MeV with a width of 0.51±0.11 MeV (FWHM). The angular distribution is consistent with a dipole excitation in the distorted-wave Born approximation framework. The observed resonance energy together with shell model calculations show the first signature that the monopole tensor interaction is important in ^{11}Li. The first ab initio calculations in the coupled cluster framework are also presented.

Power counting of contact-range currents in effective field theory

We analyze the power counting of two-body currents in nuclear effective field theories (EFTs). We find that the existence of nonperturbative physics at low energies, which is manifest in the existence of the deuteron and the ^{1}S_{0} NN virtual bound state, combined with the appearance of singular potentials in versions of nuclear EFT that incorporate chiral symmetry, modifies the renormalization-group flow of the couplings associated with contact operators that involve nucleon-nucleon pairs and external fields. The order of these couplings is thereby enhanced with respect to the naive-dimensional-analysis estimate. Consequently, short-range currents enter at a lower order in the chiral EFT than has been appreciated up until now, and their impact on low-energy observables is concomitantly larger. We illustrate the changes in the power counting with a few low-energy processes involving external probes and few-nucleon systems, including electron-deuteron elastic scattering and radiative neutron capture by protons.

Effects of three-nucleon forces and two-body currents on Gamow-Teller strengths

We optimize chiral interactions at next-to-next-to leading order to observables in two- and three-nucleon systems and compute Gamow-Teller transitions in 14C and (22,24)O using consistent two-body currents. We compute spectra of the daughter nuclei 14N and (22,24)F via an isospin-breaking coupled-cluster technique, with several predictions. The two-body currents reduce the Ikeda sum rule, corresponding to a quenching factor q2≈0.84-0.92 of the axial-vector coupling. The half-life of 14C depends on the energy of the first excited 1+ state, the three-nucleon force, and the two-body current.

Nuclear matter from effective quark-quark interaction

We study neutron matter and symmetric nuclear matter with the quark-meson model for the two-nucleon interaction. The Bethe-Bruckner-Goldstone many-body theory is used to describe the correlations up to the three hole-line approximation with no extra parameters. At variance with other nonrelativistic realistic interactions, the three hole-line contribution turns out to be non-negligible and to have a substantial saturation effect. The saturation point of nuclear matter, the compressibility, the symmetry energy, and its slope are within the phenomenological constraints. Since the interaction also reproduces fairly well the properties of the three-nucleon system, these results indicate that the explicit introduction of the quark degrees of freedom within the considered constituent quark model is expected to reduce the role of three-body forces.

Quantum Monte Carlo calculations of light nuclei using chiral potentials

We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of A=3 and A=4 systems. We also carry out the first few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon systems using quantum Monte Carlo methods with chiral 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).

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.

Quantum Monte Carlo calculations of neutron matter with nonlocal chiral interactions

We present fully nonperturbative quantum Monte Carlo calculations with nonlocal chiral effective field theory (EFT) interactions for the ground-state properties of neutron matter. The equation of state, the nucleon chemical potentials, and the momentum distribution in pure neutron matter up to one and a half times the nuclear saturation density are computed with a newly optimized chiral EFT interaction at next-to-next-to-leading order. This work opens the way to systematic order by order benchmarking of chiral EFT interactions and ab initio prediction of nuclear properties while respecting the symmetries of quantum chromodynamics.