Fact-checked by Grok 2 weeks ago

Ladder operator

In quantum mechanics, ladder operators, also referred to as raising and lowering operators or creation and annihilation operators, are linear combinations of position and momentum operators that systematically increase or decrease the eigenvalues of a given observable by discrete quanta, facilitating the algebraic solution of eigenvalue problems without direct integration of differential equations.^[1] These operators satisfy specific commutation relations, such as [a, a^\dagger] = 1 for the harmonic oscillator case, where a is the lowering operator and a^\dagger is the raising operator, enabling the construction of complete sets of eigenstates from a ground or reference state.^[2] Introduced by Paul Dirac in the foundational operator formalism of quantum mechanics during the late 1920s, ladder operators represent a cornerstone of the algebraic approach to quantum theory, shifting emphasis from wavefunctions to abstract state vectors and operator algebras.^[2] The most prominent application of ladder operators arises in the quantum harmonic oscillator, where the lowering operator a_- and raising operator a_+ are defined as a_- = \frac{1}{\sqrt{2}}(x + i p) and a_+ = \frac{1}{\sqrt{2}}(x - i p) (in units where \hbar = m = \omega = 1), acting on energy eigenstates to generate an infinite ladder of discrete levels with energies E_n = (n + \frac{1}{2}), where n = 0, 1, 2, \dots.^[3] Applying a_+ to the ground state \psi_0, which satisfies a_- \psi_0 = 0 and has the Gaussian form \psi_0(x) = \pi^{-1/4} e^{-x^2/2}, yields higher excited states \psi_n = \frac{(a_+)^n}{\sqrt{n!}} \psi_0, confirming the equidistant spectrum and non-zero zero-point energy.^[3] This framework extends to multi-particle systems, modeling phonons in solids or photons in electromagnetic fields, and forms the basis for second quantization in quantum field theory.^[1] In systems with rotational symmetry, such as angular momentum, ladder operators J_+ and J_- adjust the magnetic quantum number m by \pm \hbar while preserving the total angular momentum quantum number j, satisfying commutation relations [J_z, J_\pm] = \pm \hbar J_\pm and [J_+, J_-] = 2 \hbar J_z.^[1] These operators generate the full multiplet of states from a highest-weight state where J_+ |j, m = j\rangle = 0, applicable to orbital angular momentum in hydrogen-like atoms and spin in multi-electron systems.^[1] Beyond these, ladder operators underpin the algebraic structure of Lie groups in quantum mechanics, influencing treatments of the hydrogen atom radial equation and more complex potentials through generalized raising and lowering mechanisms.^[2]

Fundamentals

Terminology

Ladder operators refer to a pair of mathematical operators in quantum mechanics that act on the eigenstates of a given observable, shifting its eigenvalues by fixed, discrete amounts. These operators facilitate transitions between energy levels or other quantized spectra, enabling the systematic exploration of the operator's eigenspectrum.^[4] The pair typically consists of a raising operator, denoted a^\dagger, which increases the eigenvalue, and a lowering operator, denoted a, which decreases it. In the context of Hilbert space, these operators are Hermitian adjoints of each other. This relationship ensures consistency with the inner product structure, though the norms of the resulting states are scaled by factors related to the eigenvalues. This adjoint relationship underpins their role in maintaining the unitarity of quantum evolution.^[2] Common synonyms for ladder operators include raising and lowering operators, step operators, and shift operators; in quantum field theory and many-body systems, they are often termed creation and annihilation operators due to their interpretation in terms of particle number. The term "ladder" originates from the ladder-like arrangement of the discrete eigenvalues, resembling rungs on a ladder that the operators ascend or descend. These operators play a key role in generating successive eigenstates from a fundamental one.

Mathematical Motivation

Ladder operators arise in the study of eigenspectra of linear operators on vector spaces, particularly where the spectrum consists of degenerate or equally spaced eigenvalues, allowing these operators to systematically connect eigenstates differing by a fixed spectral increment. In operator theory, such operators facilitate the exploration of the structure of eigenspaces by generating chains of states from a given eigenvector, revealing the dimensionality and organization of representations associated with the operator. This approach is especially valuable in infinite-dimensional Hilbert spaces, where traditional diagonalization may be infeasible, providing a pathway to classify and navigate the spectrum without explicit solution of the characteristic equation.^[5] A profound connection exists between ladder operators and the generators of Lie algebras, notably in semisimple cases like su(2) for compact representations and su(1,1) or its real form sl(2,R) for non-compact ones, where the ladder operators serve as root vectors that shift the weights (eigenvalues of the Cartan subalgebra) by root values. In representation theory, these algebras underpin the systematic construction of irreducible representations, with ladder operators enabling the decomposition of modules into weight spaces linked by algebraic action. For instance, in su(1,1), the non-compact nature leads to unitary representations on infinite-dimensional spaces, where ladder operators generate discrete series of states with eigenvalues increasing or decreasing without bound.^[6] An illustrative example occurs in infinite-dimensional representations of sl(2,R), where repeated application of a raising ladder operator to a lowest-weight vector produces an orthonormal basis spanning the representation space, with each step incrementing the weight by a fixed amount, thus forming a "ladder" of basis vectors that diagonalizes the Cartan generator. This method not only generates the basis but also highlights the module's structure, such as its unitarity and irreducibility, essential for applications in analysis and physics.^[7] Historical mathematical precursors to ladder operators can be traced to finite-dimensional matrix theory, particularly Jordan's canonical form introduced in the late 19th century, where Jordan blocks describe chains of generalized eigenvectors connected by powers of a nilpotent operator, analogous to finite ladders terminating at the ends of the chain. These blocks, part of the broader theory of linear transformations on finite-dimensional spaces, prefigure the infinite chains in Lie algebra representations by illustrating how operator action creates structured bases for non-diagonalizable cases.^[8]

