The properties of symmetric nuclear and pure neutron matter are investigated within an extended self-consistent Greens function method that includes the effects of three-body forces. We use the ladder approximation for the study of infinite nuclear matter and incorporate the three-body interaction by means of a density-dependent two-body force. This force is obtained via a correlated average over the third particle, with an in-medium propagator consistent with the many-body calculation we perform. We analyze different prescriptions in the construction of the average and conclude that correlations provide small modifications at the level of the density-dependent force. Microscopic as well as bulk properties are studied, focusing on the changes introduced by the density dependent two-body force. The total energy of the system is obtained by means of a modified Galitskii-Migdal-Koltun sum rule. Our results validate previously used uncorrelated averages and extend the availability of chirally motivated forces to a larger density regime.
We present calculations for symmetric nuclear matter using chiral nuclear interactions within the Self-Consistent Greens Functions approach in the ladder approximation. Three-body forces are included via effective one-body and two-body interactions, computed from an uncorrelated average over a third particle. We discuss the effect of the three-body forces on the total energy, computed with an extended Galitskii-Migdal-Koltun sum-rule, as well as on single-particle properties. Saturation properties are substantially improved when three-body forces are included, but there is still some underlying dependence on the SRG evolution scale.
We have studied the scaling properties of the electromagnetic response functions of $^4$He and $^{12}$C nuclei computed by the Greens Function Monte Carlo approach, retaining only the one-body current contribution. Longitudinal and transverse scaling functions have been obtained in the relativistic and non relativistic cases and compared to experiment for various kinematics. The characteristic asymmetric shape of the scaling function exhibited by data emerges in the calculations in spite of the non relativistic nature of the model. The results are consistent with scaling of zeroth, first and second kinds. Our analysis reveals a direct correspondence between the scaling and the nucleon-density response functions.
We compute inclusive electron-nucleus cross sections using ab initio spectral functions of $^4$He and $^{16}$O obtained within the Self Consistent Greens Function approach. The formalism adopted is based on the factorization of the spectral function and the nuclear transition matrix elements. This allows to provide an accurate description of nuclear dynamics and to account for relativistic effects in the interaction vertex. Our calculations use a saturating chiral Hamiltonian in order reproduce the correct nuclear sizes. When final state interactions for the struck particle are accounted for, we find nice agreement between the data and the theory for the inclusive electron-$^{16}$O cross section. The results lay the foundations for future applications of the Self Consistent Greens Function method, in both closed and open shell nuclei, to neutrino data analysis. This work also presents results for the point-proton, charge and single-nucleon momentum distribution of the same two nuclei. The center of mass can affect these quantities for light nuclei and cannot be separated cleanly in most ab initio post-Hartree-Fock methods. In order to address this, we developed a Metropolis Monte Carlo calculation in which the center of mass coordinate can be subtracted exactly from the trial wave function and the expectation values. We gauged this effect for $^4$He by removing the center of mass effect from the Optimal Reference State wave function that is generated during the Self Consistent Greens Function calculations. Our findings clearly indicate that the residual center of mass contribution strongly modifies calculated matter distributions with respect to those obtained in the intrinsic frame. Hence, its subtraction is crucial for a correct description of light nuclei.
We present a new charge self-consistent scheme combining Density Functional and Dynamical Mean Field Theory, which uses Greens function of multiple scattering-type. In this implementation the many-body effects are incorporated into the Kohn-Sham iterative scheme without the need for the numerically ill-posed analytic continuation of the Greens function and of the self-energy. This is achieved by producing the Kohn-Sham Hamiltonian in the sub-space of correlated partial waves and allows to formulate the Greens function directly on the Matsubara axis. The spectral moments of the Matsubara Greens function enable us to put together the real space charge density, therefore the charge self-consistency can be achieved. Our results for the spectral functions (density of states) and equation of state curves for transition metal elements, Fe, Ni and FeAl compound agree very well with those of Hamiltonian based LDA+DMFT implementations. The current implementation improves on numerical accuracy, requires a minimal effort besides the multiple scattering formulation and can be generalized in several ways that are interesting for applications to real materials.
We derive from the subleading contributions to the chiral three-nucleon force (long-range terms, published in Phys.,Rev.,C,77, 064004 (2008)) a density-dependent two-nucleon interaction $V_text{med}$ in isospin-symmetric, spin-saturated nuclear matter. Following the division of the pertinent 3N-diagrams into two-pion exchange topology, two-pion-one-pion exchange topology and ring topology, we evaluate for these all self-closings and concatenations of nucleon-lines to an in-medium loop. The momentum and $k_f$-dependent potentials associated with the isospin operators ($1$ and $vectau_1!cdot!vectau_2$) and five independent spin-structures are expressed in terms of functions, which are either given in closed analytical form or require at most one numerical integration. In the same way we treat the $2pi$-exchange 3N-force up to fourth order. Our results for $V_text{med}$ are most helpful to implement the long-range subleading chiral 3N-forces into nuclear many-body calculations.