Degenerate energy levels

In quantum mechanics, an energy level is degenerate if it corresponds to two or more different measurable states of a quantum system. Conversely, two or more different states of a quantum mechanical system are said to be degenerate if they give the same value of energy upon measurement. The number of different states corresponding to a particular energy level is known as the degree of degeneracy (or simply the degeneracy) of the level. It is represented mathematically by the Hamiltonian for the system having more than one linearly independent eigenstate with the same energy eigenvalue. When this is the case, energy alone is not enough to characterize what state the system is in, and other quantum numbers are needed to characterize the exact state when distinction is desired. In classical mechanics, this can be understood in terms of different possible trajectories corresponding to the same energy.

Degeneracy plays a fundamental role in quantum statistical mechanics. For an $N$-particle system in three dimensions, a single energy level may correspond to several different wave functions or energy states. These degenerate states at the same level all have an equal probability of being filled. The number of such states gives the degeneracy of a particular energy level.



Mathematics
The possible states of a quantum mechanical system may be treated mathematically as abstract vectors in a separable, complex Hilbert space, while the observables may be represented by linear Hermitian operators acting upon them. By selecting a suitable basis, the components of these vectors and the matrix elements of the operators in that basis may be determined. If $A$ is a $N × N$ matrix, $X$ a non-zero vector, and $λ$ is a scalar, such that $$AX = \lambda X$$, then the scalar $λ$ is said to be an eigenvalue of $A$ and the vector $X$ is said to be the eigenvector corresponding to $λ$. Together with the zero vector, the set of all eigenvectors corresponding to a given eigenvalue $λ$ form a subspace of $C^{n}$, which is called the eigenspace of $λ$. An eigenvalue $λ$ which corresponds to two or more different linearly independent eigenvectors is said to be degenerate, i.e., $$A X_1 = \lambda X_1$$ and $$ A X_2 = \lambda X_2$$, where $$ X_1 $$ and $$ X_2 $$ are linearly independent eigenvectors. The dimension of the eigenspace corresponding to that eigenvalue is known as its degree of degeneracy, which can be finite or infinite. An eigenvalue is said to be non-degenerate if its eigenspace is one-dimensional.

The eigenvalues of the matrices representing physical observables in quantum mechanics give the measurable values of these observables while the eigenstates corresponding to these eigenvalues give the possible states in which the system may be found, upon measurement. The measurable values of the energy of a quantum system are given by the eigenvalues of the Hamiltonian operator, while its eigenstates give the possible energy states of the system. A value of energy is said to be degenerate if there exist at least two linearly independent energy states associated with it. Moreover, any linear combination of two or more degenerate eigenstates is also an eigenstate of the Hamiltonian operator corresponding to the same energy eigenvalue. This clearly follows from the fact that the eigenspace of the energy value eigenvalue $λ$ is a subspace (being the kernel of the Hamiltonian minus $λ$ times the identity), hence is closed under linear combinations.

$$

Effect of degeneracy on the measurement of energy
In the absence of degeneracy, if a measured value of energy of a quantum system is determined, the corresponding state of the system is assumed to be known, since only one eigenstate corresponds to each energy eigenvalue. However, if the Hamiltonian $$\hat{H}$$ has a degenerate eigenvalue $$E_n$$ of degree gn, the eigenstates associated with it form a vector subspace of dimension gn. In such a case, several final states can be possibly associated with the same result $$E_n$$, all of which are linear combinations of the gn orthonormal eigenvectors $$|E_{n,i}\rangle$$.

In this case, the probability that the energy value measured for a system in the state $$|\psi\rangle$$ will yield the value $$E_n$$ is given by the sum of the probabilities of finding the system in each of the states in this basis, i.e., $$P(E_n)=\sum_{i=1}^{g_n}|\langle E_{n,i}|\psi\rangle|^2$$

Degeneracy in different dimensions
This section intends to illustrate the existence of degenerate energy levels in quantum systems studied in different dimensions. The study of one and two-dimensional systems aids the conceptual understanding of more complex systems.