General Formulation

Definition and Construction

In quantum mechanics, ladder operators are linear operators acting on the Hilbert space of a physical system that connect successive eigenstates of a self-adjoint observable, typically the Hamiltonian H with a discrete spectrum \{\lambda_n\}_{n=0}^\infty or finite range. For an eigenstate |n\rangle satisfying H |n\rangle = \lambda_n |n\rangle, the raising operator L^+ and lowering operator L^- satisfy L^+ |n\rangle \propto |n+1\rangle and L^- |n\rangle \propto |n-1\rangle (for n \geq 1), shifting the associated eigenvalue by a fixed increment \Delta\lambda = \lambda_{n+1} - \lambda_n or its negative. Ladder operators for an observable A (often the Hamiltonian H) are defined by the commutation relations [A, L^\pm] = \pm \Delta L^\pm, where \Delta is the fixed step size, ensuring L^\pm shifts eigenvalues of A by \pm \Delta. These operators, often denoted interchangeably as raising/lowering operators, facilitate the algebraic manipulation of energy levels and quantum states without solving the full Schrödinger equation. Ladder operators play a crucial role in completing the eigenspace basis by iteratively generating states from an initial seed. Starting from a lowest-weight state |0\rangle (where L^- |0\rangle = 0), the states are constructed as |n\rangle \propto (L^+)^n |0\rangle for n = 0, 1, 2, \dots, spanning the entire irreducible representation subspace associated with the observable. This recursive generation ensures orthogonality and completeness within the relevant sector, simplifying the diagonalization of H. In cases with a highest-weight state (where L^+ |m\rangle = 0), the process reverses using powers of L^-.^[9] The nature of ladder operators distinguishes bounded from unbounded chains, reflecting the underlying algebraic structure. Bounded ladders terminate at both ends, yielding finite-dimensional representations (e.g., a finite number of states between highest and lowest weights), as seen in compact symmetry groups where repeated application eventually yields zero. Unbounded ladders extend infinitely in at least one direction, producing infinite-dimensional Hilbert spaces, typical of non-compact symmetries where no such termination occurs. This dichotomy determines the dimensionality and completeness of the state space for the system.^[9]

Algebraic Properties and Commutation Relations

In general, an observable K (such as the Hamiltonian H or J_z) and its associated ladder operators L_- and L_+ satisfy [K, L_\pm] = \pm \Delta L_\pm, where \Delta is the eigenvalue step size (often set to 1 or \hbar in appropriate units). This ensures that applying L_+ to an eigenvector | \lambda \rangle of K with eigenvalue \lambda yields K (L_+ | \lambda \rangle) = (\lambda + \Delta) (L_+ | \lambda \rangle), and similarly K (L_- | \lambda \rangle) = (\lambda - \Delta) (L_- | \lambda \rangle) for the lowering operator. For unbounded ladders starting from a lowest state where L_- | \lambda_0 \rangle = 0, the eigenvalues form the arithmetic sequence \lambda_n = \lambda_0 + n \Delta with n = 0, 1, 2, \dots, ensuring non-negative spectrum in many physical cases. In bounded representations, the sequence terminates at a highest state where L_+ | \lambda_{\max} \rangle = 0, limiting n to finite values. The full commutation relations depend on the system; for example, the quantum harmonic oscillator follows the Heisenberg-Weyl algebra with [a, a^\dagger] = 1 and [N, a^\pm] = \mp a^\pm (where N is the number operator), while angular momentum follows the su(2) algebra with [J_z, J_\pm] = \pm \hbar J_\pm and [J_+, J_-] = 2 \hbar J_z. Details for specific systems are covered in later sections.^[9] Normalization of the ladder states follows from the commutation relations, with factors determined recursively for each algebra; for instance, in unbounded cases like the harmonic oscillator, successive applications often yield factors like \sqrt{n+1}. States are chosen such that \langle n | n \rangle = 1, preserving unitarity in the representation. In cases where the ladder operators and K generate a Lie algebra (e.g., su(2) or su(1,1)), there exists a Casimir operator that commutes with all generators ([C, L_\pm] = [C, K] = 0) and takes a constant value across the ladder, labeling irreducible representations. The ladder operators generate a complete basis for the representation space, as repeated applications of L_\pm from extremal states span all eigenvectors, with orthogonality \langle n | m \rangle = \delta_{nm} following from the commutation algebra and finite dimensionality in bounded cases.

Angular Momentum

Operator Construction

In quantum mechanics, the total angular momentum operators J_x, J_y, and J_z obey the commutation relations of the su(2) Lie algebra:

[J_i, J_j] = i \hbar \epsilon_{ijk} J_k,

where i, j, k \in \{x, y, z\}, \epsilon_{ijk} is the Levi-Civita symbol, and \hbar is the reduced Planck's constant.^[10] The square of the total angular momentum J^2 = J_x^2 + J_y^2 + J_z^2 commutes with each component: [J^2, J_i] = 0.^[10] Ladder operators are constructed as linear combinations of the Cartesian components:

J_+ = J_x + i J_y, \quad J_- = J_x - i J_y.

