# 6.1: The Schrodinger Equation for the Hydrogen Atom Can Be Solved Exactly

- Page ID
- 210823

The hydrogen atom, consisting of an electron and a proton, is a two-particle system, and the internal motion of two particles around their center of mass is equivalent to the motion of a single particle with a reduced mass. This reduced particle is located at \(r\), where \(r\) is the vector specifying the position of the electron relative to the position of the proton. The length of \(r\) is the distance between the proton and the electron, and the direction of \(r\) and the direction of \(r\) is given by the orientation of the vector pointing from the proton to the electron. Since the proton is much more massive than the electron, we will assume throughout this chapter that the reduced mass equals the electron mass and the proton is located at the center of mass.

Since the internal motion of any two-particle system can be represented by the motion of a single particle with a reduced mass, the description of the hydrogen atom has much in common with the description of a diatomic molecule discussed previously. The Schrödinger Equation for the hydrogen atom

\[ \hat {H} (r , \theta , \varphi ) \psi (r , \theta , \varphi ) = E \psi ( r , \theta , \varphi) \label {6.1.1}\]

employs the same kinetic energy operator, \(\hat {T}\), written in spherical coordinates. For the hydrogen atom, however, the distance, \(r\), between the two particles can vary, unlike the diatomic molecule where the bond length was fixed, and the rigid rotor model was used. The hydrogen atom Hamiltonian also contains a potential energy term, \(\hat {V}\), to describe the attraction between the proton and the electron. This term is the Coulomb potential energy,

\[ \hat {V} (r) = - \dfrac {e^2}{4 \pi \epsilon _0 r } \label {6.1.2}\]

where r is the distance between the electron and the proton. The Coulomb potential energy depends inversely on the distance between the electron and the nucleus and does not depend on any angles. Such a potential is called a central potential.

The full expression for \(\hat {H}\) in spherical coordinates is

\[\hat {H} (r , \theta , \varphi ) = - \dfrac {\hbar ^2}{2 \mu r^2} \left [ \dfrac {\partial}{\partial r} \left (r^2 \dfrac {\partial}{\partial r} \right ) + \dfrac {1}{\sin \theta } \dfrac {\partial}{\partial \theta } \left ( \sin \theta \dfrac {\partial}{\partial \theta} \right ) + \dfrac {1}{\sin ^2 \theta} \dfrac {\partial ^2}{\partial \varphi ^2} \right ] - \dfrac {e^2}{4 \pi \epsilon _0 r } \label {6.1.3}\]

The contributions from rotational and radial components of the motion become clearer if we write out the complete Schrödinger equation,

\[ \left \{ -\dfrac {\hbar ^2}{2 \mu r^2} \left [ \dfrac {\partial}{\partial r} \left (r^2 \dfrac {\partial}{\partial r} \right ) + \dfrac {1}{\sin \theta } \dfrac {\partial}{\partial \theta } \left ( \sin \theta \dfrac {\partial}{\partial \theta} \right ) + \dfrac {1}{\sin ^2 \theta} \dfrac {\partial ^2}{\partial \varphi ^2} \right ] - \dfrac {e^2}{4 \pi \epsilon _0 r } \right \} \psi (r , \theta , \varphi ) = E \psi (r , \theta , \varphi ) \label {6.1.4}\]

multiply both sides of Equation \(\ref{6.1.4}\) by \(2 \mu r^2 \), and rearrange to obtain

\[\hbar ^2 \dfrac {\partial}{\partial r} \left ( r^2 \dfrac {\partial}{\partial r} \psi (r , \theta , \varphi ) \right ) + 2 \mu r^2 \left [ E + \dfrac {e^2}{4 \pi \epsilon _0 r } \right ] \psi (r , \theta , \varphi ) = \]

\[ - \hbar^2 \left [ \dfrac {1}{\sin \theta } \dfrac {\partial}{\partial \theta } \left ( \sin \theta \dfrac {\partial}{\partial \theta} \right ) + \dfrac {1}{\sin ^2 \theta} \dfrac {\partial ^2}{\partial \varphi ^2} \right ] \psi (r , \theta , \varphi ) \label {6.1.5}\]

Manipulating the Schrödinger equation in this way helps us recognize the square of the angular momentum operator in Equation \(\ref{6.1.5}\). The square of the angular momentum operator in Equation \(\ref{6.1.6}\).

\[\hat {M} ^2 = -\hbar ^2 \left [\dfrac {1}{\sin \theta } \dfrac {\partial}{\partial \theta } \left ( \sin \theta \dfrac {\partial}{\partial \theta} \right ) + \dfrac {1}{\sin ^2 \theta} \dfrac {\partial ^2}{\partial \varphi ^2} \right ] \label {6.1.6}\]