Degeneracy in one dimension
In several cases, analytic results can be obtained more easily in the study of one-dimensional systems. For a quantum particle with a wave function $$|\psi\rangle$$ moving in a one-dimensional potential $$V(x)$$, the time-independent Schrödinger equation can be written as $$ -\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + V\psi = E\psi$$ Since this is an ordinary differential equation, there are two independent eigenfunctions for a given energy $$E$$ at most, so that the degree of degeneracy never exceeds two. It can be proven that in one dimension, there are no degenerate bound states for normalizable wave functions. A sufficient condition on a piecewise continuous potential $$V$$ and the energy $$E$$ is the existence of two real numbers $$M,x_0$$ with $$M \neq 0$$ such that $$\forall x > x_0$$ we have $$V(x) - E \geq M^2$$. In particular, $$V$$ is bounded below in this criterion.


 * {| class="toccolours collapsible collapsed" width="60%" style="text-align:left"

!Proof of the above theorem.

$$\begin{align} - \frac{\hbar^2}{2m}\frac{d^2\psi_1}{dx^2} + V\psi_1 &= E\psi_1 \\ - \frac{\hbar^2}{2m}\frac{d^2\psi_2}{dx^2} + V\psi_2 &= E\psi_2 \end{align}$$ Multiplying the first equation by $$ \psi_2 $$ and the second by $$\psi_1$$ and subtracting one from the other, we get: $$\psi_1\frac{d^2}{dx^2}\psi_2-\psi_2\frac{d^2}{d x^2}\psi_1 = 0$$ Integrating both sides $$\psi_1\frac{d\psi_2}{dx} - \psi_2\frac{d\psi_1}{dx} = \mbox{constant}$$ In case of well-defined and normalizable wave functions, the above constant vanishes, provided both the wave functions vanish at at least one point, and we find: $$\psi_1(x) = c\psi_2(x)$$ where $$c$$ is, in general, a complex constant. For bound state eigenfunctions (which tend to zero as $$x \to \infty$$), and assuming $$V$$ and $$E$$ satisfy the condition given above, it can be shown that also the first derivative of the wave function approaches zero in the limit $$x \to \infty$$, so that the above constant is zero and we have no degeneracy.
 * Considering a one-dimensional quantum system in a potential $$V(x)$$ with degenerate states $$|\psi_1\rangle$$ and $$|\psi_2\rangle$$ corresponding to the same energy eigenvalue $$E$$, writing the time-independent Schrödinger equation for the system:
 * Considering a one-dimensional quantum system in a potential $$V(x)$$ with degenerate states $$|\psi_1\rangle$$ and $$|\psi_2\rangle$$ corresponding to the same energy eigenvalue $$E$$, writing the time-independent Schrödinger equation for the system:
 * }

Degeneracy in two-dimensional quantum systems
Two-dimensional quantum systems exist in all three states of matter and much of the variety seen in three dimensional matter can be created in two dimensions. Real two-dimensional materials are made of monoatomic layers on the surface of solids. Some examples of two-dimensional electron systems achieved experimentally include MOSFET, two-dimensional superlattices of Helium, Neon, Argon, Xenon etc. and surface of liquid Helium. The presence of degenerate energy levels is studied in the cases of Particle in a box and two-dimensional harmonic oscillator, which act as useful mathematical models for several real world systems.

Particle in a rectangular plane
Consider a free particle in a plane of dimensions $$L_x$$ and $$L_y$$ in a plane of impenetrable walls. The time-independent Schrödinger equation for this system with wave function $$|\psi\rangle$$ can be written as $$ -\frac{\hbar^2}{2m}\left(\frac{\partial^2 \psi}{{\partial x}^2} + \frac{\partial^2 \psi}{{\partial y}^2}\right) = E\psi$$ The permitted energy values are $$E_{n_x,n_y} = \frac{\pi^2 \hbar^2}{2m}\left(\frac{n_x^2}{L_x^2}+\frac{n_y^2}{L_y^2}\right)$$ The normalized wave function is $$\psi_{n_x,n_y}(x,y) = \frac 2{\sqrt{L_xL_y}} \sin\left(\frac{n_x\pi x}{L_x}\right)\sin\left(\frac{n_y\pi y}{L_y}\right)$$ where $$n_x,n_y = 1,2,3,\dots$$