These are non-Hermitian, with J_- = J_+^\dagger, and satisfy [J_z, J_\pm] = \pm \hbar J_\pm and [J^2, J_\pm] = 0.^[11] The simultaneous eigenstates of J^2 and J_z, labeled |j, m\rangle, obey J^2 |j, m\rangle = \hbar^2 j(j+1) |j, m\rangle and J_z |j, m\rangle = \hbar m |j, m\rangle, where j \geq 0 is integer or half-integer, and m is the projection quantum number.^[10] The action of the ladder operators on these states shifts m by \pm 1 while preserving j:

J_+ |j, m\rangle = \hbar \sqrt{j(j+1) - m(m+1)} \, |j, m+1\rangle,

J_- |j, m\rangle = \hbar \sqrt{j(j+1) - m(m-1)} \, |j, m-1\rangle.

These relations provide the explicit matrix elements \langle j, m' | J_\pm | j, m \rangle in the |j, m\rangle basis, which are zero unless m' = m \pm 1, with the nonzero elements given by the square roots above.^[10] For instance, in the basis ordered by decreasing m, J_z is diagonal with entries \hbar m, while J_+ and J_- are strictly super- and subdiagonal, respectively.^[11] The structure forms a finite ladder because the operators terminate the chain: J_+ |j, j\rangle = 0 and J_- |j, -j\rangle = 0, restricting m to the values -j, -j+1, \dots, j in integer steps and yielding $2j+1 orthonormal states per j.^[11] When coupling two angular momenta \mathbf{j_1} and \mathbf{j_2} to form total \mathbf{j}, the coupled states |j, m\rangle expand in the uncoupled product basis as

|j, m\rangle = \sum_{m_1, m_2} C^{j m}_{j_1 m_1, j_2 m_2} |j_1, m_1\rangle |j_2, m_2\rangle,

with Clebsch-Gordan coefficients C^{j m}_{j_1 m_1, j_2 m_2} that vanish unless m = m_1 + m_2 and |j_1 - j_2| \leq j \leq j_1 + j_2.^[12] These coefficients are constructed recursively using ladder operators: begin with the highest-weight state |j = j_1 + j_2, m = j\rangle = |j_1, j_1\rangle |j_2, j_2\rangle (where the coefficient is 1), then apply the total lowering operator J_- = J_{1-} + J_{2-} successively to both sides of the expansion, equating coefficients and normalizing at each step to ensure unitarity.^[12] This iterative process generates all coefficients for a given j_1 + j_2 without direct computation of integrals.^[11]

Applications in Atomic and Molecular Physics

In atomic physics, ladder operators J_\pm play a crucial role in determining selection rules for electric dipole transitions between states of different magnetic quantum numbers m. The matrix elements of the dipole operator involve terms proportional to J_\pm, leading to the rule \Delta m = \pm 1 for allowed transitions, while \Delta m = 0 requires the J_z component; this ensures conservation of angular momentum during photon emission or absorption.^[13] These rules govern the polarization and intensity of spectral lines in atomic spectra.^[14] In the Zeeman effect, where an external magnetic field splits degenerate energy levels, ladder operators simplify the computation of matrix elements for the perturbation Hamiltonian -\vec{\mu} \cdot \vec{B}, with \vec{\mu} being the magnetic moment operator. The splitting into $2j + 1 sublevels with energies proportional to m_j g \mu_B B (where g is the Landé g-factor) is analyzed using J_\pm to evaluate transitions between these levels, explaining observed linear and anomalous patterns in atomic spectra.^[15] For molecular rotations, the rigid rotor Hamiltonian H = \frac{\vec{J}^2}{2I} (with I the moment of inertia) uses angular momentum ladder operators to derive selection rules \Delta J = \pm 1 and \Delta m_J = 0, \pm 1 for microwave absorption spectra, predicting evenly spaced lines at frequencies \nu = 2B(J+1) (where B = h/(8\pi^2 I c)) for diatomic molecules with permanent dipoles.^[16] This framework enables precise determination of bond lengths from observed rotational constants. In multi-electron atoms, spin-orbit coupling is treated via the j-j coupling scheme, where individual electron angular momenta \vec{j}_i = \vec{l}_i + \vec{s}_i are coupled using ladder operators to form total \vec{J}; this constructs eigenstates of the coupled Hamiltonian, revealing fine structure splittings from the interaction \sum_i \vec{l}_i \cdot \vec{s}_i.^[17] These methods are experimentally verified in alkali atom spectroscopy, such as the sodium D-line doublet (arising from 3p j=3/2 and j=1/2 levels split by ~17 cm⁻¹), where transitions obey the derived selection rules and match fine structure intervals measured via laser spectroscopy.^[18]

Harmonic Oscillator

One-Dimensional Case

The one-dimensional quantum harmonic oscillator serves as a foundational model in quantum mechanics, with its Hamiltonian operator given by

H = \frac{p^2}{2m} + \frac{1}{2} m \omega^2 x^2,

where p is the momentum operator, x the position operator, m the particle mass, and \omega the angular frequency.^[19] This Hamiltonian describes systems like vibrating molecules or photons in a cavity, and its energy spectrum is discrete, as first derived using early quantum methods. To solve for the eigenstates and eigenvalues algebraically, ladder operators are introduced, building on Dirac's operator formalism.^[20] The lowering (annihilation) operator a is constructed as

a = \sqrt{\frac{m \omega}{2 \hbar}} \, x + \frac{i p}{\sqrt{2 m \omega \hbar}},

with the raising (creation) operator a^\dagger as its Hermitian adjoint:

a^\dagger = \sqrt{\frac{m \omega}{2 \hbar}} \, x - \frac{i p}{\sqrt{2 m \omega \hbar}}.

These operators satisfy the commutation relation [a, a^\dagger] = 1, which follows directly from the canonical commutation relation [x, p] = i \hbar.^[19] Using this structure, the Hamiltonian factorizes compactly as

H = \hbar \omega \left( a^\dagger a + \frac{1}{2} \right),

where a^\dagger a acts as the number operator N. This expression reveals the oscillator's spectrum without solving the differential equation. The ground state |0\rangle is defined by the condition a |0\rangle = 0, yielding the lowest energy eigenvalue E_0 = \frac{1}{2} \hbar \omega.^[19] Higher excited states are generated by successive application of the raising operator: |n\rangle = \frac{(a^\dagger)^n}{\sqrt{n!}} |0\rangle for n = 0, 1, 2, \dots, with corresponding energies E_n = \hbar \omega \left(n + \frac{1}{2}\right).^[20] These states form an orthonormal basis, and the ladder operators shift the quantum number n by \pm 1, enforcing the equidistant energy spacing characteristic of the harmonic oscillator.^[19]

Energy Levels and Coherent States

The number operator for the quantum harmonic oscillator is defined as N = a^\dagger a, where a and a^\dagger are the lowering and raising ladder operators, respectively.^[21] The eigenstates of the Hamiltonian are simultaneous eigenstates of N, denoted |n\rangle, with eigenvalues n = 0, 1, 2, \dots, corresponding to discrete energy levels E_n = \hbar \omega (n + 1/2). This structure reveals the equally spaced energy spectrum inherent to the oscillator, where ladder operators shift states between these levels by units of \hbar \omega.^[21] Coherent states, also known as Glauber states, are right eigenstates of the lowering operator satisfying a |\alpha\rangle = \alpha |\alpha\rangle, where \alpha is a complex eigenvalue. Their explicit expansion in the number basis is

|\alpha\rangle = e^{-|\alpha|^2/2} \sum_{n=0}^\infty \frac{\alpha^n}{\sqrt{n!}} |n\rangle,

a form that underscores their Poissonian photon number distribution. These states minimize the Heisenberg uncertainty relation, achieving \Delta x \Delta p = \hbar/2 with equal variances in position and momentum quadratures, and their expectation values \langle x \rangle and \langle p \rangle oscillate classically with the frequency \omega.^[22] Coherent states can be generated via the displacement operator D(\alpha) = e^{\alpha a^\dagger - \alpha^* a}, which acts on the vacuum as |\alpha\rangle = D(\alpha) |0\rangle, effectively shifting the state in phase space without altering its Gaussian shape. This operator satisfies D^\dagger(\alpha) a D(\alpha) = a + \alpha, highlighting the coherent state's role as a displaced vacuum.^[22] Squeezed states generalize coherent states by employing the non-compact su(1,1) algebra, generated by operators K_+ = (a^\dagger)^2/2, K_- = a^2/2, and K_3 = (N + 1/2)/2, to reduce uncertainty in one quadrature below the vacuum level while increasing it in the conjugate, maintaining the minimum uncertainty product.^[23] These states, produced via processes like parametric down-conversion, exhibit non-classical correlations useful in precision measurements. The set of coherent states forms an overcomplete basis, satisfying the resolution of the identity

\int \frac{d^2\alpha}{\pi} |\alpha\rangle\langle\alpha| = \hat{1},

which enables phase-space representations and quantization schemes, such as the Husimi transform, bridging quantum and classical descriptions.

Hydrogen-like Atom

Laplace–Runge–Lenz Vector

In classical mechanics, the Laplace–Runge–Lenz vector arises in the Kepler problem, describing the motion of a particle under an inverse-square central force, such as gravitational or electrostatic attraction. Defined as \vec{A} = \vec{p} \times \vec{L} - m k \hat{r}, where \vec{p} is the linear momentum, \vec{L} = \vec{r} \times \vec{p} is the angular momentum, m is the mass, k is the coupling constant (e.g., G M m for gravity or Z e^2 for the Coulomb potential), and \hat{r} = \vec{r}/r, this vector points toward the periapsis and has magnitude m k e, with e the eccentricity of the orbit.^[24] Its conservation leads to closed elliptical orbits, a hallmark of the Kepler problem's integrability, distinguishing it from generic central forces that yield rosette patterns.^[25] In quantum mechanics, Wolfgang Pauli adapted the Laplace–Runge–Lenz vector to the hydrogen atom in 1926, recognizing its role in revealing the hidden dynamical symmetry of the Coulomb potential.^[26] The quantum version is symmetrized to ensure Hermiticity: \vec{A} = \frac{1}{2m} (\vec{p} \times \vec{L} - \vec{L} \times \vec{p}) - k \frac{\vec{r}}{r}, where \vec{p} = -i \hbar \nabla and \vec{L} = -i \hbar \vec{r} \times \nabla. For bound states (E < 0), a normalized form is \vec{K} = \frac{\vec{A}}{\sqrt{-2 m E}}.^[27] This operator commutes with the hydrogen Hamiltonian H = \frac{\vec{p}^2}{2m} - \frac{k}{r}, [H, \vec{A}] = 0, ensuring conservation.^[27] The components satisfy [L_i, A_j] = i \hbar \sum_k \epsilon_{ijk} A_k, confirming that \vec{A} transforms as a vector under rotations.^[27] Together with \vec{L}, the normalized \vec{K} generates the Lie algebra so(4) for bound states via \vec{J}_\pm = \frac{1}{2} \left( \vec{L} \pm \vec{K} \right), which satisfy separate su(2) commutation relations: [J_\pm^i, J_\pm^j] = i \hbar \epsilon_{ijk} J_\pm^k and [\vec{J}_+, \vec{J}_-] = 0.^[27] This algebraic structure underlies the accidental degeneracy of hydrogen energy levels, depending only on the principal quantum number n. From the SO(4) Casimir, K^2 + L^2 + \hbar^2 = \hbar^2 n^2, or equivalently K^2 + L^2 = \hbar^2 (n^2 - 1), linking directly to the quantization condition and energy E_n = -\frac{m k^2}{2 \hbar^2 n^2}.^[27] In the classical limit, |\vec{A}| / (m k) = \sqrt{1 - l(l+1)/n^2}, analogous to the eccentricity.^[27]

