next up previous
Next: Linear Stark Effect Up: Time-Independent Perturbation Theory Previous: Quadratic Stark Effect


Degenerate Perturbation Theory

Let us now consider systems in which the eigenstates of the unperturbed Hamiltonian, $ H_0$ , possess degenerate energy levels. It is always possible to represent degenerate energy eigenstates as the simultaneous eigenstates of the Hamiltonian and some other Hermitian operator (or group of operators). Let us denote this operator (or group of operators) $ L$ . We can write

$\displaystyle H_0\, \vert n, l\rangle = E_n\, \vert n, l\rangle,$ (7.66)

and

$\displaystyle L\,\vert n,l\rangle = L_{n\,l}\, \vert n, l\rangle,$ (7.67)

where $ [H_0, L] = 0$ . Here, the $ E_n$ and the $ L_{n\,l}$ are real numbers that depend on the quantum numbers $ n$ , and $ n$ and $ l$ , respectively. It is always possible to find a sufficient number of operators that commute with the Hamiltonian in order to ensure that the $ L_{n\,l}$ are all different. In other words, we can choose $ L$ such that the quantum numbers $ n$ and $ l$ uniquely specify each eigenstate. Suppose that for each value of $ n$ there are $ N_n$ different values of $ l$ . In other words, the $ n$ th energy eigenstate is $ N_n$ -fold degenerate.

In general, $ L$ does not commute with the perturbing Hamiltonian, $ H_1$ . In this situation, we expect the perturbation to couple degenerate states with the same value of $ n$ , but different values of $ l$ . Let us naively attempt to use the standard perturbation theory of Section 7.3 to evaluate the modified energy eigenstates and energy levels. A direct generalization of Equations (7.32) and (7.33) yields

$\displaystyle E_{nl}' = E_n + e_{nlnl} + \sum_{n', l' \neq n,l} \frac{\vert e_{n'l'nl}\vert^{\,2}}{E_n - E_{n'}} + {\cal O}(\epsilon^{\,3}),$ (7.68)

and

$\displaystyle \vert n, l\rangle' = \vert n,l\rangle + \sum_{n', l'\neq n, l} \frac{e_{n'l'nl}}{E_n-E_{n'}}\,\vert n',l'\rangle + {\cal O}(\epsilon^{\,2}),$ (7.69)

where

$\displaystyle e_{n'l'nl} = \langle n',l'\vert\,H_1\,\vert n,l\rangle.$ (7.70)

It is fairly obvious that the summations in Equations (7.69) and (7.70) are not well behaved if the $ n$ th energy level is degenerate. The problem terms are those that involve coupling to unperturbed eigenstates labeled by the same value of $ n$ , but different values of $ l$ : that is, those states whose unperturbed energies are $ E_n$ . These terms give rise to singular factors $ 1/(E_n - E_n)$ in the summations. Note, however, that this problem would not exist if the matrix elements, $ e_{nl'nl}$ , of the perturbing Hamiltonian between distinct, degenerate, unperturbed energy eigenstates corresponding to the eigenvalue $ E_n$ were zero. In other words, if

$\displaystyle \langle n, l' \vert\,H_1\,\vert n, l\rangle = \lambda_{n\,l}\, \delta_{l\,l'}$ (7.71)

then all of the singular terms in Equations (7.69) and (7.70) would vanish.

In general, Equation (7.72) is not satisfied. Fortunately, we can always redefine the unperturbed energy eigenstates belonging to the eigenvalue $ E_n$ in such a manner that Equation (7.72) is satisfied. Let us define $ N_n$ new states that are linear combinations of the $ N_n$ original degenerate eigenstates corresponding to the eigenvalue $ E_n$ :

$\displaystyle \vert n,l^{\,(1)}\rangle = \sum_{k=1,N_n} \langle n,k\vert n,l^{\,(1)}\rangle \vert n,k\rangle.$ (7.72)

Note that these new states are also degenerate energy eigenstates of the unperturbed Hamiltonian corresponding to the eigenvalue $ E_n$ . The $ \vert n,l^{\,(1)}\rangle$ are chosen in such a manner that they are eigenstates of the perturbing Hamiltonian, $ H_1$ . Thus,