So, quantum numbers $$n_x$$ and $$n_y$$ are required to describe the energy eigenvalues and the lowest energy of the system is given by $$E_{1,1} = \pi^2\frac{\hbar^2}{2m}\left(\frac 1{L_x^2}+\frac 1{L_y^2}\right)$$

For some commensurate ratios of the two lengths $$L_x$$ and $$L_y$$, certain pairs of states are degenerate. If $$L_x/L_y=p/q$$, where p and q are integers, the states $$(n_x, n_y)$$ and $$(pn_y/q, qn_x/p)$$ have the same energy and so are degenerate to each other.

Particle in a square box
In this case, the dimensions of the box $$L_x = L_y = L$$ and the energy eigenvalues are given by $$E_{n_x,n_y}=\frac{\pi^2\hbar^2}{2mL^2}(n_x^2+n_y^2)$$

Since $$n_x$$ and $$n_y$$ can be interchanged without changing the energy, each energy level has a degeneracy of at least two when $$n_x$$ and $$n_y$$ are different. Degenerate states are also obtained when the sum of squares of quantum numbers corresponding to different energy levels are the same. For example, the three states (nx = 7, ny = 1), (nx = 1, ny = 7) and (nx = ny = 5) all have $$E=50 \frac{\pi^2\hbar^2}{2mL^2}$$ and constitute a degenerate set.

Degrees of degeneracy of different energy levels for a particle in a square box:

Particle in a cubic box
In this case, the dimensions of the box $$L_x = L_y =L_z= L$$ and the energy eigenvalues depend on three quantum numbers. $$E_{n_x,n_y,n_z} = \frac{\pi^2\hbar^2}{2mL^2}(n_x^2+n_y^2+n_z^2)$$

Since $$n_x$$, $$n_y$$ and $$n_z$$ can be interchanged without changing the energy, each energy level has a degeneracy of at least three when the three quantum numbers are not all equal.

Finding a unique eigenbasis in case of degeneracy
If two operators $$\hat{A}$$ and $$\hat{B}$$ commute, i.e., $$[\hat{A},\hat{B}] = 0$$, then for every eigenvector $$|\psi\rangle$$ of $$\hat{A}$$, $$\hat{B}|\psi\rang$$ is also an eigenvector of $$\hat{A}$$ with the same eigenvalue. However, if this eigenvalue, say $$\lambda$$, is degenerate, it can be said that $$\hat{B}|\psi\rangle$$ belongs to the eigenspace $$E_\lambda$$ of $$\hat{A}$$, which is said to be globally invariant under the action of $$\hat{B}$$.

For two commuting observables $A$ and $B$, one can construct an orthonormal basis of the state space with eigenvectors common to the two operators. However, $$\lambda$$ is a degenerate eigenvalue of $$\hat{A}$$, then it is an eigensubspace of $$\hat{A}$$ that is invariant under the action of $$\hat{B}$$, so the representation of $$\hat{B}$$ in the eigenbasis of $$\hat{A}$$ is not a diagonal but a block diagonal matrix, i.e. the degenerate eigenvectors of $$\hat{A}$$ are not, in general, eigenvectors of $$\hat{B}$$. However, it is always possible to choose, in every degenerate eigensubspace of $$\hat{A}$$, a basis of eigenvectors common to $$\hat{A}$$ and $$\hat{B}$$.

Choosing a complete set of commuting observables
If a given observable A is non-degenerate, there exists a unique basis formed by its eigenvectors. On the other hand, if one or several eigenvalues of $$\hat{A}$$ are degenerate, specifying an eigenvalue is not sufficient to characterize a basis vector. If, by choosing an observable $$\hat{B}$$, which commutes with $$\hat{A}$$, it is possible to construct an orthonormal basis of eigenvectors common to $$\hat{A}$$ and $$\hat{B}$$, which is unique, for each of the possible pairs of eigenvalues {a,b}, then $$\hat{A}$$ and $$\hat{B}$$ are said to form a complete set of commuting observables. However, if a unique set of eigenvectors can still not be specified, for at least one of the pairs of eigenvalues, a third observable $$\hat{C}$$, which commutes with both $$\hat{A}$$ and $$\hat{B}$$ can be found such that the three form a complete set of commuting observables.

