next up previous contents
Next: Scattering theory Up: Variational methods Previous: The helium atom   Contents

The hydrogen molecule ion

The hydrogen molecule ion consists of an electron orbiting about two protons, and is the simplest imaginable molecule. Let us investigate whether or not this molecule possesses a bound state: i.e., whether or not it possesses a ground-state whose energy is less than that of a hydrogen atom and a free proton. According to the variation principle, we can deduce that the $H_2^+$ ion has a bound state if we can find any trial wave-function for which the total Hamiltonian of the system has an expectation value less than that of a hydrogen atom and a free proton.

Figure 23: The hydrogen molecule ion.
\begin{figure}
\epsfysize =3.5in
\centerline{\epsffile{h2p.eps}}
\end{figure}

Suppose that the two protons are separated by a distance $R$. In fact, let them lie on the $z$-axis, with the first at the origin, and the second at $z=R$ (see Fig. 23). In the following, we shall treat the protons as essentially stationary. This is reasonable, since the electron moves far more rapidly than the protons.

Let us try

\begin{displaymath}
\psi({\bf r})_\pm = A\left[\psi_0({\bf r}_1) \pm \psi_0({\bf r}_2)\right]
\end{displaymath} (1202)

as our trial wave-function, where
\begin{displaymath}
\psi_0({\bf r}) = \frac{1}{\sqrt{\pi}\,a_0^{\,3/2}}\,{\rm e}^{-r/a_0}
\end{displaymath} (1203)

is a normalized hydrogen ground-state wave-function centered on the origin, and ${\bf r}_{1,2}$ are the position vectors of the electron with respect to each of the protons (see Fig. 23). Obviously, this is a very simplistic wave-function, since it is just a linear combination of hydrogen ground-state wave-functions centered on each proton. Note, however, that the wave-function respects the obvious symmetries in the problem.

Our first task is to normalize our trial wave-function. We require that

\begin{displaymath}
\int \vert\psi_\pm\vert^2\,d^3{\bf r} = 1.
\end{displaymath} (1204)