Factorization of the Hamiltonian

The Schrödinger equation for the hydrogen atom is separable in spherical coordinates, where the wave function takes the form \psi(r, \theta, \phi) = R(r) Y_{l m}(\theta, \phi), with Y_{l m} denoting spherical harmonics. The radial function R(r) then obeys the differential equation

\frac{d^2 R}{dr^2} + \left[ \frac{2 m E}{\hbar^2} + \frac{2 m k}{\hbar^2 r} - \frac{l(l+1)}{r^2} \right] R = 0,

where m is the reduced mass, k = e^2 / (4 \pi \epsilon_0), and E < 0 for bound states.^[28] To solve this algebraically, the factorization method is applied to the radial Hamiltonian H_r = -\frac{\hbar^2}{2m} \frac{d^2}{dr^2} + V_\mathrm{eff}(r), where the effective potential is V_\mathrm{eff}(r) = -\frac{k}{r} + \frac{\hbar^2 l(l+1)}{2 m r^2}. Ladder operators \eta_l^\pm are introduced such that H_r - E = \frac{\hbar^2}{2m} (\eta_l^-)^\dagger \eta_l^- (or similar paired forms for raising and lowering), with the explicit operators given by \eta_l^\pm = \frac{d}{dr} \pm \left( \frac{l+1}{r} + \frac{m k}{\hbar^2 \nu} - \frac{1}{2} \frac{\rho}{r} \right), where \nu relates to a factorization parameter and \rho to the energy scale.^[29] These operators act on the radial wave functions, shifting between solutions for successive values of the orbital quantum number l while preserving the principal quantum number n. The method proceeds by successive factorizations, starting from the highest l = n-1 and lowering to l = 0, yielding solutions in terms of associated Laguerre polynomials.^[29] The components of the Laplace–Runge–Lenz vector enable the construction of additional ladder operators that connect states within the degenerate manifold. Standard angular momentum operators L^\pm shift the magnetic quantum number m while preserving n and l; LRL-derived operators (e.g., K^\pm \propto A_x \pm i A_y) shift both l and m, generating the full set of degenerate states for a given energy level.^[26] For bound states, the ladder must terminate to ensure normalizable wave functions, which occurs after a finite number of applications equal to the principal quantum number n. This termination condition quantizes the energy as E_n = -\frac{m k^2}{2 \hbar^2 n^2}. The resulting degeneracy of each level is n^2, arising from the possible values l = 0, 1, \dots, n-1 and m = -l, \dots, l.^[26]

Relation to Group Theory

The hidden SO(4) symmetry of the hydrogen atom, generated by the angular momentum \vec{L} and the normalized Laplace–Runge–Lenz vector \vec{K}, explains the n^2 degeneracy of energy levels. The generators \vec{J}_\pm = \frac{1}{2} (\vec{L} \pm \vec{K}) form two commuting SU(2) algebras, isomorphic to SO(4) for bound states. The bound states transform under the irreducible representation of dimension n^2, labeled by the highest weight related to the principal quantum number n.^[30] This representation branches under the SO(3) subgroup (generated by \vec{L}) into irreps with l = 0, 1, \dots, n-1, each with $2l+1 states for m = -l, \dots, l, yielding the total degeneracy \sum_{l=0}^{n-1} (2l+1) = n^2. The ladder operators within SO(4), such as combinations of L^\pm and K^\pm, raise and lower both l and m, systematically generating the full multiplet from a reference state (e.g., maximum l = n-1, m = l). This algebraic framework, introduced by Pauli, provides an operator-based solution to the spectrum without solving the differential equation directly.^[26]^[30]

3D Isotropic Harmonic Oscillator

Factorization Method

The three-dimensional isotropic harmonic oscillator is described by the Hamiltonian

H = \frac{\mathbf{p}^2}{2m} + \frac{1}{2} m \omega^2 r^2,

which is separable in spherical coordinates into radial and angular components.^[31] The angular part yields the standard spherical harmonics Y_l^m(\theta, \phi) with eigenvalues l(l+1)\hbar^2 for the angular momentum operator, where l = 0, [1](/page/1), 2, \dots and m = -l, \dots, l. The radial equation for fixed l involves an effective potential consisting of the harmonic term \frac{1}{2} m \omega^2 r^2 and the centrifugal barrier \frac{\hbar^2 l(l+1)}{2m r^2}, resulting in the radial Hamiltonian

H_l = \frac{p_r^2}{2m} + \frac{1}{2} m \omega^2 r^2 + \frac{\hbar^2 l(l+1)}{2m r^2}.

