Spin projection of single-determinant wavefunctions

The components of the one- and two-particle density matrices resulting from spin projection of a single-determinant wave function are rederived by the method of expansion in terms of the natural orbitals of charge density, and the results are found to diverge from those given by Harriman and Sando for the two-particle case. The theory is generalized to include molecules with unequal numbers of electrons and basis orbitals and is applied to a number of organic molecules and ions. The correctness and internal consistency of the results argue in favor of the modification described.