It follows that the eigenfunctions of the Hamiltonian of a quantum system with a common energy value must be labelled by giving some additional information, which can be done by choosing an operator that commutes with the Hamiltonian. These additional labels required naming of a unique energy eigenfunction and are usually related to the constants of motion of the system.

Degenerate energy eigenstates and the parity operator
The parity operator is defined by its action in the $$|r\rangle$$ representation of changing r to −r, i.e. $$\langle r|P|\psi\rangle=\psi(-r)$$ The eigenvalues of P can be shown to be limited to $$\pm1$$, which are both degenerate eigenvalues in an infinite-dimensional state space. An eigenvector of P with eigenvalue +1 is said to be even, while that with eigenvalue −1 is said to be odd.

Now, an even operator $$\hat{A}$$ is one that satisfies, $$\tilde{A} = P \hat{A} P$$ $$[P,\hat{A}] = 0$$ while an odd operator $$\hat{B}$$ is one that satisfies $$P \hat{B}+\hat{B} P=0$$ Since the square of the momentum operator $$\hat{p}^2$$ is even, if the potential V(r) is even, the Hamiltonian $$\hat{H}$$ is said to be an even operator. In that case, if each of its eigenvalues are non-degenerate, each eigenvector is necessarily an eigenstate of P, and therefore it is possible to look for the eigenstates of $$\hat{H}$$ among even and odd states. However, if one of the energy eigenstates has no definite parity, it can be asserted that the corresponding eigenvalue is degenerate, and $$P|\psi\rangle$$ is an eigenvector of $$\hat{H}$$ with the same eigenvalue as $$|\psi\rangle$$.

Degeneracy and symmetry
The physical origin of degeneracy in a quantum-mechanical system is often the presence of some symmetry in the system. Studying the symmetry of a quantum system can, in some cases, enable us to find the energy levels and degeneracies without solving the Schrödinger equation, hence reducing effort.

Mathematically, the relation of degeneracy with symmetry can be clarified as follows. Consider a symmetry operation associated with a unitary operator $S$. Under such an operation, the new Hamiltonian is related to the original Hamiltonian by a similarity transformation generated by the operator $S$, such that $$H'=SHS^{-1}=SHS^\dagger$$, since $S$ is unitary. If the Hamiltonian remains unchanged under the transformation operation $S$, we have $$\begin{align} SHS^\dagger &= H \\[1ex] SHS^{-1} &= H \\[1ex] SH &= HS \\[1ex] [S,H] &= 0 \end{align}$$ Now, if $$|\alpha\rangle $$ is an energy eigenstate, $$H|\alpha\rangle = E|\alpha\rangle$$ where $E$ is the corresponding energy eigenvalue. $$HS|\alpha\rangle = SH|\alpha\rangle = SE|\alpha\rangle = ES|\alpha\rangle$$ which means that $$S|\alpha\rangle$$ is also an energy eigenstate with the same eigenvalue $E$. If the two states $$|\alpha\rangle$$ and $$S|\alpha\rangle$$ are linearly independent (i.e. physically distinct), they are therefore degenerate.

In cases where $S$ is characterized by a continuous parameter $$\epsilon$$, all states of the form $$S(\epsilon)|\alpha\rangle$$ have the same energy eigenvalue.

Symmetry group of the Hamiltonian
The set of all operators which commute with the Hamiltonian of a quantum system are said to form the symmetry group of the Hamiltonian. The commutators of the generators of this group determine the algebra of the group. An n-dimensional representation of the Symmetry group preserves the multiplication table of the symmetry operators. The possible degeneracies of the Hamiltonian with a particular symmetry group are given by the dimensionalities of the irreducible representations of the group. The eigenfunctions corresponding to a n-fold degenerate eigenvalue form a basis for a n-dimensional irreducible representation of the Symmetry group of the Hamiltonian.

Types of degeneracy
Degeneracies in a quantum system can be systematic or accidental in nature.

Systematic or essential degeneracy
This is also called a geometrical or normal degeneracy and arises due to the presence of some kind of symmetry in the system under consideration, i.e. the invariance of the Hamiltonian under a certain operation, as described above. The representation obtained from a normal degeneracy is irreducible and the corresponding eigenfunctions form a basis for this representation.

