We analyze the structure and the solutions of the irreducible k-particle Brillouin conditions (IBCk) and the irreducible contracted Schrödinger equations (ICSEk) for an n-electron system without electron interaction. This exercise is very instructive in that it gives one both the perspective and the strategies to be followed in applying the IBC and ICSE to physically realistic systems with electron interaction. The IBC1 leads to a Liouville equation for the one-particle density matrix gamma1=gamma, consistent with our earlier analysis that the IBC1 holds both for a pure and an ensemble state. The IBC1 or the ICSE1 must be solved subject to the constraints imposed by the n-representability condition, which is particularly simple for gamma. For a closed-shell state gamma is idempotent, i.e., all natural spin orbitals (NSO's) have occupation numbers 0 or 1, and all cumulants lambdak with k> or =2 vanish. For open-shell states there are NSO's with fractional occupation number, and at the same time nonvanishing elements of lambda2, which are related to spin and symmetry coupling. It is often useful to describe an open-shell state by a totally symmetric ensemble state. If one wants to treat a one-particle perturbation by means of perturbation theory, this mainly as a run-up for the study of a two-particle perturbation, one is faced with the problem that the perturbation expansion of the Liouville equation gives information only on the nondiagonal elements (in a basis of the unperturbed states) of gamma. There are essentially three possibilities to construct the diagonal elements of gamma: (i) to consider the perturbation expansion of the characteristic polynomial of gamma, especially the idempotency for closed-shell states, (ii) to rely on the ICSE1, which (at variance with the IBC1) also gives information on the diagonal elements, though not in a very efficient manner, and (iii) to formulate the perturbation theory in terms of a unitary transformation in Fock space. The latter is particularly powerful, especially, when one wishes to study realistic Hamiltonians with a two-body interaction.
We analyze the structure and the solutions of the irreducible k-particle Brillouin conditions (IBCk) and the irreducible contracted Schrödinger equations (ICSEk) for an n-electron system without electron interaction. This exercise is very instructive in that it gives one both the perspective and the strategies to be followed in applying the IBC and ICSE to physically realistic systems with electron interaction. The IBC1 leads to a Liouville equation for the one-particle density matrix gamma1=gamma, consistent with our earlier analysis that the IBC1 holds both for a pure and an ensemble state. The IBC1 or the ICSE1 must be solved subject to the constraints imposed by the n-representability condition, which is particularly simple for gamma. For a closed-shell state gamma is idempotent, i.e., all natural spin orbitals (NSO's) have occupation numbers 0 or 1, and all cumulants lambdak with k> or =2 vanish. For open-shell states there are NSO's with fractional occupation number, and at the same time nonvanishing elements of lambda2, which are related to spin and symmetry coupling. It is often useful to describe an open-shell state by a totally symmetric ensemble state. If one wants to treat a one-particle perturbation by means of perturbation theory, this mainly as a run-up for the study of a two-particle perturbation, one is faced with the problem that the perturbation expansion of the Liouville equation gives information only on the nondiagonal elements (in a basis of the unperturbed states) of gamma. There are essentially three possibilities to construct the diagonal elements of gamma: (i) to consider the perturbation expansion of the characteristic polynomial of gamma, especially the idempotency for closed-shell states, (ii) to rely on the ICSE1, which (at variance with the IBC1) also gives information on the diagonal elements, though not in a very efficient manner, and (iii) to formulate the perturbation theory in terms of a unitary transformation in Fock space. The latter is particularly powerful, especially, when one wishes to study realistic Hamiltonians with a two-body interaction.
The k-particle irreducible Brillouin conditions IBCk and the k-particle irreducible contracted Schrödinger equations ICSEk for a closed-shell state are analyzed in terms of a Møller-Plesset-type perturbation expansion. The zeroth order is Hartree-Fock. From the IBC2(1), i.e., from the two-particle IBC to first order in the perturbation parameter mu, one gets the leading correction lambda2(1) to the two-particle cumulant lambda2 correctly. However, in order to construct the second-order energy E(2), one also needs the second-order diagonal correction gammaD(2) to the one-particle density matrix gamma. This can be obtained: (i) from the idempotency of the n-particle density matrix, i.e., essentially from the requirement of n-representability; (ii) from the ICSE1(2); or (iii) by means of perturbation theory via a unitary transformation in Fock space. Method (ii) is very unsatisfactory, because one must first solve the ICSE3(2) to get lambda3(2), which is needed in the ICSE2(2) to get lambda2(2), which, in turn, is needed in the ICSE1(2) to get gamma(2). Generally the (k+1)-particle approximation is needed to obtain Ek correctly. One gains something, if one replaces the standard hierarchy, in which one solves the ICSEk, ignoring lambda(k+1) and lambda(k+2), by a renormalized hierarchy, in which only lambda(k+2) is ignored, and lambda(k+1) is expressed in terms of the lambdap of lower particle rank via the partial trace relation for lambda(k+2). Then the k-particle approximation is needed to obtain E(k) correctly. This is still poorer than coupled-cluster theory, where the k-particle approximation yields E(k+1). We also study the possibility to use some simple necessary n-representability conditions, based on the non-negativity of gamma(2) and two related matrices, in order to get estimates for gammaD(2) in terms of lambda2(1). In general these estimates are rather weak, but they can become close to the best possible bounds in special situations characterized by a very sparse structure of lambda2 in terms of a localized representation. The perturbative analysis does not encourage the use of a k-particle hierarchy based on the ICSEk (or on their reducible counterparts, the CSEk), it rather favors the approach in terms of the unitary transformation, where the k-particle approximation yields the energy correct up to E(2k-1). The problems that arise are related to the unavoidable appearance of exclusion-principle violating cumulants. The good experience with perturbation theory in terms of a unitary transformation suggests that one should abandon a linearly convergent iteration scheme based on the ICSEk hierarchy, in favor of a quadratically convergent one based on successive unitary transformations.
Recently [W. Kutzelnigg and D. Mukherjee, Chem. Phys. Lett. 317, 567 (2000); D. Mukherjee and W. Kutzelnigg, J. Chem. Phys 114, 2047 (2001)] the irreducible k-particle Brillouin conditions IBCk and the irreducible k-particle contracted Schrödinger equations ICSEk were derived. These permit the definition of a hierarchy of k-particle approximations for the direct calculation of the cumulants λk of the k-particle density matrices. Now, the spin-free form of these conditions, appropriate for a spin-free Hamiltonian, is given. This is particularly useful for open-shell states. The definition of the cumulants of the reduced densities has to be generalized for these anyway, making use of irreducible tensor operators with respect to SU2. There are two alternative definitions of spin-free cumulants, of which the one in terms of spin-free reduced density matrices appears to be preferable. Alternatively to the straight spin-free formulation, we also present a theory in terms of spin-free operators adapted to the symmetric group. Partial trace relations that relate the elements of the cumulants of different particle rank are derived. There are partial trace relations for “exchange elements,” which are determined by the total spin quantum number S. From these relations the individual exchange elements of Λ2 can be obtained in special cases. This allows a simple formulation of the stationarity conditions for open-shell states of any spin multiplicity.
Two alternative conditions for the stationarity of the energy expectation value with respect to k-particle excitations are the k-particle Brillouin conditions BCk and the k-particle contracted Schrödinger equations, CSEk. These conditions express the k-particle density matrices γk in terms of density matrices of higher particle rank. The latter can be eliminated if one expresses the γk in terms of their cumulants λk, but this is not sufficient to make the BCk or CSEk separable (extensive), i.e., they are not expressible in terms of only connected diagrams. However, in a formulation based on the recently introduced general normal ordering with respect to arbitrary wave functions, the irreducible counterparts IBCk and ICSEk of the BCk and CSEk can be defined. They are easily evaluated explicitly in terms of the generalized Wick theorem for arbitrary wave functions, and they lead to equations for the direct construction of the cumulants λk, which are additively separable quantities and which scale linearly with the system size. The IBCk or the ICSEk are necessary conditions for γ and the λk to represent an exact n-fermionic eigenstate of the given Hamiltonian. To specify the desired state, additional conditions must be satisfied as well, e.g., the partial trace relations which relate λ2 to γ and γ2. The particle number and the total spin must be specified and n-representability conditions enter implicitly. While the nondiagonal elements of γ and the λk are determined by the IBCk or the ICSEk, the additional conditions mainly serve to fix the diagonal elements. A hierarchy of k-particle approximations is defined. It is based on the fact that the expansion in terms of cumulants λk can be truncated at any particle rank, which would not be possible for the density matrices γk. For closed-shell states the one-particle approximation agrees with Hartree–Fock.
The Fock space Hamiltonian H has a simpler structure than its projection Hn to n-particle Hilbert space. It is therefore recommended to diagonalize H—to the extent that this is possible—before one specifies n. Diagonalization of H is possible if one defines the diagonal and nondiagonal parts of an operator appropriately. It is shown that the diagonalized Fock space Hamiltonian L, called an energy operator, contains all information about the eigenvalues of H in a simply coded form. For a spinfree Hamiltonian, spin can be completely eliminated and all interesting quantities are expressible in terms of spinfree excitation operators (generators of the unitary group). The construction of the wave operator W and the energy operator L is formulated in terms of perturbation theory both in the strictly degenerate and the quasidegenerate version. Three variants are discussed that differ in the normalization of W, namely, a intermediate normalization, b and c unitary normalization with two additional conditions, b: WD=WD+, c:σD=(ln W)D=0. The transformation of operator products to normal product form is achieved by means of a generalization of Wick’s theorem to spinfree quantum chemistry. The contributions to the various orders of L, W, and σ are illustrated by means of diagrams. Only variant c guarantees a connected-diagram expansion of the Fock space energy operator L. None of the three variants satisfies a generalization of the Wigner (2n+1) rule, but with a slight change in the normalization from b to b′ or c to c′ (that does, in the case of c′, not affect the connectedness of the diagram expansion) L(2n) and L(2n+1) are expressible in terms of the W(k), k⩽n. A generalization of the variation principle for the Fock space energy operator is given and its consequences are discussed. Generalized Fock space Brillouin conditions are derived and a variational scheme for the construction of the energy operator is proposed. We finally compare the present Fock space approach with the traditional degenerate perturbation theory in n-particle Hilbert space. If the diagonal part ΩD of an operator Ω is defined in the spirit of the theory of the universal wave and energy operators, then the two approaches lead to the same results (for the same normalization), but the Fock space approach looks simpler and more transparent.
The density operator (density matrix) of a quantum mechenical system can be decomposed into operators which transform as irreducible representations of the symmetry group in coordinate and spin space. Each of these components has a physical meaning connected with the expectation values of certain operators. The reduced density matrices can be decomposed in a completely analogous way. The symmetry properties of the total wave function give rise to degeneracies of the eigenvalues of the reduced density matrices. These degeneracies can be removed by requiring that the natural spin orbitals (NSO, defined as the eigenfunctions of the first-order density matrix), as well as the natural spin geminals (NSG, the eigenfunctions of the second-order density matrix) and their spinless counterparts transform as irreducible representations of the symmetry group and are eigenfunctions of S/sup 2/ and S/sub z/. In many cases this requirement is compatible with the original definition of the NSO, the NSG, etc., e.g., when there is no spatial degeneracy of the total wave function and when the Z- component of the total spin vanishes. When these conditions are not fulfilled an alternative definition of the NSO and the NSG is proposed. (auth)
The method of superposition of configurations, which provides a general solution of the quantum-mechanical many-particle problem (for fermions), is reformulated so that equations may be compared with those used in the Brueckner theory. Important differences occur in application to finite or nonuniform systems for which the Hartree-Fock or Brueckner self-consistent orbitals are not plane waves. In such cases nonvanishing single-particle matrix elements occur which cannot be described by the Brueckner formalism based on a two-particle operator. Equations for an effective two-particle operator, equivalent to the variational equations of the method of superposition of configurations, are derived for a basis of Hartree-Fock orbitals. At the expense of making the orbital basis dependent upon the effective two-particle operator, the orbital basis can be determined by a condition which is essentially that of the Brueckner method. This condition removes a class of matrix elements which do not necessarily vanish in the Hartree-Fock basis although they would be neglected to second order in a perturbation calculation starting from the Hartree-Fock wave function. The equations for the effective two-particle operator are formally the same in both cases but lead to different operators because of the different choice of basis. In neither case can the equations be written in terms of products of such operators, a formalism assumed in the Brueckner theory. It is shown that the Brueckner condition is not equivalent to the condition which would determine the best possible orbital basis for the form of wave function implied by the use of an effective operator dependent on two particles only. The argument of the present paper is limited to systems with a finite number of particles, since the variational principle used here is not applicable without modification to systems with an infinite number of particles.