Classical angular momentum is given by

$$\mathbf{L}=\mathbf{r}\times \mathbf{p}=(y{p}_{z}-z{p}_{y})\mathbf{i}+(z{p}_{x}-x{p}_{z})\mathbf{j}+(x{p}_{y}-y{p}_{x})\mathbf{k}.$$ | (259) |

The angular momentum operator $\widehat{\mathbf{L}}=\widehat{\mathbf{r}}\times \widehat{\mathbf{p}}$ has components

${\widehat{L}}_{x}=\widehat{y}{\widehat{p}}_{z}-\widehat{z}{\widehat{p}}_{y}=-i\mathrm{\hslash}(y{\displaystyle \frac{\partial}{\partial z}}-z{\displaystyle \frac{\partial}{\partial y}}),$ | (260) | ||

${\widehat{L}}_{y}=\widehat{z}{\widehat{p}}_{x}-\widehat{x}{\widehat{p}}_{z}=-i\mathrm{\hslash}(z{\displaystyle \frac{\partial}{\partial x}}-x{\displaystyle \frac{\partial}{\partial z}}),$ | (261) | ||

${\widehat{L}}_{z}=\widehat{x}{\widehat{p}}_{y}-\widehat{y}{\widehat{p}}_{x}=-i\mathrm{\hslash}(x{\displaystyle \frac{\partial}{\partial y}}-y{\displaystyle \frac{\partial}{\partial z}}).$ | (262) |

Total angular momentum is given by

$${\widehat{\mathbf{L}}}^{2}={\widehat{L}}_{x}^{2}+{\widehat{L}}_{y}^{2}+{\widehat{L}}_{z}^{2}=-{\mathrm{\hslash}}^{2}{\mathrm{\Delta}}_{\theta ,\varphi}.$$ | (263) |

The commutation relations are

$[{\widehat{L}}_{x},{\widehat{L}}_{y}]=i\mathrm{\hslash}{\widehat{L}}_{z},[{\widehat{L}}_{y},{\widehat{L}}_{z}]=i\mathrm{\hslash}{\widehat{L}}_{x},[{\widehat{L}}_{z},{\widehat{L}}_{x}]=i\mathrm{\hslash}{\widehat{L}}_{y},$ | |||

(264) | |||

$[{\widehat{L}}_{x},{\widehat{\mathbf{L}}}^{2}]=0,[{\widehat{L}}_{y},{\widehat{\mathbf{L}}}^{2}]=0,[{\widehat{L}}_{z},{\widehat{\mathbf{L}}}^{2}]=0.$ | (265) |

Heisenberg’s uncertainty principle dictates that it is not possible to know two components of angular momentum at the same time (since their commutator does not vanish). However, it is possible to determine one of the components (say ${L}_{z}$) simultaneously with the total angular momentum, since $[{\widehat{L}}_{z},{\widehat{\mathbf{L}}}^{2}]=0$.

We already determined that the spherical harmonics ${Y}_{lm}$ are joint eigenfunctions of ${\widehat{\mathbf{L}}}^{2}$ and ${\widehat{L}}_{z}$,

$${\widehat{\mathbf{L}}}^{2}{Y}_{lm}(\theta ,\varphi )={\mathrm{\hslash}}^{2}l(l+1){Y}_{lm},$$ | (266) |

$${\widehat{L}}_{z}{Y}_{lm}(\theta ,\varphi )=\mathrm{\hslash}m{Y}_{lm}.$$ | (267) |

The eigenvalues are ${\mathrm{\hslash}}^{2}l(l+1)$ and $\mathrm{\hslash}m$, respectively.
E.g., for $l=2$ the length of the angular momentum vector is $|\mathbf{L}|=\sqrt{6}\mathrm{\hslash}$ and there are 5 possible values values of ${L}_{z}$, as is illustrated in the following figure.

Instead of interpreting the angular Schrödinger equation as a differential equation, this equation can be solved very efficiently by employing the perspective of operators and vectors.

In Dirac notation, the eigenstates of the ${\widehat{L}}_{z}$ and ${\widehat{\mathbf{L}}}^{2}$ are denoted as $|lm\u27e9$, such that

$${\widehat{L}}_{z}|lm\u27e9=\mathrm{\hslash}m|lm\u27e9,{\widehat{\mathbf{L}}}^{2}|lm\u27e9={\mathrm{\hslash}}^{2}l(l+1)|lm\u27e9.$$ | (268) |