Substituting Equation \(\ref{6.1.6}\) into Equation \(\ref{6.1.5}\) produces

\[\hbar ^2 \dfrac {\partial}{\partial r } \left ( r^2 \dfrac {\partial}{\partial r} \psi (r , \theta , \varphi ) \right ) + 2 \mu r^2 [ E - \hat {V} ] \psi (r , \theta , \varphi ) = \hat {M} ^2 \psi (r, \theta , \varphi ) \label {6.1.7}\]

Since the angular momentum operator does not involve the radial variable, \(r\), we can separate variables in Equation \(\ref{6.1.7}\) by using a product wavefunction, as we did previously for rigid rotors. We know that the eigenfunctions of the angular momentum operator are the Spherical Harmonic functions (Table M4), \(Y (\theta , \varphi )\), so a good choice for a product function is

\[ \psi (r , \theta , \varphi ) = R (r) Y (\theta , \varphi ) \label {6.1.8}\]

The Spherical Harmonic functions provide information about where the electron is around the proton, and the radial function \(R(r)\) describes how far the electron is away from the proton.

To separate variables, substitute the product function, Equation \(\ref{6.1.8}\) into Equation \(\ref{6.1.7}\), evaluate partial derivatives, divide each side by R(r) \(Y (\theta, \varphi ) \), and set each side of that resulting equation equal to a constant \(\lambda\).

\[ \begin{align} \dfrac {\hbar ^2}{R (r)} \dfrac {\partial}{\partial r} r^2 \dfrac {\partial}{\partial r} R(r) + \dfrac {2 \mu r^2}{R (r)} [ E - V ] R (r) &= \lambda \label {6.1.9} \\[4pt] \dfrac {1}{Y (\theta , \varphi )} \hat {M} ^2 Y (\theta , \varphi ) &= \lambda \label {6.1.10} \end{align}\]

Equations \(\ref{6.1.9}\) and \(\ref{6.1.10}\) represent the radial differential equation and the angular differential equation, respectively. As we describe below, they are solved separately to give the \(Y (\theta , \varphi )\) angular functions and the R(r) radial functions.

Rearranging Equation \(\ref{6.1.10}\) yields

\[\hat {M} ^2 Y^{m_l}_l (\theta , \varphi ) = \lambda Y^{m_l}_l (\theta , \varphi ) \label {6.1.11}\]

where we have added the indices \(l\) and \(m_l\) to identify a particular spherical harmonic function. Note that the notation has changed from that used with the Rigid Rotor; it is customary to use \(J\) and \(m_J\) to represent the angular momentum quantum numbers for rotational states, but for electronic states, it is customary to use \(l\) and \(m_l\) to represent the same thing. Further, the electronic angular momentum is designated by \(L\) and the corresponding operator is called \(\hat {L}\). In complete electronic notation, Equation \(\ref{6.1.11}\) is

\[ \hat {L} ^2 Y^{m_l}_l (\theta , \varphi ) = \lambda Y^{m_l}_l (\theta , \varphi ) \label {6.1.12}\]

Equation \(\ref{6.1.12}\) says that \(Y^{m_l}_l (\theta , \varphi )\) must be an eigenfunction of the angular momentum operator \(\hat {L} ^2\) with eigenvalue λ. We know from the discussion of the Rigid Rotor that the eigenvalue λ is \(J(J+1)ħ^2\), or in electronic notation, \(l (l + 1) \hbar ^2\). Consequently, Equation \(\ref{6.1.12}\) becomes

\[ \hat {L} ^2 Y^{m_l}_l (\theta , \varphi ) = l (l + 1) \hbar ^2 Y^{m_l}_l (\theta , \varphi ) \label{6.1.13}\]

Using this value for λ and rearranging Equation \(\ref{6.1.9}\), we obtain

\[ - \dfrac {\hbar ^2}{2 \mu r^2} \dfrac {\partial}{\partial r} r^2 \dfrac {\partial}{\partial r} R(r) + \left [ \dfrac {l(l +1) \hbar ^2}{2 \mu r^2} + V (r) - E \right ] R (r) = 0 \label {6.1.14}\]

The details for solving Equation \(\ref{6.1.14}\) are provided elsewhere, but the procedure and consequences are similar to previously examined cases. As for the harmonic oscillator, an asymptotic solution (valid at large \(r\)) is found, and then the complete solutions are written as products of the asymptotic solution and polynomials arising from sequential truncations of a power series expansion.

## Contributors

David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski ("Quantum States of Atoms and Molecules")