Skip to main content
Physics LibreTexts

7.7: Fine Structure

Let us now consider the energy levels of hydrogen-like atoms (i.e., alkali metal atoms) in more detail. The outermost electron moves in a spherically symmetric potential $ V(r)$ due to the nuclear charge and the charges of the other electrons (which occupy spherically symmetric closed shells). The shielding effect of the inner electrons causes $ V(r)$ to depart from the pure Coulomb form. This splits the degeneracy of states characterized by the same value of $ n$ , but different values of $ l$ . In fact, higher $ l$ states have higher energies.

Let us examine a phenomenon known as fine structure, which is due to interaction between the spin and orbital angular momenta of the outermost electron. This electron experiences an electric field


$\displaystyle {\bf E} = \frac{\nabla V}{e}.$ (679)



However, a non-relativistic charge moving in an electric field also experiences an effective magnetic field


$\displaystyle {\bf B} = - \frac{{\bf v} \times {\bf E}}{c^2}.$ (680)



Now, an electron possesses a spin magnetic moment [see Equation (462)]


$\displaystyle {\mbox{\boldmath$\mu$}} = - \frac{e\, {\bf S}}{m_e}.$ (681)



We, therefore, expect a spin-orbit contribution to the Hamiltonian of the form


$\displaystyle H_{LS} = - {\mbox{\boldmath$\mu$}}\cdot {\bf B} = - \frac{e \,{\b...
...dr}\right) = \frac{1}{m_e^{\,2}\,c^2\,r} \frac{d V}{dr}\, {\bf L}\cdot {\bf S},$ (682)



where $ {\bf L} = m_e \,{\bf x}\times{\bf v}$ is the orbital angular momentum. Actually, when the above expression is compared to the observed spin-orbit interaction, it is found to be too large by a factor of two. There is a classical explanation for this, due to spin precession, which we need not go into. The correct quantum mechanical explanation requires a relativistically covariant treatment of electron dynamics (this is achieved using the so-called Dirac equation--see Chapter 11).

Let us now apply perturbation theory to a hydrogen-like atom, using $ H_{LS}$ as the perturbation (with $ H_{LS}$ taking one half of the value given above), and


$\displaystyle H_0 = \frac{{p}^2}{2\,m_e} + V(r)$ (683)



as the unperturbed Hamiltonian. We have two choices for the energy eigenstates of $ H_0$ . We can adopt the simultaneous eigenstates of $ H_0, L^2, S^2, L_z$ and $ S_z$ , or the simultaneous eigenstates of $ H_0, L^2, S^{\,2}, J^{\,2},$ and $ J_z$ , where $ {\bf J} = {\bf L} + {\bf S}$ is the total angular momentum. Although the departure of $ V(r)$ from a pure $ 1/r$ form splits the degeneracy of same $ n$ , different $ l$ , states, those states characterized by the same values of $ n$ and $ l$ , but different values of $ m_l$ , are still degenerate. (Here, $ m_l, m_s,$ and $ m_j$ are the quantum numbers corresponding to $ L_z, S_z,$ and $ J_z$ , respectively.) Moreover, with the addition of spin degrees of freedom, each state is doubly degenerate due to the two possible orientations of the electron spin (i.e., $ m_s = \pm 1/2$ ). Thus, we are still dealing with a highly degenerate system. We know, from Section 7.6, that the application of perturbation theory to a degenerate system is greatly simplified if the basis eigenstates of the unperturbed Hamiltonian are also eigenstates of the perturbing Hamiltonian. Now, the perturbing Hamiltonian, $ H_{LS}$ , is proportional to $ {\bf L}\cdot {\bf S}$ , where


$\displaystyle {\bf L} \cdot{\bf S} = \frac{J^{\,2} - L^2 - S^{\,2}}{2}.$ (684)



