\begin{equation}
\frac{1}{r_{12}} = \frac{1 - [\alpha + \beta erf (\mu r_{12})}{r_{12}}
                 + \frac{    [\alpha + \beta erf (\mu r_{12})}{r_{12}}
\end{equation}