We develop a procedure recently proposed for the treatment of correlations in finite Fermi systems in a more consistent way than the random phase approximation (RPA). The one-body density matrix in the correlated ground state is now no more assumed to be diagonal in the Hartree-Fock basis. An expression is rather provided where this quantity is expressed in terms of the X and Y amplitudes of the particle-hole phonon operators. A set of RPA-like equations depending only on the one-body density matrix is then constructed by means of the method of linearized equations of motion and solved via an iterative procedure. The procedure is employed to calculate spectra and strength distributions of several neutral and charged Na and Li metal clusters within the jellium approximation. Improvements are observed with respect to standard RPA calculations.
We present a method for the treatment of correlations in finite Fermi systems in a more consistent way than the random phase approximation (RPA). The main differences with respect to previous approaches are underlined and discussed. In particular, no use of renormalized operators is made. By means of the method of linearization of equations of motion, we derive a set of RPA-like equations depending only on the one-body density matrix. The latter is no more assumed to be diagonal and it is expressed in terms of the
X and
Y amplitudes of the particle-hole phonon operators. This set of nonlinear equations is solved via an iterative procedure, which allows us to calculate the energies and the wave functions of the excited states. After presenting our approach in its general formulation, we test its quality by using metal clusters as a good test laboratory for a generic many-body system. Comparison to RPA shows significant improvements. We discuss also how the present approach can be further extended beyond the particle-hole configuration space, getting an approximation scheme in which energy weighted sum rules are exactly preserved, thus solving a problem common to all extension of RPA proposed until now. The implementation of such approach will be afforded in a future work.
The role of the small-norm amplitudes in extended random-phase approximation (RPA) theories such as the particle-particle and hole-hole components of one-body amplitudes and the twobody amplitudes other than two-particle/two-hole components are investigated for the onedimensional Hubbard model using an extended RPA derived from the time-dependent density matrix theory.It is found that these amplitudes cannot be neglected in strongly interacting regions where the effects of ground-state correlations are significant.............