It is fairly obvious that the first group of operators ( $ H_0, L^2, S^2, L_z$ and $ S_z$ ) does not commute with $ H_{LS}$ , whereas the second group ( $ H_0, L^2, S^{\,2}, J^{\,2},$ and $ J_z$ ) does. In fact, $ {\bf L}\cdot {\bf S}$ is just a combination of operators appearing in the second group. Thus, it is advantageous to work in terms of the eigenstates of the second group of operators, rather than those of the first group.

We now need to find the simultaneous eigenstates of $ H_0, L^2, S^{\,2}, J^{\,2},$ and $ J_z$ . This is equivalent to finding the eigenstates of the total angular momentum resulting from the addition of two angular momenta: $ j_1=l$ , and $ j_2 = s = 1/2$ . According to Equation (568), the allowed values of the total angular momentum are $ j=l+1/2$ and $ j=l-1/2$ . We can write


$\displaystyle \vert l+1/2, m\rangle$ $\displaystyle = \cos\alpha\, \vert m-1/2, 1/2\rangle + \sin\alpha\, \vert m+1/2, -1/2\rangle,~~$ (685)
$\displaystyle \vert l-1/2, m\rangle$ $\displaystyle = -\sin\alpha\, \vert m-1/2, 1/2\rangle + \cos\alpha\, \vert m+1/2, -1/2\rangle.~~~~$ (686)



Here, the kets on the left-hand side are $ \vert j,m_j\rangle $ kets, whereas those on the right-hand side are $ \vert m_l, m_s\rangle$ kets (the $ j_1, j_2$ labels have been dropped, for the sake of clarity). We have made use of the fact that the Clebsch-Gordon coefficients are automatically zero unless $ m_j=m_l+m_s$ . We have also made use of the fact that both the $ \vert j,m_j\rangle $ and $ \vert m_l, m_s\rangle$ kets are orthonormal, and have unit lengths. We now need to determine


$\displaystyle \cos\alpha = \langle m-1/2,1/2\vert l+1/2, m\rangle,$ (687)



where the Clebsch-Gordon coefficient is written in $ \langle m_l, m_s\vert j, m_j\rangle$ form.

Let us now employ the recursion relation for Clebsch-Gordon coefficients, Equation (574), with $ j_1=l, j_2 = 1/2, j = l+1/2, m_1=m-1/2, m_2=1/2$ (lower sign). We obtain


  $\displaystyle [(l+1/2)\,(l+3/2)-m\,(m+1)]^{1/2} \,\langle m-1/2, 1/2\vert l+1/2, m\rangle$    
  $\displaystyle = [l\,(l+1)-(m-1/2)\,(m+1/2)]^{1/2}\, \langle m+1/2, 1/2\vert l+1/2, m+1\rangle,$ (688)



which reduces to


$\displaystyle \langle m-1/2, 1/2\vert l+1/2, m\rangle = \sqrt{\frac{l+m+1/2}{l+m+3/2}}\, \langle m+1/2, 1/2\vert l+1/2, m+1\rangle.$ (689)



We can use this formula to successively increase the value of $ m_l$ . For instance,


$\displaystyle \langle m-1/2, 1/2\vert l+1/2, m\rangle =\sqrt{\frac{l+m+1/2}{l+m+3/2}} \sqrt{\frac{l+m+3/2}{l+m+5/2}} \langle m+3/2, 1/2\vert l+1/2, m+2\rangle.$ (690)



This procedure can be continued until $ m_l$ attains its maximum possible value, $ l$ . Thus,


$\displaystyle \langle m-1/2, 1/2\vert l+1/2, m\rangle = \sqrt{\frac{l+m+1/2}{2\,l+1}}\, \langle l, 1/2\vert l+1/2, l+1/2\rangle.$ (691)



Consider the situation in which $ m_l$ and $ m$ both take their maximum values, $ l$ and $ 1/2$ , respectively. The corresponding value of $ m_j$ is $ l+1/2$ . This value is possible when $ j=l+1/2$ , but not when $ j=l-1/2$ . Thus, the $ \vert m_l, m_s\rangle$ ket $ \vert l,1/2\rangle$ must be equal to the $ \vert j,m_j\rangle $ ket $ \vert l+1/2, l+1/2\rangle$ , up to an arbitrary phase-factor. By convention, this factor is taken to be unity, giving


