Linear combination of atomic orbitals

The LCAO approximation builds a molecular orbital as a linear combination of single-centre atomic orbitals. For two equivalent normalized real $1s$ functions ${\varphi }_A$ and ${\varphi }_B$, choose their signs positive and define their orbital overlap $S$. At nonzero separation, $0<S<1$, and the normalized combinations areHere ${\sigma }_{+}={\sigma }_g$ and ${\sigma }_{-}={\sigma }_u^{\ast }$. Equal coefficients use the equivalence of the two centres; a general LCAO need not have equal coefficients. See Field, MIT 5.61, lecture 24, pp. 3–6: LCAO coefficients and overlap integrals.

For the Gaussian function toy's $a=1$, the two denominators are $\sqrt{2(1+S)}\approx 1.654013$ and $\sqrt{2(1-S)}\approx 1.124385$. To check wavefunction normalization, expand the square: $\int ({\varphi }_A\pm {\varphi }_B{)}^2=1+1\pm 2S$. Dividing by $2(1\pm S)$ gives one. The local probability density isThe cross term is the interference; adding the separate densities first would lose it. These formulas normalize the chosen linear combination exactly, while the restricted choice of atomic orbitals remains an approximation to an actual molecular orbital.