Find this document at http://www.ncsa.uiuc.edu/Classes/MATH198/kaduk
Molecular Orbital theory attempts to determine the values of the wavefunction for electrons that are bound to a molecular framework. These values are useful, because they can predict reactivity and structure of the molecular species in an empirical way. However, the Schrodenger equation for all but the most simple cases is impossible to solve explicitly; its solutions can only be approximated. Almost all MO (molecular orbital) calculations rely on the Born-Oppenheimer approximation: since the nuclei are at least 1000 times more massive than the electrons, they are relatively motionless on the timescale of electron motion. This simplifies the calculations immensely and also allows for a very simplistic approximation that ignores certain aspects of the differential equations involved, but is still reasonably accurate, and is (most notably) very simple to implement, and this is the strategy that I have chosen to implement. Since the wavefunctions for single atoms can be explicitly calculated, we assume that the wavefunction for the molecule as a whole is merely a linear combination of these atomic wavefunctions, and calculate the formula for this sum. This method can be parametrized fairly easily, allowing for user-variable parameters, such as the interatomic distance and the cutoff value used in determining the boundary of the orbital as a level set of the probability density function.