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)
For a quantum-mechanical system of
N identical fermions, the
N-representability problem is the problem of recognizing whether, for a given pth-order reduced density matrix ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})$, there exists an antisymmetric
N-particle wave function $\ensuremath{\Psi}(12\ensuremath{\cdots}N)$ such that ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})=({N}{p})\ensuremath{\int}\ensuremath{\Psi}(12\ensuremath{\cdots}N)\ifmmode\times\else\texttimes\fi{}{\ensuremath{\Psi}}^{*}({1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}},p+1\ensuremath{\cdots}N){d}_{{\ensuremath{\tau}}_{p+1}}\ensuremath{\cdots}{d}_{{\ensuremath{\tau}}_{N}}$. It is shown that if the Hamiltonian of a system is time-reversal invariant, and the number of particles,
N, is even, the necessary and sufficient condition that an approximate first-order density matrix corresponding to a nondegenerate energy eigenstate be
N-representable is that its natural spin-orbital occupation numbers be equal in pairs.
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 Tribute to A. John Coleman - The Tame Mathematician.- Reduced Density Matrices: 1929-1989.- Some Aspects on the Development of the Theory of Reduced Density Matrices and the Representability Problem.- Representability Conditions.- On the Diagonal N-Representability Problem.- Fermion N-Representability Conditions Generated by a Decomposition of the 1-Particle Identity Operator onto Mutually Orthogonal Projection Operators.- The Unitarily Invariant Decomposition of Hermitian Operators.- Building Up N-Electron States with Symplectic Symmetry.- Time Dependent Antisymmetrized Geminal Power Theory Using a Coherent State Formulation.- Griffiths Inequalities for Fermion Systems.- Entropy of Reduced Density Matrices.- A Lower Bound to the Ground State Energy of a Boson System With Fermion Source.- Reduced Density Operators, Their Related von Neumann Density Operators, Close Cousins of These, and Their Physical Interpretation.- Theory and Practice of the Spin-Adapted Reduced Hamiltonians (SRH).- Variational Principle with Built-In Pure State N-Representability Conditions. The N-Electron Case.- Wigner Distributions as Representations of the Density Matrix.- Inter-Relationships Between Various Representations of One-Matrices and Related Densities: A Road Map and an Example ..- Current Problems in Density Functional Theory.- The Interface Between Reduced Density Matrices and Density Functional Theory.- The Physics Underlying the Langreth-Mehl Scheme for Non-Uniform Systems.- Understanding Energy Differences in Density Functional Theory.- Density Functional Calculations of Molecular Bond Energies.- Non-Local Effects on Atomic and Molecular Correlation Energies Studied with a Gradient-Corrected Density Functional.- An Evaluation of Local Electron Correlation Corrections and Non-Local Exchange Corrections to the Hartree-Fock-Slater Method from Calculations on Bond Energies and Electronic Spectra of Molecular Systems.- Correlation Energy Functionals of One-Matrices and Hartree-Fock Densities.- Some Remarks on Scaling Relations in Density Functional Theory.- Deduction of Semiempirical MO Methods from Density Functional Theory.- Charge and Spin Densities in Molecular Solids: Local Density Functional Calculations Versus Experiment.- A Functional of the Two-Particle Density Matrix for the Approximate Calculation of the Electronic Correlation Energy.- Extracules, Intracules, Correlation Holes, Potentials, Coefficients and All That.- The Exact Schroedinger Equation for the Electron Density.- Adiabatic Separation, Broken Symmetries and Geometry Optimization.- Asymptotic Results for Density Matrices and Electron Density in Atoms and Nearly Spherical Molecules.- An Algorithm for Calculating Isoelectronic Changes in Energies, Densities, and One-Matrices.- Atoms and Ions in the Limit of Large Nuclear Charge.- Improved Thomas-Fermi Theory for Atoms.- A Bond Energy from Quantum Mechanics.- Measured Electron Densities and Band Structure Calculations.- X-ray Orthonormal Orbital Model for Crystallography.
Abstract A partial trace over the occupation numbers of all but k states in the density matrix of an ensemble with an arbitrary number of single‐particle states is defined as the (reduced) k‐state density matrix. This matrix is used to obtain a complete, practical solution to the problem of determining the representability of the diagonal elements of the one‐ and two‐particle (reduced) density matrices. This solution is expressed as a series of linear inequalities involving the density‐matrix elements; the inequalities are identical with those derived previously by Davidson and McCrae by a different method. In addition, our method is used to obtain nonlinear, matrix inequalities on the off‐diagonal elements of the density matrices.