^[32]^[33] The factorization method applies Schrödinger's algebraic technique to this radial Hamiltonian by constructing ladder operators that factor H_l - E into non-Hermitian factors, analogous to the one-dimensional case but modified by the centrifugal term. The lowering (annihilation) operator for fixed l is

\hat{A}_l = \frac{1}{\sqrt{2m}} \left( p_r + i \frac{\hbar (l+1)}{r} - i m \omega r \right),

which annihilates the lowest radial state for that l, and the corresponding raising (creation) operator \hat{A}_l^\dagger generates higher radial excitations. These operators factorize the Hamiltonian via relations such as \hat{A}_l H_l = (H_l - \hbar \omega) \hat{A}_l, accounting for the centrifugal barrier through l-dependent shifts, allowing the radial Hamiltonian to be expressed as H_l = \hbar \omega (\hat{A}_l^\dagger \hat{A}_l + l + 3/2).^[31]^[33] In some formulations, generalized ladder operators couple the radial and angular degrees of freedom, enabling transitions between states while preserving the factorization structure.^[32] Applying these operators iteratively from the radial ground state R_{0,l}(r) \propto r^l e^{-m \omega r^2 / 2 \hbar} yields the full set of radial eigenfunctions, which are associated Laguerre polynomials times the ground-state form. The resulting energy eigenvalues are

E_{n,l} = \hbar \omega \left( 2n + l + \frac{3}{2} \right),

where n = 0, 1, 2, \dots is the radial quantum number.^[31]^[32]^[33] For a given total quantum number N = 2n + l, all states sharing the same N are degenerate with energy E_N = \hbar \omega (N + 3/2), as the ladder operators connect configurations differing in n and l by even steps while maintaining the eigenvalue. The degeneracy for each N is (N+1)(N+2)/2, arising from the possible values of l = N, N-2, \dots, 0 or $1, each with $2l+1 magnetic substates.^[31]^[32]

Relation to Group Theory

The three-dimensional isotropic harmonic oscillator exhibits an underlying SU(3) symmetry, which arises from the invariance of its Hamiltonian under transformations generated by bilinear combinations of creation and annihilation operators. Specifically, the generators of this Lie algebra are the operators Q_{ij} = a_i^\dagger a_j for i,j = 1,2,3, where a_i^\dagger and a_i are the creation and annihilation operators for the i-th Cartesian direction. These operators satisfy the commutation relations of the su(3) algebra, [Q_{ij}, Q_{kl}] = \delta_{jk} Q_{il} - \delta_{il} Q_{kj}, ensuring that the Hamiltonian, expressed as H = \hbar \omega (N + 3/2) with total number operator N = \sum_i a_i^\dagger a_i, is invariant under SU(3) transformations. This symmetry accounts for the observed degeneracy in energy levels beyond the rotational SO(3) subgroup. Within the SU(3) framework, ladder operators emerge as specific combinations that raise or lower the oscillator quanta along different directions, facilitating transitions between states while preserving the total symmetry. These operators, derived from the Cartan-Weyl basis of su(3), include raising and lowering generators such as E_{\alpha} for positive roots α, which act to increase or decrease the occupation numbers in a manner analogous to angular momentum ladders but extended to the full three-dimensional space. The irreducible representations (irreps) of SU(3) relevant to the oscillator are labeled by the highest weight (λ, μ), where for states with total quantum number N (the principal excitation level), the irreps take the form (N, 0). This labeling explains the degeneracy of each energy level E_N = \hbar \omega (N + 3/2), given by the dimension of the (N, 0) irrep: d = \frac{(N+1)(N+2)}{2}, which matches the number of ways to distribute N indistinguishable quanta among three directions. The SU(3) symmetry further decomposes via the subgroup chain SU(3) ⊃ SO(3) ⊃ SO(2), linking the full degeneracy to angular momentum quantum numbers. Under this chain, the (N, 0) irrep branches into SO(3) representations with orbital angular momentum l ranging from N down to 0 (or 1) in steps of 2, each containing states labeled by magnetic quantum number m_l from -l to l. This structure reveals how the accidental degeneracy arises from hidden symmetries beyond the obvious rotational invariance, with SO(2) providing the final quantization along the z-axis.^[34] This SU(3) symmetry has found significant application in nuclear physics through the Elliott model, which describes the shell structure and collective rotational motion of atomic nuclei by treating nucleons in a harmonic oscillator potential with SU(3) classifications. In this model, nuclear states are classified by SU(3) irreps to capture deformation and quadrupole collectivity, providing a microscopic foundation for phenomena like rotational bands in light nuclei.

History

Early Development

The mathematical foundations of ladder operators trace back to 19th-century methods for solving differential equations through recursion relations in the theory of orthogonal polynomials. Edmond Laguerre introduced polynomials in 1879 that satisfy a specific differential equation and can be generated recursively, providing an early algebraic framework for stepping between solutions of varying degrees, akin to modern ladder techniques.^[35] These pre-quantum approaches, developed for classical problems in analysis, laid groundwork for later quantum applications without explicit operator formalism. In the formative years of quantum mechanics, Erwin Schrödinger employed a factorization method in 1926 to solve the harmonic oscillator problem, constructing operators that raise and lower energy levels by systematically relating wavefunctions. This technique, presented in his seminal paper formulating wave mechanics, predated the full integration of operator algebras but anticipated their utility in bound-state spectra. Concurrently, Werner Heisenberg and Pascual Jordan's 1925 formulation of matrix mechanics implicitly incorporated step-like transitions between energy states via non-commuting arrays representing observables, enabling calculations of spectral lines that effectively mimic ladder operations for quantized harmonic motion. Paul Dirac advanced this algebraic structure in 1927 by explicitly defining creation and annihilation operators for the quantum harmonic oscillator within his quantum theory of radiation, treating them as fundamental tools to describe emission and absorption processes.^[36] These operators facilitated the infinite ladder of equidistant energy levels, bridging matrix and wave mechanics. By the 1930s, Eugene Wigner and Carl Eckart systematized ladder operators for angular momentum in atomic spectroscopy, with Eckart's 1930 analysis of vector operators and Wigner's 1931 group-theoretic extension providing a rigorous framework for selection rules and matrix elements in multipole transitions.

