6.2: The linear Stark Effect

As a more involved example of the level degeneracy lifting by a perturbation, let us discuss the Stark effect \(^{6}\) - the atomic level splitting by an external electric field. Let us study this effect, in the linear approximation, for a hydrogen-like atom/ion. Taking the direction of the external electric field \(\mathscr{E}\) (which is practically always uniform on the atomic scale) for the \(z\)-axis, the perturbation may be represented by the following Hamiltonian: \[\hat{H}^{(1)}=-F \hat{z}=-q \mathscr{E} \hat{z}=-q \mathscr{E} r \cos \theta .\] (In the last form, the operator sign is dropped, because we will work in the coordinate representation.)

As you (should :-) remember, energy levels of a hydrogen-like atom/ion depend only on the principal quantum number \(n\) - see Eq. (3.201); hence all the states, besides the ground \(1 s\) state in which \(n=1\) and \(l=m=0\), have some orbital degeneracy, which grows rapidly with \(n\). Let us consider the lowest degenerate level with \(n=2\). Since, according to Eq. (3.203), \(0 \leq l \leq n-1\), at this level the orbital quantum number \(l\) may equal either 0 (one \(2 s\) state, with \(m=0\) ) or 1 (three \(2 p\) states, with \(m=0, \pm 1\) ). Due to this 4-fold degeneracy, \(\mathrm{H}^{(1)}\) is a \(4 \times 4\) matrix with 16 elements: \[\begin{array}{rl} & \overbrace{m=0}^{l=0} \overbrace{m=0 m=+1 m=-1}^{l=1} \\ \mathrm{H}^{(1)}= & \left(\begin{array}{llll} H_{11} & H_{12} & H_{13} & H_{14} \\ H_{21} & H_{22} & H_{23} & H_{24} \\ H_{31} & H_{32} & H_{33} & H_{34} \\ H_{41} & H_{42} & H_{43} & H_{44} \end{array}\right) m=0, \quad m=0, \quad m=0, \\ m & m=1 . \end{array}\] However, there is no need to be scared. First, due to the Hermitian nature of the operator, only ten of these matrix elements (four diagonal and six off-diagonal ones) may be substantially different from each other. Moreover, due to a high symmetry of the problem, there are a lot of zeros even among these elements. Indeed, let us have a look at the angular components \(Y_{l}^{m}\) of the corresponding wavefunctions, with \(l=0\) and \(l=1\), described by Eqs. (3.174)-(3.175). For the states with \(m=\pm 1\), the azimuthal parts of wavefunctions are proportional to \(\exp \{\pm i \varphi\}\); hence the off-diagonal elements \(H_{34}\) and \(H_{43}\) of the matrix (30), relating these functions, are proportional to \[\oint d \Omega Y_{1}^{\pm}{ }^{*} \hat{H}^{(1)} Y_{1}^{\mp} \propto \int_{0}^{2 \pi} d \varphi\left(e^{\pm i \varphi}\right)^{*}\left(e^{\mp i \varphi}\right)=0 .\] The azimuthal-angle symmetry also kills the off-diagonal elements \(H_{13}, H_{14}, H_{23}, H_{24}\) (and hence their complex conjugates \(\mathrm{H}_{31}, \mathrm{H}_{41}, \mathrm{H}_{32}\), and \(\mathrm{H}_{42}\) ), because they relate states with \(m=0\) and \(m=\pm 1\), and hence are proportional to \[\oint d \Omega Y_{1}^{0^{*}} \hat{H}^{(1)} Y_{1}^{\pm 1} \propto \int_{0}^{2 \pi} d \varphi e^{\pm i \varphi}=0\] For the diagonal matrix elements \(H_{33}\) and \(H_{44}\), corresponding to \(l=1\) and \(m=\pm 1\), the azimuthal-angle integrals do not vanish, but since the corresponding spherical functions depend on the polar angle as \(\sin \theta\), these elements are proportional to \[\oint d \Omega Y_{1}^{\pm 1^{*}} \hat{H}^{(1)} Y_{1}^{\pm 1} \propto \int_{0}^{\pi} \sin \theta d \theta \sin \theta \cos \theta \sin \theta=\int_{-1}^{+1} \cos \theta\left(1-\cos ^{2} \theta\right) d(\cos \theta),\] and hence are equal to zero - as any limit-symmetric integral of an odd function. Finally, for the states \(2 s\) and \(2 p\) with \(m=0\), the diagonal elements \(H_{11}\) and \(H_{22}\) are also killed by the polar-angle integration: \[\begin{aligned} &\oint d \Omega Y_{0}^{0^{*}} \hat{H}^{(1)} Y_{0}^{0} \propto \int_{0}^{\pi} \sin \theta d \theta \cos \theta=\int_{-1}^{1} \cos \theta d(\cos \theta)=0, \\ &\oint d \Omega Y_{0}^{1^{*}} \hat{H}^{(1)} Y_{0}^{1} \propto \int_{0}^{\pi} \sin \theta d \theta \cos ^{3} \theta=\int_{-1}^{+1} \cos ^{3} \theta d(\cos \theta)=0 . \end{aligned}\] Hence, the only non-zero elements of the matrix (30) are two off-diagonal elements \(H_{12}\) and \(H_{21}\), which relate two states with the same \(m=0\), but different \(l=\{0,1\}\), because they are proportional to \[\oint d \Omega Y_{0}^{0^{*}} \cos \theta Y_{1}^{0}=\frac{\sqrt{3}}{4 \pi} \int_{0}^{2 \pi} d \varphi \int_{0}^{\pi} \sin \theta d \theta \cos ^{2} \theta=\frac{1}{\sqrt{3}} \neq 0 .\] What remains is to use Eqs. (3.209) for the radial parts of these functions to complete the calculation of those two matrix elements: \[H_{12}=H_{21}=-\frac{q \mathscr{E}}{\sqrt{3}} \int_{0}^{\infty} r^{2} d r R_{2,0}(r) r R_{2,1}(r) .\] Due to the additive structure of the function \(R_{2,0}(r)\), the integral falls into a sum of two table integrals, both of the type MA Eq. (6.7d), finally giving \[H_{12}=H_{21}=3 q \delta r_{0},\] where \(r_{0}\) is the spatial scale (3.192); for the hydrogen atom, it is just the Bohr radius \(r_{\mathrm{B}}\) - see Eq. (1.10).

