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.
A self-consistent field Hamiltonian is constructed, for arbitrary pair interactions, which yields (in principle) an exact value for the energy of a chosen state of a system of
N particles, as well as estimates for other states. The method for extending the formalism to interactions involving arbitrary numbers of particles is indicated and is straightforward. The method is carried out within the bound-state framework, and while motivated by the recent work of Bethe, results in a Hamiltonian somewhat different from his. An argument is presented which indicates that the independent-particle eigenfunctions of this Hamiltonian are "good" representations of the original
N-body system. A discussion of the practicality of the computations is included.
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.
Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-density-matrix-functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1-RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N≤ M≤N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1-RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many- electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized effective potential (OEP) method which allows for the use of explicitly orbital-dependent functionals in the context of CSDFT. We present the derivation of corresponding equations and show results for a quantum dot in external magnetic fields in the second part of this thesis.
Many-electron theory is treated here with much attention focused on the synthesis of the density functional and density matrix theory. This synthesis shows that the exchange–correlation potential, normally obtained as the functional derivative with respect to the electron density of an assumed approximate exchange–correlation energy, can be written exactly in terms of the “near-diagonal” behavior of the fully interacting first-order density matrix and the diagonal of the second-order density matrix, thereby bypassing the functional derivative. Also involved, since the exchange–correlation potential is an “artifact” of reducing the many-electron Schrödinger equation to one-body Slater–Kohn–Sham equations, are the one-body orbitals derived from these latter equations, not just the density thereby obtained. This situation has then motivated a discussion of what constraints are placed on the low-order density matrices of atoms and molecules, described by an assumed spin-independent Hamiltonian, by the requirement that these matrices are constructed from an N-electron wave function which is an eigenfunction of the spin operators S2 and Sz. Some briefer discussion is added on the optimally convergent determinantal expansion of such an N-electron wave function, the relation between Brueckner and Löwdin natural orbitals being commented on as a by-product. In this latter context, it is noted that the exchange-only potential of density functional theory can be represented as a sum of two line integrals, each of which is path-dependent, the sum being, of course, path-independent. But in addition to knowledge of Slater–Kohn–Sham orbitals appearing as in the Harbola–Sahni exchange-only potential (which is one of the path-dependent terms), there is a kinetic correction. This involves Löwdin natural orbitals and their occupation numbers, but only now in the weak coupling limit of the electron–electron interaction, as emphasized in the kinetic “correction” to the Harbola–Sahni formula exhibited by Levy and the writer. Some directions for further work are briefly set out. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 163–176, 1999
We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.