Skip to content

Wigner D-Matrices

Wigner D-matrices are the matrix elements of rotation operators in an angular momentum basis. They make the representation-theoretic statement

SO(3) or SU(2)⟶unitary matrices on VjSO(3)\ \text{or}\ SU(2) \quad\longrightarrow\quad \text{unitary matrices on }V_j

explicit in the basis {∣j,m⟩}m=−jj\{\lvert j,m\rangle\}_{m=-j}^{j}.

This page records the mathematical convention used across angular momentum calculations. The underlying Lie algebra is Angular Momentum Algebra, the spinor double-cover issue is SU(2) versus SO(3), and the physical rotation convention is summarized in Notation and Conventions.

Let VjV_j be the irreducible spin-jj representation space with basis vectors

J2∣j,m⟩=ℏ2j(j+1)∣j,m⟩,Jz∣j,m⟩=ℏm∣j,m⟩.J^2\lvert j,m\rangle = \hbar^2 j(j+1)\lvert j,m\rangle, \qquad J_z\lvert j,m\rangle = \hbar m\lvert j,m\rangle.

For an active rotation RR, let Uj(R)U_j(R) be the corresponding unitary operator on VjV_j. The Wigner D-matrix entries are

Dm′mj(R)=⟨j,m′∣Uj(R)∣j,m⟩.D^j_{m'm}(R) = \langle j,m'\vert U_j(R)\vert j,m\rangle.

The row label is the final JzJ_z label m′m' and the column label is the initial JzJ_z label mm. With this convention,

Uj(R)∣j,m⟩=∑m′=−jj∣j,m′⟩Dm′mj(R).U_j(R)\lvert j,m\rangle = \sum_{m'=-j}^{j} \lvert j,m'\rangle D^j_{m'm}(R).

Because Uj(R)U_j(R) is unitary, Dj(R)D^j(R) is a unitary matrix:

∑m′′=−jjDm′′mj(R)∗Dm′′nj(R)=δmn.\sum_{m''=-j}^{j} D^j_{m''m}(R)^*D^j_{m''n}(R) = \delta_{mn}.

The most common angular momentum convention uses a zz-yy-zz Euler-angle decomposition:

R(α,β,γ)=Rz(α)Ry(β)Rz(γ).R(\alpha,\beta,\gamma) = R_z(\alpha)R_y(\beta)R_z(\gamma).

For active quantum rotations this page uses

Uj(α,β,γ)=exp⁡(−iℏαJz)exp⁡(−iℏβJy)exp⁡(−iℏγJz).U_j(\alpha,\beta,\gamma) = \exp\left(-\frac{i}{\hbar}\alpha J_z\right) \exp\left(-\frac{i}{\hbar}\beta J_y\right) \exp\left(-\frac{i}{\hbar}\gamma J_z\right).

Since the first and last factors are diagonal in the ∣j,m⟩\lvert j,m\rangle basis,

Dm′mj(α,β,γ)=e−im′αdm′mj(β)e−imγ.D^j_{m'm}(\alpha,\beta,\gamma) = e^{-im'\alpha} d^j_{m'm}(\beta) e^{-im\gamma}.

The reduced Wigner matrix, or small d-matrix, is

dm′mj(β)=⟨j,m′∣exp⁡(−iℏβJy)∣j,m⟩.d^j_{m'm}(\beta) = \langle j,m'\vert \exp\left(-\frac{i}{\hbar}\beta J_y\right) \vert j,m\rangle.

Different books sometimes reverse the order of active rotations, use passive coordinate rotations, or place the complex conjugate on the spherical harmonics relation. Always compare the definition of Dm′mjD^j_{m'm} before importing signs from a table.

With the Condon–Shortley phase convention, a standard finite-sum expression is

dm′mj(β)=(j+m′)!(j−m′)!(j+m)!(j−m)!×∑k(−1)k−m′+m(cos⁡β2)2j+m−m′−2k(sin⁡β2)m′−m+2k(j+m−k)! k! (m′−m+k)! (j−m′−k)!.\begin{aligned} d^j_{m'm}(\beta) &= \sqrt{ (j+m')!(j-m')!(j+m)!(j-m)! } \\ &\quad\times \sum_k \frac{ (-1)^{k-m'+m} \left(\cos\frac{\beta}{2}\right)^{2j+m-m'-2k} \left(\sin\frac{\beta}{2}\right)^{m'-m+2k} }{ (j+m-k)!\,k!\,(m'-m+k)!\,(j-m'-k)! }. \end{aligned}

The sum runs over the integers kk for which every factorial argument is a nonnegative integer. This condition is the compact way to encode all allowed terms.

For many calculations, the defining matrix element or a recurrence relation is safer than memorizing the finite sum. The finite sum is best treated as a convention-sensitive table generator.

For j=0j=0,

D0(R)=1.D^0(R)=1.

For j=1/2j=1/2, use the ordered basis

(∣12,12⟩,∣12,−12⟩).\left( \lvert\tfrac12,\tfrac12\rangle, \lvert\tfrac12,-\tfrac12\rangle \right).

The small d-matrix is

d1/2(β)=(cos⁡(β/2)−sin⁡(β/2)sin⁡(β/2)cos⁡(β/2)).d^{1/2}(\beta) = \begin{pmatrix} \cos(\beta/2) & -\sin(\beta/2)\\ \sin(\beta/2) & \cos(\beta/2) \end{pmatrix}.

The full Wigner matrix is obtained by multiplying the first row by e−iα/2e^{-i\alpha/2} or eiα/2e^{i\alpha/2} according to m′m', and the first or second column by e−iγ/2e^{-i\gamma/2} or eiγ/2e^{i\gamma/2} according to mm.

For j=1j=1, use the ordered basis

(∣1,1⟩,∣1,0⟩,∣1,−1⟩).\left( \lvert1,1\rangle, \lvert1,0\rangle, \lvert1,-1\rangle \right).

Then

d1(β)=(1+cos⁡β2−sin⁡β21−cos⁡β2sin⁡β2cos⁡β−sin⁡β21−cos⁡β2sin⁡β21+cos⁡β2).d^1(\beta) = \begin{pmatrix} \frac{1+\cos\beta}{2} & -\frac{\sin\beta}{\sqrt2} & \frac{1-\cos\beta}{2} \\ \frac{\sin\beta}{\sqrt2} & \cos\beta & -\frac{\sin\beta}{\sqrt2} \\ \frac{1-\cos\beta}{2} & \frac{\sin\beta}{\sqrt2} & \frac{1+\cos\beta}{2} \end{pmatrix}.

This is the spin-11 representation of the same rotation, expressed in the spherical basis rather than the Cartesian basis.

Wigner D-matrices are representation matrices. If R1R_1 and R2R_2 are rotations, then

Dm′mj(R1R2)=∑n=−jjDm′nj(R1)Dnmj(R2).D^j_{m'm}(R_1R_2) = \sum_{n=-j}^{j} D^j_{m'n}(R_1)D^j_{nm}(R_2).

The inverse is the Hermitian adjoint:

Dm′mj(R−1)=Dmm′j(R)∗.D^j_{m'm}(R^{-1}) = D^j_{mm'}(R)^*.

The character of the spin-jj representation depends only on the rotation angle ω\omega:

χj(R)=∑m=−jjDmmj(R)=sin⁡((j+12)ω)sin⁡(ω/2).\chi_j(R) = \sum_{m=-j}^{j}D^j_{mm}(R) = \frac{ \sin\left((j+\frac12)\omega\right) }{ \sin(\omega/2) }.

This character is useful when decomposing representations, checking degeneracies, or comparing SU(2)SU(2) and SO(3)SO(3) representations.

Spherical harmonics are closely related to Wigner D-matrices. With the Condon–Shortley convention used here,

Yℓm(θ,ϕ)=2ℓ+14π Dm0ℓ(ϕ,θ,0)∗.Y_\ell^m(\theta,\phi) = \sqrt{\frac{2\ell+1}{4\pi}}\, D^{\ell}_{m0}(\phi,\theta,0)^*.

Equivalently,

Dm0ℓ(ϕ,θ,0)=4π2ℓ+1 Yℓm(θ,ϕ)∗.D^{\ell}_{m0}(\phi,\theta,0) = \sqrt{\frac{4\pi}{2\ell+1}}\, Y_\ell^m(\theta,\phi)^*.

This identity is one reason D-matrices appear in orbital angular momentum, rigid rotor problems, molecular rotations, multipole expansions, and angular-momentum coupling.

The same matrices control spherical tensor operators. If TqkT^k_q is a rank-kk spherical tensor component, then

U(R)TqkU(R)†=∑q′=−kkDq′qk(R)Tq′k.U(R)T^k_qU(R)^\dagger = \sum_{q'=-k}^{k} D^k_{q'q}(R)T^k_{q'}.

This transformation law is the representation-theoretic input behind the Wigner–Eckart theorem and selection rules. The matrix elements of TqkT^k_q are constrained because states and operators transform in compatible SU(2)SU(2) representations.

Wigner D-matrices appear whenever a calculation needs the explicit action of rotations on angular momentum states:

  • rotating spin states or preparing spin coherent states;
  • transforming spherical harmonics under a rotated coordinate axis;
  • solving rigid rotor and molecular spectroscopy problems;
  • expressing angular distributions in scattering and decay;
  • deriving selection rules for tensor operators;
  • relating laboratory-frame and body-frame angular momentum bases;
  • building numerical rotation matrices for finite-dimensional spin systems.

They are not new dynamical laws. They are the representation matrices for rotations in a chosen angular momentum basis.

  • Mixing active and passive rotation conventions without taking an inverse or complex conjugate.
  • Reading Dm′mjD^j_{m'm} with the row and column labels reversed.
  • Forgetting the phase factors e−im′αe^{-im'\alpha} and e−imγe^{-im\gamma} around the small d-matrix.
  • Importing a table that uses a different Euler-angle order.
  • Treating the j=1j=1 spherical-basis matrix as the same array as the Cartesian 3×33\times3 rotation matrix.
  • Forgetting that half-integer jj is a representation of SU(2)SU(2) and only projectively represents SO(3)SO(3).
  • A. R. Edmonds, Angular Momentum in Quantum Mechanics, Princeton University Press, 1957.
  • D. M. Brink and G. R. Satchler, Angular Momentum, 3rd ed., Oxford University Press, 1993.
  • D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum Theory of Angular Momentum, World Scientific, 1988.
  • E. P. Wigner, Group Theory and Its Application to the Quantum Mechanics of Atomic Spectra, Academic Press, 1959.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  1. Show that the Euler-angle factorization gives the phase form of Dm′mj(α,β,γ)D^j_{m'm}(\alpha,\beta,\gamma).
Solution

Insert the identity in the ∣j,n⟩\lvert j,n\rangle basis between the three exponential factors:

Dm′mj=⟨j,m′∣e−iαJz/ℏe−iβJy/ℏe−iγJz/ℏ∣j,m⟩=e−im′α⟨j,m′∣e−iβJy/ℏ∣j,m⟩e−imγ.\begin{aligned} D^j_{m'm} &= \langle j,m'\vert e^{-i\alpha J_z/\hbar} e^{-i\beta J_y/\hbar} e^{-i\gamma J_z/\hbar} \vert j,m\rangle \\ &= e^{-im'\alpha} \langle j,m'\vert e^{-i\beta J_y/\hbar}\vert j,m\rangle e^{-im\gamma}. \end{aligned}

The middle matrix element is dm′mj(β)d^j_{m'm}(\beta).

  1. Verify that the j=1/2j=1/2 small d-matrix is unitary.
Solution

Let c=cos⁡(β/2)c=\cos(\beta/2) and s=sin⁡(β/2)s=\sin(\beta/2). Then

d1/2(β)=(c−ssc).d^{1/2}(\beta) = \begin{pmatrix} c & -s\\ s & c \end{pmatrix}.

The columns have norms c2+s2=1c^2+s^2=1 and their inner product is −cs+sc=0-cs+sc=0. Hence the matrix is orthogonal and therefore unitary.

  1. Use the spherical harmonic relation to express D00ℓ(ϕ,θ,0)D^\ell_{00}(\phi,\theta,0).
Solution

For m=0m=0,

Yℓ0(θ,ϕ)=2ℓ+14π D00ℓ(ϕ,θ,0)∗.Y_\ell^0(\theta,\phi) = \sqrt{\frac{2\ell+1}{4\pi}}\, D^\ell_{00}(\phi,\theta,0)^*.

The m=0m=0 spherical harmonic is real in the standard convention, so

D00ℓ(ϕ,θ,0)=4π2ℓ+1 Yℓ0(θ,ϕ)=Pℓ(cos⁡θ).D^\ell_{00}(\phi,\theta,0) = \sqrt{\frac{4\pi}{2\ell+1}}\, Y_\ell^0(\theta,\phi) = P_\ell(\cos\theta).
  1. Explain why a spin-1/21/2 D-matrix changes sign under a 2π2\pi rotation while a spin-11 D-matrix does not.
Solution

For a rotation about zz,

Dm′mj(0,0,γ)=δm′me−imγ.D^j_{m'm}(0,0,\gamma) = \delta_{m'm}e^{-im\gamma}.

If γ=2π\gamma=2\pi and j=1/2j=1/2, then m=±1/2m=\pm1/2 and both diagonal entries equal −1-1. If j=1j=1, then m=1,0,−1m=1,0,-1 and all diagonal entries equal 11. This reflects the fact that half-integer spin is a true representation of SU(2)SU(2) but only a projective representation of SO(3)SO(3).