Accidental degeneracy
It is a type of degeneracy resulting from some special features of the system or the functional form of the potential under consideration, and is related possibly to a hidden dynamical symmetry in the system. It also results in conserved quantities, which are often not easy to identify. Accidental symmetries lead to these additional degeneracies in the discrete energy spectrum. An accidental degeneracy can be due to the fact that the group of the Hamiltonian is not complete. These degeneracies are connected to the existence of bound orbits in classical Physics.

Examples: Coulomb and Harmonic Oscillator potentials
For a particle in a central $1/r$ potential, the Laplace–Runge–Lenz vector is a conserved quantity resulting from an accidental degeneracy, in addition to the conservation of angular momentum due to rotational invariance.

For a particle moving on a cone under the influence of $1/r$ and $r^{2}$ potentials, centred at the tip of the cone, the conserved quantities corresponding to accidental symmetry will be two components of an equivalent of the Runge-Lenz vector, in addition to one component of the angular momentum vector. These quantities generate SU(2) symmetry for both potentials.

Example: Particle in a constant magnetic field
A particle moving under the influence of a constant magnetic field, undergoing cyclotron motion on a circular orbit is another important example of an accidental symmetry. The symmetry multiplets in this case are the Landau levels which are infinitely degenerate.

The hydrogen atom
In atomic physics, the bound states of an electron in a hydrogen atom show us useful examples of degeneracy. In this case, the Hamiltonian commutes with the total orbital angular momentum $$\hat{L}^2$$, its component along the z-direction, $$\hat{L}_z$$, total spin angular momentum $$\hat{S}^2$$ and its z-component $$\hat{S}_z$$. The quantum numbers corresponding to these operators are $$\ell$$, $$m_\ell$$, $$s$$ (always 1/2 for an electron) and $$m_s$$ respectively.

The energy levels in the hydrogen atom depend only on the principal quantum number $n$. For a given $n$, all the states corresponding to $$\ell = 0, \ldots, n-1$$ have the same energy and are degenerate. Similarly for given values of $n$ and $\ell$, the $$(2\ell+1)$$, states with $$m_\ell = -\ell, \ldots, \ell$$ are degenerate. The degree of degeneracy of the energy level En is therefore $$\sum_{\ell \mathop = 0}^{n-1} (2\ell+1) = n^2,$$ which is doubled if the spin degeneracy is included.

The degeneracy with respect to $$m_\ell$$ is an essential degeneracy which is present for any central potential, and arises from the absence of a preferred spatial direction. The degeneracy with respect to $$\ell$$ is often described as an accidental degeneracy, but it can be explained in terms of special symmetries of the Schrödinger equation which are only valid for the hydrogen atom in which the potential energy is given by Coulomb's law.

Isotropic three-dimensional harmonic oscillator
It is a spinless particle of mass m moving in three-dimensional space, subject to a central force whose absolute value is proportional to the distance of the particle from the centre of force. $$F = -kr$$ It is said to be isotropic since the potential $$V(r)$$ acting on it is rotationally invariant, i.e., $$V(r) = \tfrac{1}{2} m \omega^2 r^2$$ where $$\omega$$ is the angular frequency given by $\sqrt{k/m}$.

Since the state space of such a particle is the tensor product of the state spaces associated with the individual one-dimensional wave functions, the time-independent Schrödinger equation for such a system is given by- $$-\frac{\hbar^2}{2m} \left(\frac{\partial^2 \psi}{\partial x^2}+ \frac{\partial^2 \psi}{\partial y^2}+ \frac{\partial^2 \psi}{\partial z^2}\right) +\frac{1}{2}{m \omega^2 \left(x^2+y^2+z^2\right)\psi}=E\psi$$

