Now we would like to discuss how stationary states of a matter wave look like in a potential that exhibits sperical symmetry (\(E_{\mathrm{pot}} = E_{\mathrm{pot}} \left( r \right) = V \left( r \right)\)). In order to do so, we state the three-dimensional stationary Schrödinger equation
\[\begin{split}\begin{eqnarray}
H \psi \left( x,y,z \right) & = & E \psi \left( x,y,z \right)\\
- \frac{\hbar^2}{2m} \Delta \psi \left( x,y,z \right) + V \left( r \right) \psi \left( x,y,z \right) & = & E \psi \left( x,y,z \right)
\end{eqnarray}\end{split}\]
The solution of the Schrödinger equation in the case of spherically symmetric potentials is easier, if we express this equation in sperical coordinates,
\[\begin{split}\begin{eqnarray}
x & = & r \sin \left( \vartheta \right) \cos \left( \varphi \right)\\
y & = & r \sin \left( \vartheta \right) \sin \left( \varphi \right)\\
z & = & r \cos \left( \vartheta \right) \mathrm{.}
\end{eqnarray}\end{split}\]
The Laplace operator in sperical coordinates then reads as
In order to solve the time-indepedent Schrödinger equation we make use of separation of variables as we did before and use the ansatz
\[\psi \left( r,\vartheta,\varphi \right) = R \left( r \right) \cdot \Theta \left(\vartheta\right) \cdot \Phi \left(\varphi \right) \mathrm{.}\]
Using this ansatz in the last form of the Schrödinger equation results in
\[\begin{eqnarray}
0 & = &
\frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial R \left( r \right)}{\partial r} \right) \Theta \left(\vartheta\right) \Phi \left(\varphi \right)
+ \frac{1}{r^2 \sin \left( \vartheta \right)} \frac{\partial}{\partial \vartheta} \left( \sin \left( \vartheta \right) \frac{\partial \Theta \left(\vartheta\right)}{\partial \vartheta} \right) R \left( r \right) \Phi \left(\varphi \right)
+ \frac{1}{r^2 \sin^2 \left( \vartheta \right)} \frac{\partial^2 \Phi \left(\varphi \right)}{\partial \varphi^2} R \left( r \right) \Theta \left(\vartheta\right)
+ \frac{2m}{\hbar^2} \left( E - V\left( r \right) \right) R \left( r \right) \Theta \left(\vartheta\right) \Phi \left(\varphi \right) \mathrm{,}
\end{eqnarray}\]
which we multiply by \(r^2 \sin^2 \left( \vartheta \right) / \left[ R \left( r \right) \Theta \left(\vartheta\right) \Phi \left(\varphi \right)\right]\) in order to obtain
which allows us to compared both sides. The left hand side of the last equation depends solely on \(\varphi\), whereas the right hand side of the very same equation depends solely on \(r\) and \(\vartheta\).
Solving the solution function \(\Phi \left( \varphi \right)\)¶
So far we have used the ansatz of separation of variables and reshaped the Schrödinger equation leading to the situation that both sides of the same equation depend on different variables. Because this equation is supposed to be valid for every values of \(r\), \(\vartheta\), and \(\varphi\), both sides are identical to a constant value\(C_1\). Thus, for the left hand side it immediately follows
and deivide it through \(\sin^2 \left( \vartheta \right)\) before we sort the terms with respect to the parameters \(r\) and \(\vartheta\) resulting in
Similar to the case of the function \(\Phi \left( \varphi \right)\), here the left hand side of the last equation does solely depend on the parameter \(\vartheta\), whereas the right hand side of the very same equation solely depends on \(r\). As a consequence both sides of the same equation have to be equal to a constant \(C_2\). Concerning the funcion \(\Theta \left( \vartheta \right)\) we get
The solution of this differential equation we set in form of a power series
\[\theta = a_0 + a_1 \xi + a_2 \xi^2 + \ldots\]
Moreover, in order to obtain finite values for the series especially in the case of \(\xi = 1\) (\(\vartheta = 0^{\circ}\) or \(\vartheta = 180^{\circ}\)), it is mandatory that the series has only a finite number of elements. If we use this series in the Legendre differential equation and sort with respect to coefficients of the same power \(\xi^{k}\) one obtains the recursive formula
Since the series has to be finite, we set the \(l\)-th element which is \(a_l \cdot \xi^l\) as the last element with \(a_l \ne 0\) but \(a_{l+2}=0\). Thus, it direcly follows
The solutions of the Legendre differential equation are the Legendre polynomials (similar to the Hermite differential equation and Hermite polynomials),
Since the propability density \(\left| \psi \left( r,\vartheta,\varphi \right) \right|^2\) has to be unique at the same position \(\vartheta\) and \(\vartheta + \phi\) position, it follows \(\theta^2 \left( \vartheta \right) = \theta^2 \left( \vartheta + \pi\right)\) and \(\theta \left( \vartheta \right) = \pm \theta \left( \vartheta + \pi\right)\). As a consequence, every function represented through the power series of \(\xi^{k}\) contains only even or only odd
powers of \(\xi\).
In the case \(m\ne0\), we can solve the above differential equation (equal to \(-C_2\)) with the aid of the associated Legendre polynomials\(P_l^m \left( \cos \left( \vartheta \right) \right)\), that can be determined on the basis of the Legendre function\(\theta_l \left( \xi \right) = \cdot P_l \left( \cos \left( \vartheta \right) \right)\) and the conditional equation
Since \(P_l \left( \cos \left( \vartheta \right) \right) = P_l \left( \xi \right)\) is a power series of \(\xi\) with \(\xi^l\) as the highest power, the condition
\[\left| m \right| \le l\]
follows immediately. The constant prefactor in \(P_l^m \left( \cos \left( \vartheta \right) \right)\) is determined by means of the normalization condition
The square of the product function \(\left| Y_l^m \left( \vartheta, \varphi \right) \right|^2\) represents the \(\vartheta\) and \(\varphi\) dependence of the propability denisty of a quantum particle in a sperically symmetric potential. It is evident that for a given energy \(E\) and quantum number \(l\) there are \(2l+1\) distinct spherial hamonics \(Y_l^m\) as a consequnce of $ m \in `:nbsphinx-math:left[ -l ; +l right]`$.
The normalized spherical harmonics up to \(l=2\) are listed in the following table.
Solving the solution function \(R \left( r \right)\)¶
Previously we have derived the sperical harmonics and their squares as angular dependence of the probability density of a perticle within a spericaly symmetric potential. Now, we would like to solve the last dependence, namely the solution for the radial function \(R \left( r \right)\). Above we have derived the constant \(C_2\) with
\[\begin{split}\begin{eqnarray}
C_2 & = & l \left( l+1 \right) \\
{} & = & \frac{1}{R \left( r \right)} \frac{\mathrm{d}}{\mathrm{d} r} \left( r^2 \frac{\mathrm{d} R \left( r \right)}{\mathrm{d} r} \right)
+ \frac{2m}{\hbar^2} \left( E - V\left( r \right) \right) r^2
\end{eqnarray} \mathrm{.}\end{split}\]
If we multiply this equation with \(R \left( r \right)\) and divide through \(r^2\), we obtain
\[\begin{eqnarray}
\frac{l \left( l+1 \right)}{r^2} R \left( r \right) & = &
\frac{1}{r^2} \frac{\mathrm{d}}{\mathrm{d} r} \left( r^2 \frac{\mathrm{d} R \left( r \right)}{\mathrm{d} r} \right)
+ \frac{2m}{\hbar^2} \left( E - V\left( r \right) \right) R \left( r \right)
\end{eqnarray} \mathrm{,}\]
where the qunatum number \(l\) denotes the angular momentum \(\left| \vec{L} \right| = \sqrt{l \left( l+1 \right)} \cdot \hbar\) of our quantum particle. If we now make use of the Coulomb potential
as an example for a sperically symmetric potential and replace the mass \(m\) by the reduced mass \(\mu\) of an electron in the Coulomb potential of a nucleus with charge \(Z \cdot e\) and mass \(m_n \gg m_e\), as well as calculate the derivative of the first term, we obtain
\[\begin{eqnarray}
0 & = &
\frac{\mathrm{d}^2 R \left( r \right)}{\mathrm{d} r^2} + \frac{2}{r} \frac{\mathrm{d} R \left( r \right)}{\mathrm{d} r}
+ \left( \frac{2\mu}{\hbar^2} \left( E + \frac{Z e^2}{4 \pi \varepsilon_0 r} \right) - \frac{l \left( l+1 \right)}{r^2} \right) R \left( r \right)
\end{eqnarray} \mathrm{.}\]
In order to solve this equation for \(E< 0\), we use the ansatz
\[R \left( r \right) = u \left( r \right) \cdot \mathrm{e}^{- \kappa r} \mathrm{,}\]
were we made use of the abbreviation
\[\kappa = \frac{\sqrt{-2E}}{\hbar} \mathrm{.}\]
With the substtution
\[a = \frac{\mu Z e^2}{4 \pi \varepsilon_0 \hbar^2}\]
we obtain
\[\begin{eqnarray}
0 & = &
\frac{\mathrm{d}^2 u \left( r \right)}{\mathrm{d} r^2}
+ 2 \left( \frac{1}{r} - \kappa\right) \frac{\mathrm{d} u \left( r \right)}{\mathrm{d} r}
+ \left(2 \frac{a-\kappa}{r} - \frac{l \left( l+1 \right)}{r^2} \right) u \left( r \right)
\end{eqnarray} \mathrm{.}\]
In order to solve this problem we use the ansatz of a power series
\[u \left(r\right) = \sum_j b_j r^j \mathrm{.}\]
If we use this series in the differential equation and compare the coefficients, we obtain the recursive formula
In order to normalize \(R\left( r\right)\), the series has to have a finite number of elements. Furthermore, from the denomonator we get the condition \(j < l\). Let us assume the element \(j=n-1\) is the last non-zero element (\(b_{n-1} \ne 0\)) of this series, but the element \(j = n\) ist the first element that is zero (\(b_{n} = 0\)). Thus, we introduced the condition
\[j < n \mathrm{.}\]
Furthermore, on the basis of the recursive formula we obtain
The qunatization of the energy eigenstates arise from the condition \(\psi \left( r\longrightarrow\infty \right) \longrightarrow 0\). Furthermore, from the denominator of the recursive formula we derive
\[l \le j \le n-1\]
which leads to the condition for the quantum number of the angular momentum
\[l \le n-1 \mathrm{.}\]
On teh basis of teh recursicve formula w can calculate the function \(u \left( r \right)\) and thus \(R \left( r \right)\). This function \(R \left( r \right) = R_{n,l} \left( r \right)\) depends on the principal quantum number\(n\) (due to \(j < n\) in the sum formula) and on the quantum number of the angular momentum\(l\) (due to the recursive formula).
The first three normalized radial functions \(R_{n,l} \left( r \right)\) are listed in the following table (\(N = \left(Z/\left( a_0 n\right)\right)^{3/2}\), \(x = Z r / \left(a_0 n\right)\), \(a_0 = 4 \pi \varepsilon_0 \hbar^2 / \left(\mu e^2 \right)\)).
We can see that the energy of a quantum state depends only on the principal quantum number \(n\), but not on the other quantum numbers \(l\) and \(m\). Because \(-l \le m \le +l\), there exists \(2l + 1\) energeticall equal (degenerated) states for every particular value of \(l\). Furthermore, for every particular values of \(n\), there exists
\[\sum_{l=0}^{n-1} \left(2l+1\right) = n^2\]
distinct states \(\left( n,l,m \right)\) with \(n^2\) distinct wave functions