Evolution of nuclear spectra with nuclear forces

Direct Observation of Competing M1 and M3 Transitions in ^{10}B

Excited states in ^{10}B were populated with the ^{10}B(p,p^{'}γ)^{10}B^{*} reaction at 8.5 MeV and their γ decay was investigated via coincidence γ-ray spectroscopy. The emitted γ rays were measured using large-volume LaBr_{3}:Ce and CeBr_{3} detectors placed in anti-Compton shields. This allowed the observation of weak γ-ray transitions, such as the M3 transition between the J^{π},T=0^{+},1 isobaric analog state (IAS) and the J^{π},T=3^{+},0 ground state and the E2 transition between the J^{π},T=2_{1}^{+},0 state and the IAS, i.e., performing measurements of branching ratios at the level of λ≥10^{-4}. For the first time in ^{10}B, the competing M1 and M3 transitions from the decay of the IAS have been observed in a γ spectroscopy experiment. The experimental results are compared with ab initio no-core shell model calculation using the newest version of the local position-space chiral N^{3}LO nucleon-nucleon interaction. The calculations reproduce correctly the ordering of the bound states in ^{10}B, and are in reasonable agreement with the observed branching ratios and reduced transition probabilities.

Perturbative Quantum Monte Carlo Method for Nuclear Physics

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. We demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to ^{16}O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second-order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

The 1S0 Pairing Gap in Neutron Matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

Nonperturbative Extraction of the Effective Mass in Neutron Matter

We carry out nonperturbative calculations of the single-particle excitation spectrum in strongly interacting neutron matter. These are microscopic quantum Monte Carlo computations of many-neutron energies at different densities as well as several distinct excited states. As input, we employ both phenomenological and chiral two- and three-nucleon interactions. We use the single-particle spectrum to extract the effective mass in neutron matter. With a view to systematizing the error involved in this extraction, we carefully assess the impact of finite-size effects on the quasiparticle dispersion relation. We find an effective-mass ratio that drops from 1 as the density is increased. We conclude by connecting our results with the physics of ultracold gases as well as with energy-density functional theories of nuclei and neutron-star matter.

Exploring finite-size effects in strongly correlated systems

Complexities greatly limit any study of strongly correlated systems to a small number of particles. Thus, any attempt at understanding infinite systems such as those arising from neutron matter (NM) must consider finite-size (FS) effects at play when below the thermodynamic limit (TL). In these conference proceedings we provide some examples of FS effects at work and discuss our prescription for extrapolating the physics of extended systems. We present our methodology and calculations performed for an assortment of strongly correlated (SC) systems. Ab initio, non-perturbative Quantum Monte Carlo (QMC) methods can be employed to accurately compute ground-state energies and finite-temperature properties. We apply these to periodically modulated NM and use our results to constrain phenomenological theories of nuclei and study the static response of NM.

Onset of η nuclear binding

Recent studies of η nuclear quasibound states by the Jerusalem-Prague Collaboration are reviewed, focusing on stochastic variational method self consistent calculations of η few-nucleon systems. These calculations suggest that a minimum value Re aηN ≈ 1 fm (0.7 fm) is needed to bind η3He (η4He).

Static Response of Neutron Matter

We generalize the problem of strongly interacting neutron matter by adding a periodic external modulation. This allows us to study from first principles a neutron system that is extended and inhomogeneous, with connections to the physics of both neutron-star crusts and neutron-rich nuclei. We carry out fully nonperturbative microscopic quantum Monte Carlo calculations of the energy of neutron matter at different densities, as well as different strengths and periodicities of the external potential. In order to remove systematic errors, we examine finite-size effects and the impact of the wave function ansatz. We also make contact with energy-density functional theories of nuclei and disentangle isovector gradient contributions from bulk properties. Finally, we calculate the static density-density linear response function of neutron matter and compare it with the response of other physical systems.

Hyperon puzzle: hints from quantum Monte Carlo calculations

The onset of hyperons in the core of neutron stars and the consequent softening of the equation of state have been questioned for a long time. Controversial theoretical predictions and recent astrophysical observations of neutron stars are the grounds for the so-called hyperon puzzle. We calculate the equation of state and the neutron star mass-radius relation of an infinite systems of neutrons and Λ particles by using the auxiliary field diffusion Monte Carlo algorithm. We find that the three-body hyperon-nucleon interaction plays a fundamental role in the softening of the equation of state and for the consequent reduction of the predicted maximum mass. We have considered two different models of three-body force that successfully describe the binding energy of medium mass hypernuclei. Our results indicate that they give dramatically different results on the maximum mass of neutron stars, not necessarily incompatible with the recent observation of very massive neutron stars. We conclude that stronger constraints on the hyperon-neutron force are necessary in order to properly assess the role of hyperons in neutron stars.

Collective modes in light nuclei from first principles

Results for ab initio no-core shell model calculations in a symmetry-adapted SU(3)-based coupling scheme demonstrate that collective modes in light nuclei emerge from first principles. The low-lying states of 6Li, 8Be, and 6He are shown to exhibit orderly patterns that favor spatial configurations with strong quadrupole deformation and complementary low intrinsic spin values, a picture that is consistent with the nuclear symplectic model. The results also suggest a pragmatic path forward to accommodate deformation-driven collective features in ab initio analyses when they dominate the nuclear landscape.

New isomers in the full seniority scheme of neutron-rich lead isotopes: the role of effective three-body forces

The neutron-rich lead isotopes, up to (216)Pb, have been studied for the first time, exploiting the fragmentation of a primary uranium beam at the FRS-RISING setup at GSI. The observed isomeric states exhibit electromagnetic transition strengths which deviate from state-of-the-art shell-model calculations. It is shown that their complete description demands the introduction of effective three-body interactions and two-body transition operators in the conventional neutron valence space beyond (208)Pb.

Continuum effects and three-nucleon forces in neutron-rich oxygen isotopes

We employ interactions from chiral effective field theory and compute binding energies, excited states, and radii for isotopes of oxygen with the coupled-cluster method. Our calculation includes the effects of three-nucleon forces and of the particle continuum, both of which are important for the description of neutron-rich isotopes in the vicinity of the nucleus 24O. Our main results are the placement of the neutron drip line at 24O, the assignment of spins, parities and resonance widths for several low-lying states of the drip line nucleus, and an efficient approximation that incorporates the effects of three-body interactions.

Tensor force manifestations in ab initio study of the 2H(d,γ)4He, 2H(d,p)3H, and 2H(d,n)3He reactions

The (2)H(d,p)(3)H, (2)H(d,n)(3)He, and (2)H(d,γ)(4)He reactions are studied at low energies in a multichannel ab initio model that takes into account the distortions of the nuclei. The internal wave functions of these nuclei are given by the stochastic variational method with the AV8' realistic interaction and a phenomenological three-body force included to reproduce the two-body thresholds. The obtained astrophysical S factors are all in very good agreement with the experiment. The most important channels for both transfer and radiative capture are identified by comparing to calculations with an effective central force. They are all found to dominate thanks to the tensor force.

Equation of state of superfluid neutron matter and the calculation of the 1S0 pairing gap

We present a quantum Monte Carlo study of the zero-temperature equation of state of neutron matter and the computation of the 1S0 pairing gap in the low-density regime with rho < 0.04 fm(-3). The system is described by a nonrelativistic nuclear Hamiltonian including both two- and three-nucleon interactions of the Argonne and Urbana type. This model interaction provides very accurate results in the calculation of the binding energy of light nuclei. A suppression of the gap with respect to the pure BCS theory is found, but sensibly weaker than in other works that attempt to include polarization effects in an approximate way.

Proton-neutron and proton-proton correlations in medium-weight nuclei and the role of the tensor force

A detailed analysis of the effect of short-range and tensor correlations on one- and two-nucleon momentum distributions of medium-weight nuclei (12<or=A<or=40) is carried out. Although our Letter is primarily aimed at understanding the role of the tensor force on nucleon momentum distributions of medium-weight nuclei, the possible relevance of our results for the interpretation of (e, e'N) and (e, e'2N) processes at high Q2, aimed at investigating nucleon-nucleon correlations, is discussed.

Auxiliary field diffusion Monte Carlo calculation of nuclei with A < or = 40 with tensor interactions

We calculate the ground-state energy of (4)He, (8)He, (16)O, and (40)Ca using the auxiliary field diffusion Monte Carlo method in the fixed-phase approximation and the Argonne v(6)' interaction which includes a tensor force. Comparison of our light nuclei results to those of Green's function Monte Carlo calculations shows the accuracy of our method for both open and closed-shell nuclei. We also apply it to (16)O and (40)Ca to show that quantum Monte Carlo methods are now applicable to larger nuclei.

Tensor forces and the ground-state structure of nuclei

Two-nucleon momentum distributions are calculated for the ground states of nuclei with mass number A< or =8, using variational Monte Carlo wave functions derived from a realistic Hamiltonian with two- and three-nucleon potentials. The momentum distribution of np pairs is found to be much larger than that of pp pairs for values of the relative momentum in the range (300-600) MeV/c and vanishing total momentum. This order of magnitude difference is seen in all nuclei considered and has a universal character originating from the tensor components present in any realistic nucleon-nucleon potential. The correlations induced by the tensor force strongly influence the structure of np pairs, which are predominantly in deuteronlike states, while they are ineffective for pp pairs, which are mostly in 1S0 states. These features should be easily observable in two-nucleon knockout processes, such as A(e,e'np) and A(e,e'pp).

Quantum Monte Carlo calculations of symmetric nuclear matter

We present an accurate numerical study of the equation of state of nuclear matter based on realistic nucleon-nucleon interactions by means of auxiliary field diffusion Monte Carlo (AFDMC) calculations. The AFDMC method samples the spin and isospin degrees of freedom allowing for quantum simulations of large nucleonic systems and represents an important step forward towards a quantitative understanding of problems in nuclear structure and astrophysics.

Ab-initio coupled-cluster study of 16O

We report converged results for the ground and excited states and matter density of 16O using realistic two-body nucleon-nucleon interactions and coupled-cluster methods and algorithms developed in quantum chemistry. Most of the binding is obtained with the coupled-cluster singles and doubles approach. Additional binding due to three-body clusters (triples) is minimal. The coupled-cluster method with singles and doubles provides a good description of the matter density, charge radius, charge form factor, and excited states of a one-particle, one-hole nature, but it cannot describe the first-excited 0(+) state. Incorporation of triples has no effect on the latter finding.

Neutrino-12C scattering in the Ab initio shell model with a realistic three-body interaction

We investigate cross sections for neutrino-12C exclusive scattering and for muon capture on 12C using wave functions obtained in the ab initio no-core shell model. In our parameter-free calculations with basis spaces up to the 6 variant Planck's over 2pi Omega we show that realistic nucleon-nucleon interactions, like, e.g., the CD-Bonn, underpredict the experimental cross sections by more than a factor of 2. By including a realistic three-body interaction, Tucson-Melbourne TM'(99), the cross sections are enhanced significantly and a much better agreement with experiment is achieved. At the same time, the TM'(99) interaction improves the calculated level ordering in 12C. The comparison between the CD-Bonn and the three-body calculations provides strong confirmation for the need to include a realistic three-body interaction to account for the spin-orbit strength in p-shell nuclei.

Can modern nuclear Hamiltonians tolerate a bound tetraneutron?

I show that it does not seem possible to change modern nuclear Hamiltonians to bind a tetraneutron without destroying many other successful predictions of those Hamiltonians. This means that, should a recent experimental claim of a bound tetraneutron be confirmed, our understanding of nuclear forces will have to be significantly changed. I also point out some errors in previous theoretical studies of this problem.