Molecular Orbital Studies of Diatomic Molecules. I. Method of Computation for Single Configurations of Heteronuclear Systems

This paper describes a method of molecular orbital calculation for heteronuclear diatomic molecules, in the approximation of a single spatial configuration which is a product of one‐electron functions. The method is applicable to systems in which the one‐electron orbitals are not required to be either orthogonal or filled in pairs. A practical procedure is presented for treating permutational symmetry and spin for an arbitrary spatial wave function, and using this procedure we consider simultaneously all the spin configurations associated with the chosen spatial wave function. The wave functions used are parametric expressions in spheroidal coordinates, and include the usual Slater orbitals as special cases. All necessary integrals are evaluated in a manner suitable for machine computation. A computer program to implement this method of calculation has been prepared, and with it an approximate energy level for a 6‐electron molecule can be obtained in about 8 min. Provision was made to program the computer to vary the parameters of the wave functions so as to find automatically the wave function of minimum energy of the parametric form chosen. Only a part of the calculation need be repeated after varying parameters, and the optimum retention of computed data for re‐use is discussed.