$\displaystyle \langle l, 1/2\vert l+1/2, l+1/2\rangle = 1.$ (692)



It follows from Equation (691) that


$\displaystyle \cos\alpha=\langle m-1/2, 1/2\vert l+1/2, m\rangle = \sqrt{\frac{l+m+1/2}{2\,l+1}}.$ (693)





$\displaystyle \sin^2\alpha = 1 - \frac{l+m+1/2}{2\,l+1} = \frac{l-m+1/2}{2\,l+1}.$ (694)



We now need to determine the sign of $ \sin\alpha$ . A careful examination of the recursion relation, Equation (574), shows that the plus sign is appropriate. Thus,


$\displaystyle \vert l+1/2, m\rangle$ $\displaystyle =\sqrt{\frac{l+m+1/2}{2\,l+1}}\,\vert m-1/2, 1/2\rangle+\sqrt{\frac{l-m+1/2}{2\,l+1}}\,\vert m+1/2, -1/2\rangle,$ (695)
$\displaystyle \vert l-1/2, m\rangle$ $\displaystyle = - \sqrt{\frac{l-m+1/2}{2\,l+1}} \,\vert m-1/2,1/2\rangle + \sqrt{\frac{l+m+1/2}{2\,l+1}} \,\vert m+1/2, -1/2\rangle.$ (696)



It is convenient to define so called spin-angular functions using the Pauli two-component formalism:


