A general formulation of multiconfigurational self-consistent-field theory is given. The fundamental condition to be satisfied by the
M occupied orbitals ${\ensuremath{\psi}}_{\ensuremath{\mu}}$ of an $N(\ensuremath{\le}M)$-body system, in order that the total energy of one state be extremalized, is expressed in terms of the fundamental invariant $\mathbf{\ensuremath{\varrho}}=\ensuremath{\Sigma}{\ensuremath{\mu}=1}^{M}|{\ensuremath{\psi}}_{\ensuremath{\mu}}〉〈{\ensuremath{\psi}}_{\ensuremath{\mu}}|$. The condition placed on $\mathbf{\ensuremath{\varrho}}$ is of the same form as the condition on the fundamental invariant of the Hartree-Fock theory. Setting
M=N, we derive the unrestricted Hartree-Fock equations from the condition on $\mathbf{\ensuremath{\varrho}}$. For
M>N there are three freedoms one may take with the fundamental condition on $\mathbf{\ensuremath{\varrho}}$. Exploiting any or all of these freedoms yields alternative forms of the fundamental condition. This enables us to derive an effective one-body Hamiltonian which is the sum of a Hartree-type Hamiltonian and a correlation and exchange operator. For finite
M>N, the one-body Hamiltonian contains a nonlocal exchange and correlation operator. This operator is defined in terms of the one- and two-body density matrices. The connection between the orbitals of this theory and those of the Hartree-Fock theory is explored. The theory as outlined here is applicable to any system of
N identical particles, but our discussion is oriented towards electronic systems. The theory contains most self-consistent-field theories as special cases, and gives a basis for the self-consistent-field formulation of others.
The physical origin of the feature of strong correlation between electrons is mainly ascribed to doubly occupied atomic orbitals. Two-particle wave functions are believed necessary to describe strongly correlated systems. We emphasize that the on-site Coulomb interaction, which is important for strong correlation, should appear in equations as a natural term. Hartree-Fock-like self-consistent equations are derived with wave functions including two-particle wave functions for such systems. In the self-consistent equations, a strong correlation term is separated from the mean field. To our knowledge, this is the first attempt to formulate self-consistent equations in which the strong correlation effect is directly embodied. Some features of the equations are discussed. Koopmans' theorem and Brillouin's theorem are applied to the equations. The former is valid for singly occupied orbitals, and the latter is entirely available. By separating an angular correlation factor which is believed important for a strong correlation effect, an approximate form of a two-particle wave function is suggested. This kind of wave function causes the six multiple integrals to degrade to triple integrals. Under such forms of two-particle wave functions, corresponding Roothaan equations are established. Through these steps, it is believed that the equations become solvable within present computational limits.
A system of
N antisymmetric particles, moving under the influence of a fixed potential and their mutual many-particle interactions, is investigated in the ordinary Hartree-Fock scheme, having the total wave function approximated by a single Slater determinant. It is shown that all the density matrices of various orders, the wave function, and the entire physical situation depends only on a fundamental invariant $\ensuremath{\rho}({x}_{1}, {x}_{2})$, which is identical with the first-order density matrix. The Hartree-Fock equations are expressed in terms of this quantity.The Hartree-Fock equations are also solved by expanding the eigenfunctions in a given complete set, and applications to the MO-LCAO theory of the electronic structure of molecules, and crystals are given. It is shown that, in this scheme, the entire physical situation depends on a charge- and bond-order matrix $R(\ensuremath{\nu}\ensuremath{\mu})$ with respect to the ordinary atomic spin-orbitals involved. The Hartree-Fock equations for this matrix are investigated.Finally, the ionized and excited states are investigated, and it is shown that the Hartree-Fock scheme has a high degree of physical visuality also in case of many-particle interactions. The excitation energy of the system is the difference (${{\ensuremath{\omega}}_{j}}^{\ensuremath{'}}\ensuremath{-}{\ensuremath{\omega}}_{i}$) between two "spin-orbital energies," being eigenvalues to the effective Hamiltonians associated with the two states under consideration.
In order to calculate the average value of a physical quantity containing also many-particle interactions in a system of
N antisymmetric particles, a set of generalized density matrices are defined. In order to permit the investigation of the same physical situation in two complementary spaces, the Hermitean density matrix of order
k has two sets of indices of each
k variables, and it is further antisymmetric in each set of these indices.Every normalizable antisymmetric wave function may be expanded in a series of determinants of order
N over all ordered configurations formed from a basic complete set of one-particle functions ${\ensuremath{\psi}}_{k}$, which gives a representation of the wave function and its density matrices also in the discrete
k-space. The coefficients in an expansion of an eigenfunction to a particular operator may be determined by the variation principle, leading to the ordinary secular equation of the method of configurational interaction. It is shown that the first-order density matrix may be brought to diagonal form, which defines the "natural spin-orbitals" associated with the system. The situation is then partly characterized by the corresponding occupation numbers, which are shown to lie between 0 and 1 and to assume the value 1, only if the corresponding spin-orbital occurs in all configurations necessary for describing the situation. If the system has exactly
N spin-orbitals which are fully occupied, the total wave function may be reduced to a single Slater determinant. However, due to the mutual interaction between the particles, this limiting case is never physically realized, but the introduction of natural spin-orbitals leads then instead to a configurational expansion of most rapid convergence.In case the basic set is of finite order
M, the best choice of this set is determined by a form of extended Hartree-Fock equations. It is shown that, in this case, the natural spin-orbitals approximately fulfill some equations previously proposed by Slater.
A method for the calculation of translationally invariant wave functions for systems of identical fermions with arbitrary potential of pair interaction is developed.It is based on the well-known result that the essential dynamic part of Hamiltonian for the system of identical particles is the reduced Hamiltonian operator describing a relative movement of two particles inside the system.The eigenfunctions of this operator take into account all correlations caused by interaction.These eigenfunctions are basic elements for building the components (i.e., the functions with a lower degree of antisymmetry) used to construct the total antisymmetric wave function of the system.The main problem of this approach appears to be the antisymmetrization of the components.The developed universal algorithm for antisymmetrization gives a possibility to carry out this operation in a simple way and to keep numerical approximations under control.
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid
He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$.
F is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the
f(rij) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements
Hmn,
Nmn of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the
Hmn are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (
Hmn) and (
Nmn), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (
Hmn) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (
Hmn). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (
Hmn). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any
F in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{},
Q-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid
He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.