This page was generated from `/home/lectures/exp3/source/notebooks/L28_AMA/L28_spherical_potential.ipynb`_.
A particle in a spherically symmetric potential¶
The sperically symmetric Schrödinger eqution¶
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
with the Laplace operator
The solution of the Schrödinger equation in the case of spherically symmetric potentials is easier, if we express this equation in sperical coordinates,
The Laplace operator in sperical coordinates then reads as
and we can reform the Schrödinger equation resulting in
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
Using this ansatz in the last form of the Schrödinger equation results in
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
The very last equation can be re-written in the form
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
which has the solution
Since \(\Phi\) has to be unique for the whole space, it follows
Thus, \(\sqrt{C_1}\) has to be an integer and we can refine the solution
Moreover, we can normalize the solution in accord to
with the consequence for the normalization constant
Now we can state the normalized function
In addition, these functions are normalized, namely
Solving the solution function \(\Theta \left( \vartheta \right)\)¶
In order to solve the solution for \(\Theta \left( \vartheta \right)\) we again use the reshaped Schrödinger equation,
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
In the case \(C_1 = 0\) (\(m = 0\)) and making use of \(\xi = \cos \left( \vartheta \right)\) we result in the Legendre differential equation
The solution of this differential equation we set in form of a power series
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
follows immediately. The constant prefactor in \(P_l^m \left( \cos \left( \vartheta \right) \right)\) is determined by means of the normalization condition
The product functions
are called sperical harmonics. Since both factors are allready normalized it immediately follows for the product
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.
\(l\) |
\(m\) |
\(Y_l^m \left( \vartheta,\varphi \right)\) |
---|---|---|
0 |
0 |
\(\frac{1}{2\sqrt{\pi}}\) |
1 |
\(\pm 1\) |
\(\mp \sqrt{\frac{3}{8\pi}} \sin \left( \vartheta\right) \mathrm{e}^{\pm i\varphi}\) |
1 |
0 |
\(\sqrt{\frac{3}{4\pi}} \cos \left( \vartheta\right)\) |
2 |
\(\pm 2\) |
\(\sqrt{\frac{15}{32\pi}} \sin^2 \left( \vartheta\right) \mathrm{e}^{\pm 2 i\varphi}\) |
2 |
\(\pm 1\) |
\(\mp \sqrt{\frac{15}{8\pi}} \cos \left( \vartheta\right) \sin \left( \vartheta\right) \mathrm{e}^{\pm i\varphi}\) |
2 |
0 |
\(\sqrt{\frac{5}{16\pi}} \left( 3 \cos^2 \left( \vartheta\right) -1 \right)\) |
Fig.: Spherical harmonics :math:`Y_l^m left(vartheta, varphiright)` from left to the right: :math:`Y_{0}^{0}`, :math:`Y_{1}^{-1}`, :math:`Y_{1}^{0}`, and :math:`Y_{1}^{1}`. Note the swaped colors of the dumbbells, if changing :math:`m = -1 rightarrow +1`.
Fig.: Spherical harmonics :math:`Y_l^m left(vartheta, varphiright)` from left to the right: :math:`Y_{2}^{-2}`, :math:`Y_{2}^{-1}`, :math:`Y_{2}^{0}`, :math:`Y_{2}^{1}`, and :math:`Y_{2}^{2}`. Note the swaped colors of the dumbbells, if changing :math:`m = -1 rightarrow +1`. Changing :math:`m = -2 rightarrow +2` only affects the complex exponential function and is not visualized.
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
If we multiply this equation with \(R \left( r \right)\) and divide through \(r^2\), we obtain
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
In order to solve this equation for \(E< 0\), we use the ansatz
were we made use of the abbreviation
With the substtution
we obtain
In order to solve this problem we use the ansatz of a power series
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
Furthermore, on the basis of the recursive formula we obtain
which leads to
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
which leads to the condition for the quantum number of the angular momentum
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)\)).
\(n\) |
\(l\) |
\(R_{n,l} \left( r \right)\) |
---|---|---|
1 |
0 |
\(2N \mathrm{e}^{-x}\) |
2 |
0 |
\(2N \mathrm{e}^{-x} \left(1-x\right)\) |
2 |
1 |
\(\frac{1}{\sqrt{3}}2N \mathrm{e}^{-x} x\) |
3 |
0 |
\(2N \mathrm{e}^{-x} \left(1-2x+\frac{2}{3}x^2\right)\) |
3 |
1 |
\(\frac{\sqrt{2}}{3}2N \mathrm{e}^{-x} x \left(2-x\right)\) |
3 |
2 |
\(\frac{\sqrt{2}}{3\sqrt{10}}2N \mathrm{e}^{-x} x^2\) |
Fig.: (left) The radial function :math:`R_{n,l} left(rright)` and (right) the corresponding probability density :math:`left| R_{n,l} left(rright) right|^2` for :math:`n=1`.
Fig.: (left) The radial function :math:`R_{n,l} left(rright)` and (right) the corresponding probability density :math:`left| R_{n,l} left(rright) right|^2` for :math:`n=2`.
Fig.: (left) The radial function :math:`R_{n,l} left(rright)` and (right) the corresponding probability density :math:`left| R_{n,l} left(rright) right|^2` for :math:`n=3`.
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
distinct states \(\left( n,l,m \right)\) with \(n^2\) distinct wave functions
.