We now introduce the *ladder operators* ${\widehat{L}}_{\pm}={\widehat{L}}_{x}\pm i{\widehat{L}}_{y}$, which are *non-hermitian* and related by ${\widehat{L}}_{+}^{\u2020}={\widehat{L}}_{-}$. These operators fulfill the commutation relation

$$[{\widehat{L}}_{z},{\widehat{L}}_{\pm}]=[{\widehat{L}}_{z},{\widehat{L}}_{x}]\pm i[{\widehat{L}}_{z},{\widehat{L}}_{y}]=i\mathrm{\hslash}{\widehat{L}}_{y}\pm \mathrm{\hslash}{L}_{x}=\pm \mathrm{\hslash}{L}_{\pm}.$$ | (269) |

The utility of these operators arises from the fact that they relate the eigenstates $|lm\u27e9$ of fixed $l$ but different $m$ in a systematic way. Indeed, the following calculation shows that ${\widehat{L}}_{\pm}|lm\u27e9$ is still an eigenstate of ${L}_{z}$,

$${\widehat{L}}_{z}({\widehat{L}}_{\pm}|lm\u27e9)={\widehat{L}}_{\pm}{\widehat{L}}_{z}|lm\u27e9+[{\widehat{L}}_{z},{\widehat{L}}_{\pm}]|lm\u27e9=\mathrm{\hslash}(m\pm 1){\widehat{L}}_{\pm}|lm\u27e9,$$ | (270) |

with eigenvalue $\mathrm{\hslash}(m\pm 1)$. Furthermore, because $[{\widehat{\mathbf{L}}}^{2},{\widehat{L}}_{\pm}]=0$ this state is also an eigenstate of ${\widehat{\mathbf{L}}}^{2}$, with eigenvalue ${\mathrm{\hslash}}^{2}l(l+1)$. Thus,

$${\widehat{L}}_{\pm}|l,m\u27e9={c}_{lm}^{\pm}|l,m\pm 1\u27e9,$$ | (271) |

where the normalisation constants work out as

$${c}_{lm}^{+}=\mathrm{\hslash}\sqrt{(l-m)(l+m+1)},{c}_{lm}^{-}={c}_{l,-m}^{+}.$$ | (272) |

Since ${c}_{ll}^{+}=0$, the state $|ll\u27e9$ is determined by the condition

$${\widehat{L}}_{+}|ll\u27e9=|\mathbb{\U0001d7d8}\u27e9,$$ | (273) |

while all other states follow by repeated application of ${\widehat{L}}_{-}$.

This construction of the angular momentum eigenstates and determination of their eigenvalues is purely algebraic, and sidesteps any explicit reference to spherical harmonics. Their explicit form can be recovered by writing $\u27e8\theta ,\varphi |lm\u27e9={Y}_{lm}(\theta ,\varphi )$, where $|\theta ,\varphi \u27e9$ denotes position basis states on the unit sphere, parameterised in spherical polar coordinates. ${Y}_{ll}$ is then obtained from Eq. (273), where ${\widehat{L}}_{+}$ is expressed as the associated differential operator, and the other spherical harmonics follow by successive application of the differential operator associated with ${\widehat{L}}_{-}$.

The orbital motion of an electron with angular momentum $\mathbf{L}$ gives rise to a magnetic moment

$$\mathbf{m}=-\frac{e}{2{m}_{e}}\mathbf{L}.$$ | (274) |

When a magnetic field $\mathbf{B}$ is applied, this gives rise to an interaction energy

$${V}_{B}=-\mathbf{m}\cdot \mathbf{B},$$ | (275) |

which can be used to measure the angular momentum. We assume that the magnetic field is applied in $z$ direction, $\mathbf{B}={B}_{z}\mathbf{k}$, therefore

$${V}_{B}=\frac{e{B}_{z}}{2{m}_{e}}{L}_{z}.$$ | (276) |

Quantum mechanically the interaction energy is represented by the operator

$${\widehat{V}}_{B}=\frac{e{B}_{z}}{2{m}_{e}}{\widehat{L}}_{z}.$$ | (277) |

In a central potential a magnetic field hence lifts the $2l+1$ fold degeneracy of the angular part of the wavefunction: For each value of m, the energies of the orbitals are shifted by an amount

$$\delta {E}_{m}={\mu}_{B}{B}_{z}m,$$ | (278) |

which depends on the magnetic quantum number $m$ (hence the name). Here ${\mu}_{B}=\frac{e\mathrm{\hslash}}{2{m}_{e}}$ is the so-called Bohr magneton.