$\displaystyle {\cal Y}_{l\,m}^{j=l\pm 1/2}\equiv {\cal Y}_{l\,m}^{\pm}$ $\displaystyle = \pm \sqrt{\frac{l\pm m+1/2}{2\,l+1}}\, Y_{l\,\,m-1/2}(\theta, \...
...i_+ +\sqrt{\frac{l\mp m+1/2}{2\,l+1}} \,Y_{l\,\,m+1/2}(\theta,\varphi)\, \chi_-$    
  $\displaystyle = \frac{1}{\sqrt{2\,l+1}}\left( \begin{array}{c} \pm \sqrt{l\pm m...
....5ex] \sqrt{l\mp m+1/2}\,\,Y_{l\,\,m+1/2}(\theta, \varphi) \end{array} \right).$ (697)



These functions are eigenfunctions of the total angular momentum for spin one-half particles, just as the spherical harmonics are eigenfunctions of the orbital angular momentum. A general wavefunction for an energy eigenstate in a hydrogen-like atom is written


$\displaystyle \psi_{nlm\pm} = R_{n\,l}(r)\, {\cal Y}_{l\,m}^\pm.$ (698)



The radial part of the wavefunction, $ R_{n\,l}(r)$ , depends on the radial quantum number $ n$ and the angular quantum number $ l$ . The wavefunction is also labeled by $ m$ , which is the quantum number associated with $ J_z$ . For a given choice of $ l$ , the quantum number $ j$ (i.e., the quantum number associated with $ J^{\,2}$ ) can take the values $ l\pm 1/2$ .

The $ \vert l\pm 1/2, m\rangle$ kets are eigenstates of $ {\bf L}\cdot {\bf S}$ , according to Equation (684). Thus,


$\displaystyle {\bf L} \cdot{\bf S}\, \vert j=l\pm 1/2,m_j= m\rangle = \frac{\hbar^2}{2} \left[ j\,(j+1) - l\,(l+1) - 3/4\right]\vert j,m\rangle,$ (699)





$\displaystyle {\bf L} \cdot{\bf S}\, \vert l+ 1/2, m\rangle$ $\displaystyle = \frac{l \,\hbar^2}{2}\, \vert l+ 1/2, m\rangle,$ (700)
$\displaystyle {\bf L} \cdot{\bf S}\, \vert l- 1/2, m\rangle$ $\displaystyle = -\frac{(l+1)\, \hbar^2}{2}\, \vert l- 1/2, m\rangle.$ (701)



It follows that


$\displaystyle \oint d{\mit\Omega}\,({\cal Y}_{l\,m}^+)^\dagger \,{\bf L} \cdot{\bf S}\, {\cal Y}_{l \,m}^+$ $\displaystyle = \frac{l \,\hbar^2}{2},$ (702)
$\displaystyle \oint d{\mit\Omega}\,({\cal Y}_{l \,m}^-)^\dagger \,{\bf L} \cdot{\bf S}\, {\cal Y}_{l\,m}^-$ $\displaystyle = -\frac{(l+1) \,\hbar^2}{2},$ (703)



where the integrals are over all solid angle, $ d{\mit\Omega} = \sin\theta\,d\theta\,d\varphi$ .

Let us now apply degenerate perturbation theory to evaluate the shift in energy of a state whose wavefunction is $ \psi_{nlm\pm}$ due to the spin-orbit Hamiltonian, $ H_{LS}$ . To first order, the energy-shift is given by


$\displaystyle {\mit\Delta} E_{nlm\pm} = \int dV\,(\psi_{nlm\pm})^\dagger\, H_{LS}\,\psi_{nlm\pm},$ (704)



where the integral is over all space, $ dV = r^2\,d{\mit\Omega}$ . Equations (682) (remembering the factor of two), (698), and (702)-(703) yield


$\displaystyle {\mit\Delta} E_{nlm+}$ $\displaystyle = +\frac{1}{2\,m_e^{\,2}\,c^2} \left\langle \frac{1}{r}\frac{dV}{dr} \right\rangle \frac{l\,\hbar^2}{2},$ (705)
$\displaystyle {\mit\Delta} E_{nlm-}$ $\displaystyle =- \frac{1}{2\,m_e^{\,2}\,c^2} \left\langle \frac{1}{r}\frac{dV}{dr} \right\rangle \frac{(l+1)\,\hbar^2}{2},$ (706)





$\displaystyle \left\langle \frac{1}{r}\frac{dV}{dr} \right\rangle = \int_0^\infty dr\, r^2\,(R_{n\,l})^\ast \,\frac{1}{r}\frac{dV}{dr}\, R_{n\,l}.$ (707)



Let us now apply the above result to the case of a sodium atom. In chemist's notation, the ground state is written


$\displaystyle (1s)^2 (2s)^2(2p)^6(3s).$ (708)



The inner ten electrons effectively form a spherically symmetric electron cloud. We are interested in the excitation of the eleventh electron from $ 3s$ to some higher energy state. The closest (in energy) unoccupied state is $ 3p$ . This state has a higher energy than $ 3s$ due to the deviations of the potential from the pure Coulomb form. In the absence of spin-orbit interaction, there are six degenerate $ 3p$ states. The spin-orbit interaction breaks the degeneracy of these states. The modified states are labeled $ (3p)_{1/2}$ and $ (3p)_{3/2}$ , where the subscript refers to the value of $ j$ . The four $ (3p)_{3/2}$ states lie at a slightly higher energy level than the two $ (3p)_{1/2}$ states, because the radial integral (707) is positive. The splitting of the $ (3p)$ energy levels of the sodium atom can be observed using a spectroscope. The well-known sodium D line is associated with transitions between the $ 3p$ and $ 3s$ states. The fact that there are two slightly different $ 3p$ energy levels (note that spin-orbit coupling does not split the $ 3s$ energy levels) means that the sodium D line actually consists of two very closely spaced spectroscopic lines. It is easily demonstrated that the ratio of the typical spacing of Balmer lines to the splitting brought about by spin-orbit interaction is about $ 1 : \alpha^2$ , where


$\displaystyle \alpha = \frac{e^2}{2\,\epsilon_0\, h\, c} = \frac{1}{137}$ (709)


is the fine structure constant. Actually, Equations (705)-(706) are not entirely correct, because we have neglected an effect (namely, the relativistic mass correction of the electron) that is the same order of magnitude as spin-orbit coupling. (See Exercises 3 and 4.)