When we treat an atom containing a number n of electrons by the method of the self-consistent field, we assume that each electron has its own particular “orbit,” specified by a wave function ( q | r ) in four variables q . These four variables are usually taken to be the three coordinates of the electron together with a variable describing the spin, but according to the transformation theory of quantum mechanics, they may be any four independent commuting functions of the coordinates, momenta and spin variables.
When we treat an atom containing a number n of electrons by the method of the self-consistent field, we assume that each electron has its own particular “orbit,” specified by a wave function ( q | r ) in four variables q . These four variables are usually taken to be the three coordinates of the electron together with a variable describing the spin, but according to the transformation theory of quantum mechanics, they may be any four independent commuting functions of the coordinates, momenta and spin variables.
Abstract The spin-independent properties of an electronic system are completely determined by spinless density functions, P1(1; 1') and P 2( 1,2; 1',2 '), derived from the 1-and 2-electron density matrices, p1(1; 1') and p2(1,2 ; 1',2'). These functions, and their physical inter-pretation, have been dealt with in part I: the present paper contains the generalizations necessary when the Hamiltonian contains spin terms (as in discussions of electron and nuclear spin resonance). The density matrices, and more generally the transition density matrices p1(kλ|1; 1') and p2(kλ|l,2; 1',2') associated with states and Φkand Φλ are resolved into spinless components associated with different spin situations: p1 has four such components (one giving the density of a-spin electrons, a second the density of β-spin electrons, and two with less direct physical significance), while p2 has sixteen components. It is shown that all components, for states k=(S, M),λ= (S',M') of given S, S', are determined by three 2-electron functions.For k=λ these reduce to p2, the spinless pair function, from which the electron density p1 may be derived by integration; Q2a ‘conditional' spin density, from which the (1-electron) spin density Q1 may be derived (as in the case of p1 and P2); Q c. a., a ‘coupling anisotropy function', which describes the anisotropy of spin-spin coupling associated with quantization of the spin z-component. It is also useful to introduce (though it is not an independent function). Qc, a ‘coupling density function', which measures the coupling between the spins of two electrons at different points in space. All spin-dependent effects, including the usual spin-spin and spin-orbit couplings and also (indirect) nuclear and electron-nuclear spin-spin couplings, may be related to the many- electron wave function through the above density functions. All the densities have a physical interpretation, independent of the precise nature of the wave functions (which may, for example, be the exact eigenfunctions of a spinless Hamiltonian or merely approximations of orbital form). The derivation of the density functions from approximate wave functions is considered in some detail. For functions describing weakly coupled electron groups (e. g. inner-shell and valence electrons) the densities may be expressed in terms of those of the separate groups. This reduction may be followed by an orbital description of each group; and subsequently, by a linear-combination-of-atomic-orbitals (l. c. a. o.) approximation. In the simpler approximations, all the coupling effects discussed are determined by electron and spin ‘populations'. Preliminary calculations suggest that such approximations may in many cases be adequate.
The vector model is examined in the light of the quantum theory of wave fields. First the distributions of electrons among the electronic states are limited to correspond to the atomic states arising from one electronic configuration. A set of independent angular momentum vectors 11, 12,...S1, S2... is found, in which one orbital and one spin vector are associated with each group of electronic states, rather than individual electrons, having given values of the quantum numbers n and l. Any other vector T may be decomposed into a sum of independent vectors ti, one associated with each of the same group of electronic states. The usual commutation relations obtain among the vectors ti, li and si. On further limiting the distributions to correspond to atomic states arising from electronic configurations containing no equivalent electrons, it is shown that the matrix of ${\ensuremath{\Sigma}}_{k}f({{\mathrm{l}}_{k}}^{0}, {{\mathrm{s}}_{k}}^{0})$ (an operator in coordinate space in which f is any algebraic function of the components of 1 and s) can be constructed from the vectors 11, 12...S1, S2,... according to the usual laws of matrix mechanics.
For systems which obey quantum mechanics, the formulation of the problem of treating large numbers of particles is, of course, somewhat different than it is for classical systems. The microscopic description of the system is provided by a wave function which (in the absence of spin) is a function of the classical coordinates { q i }. The mathematical model is provided by a Hamiltonian operator H which is often obtained from the corresponding classical Hamiltonian by the replacement p i → ( ħ /i)(∂/∂ q i ). In other cases the form of the Hamiltonian operator is simply postulated. The microscopic dynamics are provided by the Schrödinger equation i ħ (∂Ψ/∂ t ) = H Ψ which requires as initial condition the knowledge of the wave function Ψ({ q i }, t ) at some initial time t 0 . (Boundary conditions on Ψ({ q i }, t ) must be specified as part of the description of the model as well.) The results of experiments in quantum mechanics are characterized by operators, usually obtained, like the Hamiltonian, from their classical forms and termed observables. Operators associated with observables must be Hermitian. In general, the various operators corresponding to observables do not commute with one another. It is possible to find sets of commuting operators whose mutual eigenstates span the Hilbert space in which the wave function is confined by the Schrödinger equation and the boundary conditions. A set of such (time independent) eigenstates, termed ψ ν ( q ), is a basis for the Hilbert space.
In two recent papers Dirac has shown how the “duplexity” phenomena of the atom can be accounted for without recourse to the hypothesis of the spinning electron. The investigation is carried out by the methods of non-commutative algebra, the wave function ψ being a matrix of the fourth order. An alternative presentation of the theory, using the methods of wave mechanics, has been given by Darwin. The four-rowed matrix ψ is replaced by four wave functions ψ1, ψ2, ψ3, ψ4 satisfying four linear differential equations of the first order. These functions are related to one particular direction, and the work can only be given invariance of form at the expense of much additional complication, the four wave functions being replaced by sixteen.
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.