Key Contributions

In the post-war era of the 1940s and 1950s, ladder operator techniques saw key refinements through group-theoretic interpretations. Valentin Bargmann's 1954 analysis of irreducible unitary representations of the Lorentz group introduced a realization of the su(1,1) Lie algebra applicable to the harmonic oscillator, enabling algebraic generation of its energy spectrum via non-compact symmetry transformations. Concurrently, Vladimir Fock's work from the 1930s recognized the SO(4) dynamical symmetry of the hydrogen atom, integrating the Laplace-Runge-Lenz vector to explain spectral degeneracies.^[37] Fock's ongoing contributions through the 1950s also laid foundational groundwork for coherent states via the Fock space construction, which provided a multi-particle representation framework, and extended to group-theoretic analyses of quantum systems, including representations that anticipated modern uses of ladder operators in state generation. By the 1960s, Albert Messiah's comprehensive textbook formalized ladder operator methods across quantum mechanics, standardizing their use for solving eigenvalue problems in systems like the harmonic oscillator and angular momentum, thereby influencing pedagogical and research applications. In nuclear physics, John P. Elliott's 1958 SU(3) model applied ladder operators to classify collective excitations in the shell model, capturing rotational and vibrational spectra of nuclei through symmetry-adapted states. Subsequent decades expanded ladder operators into specialized domains. The 1963 Jaynes-Cummings model in quantum optics utilized them to describe resonant interactions between a two-level system and a quantized field mode, forming the basis for cavity quantum electrodynamics. In the 1980s, Edward Witten's formulation of supersymmetric quantum mechanics incorporated supercharges as ladder operators bridging bosonic and fermionic Hilbert spaces, revealing exact solvability and ground-state properties in paired potentials. Ladder operator concepts further evolved into quantum field theory during this period, where mode-specific raising and lowering operators in second quantization underpin the Fock space description of multi-particle states and field excitations.

References

