For the first time the explicit form of the effective nonlocal potential for the natural spin orbitals is derived and analyzed. It is shown that in the case of the degenerate one-electron reduced density matrix the potential is not unique. The knowledge of the effective potential allows one to establish one-electron equations for the natural spin orbitals that may be of great value for efficient density matrix functional theory calculations.
For the first time the explicit form of the effective nonlocal potential for the natural spin orbitals is derived and analyzed. It is shown that in the case of the degenerate one-electron reduced density matrix the potential is not unique. The knowledge of the effective potential allows one to establish one-electron equations for the natural spin orbitals that may be of great value for efficient density matrix functional theory calculations.
We derive a variation equation for a time-dependent effective potential that is local both in time and in space and show that this potential generates and propagates spin orbitals with which a many-particle, time-dependent, multiconfigurational wave function $\mathrm{\ensuremath{\Psi}}(t)$ is constructed. Within this approach, the wave function and the effective potential are determined simultaneously in a self-consistent manner both in time-independent and in time-dependent cases. We also derive an equation that determines a real-valued effective potential and show that the equation establishes a relation among the first-order and second-order reduced density matrices and the electron repulsion integrals in the spin-orbital representation. By introducing the first-order density equation, we show that the variation equation for the effective potential can be simplified for an exact wave function. We also show that we can derive an effective potential with which a ground-state wave function that fulfills the Brillouin-Brueckner condition is constructed and that we can derive the effective potential proposed by Slater, with which a ground-state wave function represented by the multiconfiguration expansion is calculated, when an additional constraint is imposed on the varitation equation.
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.
The formulation of the density matrix functional theory (DMFT) in which the correlation component U of the electron–electron repulsion energy is expressed in terms of a model two-electron density cumulant matrix (a.k.a. the 2-cumulant) and the tensor of two-electron integrals g is critically analyzed. The dependence of on both the vector of occupation numbers n and g derived from the corresponding vector (1) of natural spinorbitals, and its concomitant invariance to the nature of the spin-independent interparticle interaction potential that enters the definition of g are emphasized. In the case of spin-polarized systems, (and thus U ) is found to be a function of not only n and g, but also of the vector and the matrix, where is the spin-flip operator. The presence of spin polarization imposes additional constraints upon, including a simple condition that, when satisfied, assures the underlying wavefunction being an eigenstate of both the and operators with the eigenvalues and, respectively. This feature allows targeting electronic states with definite multiplicities, which is virtually impossible in the case of the Kohn–Sham implementation of density functional theory. Among the four possible pairing schemes for the natural spinorbitals that give rise to approximations employing with only two independent indices, three are found to result in unphysical constraints even for spin-unpolarized systems, whereas the failure of the fourth one turns out to be precipitated by the presence of spin polarization. Consequently, any implementation of DMFT based upon “two-index” is shown to be generally unsuitable for spin-polarized systems (and incapable of yielding the spin-parallel components of U for the spin-unpolarized ones). A clear distinction is made between the genuine 1-matrix functionals that are defined for arbitrary N -representable 1-matrices and general energy expressions that depend on auxiliary quantities playing the role of fictitious 1-matrices subject to additional (often unphysical) constraints.
The Hohenberg–Kohn theorem implies the existence of an energy functional based solely on the first-order reduced density matrix of the ground state of an atomic or molecular system. Application of the variational principle for the functional EvM[γ] generates a set of coupled Euler equations for the representation coefficients and spin orbitals of a rank-M approximation to the exact ground-state density matrix. Defining the (assumed Hermitian) kernel FM[γ;x′,x]≡δEvM/δγ (x,x′), the equations in an arbitrary representation for the approximate density matrix, γ (x,x′) =ΣijMψ1(x) γijψj* (x′) are the following: ℱdξFM[γ;x′,ξ]ψi(ξ) =μψi(x′); i=1, 2,...,M; ℱdξFM[γ;x′,ξ] ΣMiψi(ξ) γ ij=ΣMiψi(x′) λij; j=1, 2, ...,M. The quantity μ is the chemical potential of the system of interest and the λij are a set of M2 Lagrange multipliers constraining the orthonormality of the spin orbital basis {ψ}. The coefficients γij must be chosen such that FM has the degenerate eigenvalue spectrum, FiiM=μ, i=1, 2,...,M, for all partially occupied orbitals. The spin orbitals, for a fixed set of coefficients, must be simultaneously determined so that λ is Hermitian. Provided all occupation numbers lie on the open interval (0,1), the stationary density matrix itself obeys the eigenvalue equation ℱdξFM[γ;x,ξ]γ (ξ,x′) =μγ (x,x′), and for any stationary density matrix the following commutation rules are valid: [FM[γ], γ]=0; [λ, γ]=0. The matrices γ and λ consequently can be brought simultaneously to diagonal form, and the canonical representation of the energy functional provides an eigenvalue equation determining the natural spin orbitals. The value of the chemical potential is μ=εi/ni, i=1, 2, ...,M; ε=diag (λ); from which follows the distribution function governing the occupation numbers of a stationary density matrix. Two limiting forms of the variational principle are examined, the exact and Hartree–Fock functionals, and related previous work by Gilbert is discussed. The physical content of the equations is illuminated by identification of the chemical potential as the negative of the electronegativity; Sanderson’s Principle of Electronegativity Equalization follows.