No Arabic abstract
We present a computational approach to infinite nuclear matter employing Hartree-Fock theory, many-body perturbation theory and coupled cluster theory. These lectures are closely linked with those of chapters 9, 10 and 11 and serve as input for the correlation functions employed in Monte Carlo calculations in chapter 9, the in-medium similarity renormalization group theory of dense fermionic systems of chapter 10 and the Greens function approach in chapter 11. We provide extensive code examples and benchmark calculations, allowing thereby an eventual reader to start writing her/his own codes. We start with an object-oriented serial code and end with discussions on strategies for porting the code to present and planned high-performance computing facilities.
We present a simplified method to generate the Hartree-Fock Gamow basis from realistic nuclear forces. The Hartree-Fock iteration in the harmonic-oscillator basis is first performed, and then the obtained HF potential is analytically continued to the complex-k plane, finally by solving the Schrodinger equation in the complex-k plane the Gamow basis is obtained. As examples, the method is applied to 4He and 22O with the renormalized chiral N3LO potential. The basis obtained which includes bound, resonant and scattering states can be further used in many-body calculations to study weakly bound nuclei.
On the way of a microscopic derivation of covariant density functionals, the first complete solution of the relativistic Brueckner-Hartree-Fock (RBHF) equations is presented for symmetric nuclear matter. In most of the earlier investigations, the $G$-matrix is calculated only in the space of positive energy solutions. On the other side, for the solution of the relativistic Hartree-Fock (RHF) equations, also the elements of this matrix connecting positive and negative energy solutions are required. So far, in the literature, these matrix elements are derived in various approximations. We discuss solutions of the Thompson equation for the full Dirac space and compare the resulting equation of state with those of earlier attempts in this direction.
We investigate the order-by-order convergence behavior of many-body perturbation theory (MBPT) as a simple and efficient tool to approximate the ground-state energy of closed-shell nuclei. To address the convergence properties directly, we explore perturbative corrections up to 30th order and highlight the role of the partitioning for convergence. The use of a simple Hartree-Fock solution to construct the unperturbed basis leads to a convergent MBPT series for soft interactions, in contrast to, e.g., a harmonic oscillator basis. For larger model spaces and heavier nuclei, where a direct high-order MBPT calculation in not feasible, we perform third-order calculation and compare to advanced ab initio coupled-cluster calculations for the same interactions and model spaces. We demonstrate that third-order MBPT provides ground-state energies for nuclei up into tin isotopic chain that are in excellent agreement with the best available coupled-cluster results at a fraction of the computational cost.
Brueckner-Hartree-Fock theory allows to derive the $G$-matrix as an effective interaction between nucleons in the nuclear medium. It depends on the center of mass momentum $bm{P}$ of the two particles and on the two relative momenta $bm{q}$ and $bm{q}$ before and after the scattering process. In the evaluation of the total energy per particle in nuclear matter usually the angle averaged center of mass momentum approximation has been used. We derive in detail the exact expressions of the angular integrations of the momentum $bm{P}$ within relativistic Brueckner-Hartree-Fock (RBHF) theory, especially for the case of asymmetric nuclear matter. In order to assess the reliability of the conventional average momentum approximation for the binding energy, the saturation properties of symmetric and asymmetric nuclear matter are systematically investigated based on the realistic Bonn nucleon-nucleon potential. It is found that the exact treatment of the center of mass momentum leads to non-negligible contributions to the higher order physical quantities. The correlation between the symmetry energy $E_{mathrm{sym}}$, the slope parameter $L$, and the curvature $K_{mathrm{sym}}$ of the symmetry energy are investigated. The results of our RBHF calculations for the bulk parameters characterizing the equation of state are compared with recent constraints extracted from giant monopole resonance and isospin diffusion experiments.
With the relativistic representation of the nuclear tensor force that is included automatically by the Fock diagrams, we explored the self-consistent tensor effects on the properties of nuclear matter system. The analysis were performed within the density-dependent relativistic Hartree-Fock (DDRHF) theory. The tensor force is found to notably influence the saturation mechanism, the equation of state and the symmetry energy of nuclear matter, as well as the neutron star properties. Without introducing any additional free parameters, the DDRHF approach paves a natural way to reveal the tensor effects on the nuclear matter system.