Thus, the perturbation matrix (30) is reduced to \[\mathrm{H}^{(1)}=\left(\begin{array}{cccc} 0 & 3 q \mathscr{E} r_{0} & 0 & 0 \\ 3 q \mathscr{E} r_{0} & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{array}\right),\] so that the condition (26) of self-consistency of the system (25), \[\left|\begin{array}{cccc} -E_{2}^{(1)} & 3 q \delta r_{0} & 0 & 0 \\ 3 q \delta r_{0} & -E_{2}^{(1)} & 0 & 0 \\ 0 & 0 & -E_{2}^{(1)} & 0 \\ 0 & 0 & 0 & -E_{2}^{(1)} \end{array}\right|=0\] gives a very simple characteristic equation \[\left(E_{2}^{(1)}\right)^{2}\left[\left(E_{2}^{(1)}\right)^{2}-\left(3 q \mathscr{E} r_{0}\right)^{2}\right]=0 .\] with four roots: \[\left(E_{2}^{(1)}\right)_{1,2}=0, \quad\left(E_{2}^{(1)}\right)_{3,4}=\pm 3 q \mathscr{E} r_{0} .\] so that the degeneracy is only partly lifted - see the levels in Fig. 3. 7

Fig. 6.3. The linear Stark effect for the level \(n=2\) of a hydrogen-like atom.

Generally, in order to understand the nature of states corresponding to these levels, we should return to Eq. (25) with each calculated value of \(E_{2}{ }^{(1)}\), and find the corresponding expansion coefficients \(\left\langle n^{\prime \prime(0)} \mid n^{\prime}\right\rangle\) that describe the perturbed states. However, in our simple case, the outcome of this procedure is clear in advance. Indeed, since the states with \(\{l=1, m=\pm 1\}\) are not affected by the perturbation at all (in the linear approximation in the electric field), their degeneracy is not lifted, and energy is not affected - see the middle line in Fig. 3. On the other hand, the partial perturbation matrix connecting the states \(2 s\) and \(2 p\), i.e. the top left \(2 \times 2\) part of the full matrix (39), is proportional to the Pauli matrix \(\sigma_{x}\), and we already know the result of its diagonalization - see Eqs. (4.113)-(4.114). This means that the upper and lower split levels correspond to very simple linear combinations of the previously degenerate states with \(m=0\), \[|\pm\rangle=\frac{1}{\sqrt{2}}(|2 s\rangle \pm|2 p\rangle) .\] Finally, let us estimate the magnitude of the linear Stark effect for a hydrogen atom. For a very high electric field of \(\mathscr{E}=3 \times 10^{6} \mathrm{~V} / \mathrm{m},{ }^{8}|q|=e \approx 1.6 \times 10^{-19} \mathrm{C}\), and \(r_{0}=r_{\mathrm{B}} \approx 0.5 \times 10^{-10} \mathrm{~m}\), we get a level splitting of \(3 q \delta r_{0} \approx 0.8 \times 10^{-22} \mathrm{~J} \approx 0.5 \mathrm{meV}\). This number is much lower than the unperturbed energy of the level, \(E_{2}=-E_{\mathrm{H}} /\left(2 \times 2^{2}\right) \approx-3.4 \mathrm{eV}\), so that the perturbative result is quite applicable. On the other hand, the calculated splitting is much larger than the resolution limit imposed by the line’s natural width \(\left(\sim 10^{-7} E_{2}\right.\), see Chapter 9), so that the effect is quite observable even in substantially lower electric fields. Note, however, that our simple results are quantitatively correct only when the Stark splitting (42) is much larger that the fine-structure splitting of the same level in the absence of the field- see the next section.

---

\({ }^{6}\) This effect was discovered experimentally in 1913 by Johannes Stark and independently by Antonio Lo Surdo, so it is sometimes (and more fairly) called the "Stark - Lo Surdo effect". Sometimes this name is used with the qualifier "dc" to distinguish it from the \(a c\) Stark effect - the energy level shift under the effect of an ac field - see Sec. 5 below.

\({ }^{7}\) The proportionality of this splitting to the small field is responsible for the qualifier "linear" in the name of this effect. If observable effects grow only as \(\varepsilon^{2}\) (see, e.g., Problem 9), the term quadratic Stark effect is used instead.

\({ }^{8}\) This value approximately corresponds to the threshold of electric breakdown in air at ambient conditions, due to the impact ionization. As a result, experiments with higher dc fields are rather difficult.