Hence, from (1202), $A = I^{-1/2}$, where
\begin{displaymath}
I = \int\left[\vert\psi_0({\bf r}_1)\vert^2 + \vert\psi_0({\...
...2 \pm
2\,\psi_0({\bf r}_1)\,\psi({\bf r}_2)\right] d^3{\bf r}.
\end{displaymath} (1205)

It follows that
\begin{displaymath}
I = 2\,(1\pm J),
\end{displaymath} (1206)

with
\begin{displaymath}
J = \int \psi_0({\bf r}_1)\,\psi_0({\bf r}_2)\,d^3{\bf r}.
\end{displaymath} (1207)

Let us employ the standard spherical polar coordinates ($r$, $\theta $, $\phi$). Now, it is easily seen that $r_1=r$ and $r_2=(r^2+R^2-2\,r\,R\,\cos\theta)^{1/2}$. Hence,

\begin{displaymath}
J = 2\int_0^\infty \int_0^\pi \exp\left[-x-(x^2+X^2-2\,x\,X\,\cos\theta)^{1/2}\right]\,x^2\,dx\,\sin\theta\,d\theta,
\end{displaymath} (1208)

where $X=R/a_0$. Here, we have already performed the trivial $\phi$ integral. Let $y=(x^2+X^2-2\,x\,X\,\cos\theta)^{1/2}$. It follows that $d(y^2)=2\,y\,dy = 2\,x\,X\,\sin\theta\,d\theta$, giving
$\displaystyle \int_0^\pi {\rm e}^{\,(x^2+X^2-2\,x\,X\,\cos\theta)^{1/2}}\,\sin\theta\,d\theta$ $\textstyle =$ $\displaystyle \frac{1}{x\,X}\int_{\vert x-X\vert}^{x+X}
{\rm e}^{-y}\,y\,dy$ (1209)
  $\textstyle =$ $\displaystyle - \frac{1}{x\,X}\left[{\rm e}^{-(x+X)}\,(1+x+X) - {\rm e}^{-\vert x-X\vert}\,(1+\vert x-X\vert)\right].$  

Thus,
$\displaystyle J$ $\textstyle =$ $\displaystyle - \frac{2}{X}\,{\rm e}^{-X}\int_0^X \left[{\rm e}^{-2\,x}\,(1+X+x)-
(1+X-x)\right]x\,dx$  
    $\displaystyle -\frac{2}{X}\int_X^\infty {\rm e}^{-2\,x}\left[{\rm e}^{-X}\,(1+X+x)-
{\rm e}^X\,(1-X+x)\right] x\,dx,$ (1210)

which evaluates to
\begin{displaymath}
J = {\rm e}^{-X}\left(1+X+\frac{X^3}{3}\right).
\end{displaymath} (1211)

Now, the Hamiltonian of the electron is written

\begin{displaymath}
H = -\frac{\hbar^2}{2\,m_e}\,\nabla^2 - \frac{e^2}{4\pi\,\epsilon_0}\left(\frac{1}{r_1}+\frac{1}{r_2}\right).
\end{displaymath} (1212)

Note, however, that
\begin{displaymath}
\left(-\frac{\hbar^2}{2\,m_e}\,\nabla^2 - \frac{e^2}{4\pi\,\...
...,2}}\right)\psi_0({\bf r}_{1,2}) = E_0\,\psi_0({\bf r}_{1,2}),
\end{displaymath} (1213)

since $\psi_0({\bf r}_{1,2})$ are hydrogen ground-state wave-functions. It follows that
$\displaystyle H\,\psi_\pm$ $\textstyle =$ $\displaystyle A\left[-\frac{\hbar^2}{2\,m_e}\,\nabla^2 - \frac{e^2}{4\pi\,\epsi...
...frac{1}{r_2}\right)\right]
\left[\psi_0({\bf r}_1) \pm \psi_0({\bf r}_2)\right]$  
  $\textstyle =$ $\displaystyle E_0\,\psi - A\,\left(\frac{e^2}{4\pi\,\epsilon_0}\right)\left[
\frac{\psi_0({\bf r}_1)}{r_2}\pm \frac{\psi_0({\bf r}_2)}{r_1}\right].$ (1214)

Hence,
\begin{displaymath}
\langle H\rangle = E_0 + 4\,A^2\,(D\pm E)\,E_0,
\end{displaymath} (1215)

where
$\displaystyle D$ $\textstyle =$ $\displaystyle \left\langle \psi_0({\bf r}_1)\left\vert\frac{a_0}{r_2}\right\vert\psi_0({\bf r}_1)\right\rangle,$ (1216)
$\displaystyle E$ $\textstyle =$ $\displaystyle \left\langle \psi_0({\bf r}_1)\left\vert\frac{a_0}{r_1}\right\vert\psi_0({\bf r}_2)\right\rangle.$ (1217)

Now,

\begin{displaymath}
D = 2\int_0^\infty\int_0^\pi \frac{{\rm e}^{-2\,x}}{(x^2+X^2-2\,x\,X\,\cos\theta)^{1/2}}\,x^2\,dx\,\sin\theta\,d\theta,
\end{displaymath} (1218)

which reduces to
\begin{displaymath}
D =\frac{4}{X}\,\int_0^X {\rm e}^{-2\,x}\,x^2\,dx
+ 4\int_X^\infty {\rm e}^{-2\,x}\,x\,dx,
\end{displaymath} (1219)

giving
\begin{displaymath}
D = \frac{1}{X} \left( 1-[1+X]\,{\rm e}^{-2\,X}\right).
\end{displaymath} (1220)

Furthermore,
\begin{displaymath}
E = 2\int_0^\infty \int_0^\pi \exp\left[-x-(x^2+X^2-2\,x\,X\,\cos\theta)^{1/2}\right]\,x\,dx\,\sin\theta\,d\theta,
\end{displaymath} (1221)

which reduces to
$\displaystyle E$ $\textstyle =$ $\displaystyle - \frac{2}{X}\,{\rm e}^{-X}\int_0^X \left[{\rm e}^{-2\,x}\,(1+X+x)-
(1+X-x)\right]dx$  
    $\displaystyle -\frac{2}{X}\int_X^\infty {\rm e}^{-2\,x}\left[{\rm e}^{-X}\,(1+X+x)-
{\rm e}^X\,(1-X+x)\right] dx,$ (1222)

yielding
\begin{displaymath}
E = (1+X)\,{\rm e}^{-X}.
\end{displaymath} (1223)

Our expression for the expectation value of the electron Hamiltonian is

\begin{displaymath}
\langle H\rangle = \left[1+ 2\,\frac{(D\pm E)}{(1\pm J)}\right] E_0,
\end{displaymath} (1224)

where $J$, $D$, and $E$ are specified as functions of $X=R/a_0$ in Eqs. (1211), (1220), and (1223), respectively. In order to obtain the total energy of the molecule, we must add to this the potential energy of the two protons. Thus,
\begin{displaymath}
E_{total} = \langle H\rangle + \frac{e^2}{4\pi\,\epsilon_0\,R} = \langle H\rangle - \frac{2}{X}\,E_0,
\end{displaymath} (1225)

since $E_0=-e^2/(8\pi\,\epsilon_0\,a_0)$. Hence, we can write
\begin{displaymath}
E_{total} = - F_\pm(R/a_0)\,E_0,
\end{displaymath} (1226)

where $E_0$ is the hydrogen ground-state energy, and
\begin{displaymath}
F_\pm(X) = -1 + \frac{2}{X}\left[\frac{(1+X)\,{\rm e}^{-2\,X...
...X^2/3)
\,{\rm e}^{-X}}{1\pm (1+X+X^2/3)\,{\rm e}^{-X}}\right].
\end{displaymath} (1227)

The functions $F_+(X)$ and $F_-(X)$ are both plotted in Fig. 24. Recall that in order for the $H_2^+$ ion to be in a bound state it must have a lower energy than a hydrogen atom and a free proton: i.e., $E_{total}< E_0$. It follows from Eq. (1226) that a bound state corresponds to $F_\pm < -1$. Clearly, the even trial wave-function $\psi_+$ possesses a bound state, whereas the odd trial wave-function $\psi_-$ does not [see Eq. (1202)]. This is hardly surprising, since the even wave-function maximizes the electron probability density between the two protons, thereby reducing their mutual electrostatic repulsion. On the other hand, the odd wave-function does exactly the opposite. The binding energy of the $H_2^+$ ion is defined as the difference between its energy and that of a hydrogen atom and a free proton: i.e.,
\begin{displaymath}
E_{bind} = E_{total} - E_0 = - (F_+ +1)\,E_0.
\end{displaymath} (1228)

According to the variational principle, the binding energy is less than or equal to the minimum binding energy which can be inferred from Fig. 24. This minimum occurs when $X\simeq 2.5$ and $F_+\simeq -1.13$. Thus, our estimates for the separation between the two protons, and the binding energy, for the $H_2^+$ ion are $R = 2.5\,a_0 = 1.33\times 10^{-10}\,{\rm m}$ and $E_{bind} = 0.13 \,E_0 = -1.77$ eV, respectively. The experimentally determined values are $R=1.06\times 10^{-10}$ m, and $E_{bind}=-2.8$ eV, respectively. Clearly, our estimates are not particularly accurate. However, our calculation does establish, beyond any doubt, the existence of a bound state of the $H_2^+$ ion, which is all that we set out to achieve.

Figure 24: The functions $F_+(X)$ (solid curve) and $F_-(X)$ (dashed curve).
\begin{figure}
\epsfysize =3.5in
\centerline{\epsffile{h2pa.eps}}
\end{figure}

next up previous contents
Next: Scattering theory Up: Variational methods Previous: The helium atom   Contents
Richard Fitzpatrick 2006-12-12