So, the energy eigenvalues are $$E_{n_x,n_y,n_z} = \left(n_x + n_y + n_z + \tfrac{3}{2}\right) \hbar\omega$$ or, $$E_n = \left(n + \tfrac{3}{2}\right) \hbar\omega$$ where n is a non-negative integer. So, the energy levels are degenerate and the degree of degeneracy is equal to the number of different sets $$\{n_x, n_y, n_z\}$$ satisfying $$n_x + n_y + n_z = n$$ The degeneracy of the $$n$$-th state can be found by considering the distribution of $$n$$ quanta across $$n_x$$, $$n_y$$ and $$n_z$$. Having 0 in $$n_x$$ gives $$n + 1$$ possibilities for distribution across $$n_y$$ and $$n_z$$. Having 1 quanta in $$n_x$$ gives $$n$$ possibilities across $$n_y$$ and $$n_z$$ and so on. This leads to the general result of $$n - n_x + 1$$ and summing over all $$n$$ leads to the degeneracy of the $$n$$-th state, $$\sum_{n_x=0}^n (n-n_x+1)=\frac{(n+1)(n+2)}{2}$$ For the ground state $$n = 0$$, the degeneracy is $$1$$ so the state is non-degenerate. For all higher states, the degeneracy is greater than 1 so the state is degenerate.

Removing degeneracy
The degeneracy in a quantum mechanical system may be removed if the underlying symmetry is broken by an external perturbation. This causes splitting in the degenerate energy levels. This is essentially a splitting of the original irreducible representations into lower-dimensional such representations of the perturbed system.

Mathematically, the splitting due to the application of a small perturbation potential can be calculated using time-independent degenerate perturbation theory. This is an approximation scheme that can be applied to find the solution to the eigenvalue equation for the Hamiltonian H of a quantum system with an applied perturbation, given the solution for the Hamiltonian H0 for the unperturbed system. It involves expanding the eigenvalues and eigenkets of the Hamiltonian H in a perturbation series. The degenerate eigenstates with a given energy eigenvalue form a vector subspace, but not every basis of eigenstates of this space is a good starting point for perturbation theory, because typically there would not be any eigenstates of the perturbed system near them. The correct basis to choose is one that diagonalizes the perturbation Hamiltonian within the degenerate subspace.

Physical examples of removal of degeneracy by a perturbation
Some important examples of physical situations where degenerate energy levels of a quantum system are split by the application of an external perturbation are given below.

Symmetry breaking in two-level systems
A two-level system essentially refers to a physical system having two states whose energies are close together and very different from those of the other states of the system. All calculations for such a system are performed on a two-dimensional subspace of the state space.

If the ground state of a physical system is two-fold degenerate, any coupling between the two corresponding states lowers the energy of the ground state of the system, and makes it more stable.

If $$E_1$$ and $$E_2$$ are the energy levels of the system, such that $$E_1=E_2=E$$, and the perturbation $$W$$ is represented in the two-dimensional subspace as the following 2×2 matrix $$\mathbf{W} = \begin{bmatrix} 0 & W_{12} \\[1ex] W_{12}^* & 0 \end{bmatrix}. $$ then the perturbed energies are $$\begin{align} E_{+} &= E+|W_{12}| \\ E_{-} &= E-|W_{12}| \end{align}$$

Examples of two-state systems in which the degeneracy in energy states is broken by the presence of off-diagonal terms in the Hamiltonian resulting from an internal interaction due to an inherent property of the system include:
 * Benzene, with two possible dispositions of the three double bonds between neighbouring Carbon atoms.
 * Ammonia molecule, where the Nitrogen atom can be either above or below the plane defined by the three Hydrogen atoms.
 * molecule, in which the electron may be localized around either of the two nuclei.

Fine-structure splitting
The corrections to the Coulomb interaction between the electron and the proton in a Hydrogen atom due to relativistic motion and spin–orbit coupling result in breaking the degeneracy in energy levels for different values of l corresponding to a single principal quantum number n.

The perturbation Hamiltonian due to relativistic correction is given by $$H_r = -p^4 / 8m^3c^2$$ where $$p$$ is the momentum operator and $$m$$ is the mass of the electron. The first-order relativistic energy correction in the $$|nlm\rangle$$ basis is given by $$E_r = \left(-1/8m^3c^2\right) \left\langle n\ell m\right| p^4 \left|n\ell m \right\rangle$$

Now $$p^4 = 4m^2(H^0+e^2/r)^2$$ $$\begin{aligned} E_r &= -\frac{1}{2mc^2} \left[E_n^2 + 2 E_n e^2 \left\langle \frac{1}{r} \right\rangle + e^4 \left\langle \frac{1}{r^2} \right\rangle \right] \\ &= -\frac{1}{2} mc^2 \alpha^4 \left[-3/(4n^4)+1/{n^3(\ell+1/2)}\right] \end{aligned}$$ where $$\alpha$$ is the fine structure constant.

