No Arabic abstract
We present a very brief description of the Hartree-Fock method in nuclear structure physics, discuss the numerical methods used to solve the self-consistent equations, and analyze the precision and convergence properties of solutions. As an application we present results pertaining to quadrupole moments and single-particle quadrupole polarizations in superdeformed nuclei with A~60.
We study relativistic nuclear matter in the $sigma - omega$ model including the ring-sum correlation energy. The model parameters are adjusted self-consistently to give the canonical saturation density and binding energy per nucleon with the ring energy included. Two models are considered, mean-field-theory where we neglect vacuum effects, and the relativistic Hartree approximation where such effects are included but in an approximate way. In both cases we find self-consistent solutions and present equations of state. In the mean-field case the ring energy completely dominates the attractive part of the energy density and the elegant saturation mechanism of the standard approach is lost, namely relativistic quenching of the scalar attraction. In the relativistic Hartree approach the vacuum effects are included in an approximate manner using vertex form factors with a cutoff of 1 - 2 GeV, the range expected from QCD. Due to the cutoff, the ring energy for this case is significantlysmaller, and we obtain self-consistent solutions which preserve the basic saturation mechanism of the standard relativistic approach.
In recent years, the combination of precise quantum Monte Carlo (QMC) methods with realistic nuclear interactions and consistent electroweak currents, in particular those constructed within effective field theories (EFTs), has lead to new insights in light and medium-mass nuclei, neutron matter, and electroweak reactions. This compelling new body of work has been made possible both by advances in QMC methods for nuclear physics, which push the bounds of applicability to heavier nuclei and to asymmetric nuclear matter and by the development of local chiral EFT interactions up to next-to-next-to-leading order and minimally nonlocal interactions including $Delta$ degrees of freedom. In this review, we discuss these recent developments and give an overview of the exciting results for nuclei, neutron matter and neutron stars, and electroweak reactions.
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.
The prospects of extracting new physics signals in a coherent elastic neutrino-nucleus scattering (CE$ u$NS) process are limited by the precision with which the underlying nuclear structure physics, embedded in the weak nuclear form factor, is known. We present microscopic nuclear structure physics calculations of charge and weak nuclear form factors and CE$ u$NS cross sections on $^{12}$C, $^{16}$O, $^{40}$Ar, $^{56}$Fe and $^{208}$Pb nuclei. We obtain the proton and neutron densities, and charge and weak form factors by solving Hartree-Fock equations with a Skyrme (SkE2) nuclear potential. We validate our approach by comparing $^{208}$Pb and $^{40}$Ar charge form factor predictions with elastic electron scattering data. In view of the worldwide interest in liquid-argon based neutrino and dark matter experiments, we pay special attention to the $^{40}$Ar nucleus and make predictions for the $^{40}$Ar weak form factor and the CE$ u$NS cross sections. Furthermore, we attempt to gauge the level of theoretical uncertainty pertaining to the description of the $^{40}$Ar form factor and CE$ u$NS cross sections by comparing relative differences between recent microscopic nuclear theory and widely-used phenomenological form factor predictions. Future precision measurements of CE$ u$NS will potentially help in constraining these nuclear structure details that will in turn improve prospects of extracting new physics.
Microscopic calculations of the electromagnetic response of medium-mass nuclei are now feasible thanks to the availability of realistic nuclear interactions with accurate saturation and spectroscopic properties, and the development of large-scale computing methods for many-body physics. The purpose is to compute isovector dipole electromagnetic (E1) response and related quantities, i.e. integrated dipole cross section and polarizability, and compare with data from photoabsorption and Coulomb excitation experiments. The single-particle propagator is obtained by solving the Dyson equation, where the self-energy includes correlations non-perturbatively through the Algebraic Diagrammatic Construction (ADC) method. The particle-hole ($ph$) polarization propagator is treated in the Dressed Random Phase Approximation (DRPA), based on an effective correlated propagator that includes some $2p2h$ effects but keeps the same computation scaling as the standard Hartree-Fock propagator. The E1 responses for $^{14,16,22,24}$O, $^{36,40,48,52,54,70}$Ca and $^{68}$Ni have been computed: the presence of a soft dipole mode of excitation for neutron-rich nuclei is found, and there is a fair reproduction of the low-energy part of the experimental excitation spectrum. This is reflected in a good agreement with the empirical dipole polarizability values. For a realistic interaction with an accurate reproduction of masses and radii up to medium-mass nuclei, the Self-Consistent Greens Function method provides a good description of the E1 response, especially in the part of the excitation spectrum below the Giant Dipole Resonance. The dipole polarizability is largely independent from the strategy of mapping the dressed propagator to a simplified one that is computationally manageable