Spectra and binding energy predictions of chiral interactions for 7Li

Using the no-core shell model approach, we report on the first results for 7Li based on the next-to-next-to-leading order chiral nuclear interaction. Both, two-nucleon and three-nucleon interactions are taken into account. We show that the p-shell nuclei are sensitive to the subleading parts of the chiral interactions including three-nucleon forces. Though chiral interactions are soft, we do not observe overbinding for this p-shell nucleus and find a realistic description for the binding energy, excitation spectrum and radius.Comment: 12 pages, 12 figure

From non-Hermitian effective operators to large-scale no-core shell model calculations for light nuclei

No-core shell model (NCSM) calculations using ab initio effective interactions are very successful in reproducing experimental nuclear spectra. The main theoretical approach is the use of effective operators, which include correlations left out by the truncation of the model space to a numerically tractable size. We review recent applications of the effective operator approach, within a NCSM framework, to the renormalization of the nucleon-nucleon interaction, as well as scalar and tensor operators.Comment: To be submited to J. Phys. A, special issue on "The Physics of Non-Hermitian Operators

Structure of A=10-13 nuclei with two- plus three-nucleon interactions from chiral effective field theory

Properties of finite nuclei are evaluated with two-nucleon (NN) and three-nucleon (NNN) interactions derived within chiral effective field theory (EFT). The nuclear Hamiltonian is fixed by properties of the A=2 system, except for two low-energy constants (LECs) that parameterize the short range NNN interaction. We constrain those two LECs by a fit to the A=3 system binding energy and investigate sensitivity of 4He, 6Li, 10,11B and 12,13C properties to the variation of the constrained LECs. We identify a preferred choice that gives globally the best description. We demonstrate that the NNN interaction terms significantly improve the binding energies and spectra of mid-p-shell nuclei not just with the preferred choice of the LECs but even within a wide range of the constrained LECs. At the same time, we find that a very high quality description of these nuclei requires further improvements to the chiral Hamiltonian.Comment: 4 pages, 4 figure

Ab initio Translationally Invariant Nonlocal One-body Densities from No-core Shell-model Theory

[Background:] It is well known that effective nuclear interactions are in general nonlocal. Thus if nuclear densities obtained from {\it ab initio} no-core-shell-model (NCSM) calculations are to be used in reaction calculations, translationally invariant nonlocal densities must be available. [Purpose:] Though it is standard to extract translationally invariant one-body local densities from NCSM calculations to calculate local nuclear observables like radii and transition amplitudes, the corresponding nonlocal one-body densities have not been considered so far. A major reason for this is that the procedure for removing the center-of-mass component from NCSM wavefunctions up to now has only been developed for local densities. [Results:] A formulation for removing center-of-mass contributions from nonlocal one-body densities obtained from NCSM and symmetry-adapted NCSM (SA-NCSM) calculations is derived, and applied to the ground state densities of $^4$He, $^6$Li, $^{12}$C, and $^{16}$O. The nonlocality is studied as a function of angular momentum components in momentum as well as coordinate space [Conclusions:] We find that the nonlocality for the ground state densities of the nuclei under consideration increases as a function of the angular momentum. The relative magnitude of those contributions decreases with increasing angular momentum. In general, the nonlocal structure of the one-body density matrices we studied is given by the shell structure of the nucleus, and can not be described with simple functional forms.Comment: 13 pages, 11 Figure

Ab initio Folding Potentials for Nucleon-Nucleus Scattering based on NCSM One-Body Densities

Calculating microscopic optical potentials for elastic nucleon-nucleus scattering has already led to large body of work in the past. For folding first-order calculations the nucleon-nucleon (NN) interaction and the one-body density of the nucleus were taken as input to rigorous calculations in a spectator expansion of the multiple scattering series. Based on the Watson expansion of the multiple scattering series we employ a nonlocal translationally invariant nuclear density derived from a chiral next-to-next-to-leading order (NNLO) and the very same interaction for consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for light nuclei. We calculate scattering observables, such as total, reaction, and differential cross sections as well as the analyzing power and the spin-rotation parameter, for elastic scattering of protons and neutrons from $^4$He, $^{6}$He, $^{12}$C, and $^{16}$O, in the energy regime between 100 and 200~MeV projectile kinetic energy, and compare to available data. Our calculations show that the effective nucleon-nucleus potential obtained from the first-order term in the spectator expansion of the multiple scattering expansion describes experiments very well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the NN amplitudes and the one-body nuclear density.Comment: 10 pages, 14 figures, 1 tabl

Extrapolation Method for the No-Core Shell Model

Nuclear many-body calculations are computationally demanding. An estimate of their accuracy is often hampered by the limited amount of computational resources even on present-day supercomputers. We provide an extrapolation method based on perturbation theory, so that the binding energy of a large basis-space calculation can be estimated without diagonalizing the Hamiltonian in this space. The extrapolation method is tested for 3H and 6Li nuclei. It will extend our computational abilities significantly and allow for reliable error estimates.Comment: 8 pages, 7 figures, PRC accepte

Nuclear Structure in the UCOM Framework: From Realistic Interactions to Collective Excitations

The Unitary Correlation Operator Method (UCOM) provides a means for nuclear structure calculations starting from realistic NN potentials. The dominant short-range central and tensor correlations are described explicitly by a unitary transformation. The application of UCOM in the context of the no-core shell model provides insight into the interplay between dominant short-range and residual long-range correlations in the nuclear many-body problem. The use of the correlated interaction within Hartree-Fock, many-body perturbation theory, and Random Phase Approximation gives access to various nuclear structure observables throughout the nuclear chart.Comment: 9 pages, 3 figures, invited talk at the 2nd Int. Conf. on "Collective Motion in Nuclei under Extreme Conditions" (COMEX 2), Sankt Goar, German

Faddeev calculations of break-up reactions with realistic experimental constraints

We present a method to integrate predictions from a theoretical model of a reaction with three bodies in the final state over the region of phase space covered by a given experiment. The method takes into account the true experimental acceptance, as well as variations of detector efficiency, and eliminates the need for a Monte-Carlo simulation of the detector setup. The method is applicable to kinematically complete experiments. Examples for the use of this method include several polarization observables in dp break-up at 270 MeV. The calculations are carried out in the Faddeev framework with the CD Bonn nucleon-nucleon interaction, with or without the inclusion of an additional three-nucleon force.Comment: 18 pages, 9 figure

Exact calculation of three-body contact interaction to second order

For a system of fermions with a three-body contact interaction the second-order contributions to the energy per particle $\bar E(k_f)$ are calculated exactly. The three-particle scattering amplitude in the medium is derived in closed analytical form from the corresponding two-loop rescattering diagram. We compare the (genuine) second-order three-body contribution to $\bar E(k_f)\sim k_f^{10}$ with the second-order term due to the density-dependent effective two-body interaction, and find that the latter term dominates. The results of the present study are of interest for nuclear many-body calculations where chiral three-nucleon forces are treated beyond leading order via a density-dependent effective two-body interaction.Comment: 9 pages, 6 figures, to be published in European Journal