The spin–orbit interaction refers to the interaction between the intrinsic magnetic moment of the electron with the magnetic field experienced by it due to the relative motion with the proton. The interaction Hamiltonian is $$H_{so} = - \frac{e}{mc}\frac{\mathbf{m}\cdot\mathbf{L}}{r^3} = \frac{e^2}{m^2 c^2 r^3} \mathbf{S}\cdot\mathbf{L} $$ which may be written as $$H_{so} = \frac{e^2}{4m^2c^2r^3} \left[J^2 - L^2 - S^2\right]$$

The first order energy correction in the $$|j,m,\ell,1/2\rangle$$ basis where the perturbation Hamiltonian is diagonal, is given by $$E_{so} = \frac{\hbar^2 e^2}{4 m^2 c^2} \frac{j(j+1)-\ell(\ell+1)-\frac{3}{4}}{a_0^3 n^3 \ell(\ell+\frac{1}{2})(\ell+1)}$$ where $$a_0$$ is the Bohr radius. The total fine-structure energy shift is given by $$E_{fs} = - \frac{mc^2\alpha^4}{2n^3} \left[1/(j+1/2) - 3/4n\right]$$ for $j = \ell \pm \tfrac{1}{2}$.

Zeeman effect
The splitting of the energy levels of an atom when placed in an external magnetic field because of the interaction of the magnetic moment $$\vec{m}$$ of the atom with the applied field is known as the Zeeman effect.

Taking into consideration the orbital and spin angular momenta, $$\mathbf{L}$$ and $$\mathbf{S}$$, respectively, of a single electron in the Hydrogen atom, the perturbation Hamiltonian is given by $$\hat{V} = -(\mathbf{m}_\ell + \mathbf{m}_s)\cdot \mathbf{B}$$ where $$\mathbf{m}_\ell = -e \mathbf{L}/2m$$ and $$\mathbf{m}_s = -e \mathbf{S}/m$$. Thus, $$\hat{V} = \frac{e}{2m} (\mathbf{L}+2\mathbf{S})\cdot\mathbf{B}$$ Now, in case of the weak-field Zeeman effect, when the applied field is weak compared to the internal field, the spin–orbit coupling dominates and $\mathbf{L}$ and $\mathbf{S}$  are not separately conserved. The good quantum numbers are n, $\ell$, j and mj, and in this basis, the first order energy correction can be shown to be given by $$E_z = -\mu_B g_j B m_j,$$ where $$\mu_B = {e\hbar}/2m$$ is called the Bohr Magneton. Thus, depending on the value of $$m_j$$, each degenerate energy level splits into several levels.



In case of the strong-field Zeeman effect, when the applied field is strong enough, so that the orbital and spin angular momenta decouple, the good quantum numbers are now n, l, ml, and ms. Here, Lz and Sz are conserved, so the perturbation Hamiltonian is given by- $$\hat{V} = eB (L_z + 2S_z)/2m$$ assuming the magnetic field to be along the z-direction. So, $$\hat{V} = eB(m_\ell+2m_s)/2m$$ For each value of mℓ, there are two possible values of ms, $$\pm 1/2$$.

Stark effect
The splitting of the energy levels of an atom or molecule when subjected to an external electric field is known as the Stark effect.

For the hydrogen atom, the perturbation Hamiltonian is $$\hat{H}_{s}=-|e|Ez$$ if the electric field is chosen along the z-direction.

The energy corrections due to the applied field are given by the expectation value of $$\hat{H}_{s}$$ in the $$|n \ell m\rangle$$ basis. It can be shown by the selection rules that $$\langle n \ell m_\ell| z |n_1 \ell_1 m_{\ell1}\rangle \ne 0$$ when $$\ell = \ell_1 \pm 1$$ and $$m_\ell = m_{\ell1}$$.

The degeneracy is lifted only for certain states obeying the selection rules, in the first order. The first-order splitting in the energy levels for the degenerate states $$|2,0,0\rangle$$ and $$|2,1,0\rangle$$, both corresponding to n = 2, is given by $$\Delta E_{2,1,m_\ell} = \pm|e| \hbar^2 /(m_e e^2)E$$.