Modern valence bond theory is the application of valence bond theory, with computer programs that are competitive in accuracy and economy with programs for the Hartree–Fock method and other molecular orbital based methods. The latter methods dominated quantum chemistry from the advent of digital computers because they were easier to program. The early popularity of valence bond methods thus declined. It is only recently that the programming of valence bond methods has improved. These developments are due to and described by Gerratt, Cooper, Karadakov and Raimondi (1997); Li and McWeeny (2002); Joop H. van Lenthe and co-workers (2002); Song, Mo, Zhang and Wu (2005); and Shaik and Hiberty (2004).
In its simplest form the overlapping atomic orbitals are replaced by orbitals which are expanded as linear combinations of the atom-based basis functions, forming linear combinations of atomic orbitals (LCAO). This expansion is optimized to give the lowest energy. This procedure gives good energies without including ionic structures.
For example, in the hydrogen molecule, classic valence bond theory uses two 1s atomic orbitals (a and b) on the two hydrogen atoms respectively and then constructs a covalent structure:-
and then an ionic structure:-
The final wave function is a linear combination of these two functions. Coulson and Fischer pointed out that a completely equivalent function is:-
as expanding this out gives a linear combination of the covalent and ionic structures. Modern valence bond theory replaces the simple linear combination of the two atomic orbitals with a linear combination of all orbitals in a larger basis set. The two resulting valence bond orbitals look like an atomic orbital on one hydrogen atom slightly distorted towards the other hydrogen atom. Modern valence bond theory is thus an extension of this Coulson–Fischer method.