$\displaystyle H_1\, \vert n, l^{\,(1)}\rangle = \lambda_{n\,l} \,\vert n, l^{\,(1)}\rangle.$ (7.73)

The $ \vert n,l^{\,(1)}\rangle$ are also chosen so that they are orthonormal. It follows that

$\displaystyle \langle n, l'^{\,(1)} \vert \,H_1\,\vert n, l^{\,(1)}\rangle = \lambda_{n\,l}\,\delta_{l\,l'}.$ (7.74)

Thus, if we use the new eigenstates, instead of the old ones, then we can employ Equations (7.69) and (7.70) directly, because all of the singular terms vanish. The only remaining difficulty is to determine the new eigenstates in terms of the original ones.

Now,

$\displaystyle \sum_{l=1,N_n} \vert n,l\rangle \langle n,l\vert = 1,$ (7.75)

where 1 denotes the identity operator in the sub-space of all unperturbed energy eigenkets corresponding to the eigenvalue $ E_n$ . Using this completeness relation, the operator eigenvalue equation (7.74) can be transformed into a straightforward matrix eigenvalue equation:

$\displaystyle \sum_{l''=1,N_n}\langle n, l'\vert H_1\vert n, l''\rangle \langle...
... n, l^{\,(1)}\rangle = \lambda_{n\,l}\, \langle n, l'\vert n, l^{\,(1)}\rangle.$ (7.76)

This can be written more transparently as

$\displaystyle {\bf U} \,{\bf x} = \lambda \,{\bf x},$ (7.77)

where the elements of the $ N_n\times N_n$ Hermitian matrix $ {\bf U}$ are

$\displaystyle U_{j\,k} = \langle n, j\vert\, H_1\,\vert n, k\rangle.$ (7.78)

Provided that the determinant of $ {\bf U}$ is non-zero, Equation (7.78) can always be solved to give $ N_n$ eigenvalues $ \lambda_{n\,l}$ (for $ l=1$ to $ N_n$ ), with $ N_n$ corresponding eigenvectors $ {\bf x}_{n\,l}$ [92]. The eigenvectors specify the weights of the new eigenstates in terms of the original eigenstates: that is,

$\displaystyle ({\bf x}_{n\,l})_k = \langle n, k\vert n, l^{\,(1)}\rangle,$ (7.79)

for $ k=1$ to $ N_n$ . In our new scheme, Equations (7.69) and (7.70) yield

$\displaystyle E_{nl}' = E_n + \lambda_{n\,l} + \sum_{n'\neq n, l'} \frac{\vert e_{n'l'nl}\vert^{\,2}}{E_n - E_{n'}} + {\cal O}(\epsilon^{\,3}),$ (7.80)

and

$\displaystyle \vert n, l^{\,(1)}\rangle' = \vert n,l^{\,(1)}\rangle + \sum_{n'\...
...e_{n'l'nl}}{E_n-E_{n'}}\,\vert n',l'^{\,(1)}\rangle + {\cal O}(\epsilon^{\,2}).$ (7.81)

There are no singular terms in these expressions, because the summations are over $ n'\neq n$ . In other words, they specifically exclude the problematic, degenerate, unperturbed energy eigenstates corresponding to the eigenvalue $ E_n$ . Note that the first-order energy-shifts are equivalent to the eigenvalues of the matrix equation (7.78).

Incidentally, it is clear, from the previous analysis, that if the perturbing Hamiltonian, $ H_1$ , commutes with $ H_0$ and the operator (or group of operators) $ L$ then there is no danger of singular terms appearing in the perturbation expansion to second order, because $ e_{n'l'nl}=e_{n'l'nl}\,\delta_{ll'}$ . Another way of saying this is that there are no singular terms if the simultaneous eigenstates of $ H_0$ and $ L$ are also eigenstates of $ H_1$ .


next up previous
Next: Linear Stark Effect Up: Time-Independent Perturbation Theory Previous: Quadratic Stark Effect
Richard Fitzpatrick 2016-01-22