Discrete wave-packet representation in nuclear matter calculations

Abstract

The Lippmann-Schwinger equation for the nucleon-nucleon $t$-matrix as well as the corresponding Bethe-Goldstone equation to determine the Brueckner reaction matrix in nuclear matter are reformulated in terms of the resolvents for the total two-nucleon Hamiltonians defined in free space and in medium correspondingly. This allows to find solutions at many energies simultaneously by using the respective Hamiltonian matrix diagonalization in the stationary wave packet basis. Among other important advantages, this approach simplifies greatly the whole computation procedures both for coupled-channel $t$-matrix and the Brueckner reaction matrix. Therefore this principally novel scheme is expected to be especially useful for self-consistent nuclear matter calculations because it allows to accelerate in a high degree single-particle potential iterations. Furthermore the method provides direct access to the properties of possible two-nucleon bound states in the nuclear medium. The comparison between reaction matrices found via the numerical solution of the Bethe-Goldstone integral equation and the straightforward Hamiltonian diagonalization shows a high accuracy of the method suggested. The proposed fully discrete approach opens a new way to an accurate treatment of two- and three-particle correlations in nuclear matter on the basis of three-particle Bethe-Faddeev equation by an effective Hamiltonian diagonalization procedure.