No Arabic abstract
We develop the variational/parquet diagram approach to the structure of nuclear systems with strongly state-dependent interactions. For that purpose, we combine ideas of the general Jastrow-Feenberg variational method and the local parquet-diagram theory for bosons with state-dependent interactions (R. A. Smith and A. D. Jackson, Nucl. Phys. {bf 476}, 448 (1988)). The most tedious aspect of variational approaches, namely the symmetrization of an operator dependent variational wave function, is thereby avoided. We carry out calculations for neutron matter interacting via the Reid and Argonne $v_6$ models of the nucleon-nucleon interaction. While the equation of state is a rather robust quantity that comes out reasonably well even in very simplistic approaches, we show that effective interactions, which are the essential input for calculating dynamic properties, depend sensitively on the quality of the treatment of the many-body problem.
We apply parquet-diagram summation methods for the calculation of the superfluid gap in $S$-wave pairing in neutron matter for realistic nucleon-nucleon interactions such as the Argonne $v_6$ and the Reid $v_6$ potentials. It is shown that diagrammatic contributions that are outside the parquet class play an important role. These are, in variational theories, identified as so-called commutator contributions. Moreover, using a particle-hole propagator appropriate for a superfluid system results in the suppression of the spin-channel contribution to the induced interaction. Applying these corrections to the pairing interaction, our results agree quite well with Quantum Monte Carlo data.
We develop a manifestly microscopic method to deal with strongly interacting nuclear systems that have different interactions in spin-singlet and spin-triplet states. In a first step we analyze variational wave functions that have been suggested to describe such systems, and demonstrate that the so-called commutator contributions can have important effects whenever the interactions in the spin-singlet and the spin-triplet states are very different. We then identify these contributions as terms that correspond, in the language of perturbation theory, to non-parquet diagrams. We include these diagrams in a way that is suggested by the Jastrow-Feenberg approach and show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects.
We develop the variational and correlated basis functions/parquet-diagram theory of strongly interacting normal and superfluid systems. The first part of this contribution is devoted to highlight the connections between the Euler equations for the Jastrow-Feenberg wave function on the one hand side, and the ring, ladder, and self-energy diagrams of parquet-diagram theory on the other side. We will show that these subsets of Feynman diagrams are contained, in a local approximation, in the variational wave function. In the second part of this work, we derive the fully optimized Fermi-Hypernetted Chain (FHNC-EL) equations for a superfluid system. Close examination of the procedure reveals that the naive application of these equations exhibits spurious unphysical properties for even an infinitesimal superfluid gap. We will conclude that it is essential to go {em beyond/} the usual Jastrow-Feenberg approximation and to include the exact particle-hole propagator to guarantee a physically meaningful theory and the correct stability range. We will then implement this method and apply it to neutron matter and low density Fermi liquids interacting via the Lennard-Jones model interaction and the Poschl-Teller interaction. While the quantitative changes in the magnitude of the superfluid gap are relatively small, we see a significant difference between applications for neutron matter and the Lennard-Jones and Poschl-Teller systems. Despite the fact that the gap in neutron matter can be as large as half the Fermi energy, the corrections to the gap are relatively small. In the Lennard-Jones and Poschl-Teller models, the most visible consequence of the self-consistent calculation is the change in stability range of the system.
The equation of state (EOS) for neutron star (NS) crusts is studied in the Thomas-Fermi (TF) approximation using the EOS for uniform nuclear matter obtained by the variational method with the realistic nuclear Hamiltonian. The parameters associated with the nuclear three-body force, which are introduced to describe the saturation properties, are finely adjusted so that the TF calculations for isolated atomic nuclei reproduce the experimental data on masses and charge distributions satisfactorily. The resulting root-mean-square deviation of the masses from the experimental data for mass-measured nuclei is about 3 MeV. With use of the nuclear EOS thus determined, the nuclei in the crust of NS at zero temperature are calculated. The predicted proton numbers of the nuclei in the crust of NS are close to the gross behavior of the results by Negele and Vautherin, while they are larger than those for the EOS by Shen et al. due to the difference in the symmetry energy. The density profile of NS is calculated with the constructed EOS.
We discuss a variational calculation for nuclear shell-model calculations and propose a new procedure for the energy-variance extrapolation (EVE) method using a sequence of the approximated wave functions obtained by the variational calculation. The wave functions are described as linear combinations of the parity, angular-momentum projected Slater determinants, the energy of which is minimized by the conjugate gradient method obeying the variational principle. The EVE generally works well using the wave functions, but we found some difficult cases where the EVE gives a poor estimation. We discuss the origin of the poor estimation concerning shape coexistence. We found that the appropriate reordering of the Slater determinants allows us to overcome this difficulty and to reduce the uncertainty of the extrapolation.