[1]
Ladder Operators (Creation/Annihilation Operators)
Jan 29, 2023 · Ladder Operators are operators that increase or decrease eigenvalue of another operator. There are two types; raising operators and lowering ...What are Ladder Operators · Quantum Harmonic Oscillator...
[2]
[PDF] Operator methods in quantum mechanics - TCM
Such a representation was developed by. Dirac early in the formulation of quantum mechanics. ... The operators a and a† represent ladder operators and have the ...
[3]
[PDF] 17. Ladder Operators - Physics
The study of these continuous vibrating quantum systems is called quantum field theory, and ladder operators are a fundamental tool of quantum field theorists.
[4]
[PDF] An explicit basis of lowering operators for irreducible ... - Math@LSU
One systematic way of obtaining bases is to generalize the ladder operator approach to the representations of SU(2) by using the formalism of lowering operators ...
[5]
Renormalisation with SU(1, 1) coherent states on the LQC Hilbert ...
Aug 12, 2020 · The SU(1, 1) group structure enters our analysis by considering a classical set of phase space functions that generates the Lie algebra su(1, 1) ...
[6]
[PDF] The su(1,1) dynamical algebra from the - arXiv
May 21, 2010 · Abstract. We apply the Schrödinger factorization to construct the ladder operators for hydrogen atom, Mie-type potential, harmonic oscillator.
[7]
Jordan chain - StatLect
A Jordan chain is a set of generalized eigenvectors that are obtained by repeatedly applying a nilpotent operator to the same vector.
[8]
[PDF] Abstract ladder operators and their applications - arXiv
Sep 21, 2021 · In this paper, other than considering specific quantum systems, we will also propose an abstract settings for operators of this kind, and we ...
[9]
[PDF] 3.1. Rotations and Angular Momentum Commutation Relations
Jul 3, 2010 · This chapter is concerned with a systematic treatment of angular momen- tum and related topics. The importance of angular momentum in modern.
[10]
[PDF] Ladder operators approach to representation classification ... - arXiv
Jul 5, 2022 · This approach is based on commutative algebra relations and can be applied to the analysis of irreducible representations of various Lie ...
[11]
[PDF] Chapter 13. Angular Momentum: General Theory
(j + 1). The operators ˆL+ and ˆL− are called ladder operators; ˆL+ is called a raising operator and ˆL−. , a lowering operator. When ˆL+ acts on |j, mi ...
[12]
[PDF] Angular Momentum 1 Angular momentum in Quantum Mechanics
which is the same as we had in (1.16). Using these results, it is now easy to write the ladder operators L±. = Lx ± iLy in spherical coordinates:.
[13]
[PDF] chapter 12 total angular momentum
their construction, so will jump to the matrix without deriving all the zeros. ... 2 Schiff Quantum Mechanics (McGraw–Hill, New York, 1968), 3rd ed., p 218 ...
[14]
None
**Summary of Ladder Operators for Clebsch-Gordan Coefficients (University of Kentucky Physics Lecture Notes)**
[15]
[PDF] 7.1 Dipole Transitions
The ladder operators change the z-component of their corresponding angular momentum by 1, and therefore we know that ∆l = 0, ∆ml = ±1, and ∆ms = ±1. The ...
[16]
[PDF] The Wigner-Eckart Theorem - University of California, Berkeley
The selection rule m = q + m′ in the Clebsch-Gordan coefficient (16) means that many of the integrals vanish, so we have exaggerated the total number of ...Missing: ladder | Show results with:ladder
[17]
[PDF] HSO + HZeeman: Coupled vs. Uncoupled Basis Sets
ladder operators and orthogonality for transformation between basis sets. Necessary to be able to evaluate matrix elements of HZeeman in “coupled basis”.
[18]
[PDF] 5.61_F17 Lecture 18: Rigid Rotor II - MIT OpenCourseWare
Jul 8, 2017 · Selection rule for rotational transitions ... We want to derive all of the properties of all angular momentum operators, operating on angular.
[19]
[PDF] Multiplets and Spin-Orbit Coupling - cond-mat.de
The key tools for constructing multiplet state are the ladder operators of total orbital momentum ... Because of the product rule, the angular momentum operators ...
[20]
[PDF] Fine Structure in Hydrogen and Alkali Atoms
The standard spectroscopic notation for the states of hydrogen-like or alkali atoms is. nℓj, where ℓ is represented by one of the standard symbols, s, p, d ...
[21]
[PDF] A completely algebraic solution of the simple harmonic oscillator
20. Next, Born and. Jordan's 1930 textbook14 was the first to represent the ladder operators in the matrix mechanics formalism, but Dirac's. 1930 textbook7 ...Missing: primary source
[22]
The Principals Of Quantum Mechanics : Dirac. P.a.m - Internet Archive
Oct 11, 2020 · The Principals Of Quantum Mechanics. by: Dirac. P.a.m. Publication date: 1930 ... PDF WITH TEXT download · download 1 file · SINGLE PAGE PROCESSED ...Missing: online | Show results with:online
[23]
Coherent and Incoherent States of the Radiation Field | Phys. Rev.
Glauber, Phys. Rev. Letters 10, 84 (1963); R. J. Glauber, in Proceedings of the Third International Conference on Quantum Electronics, Paris, France, 1963 ...
[24]
[PDF] Roy J. Glauber - Nobel Lecture
The coherent states form a complete set of states in the sense that any state of the mode can be expressed as a suitable sum taken over them. As we have ...<|control11|><|separator|>
[25]
Multimode minimum uncertainty squeezed states - IOPscience
The concept of the single mode 'squeezed state', a minimum uncertainty state with reduced fluctuations in either of the two canonically conjugate variables, is ...Missing: original | Show results with:original
[26]
[PDF] the laplace runge lenz vector - IFSC/USP
Runge repeated the derivation in a popular German text on vector analysis (1919) and was quoted as a reference by W. Lenz in a 1924 paper on quantum.
[27]
The Laplace–Runge–Lenz vector revisited - AIP Publishing
Mar 1, 2003 · We define a conserved dynamical variable α that characterizes the orientation of the orbit in two-dimensional configuration space for the Kepler ...<|control11|><|separator|>
[28]
[PDF] Pauli-1926-On-the-spectrum-of-hydrogen-atom--English-translation ...
The present paper sets out to apply the new theory to an atom having a ... Lenz, Zs. f. Phys. 24 (1924) 197. The numbers designated as ni and ng in ...
[29]
[PDF] Laplace-Runge-Lenz vector in quantum mechanics in ... - arXiv
Sep 18, 2013 · Abstract. The main point of this paper is to examine a ”hidden” dynami- cal symmetry connected with the conservation of Laplace-Runge-Lenz.
[30]
Hydrogen Atom - Galileo
Hence u(ρ)=ρe−ρw0 and the actual radial wave function is this divided by r, and of course suitably normalized. To write the wave function in terms of r, we need ...
[31]
[PDF] Three-Dimensional Isotropic Harmonic Oscillator and SU3
Creation and annihilation operators, symmetry and supersymmetry of the 3D isotropic harmonic oscillator ... ladder operators of the three-dimensional isotropic ...
[32]
[PDF] Employing an operator form of the Rodrigues formula to calculate ...
Jul 23, 2025 · OSCILLATOR. The simplest way to solve the three-dimensional isotropic harmonic oscillator using the factorization method is to use separation ...
[33]
[PDF] Ladder Operators for the Spherical 3D Harmonic Oscillator - SciELO
Nov 9, 2020 · Ladder Operators for the Spherical 3D Harmonic Oscillator. Table 1: Radial eigenfunctions and energy eigenvalues for the spherical harmonic ...
[34]
https://www.semanticscholar.org/paper/Three-Dimensional-Isotropic-Harmonic-Oscillator-and-Fradkin/dec5cb786a2e8de58d8f371176505462b6c620b3
[35]
Edmond Laguerre (1834 - 1886) - Biography - MacTutor
... Laguerre polynomials which are solutions of the Laguerre differential equations. This work came out of his paper published in 1879 which examined. ∫ 0 ∞ e ...
[36]
The quantum theory of the emission and absorption of radiation
The new quantum theory, based on the assumption that the dynamical variables do not obey the commutative law of multiplication, has by now been developed ...