Coupled-cluster method in Fock space. I. General formalism

Abstract

The problem of finding the spectrum of the Fock-space Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{} for a system of many fermions is analyzed. The quasiparticle formalism is employed, with ``holes'' and ``particles'' defined with respect to some N-particle determinantal wave function (the model vacuum). The basic idea behind the proposed approach is to perform a similarity transformation of Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{}, such that the resulting effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{} is, unlike \ifmmode \hat{H}\else \^{H}\fi{}, a quasiparticle-number--conserving operator. It is shown that eigenvalues of \ifmmode \hat{G}\else \^{G}\fi{}, corresponding to small numbers of quasiparticles (0,1,2) can be easily calculated. This is equivalent to finding eigenvalues of \ifmmode \hat{H}\else \^{H}\fi{} for certain states of N, N\ifmmode\pm\else\textpm\fi{}1, and N\ifmmode\pm\else\textpm\fi{}2 particles. The construction of the operator transforming \ifmmode \hat{H}\else \^{H}\fi{} into \ifmmode \hat{G}\else \^{G}\fi{} (the wave operator) stems from an analysis of the structure of the algebra of linear operators acting in a (finite-dimensional) Fock space. The exponential Ansatz for the wave operator is used, resulting in a generalization of the coupled-cluster (CC) method of Coester [Nucl. Phys. 7, 421 (1958)]. The generalized CC equations determining the wave operator, and equations determining the effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{}, are presented in a diagrammatic form. An effort has been made to obtain a concise notation for expressing these equations in an algebraic form. Approximation schemes, necessary for practical applications of the proposed method, are also studied.