Physical Applications of Group Theory

Group theory earns its place in physics the moment you notice that a Hamiltonian commutes with a set of symmetry operations. If H^\hat{H} is invariant under a group GG of unitary operators {U(g)}\{U(g)\}, then [H^,U(g)]=0[\hat{H},U(g)]=0 for every gg, and each U(g)U(g) carries an eigenspace of H^\hat{H} into itself. Those eigenspaces therefore carry representations of GG, and — barring accidental degeneracy — the representations are irreducible. That single sentence converts the vague word "symmetry" into hard, countable predictions: which levels are degenerate, how a perturbation splits them, which spectral lines are allowed, and how many molecular vibrations are infrared or Raman active. No wavefunction is solved along the way; the answers are read off character tables. The engine that does the reading is the character — the trace of a representation matrix. Because traces are invariant under change of basis and constant on conjugacy classes, decomposing any reducible representation into its irreducible pieces reduces to a single weighted inner product of two short lists of numbers. This "magic formula" splits a molecule's 3N3N displacement representation into translations, rotations, and genuine vibrations, and splits an atomic multiplet dropped into a crystal into the irreps of the crystal site. Selection rules follow from the same arithmetic: a transition amplitude is an integral over all space, and by Schur's lemma it vanishes unless its integrand contains the totally symmetric representation. The framework is uniform across the discrete point groups of spectroscopy and the continuous rotation group SO(3) of angular momentum — the same characters, orthogonality, and decomposition — which is why crystal-field splittings are computed by subducing SO(3) irreps to a finite point group.

Key takeaways
  • A symmetry group GG of H^\hat{H} organizes the Hilbert space into irreducible representations: energy levels are degenerate precisely because they span an irrep, and the degeneracy equals the irrep's dimension.
  • Any reducible representation is decomposed into irreps by one weighted character inner product — the magic formula nμ=1∣G∣∑CNC χ(μ)(C)∗χ(C)n_\mu=\frac{1}{|G|}\sum_C N_C\,\chi^{(\mu)}(C)^\ast\chi(C) — with no diagonalization required.
  • A matrix element ⟨f∣O^∣i⟩\langle f|\hat{O}|i\rangle is forced to zero unless the trivial representation appears in Γf∗⊗ΓO⊗Γi\Gamma_f^\ast\otimes\Gamma_O\otimes\Gamma_i — the origin of every optical selection rule.
  • Molecular normal modes and their IR/Raman activity come from building Γ3N\Gamma_{3N} from unmoved-atom counts and subtracting the 3+33+3 translations and rotations; H₂O gives Γvib=2A1⊕B1\Gamma_{\text{vib}}=2A_1\oplus B_1.
  • The same machinery handles continuous groups: SO(3) has irreps D(ℓ)D^{(\ell)} of dimension 2ℓ+12\ell+1, and subducing D(2)D^{(2)} to the octahedral group gives the ligand-field splitting d→eg⊕t2gd\to e_g\oplus t_{2g}.
  • Key idea: A symmetry GG of H^\hat{H} block-diagonalizes the problem by irrep — degeneracies are irrep dimensions, splittings are decompositions of subduced or perturbed representations, and a matrix element survives only if its integrand contains the trivial representation. Everything reduces to character arithmetic.

Where this lives today

The character-table calculations in this article are run, in effect, every day. Vibrational spectroscopists still build Γvib\Gamma_{\text{vib}} to assign infrared and Raman bands, and the mutual exclusion rule remains a first-line diagnostic for a centre of symmetry in an unknown compound. Ligand-field theory — the subduction of D(ℓ)D^{(\ell)} to a site point group — underlies the colours and magnetic moments of transition-metal and lanthanide complexes. In the solid state, the same little-group analysis at high-symmetry points of the Brillouin zone classifies electronic bands and phonons by their irreps; multidimensional irreps force the degeneracies and symmetry-protected crossings that define topological insulators and Dirac and Weyl semimetals, whose entire classification is a group-theory exercise. The Wigner–Eckart theorem is the backbone of atomic-clock and trapped-ion quantum-computing calculations, where transition strengths between sublevels reduce to Clebsch–Gordan coefficients. And in particle physics the very notion of a particle is an irreducible representation of a symmetry group — the idea, scaled from SO(3) to the Poincaré and gauge groups, that Wigner set in motion.

The physical picture

Picture the Hamiltonian as a big matrix that must look identical after any operation of the symmetry group. Change to a basis of functions that transform "cleanly" — each set spanning an irrep — and that matrix cannot connect functions belonging to different irreps, because a symmetry-invariant operator has no invariant way to link them. The matrix therefore falls apart into independent blocks, one per irrep, and the states inside a single block are automatically degenerate. Symmetry performs this block-diagonalization for free, before any dynamics: the block sizes are fixed by group theory alone, and only each block's energy is left to the physics. A doubly degenerate EE level or a triply degenerate TT level exists because the group offers no operation that singles out one partner state from the others.

Figure
A symmetry group block-diagonalizes the Hamiltonian: in an irrep-adapted basis a fully coupled 7×7 matrix decouples into blocks 2A₁ ⊕ E ⊕ T₂ of sizes 1, 1, 2 and 3, so seven states collapse to four distinct energy levels whose degeneracies are exactly the irrep dimensions
Symmetry does the block-diagonalization for free: degeneracy = irrep dimension Ĥ, generic basis all elements coupled irrep-adapted basis Ĥ, symmetry basis A₁ · A₁ · E · T₂ → degeneracies 1, 1, 2, 3 7 states = 1 + 1 + 2 + 3 → 4 distinct energy levels off-diagonal blocks vanish because Ĥ commutes with every U(g)

A matrix element tells the complementary story. It is an integral over all space, and if you rotate the entire picture by a symmetry of the system the integral cannot change its value. But the only piece of an integrand that survives averaging over the whole group is the piece that is itself invariant — the totally symmetric representation. So a transition is allowed exactly when the product of "how the bra, the operator, and the ket transform" contains that invariant component; otherwise the integral is not merely small, it is identically zero. Lower the symmetry with a perturbation, and a level splits by exactly the way its irrep decomposes when restricted to the smaller group — which is what the crystal-field diagram below computes for a dd electron placed at the centre of an octahedron.

Figure
A free-ion d level with ℓ = 2 is five-fold degenerate; an octahedral field splits it into a lower triplet t₂g and an upper doublet eg separated by 10Dq, and because the barycentre of the level is preserved t₂g drops by 4Dq while eg rises by 6Dq so that 3×(−4Dq) + 2×(+6Dq) = 0
Subduce D⁽²⁾ from SO(3) to the octahedral group: the ligand-field splitting d, ℓ = 2 D⁽²⁾, 5-fold barycentre eg (2-fold) +6Dq { d_z² , d_x²−y² } t₂g (3-fold) −4Dq { d_xy , d_yz , d_zx } 10Dq 2 + 3 = 5 and 2×(+6Dq) + 3×(−4Dq) = 0

The core equations

The starting point is the defining statement of a symmetry: the Hamiltonian is unchanged by every group operation, so it commutes with each unitary representative. This is what forces energy eigenspaces to be invariant subspaces, and hence to carry representations of GG.

[H^,U(g)]=0∀ g∈G(1)[\hat{H},U(g)]=0\quad\forall\,g\in G \qquad (1)

where: H^\hat{H} — the Hamiltonian; U(g)U(g) — the unitary operator representing group element gg on the Hilbert space; GG — the symmetry group. Reading it: if H^ψ=Eψ\hat{H}\psi=E\psi then H^ U(g)ψ=U(g)H^ψ=E U(g)ψ\hat{H}\,U(g)\psi=U(g)\hat{H}\psi=E\,U(g)\psi, so U(g)ψU(g)\psi is another eigenstate at the same energy — the eigenspace is closed under GG. Everything computational rests on decomposing a reducible representation Γ\Gamma into irreps. Because characters are basis-independent and constant on conjugacy classes, the multiplicity of each irrep is a single inner product, best summed class-by-class with the class sizes as weights.

nμ=1∣G∣∑g∈Gχ(μ)(g)∗ χ(g)=1∣G∣∑CNC χ(μ)(C)∗ χ(C)(2)n_\mu=\frac{1}{|G|}\sum_{g\in G}\chi^{(\mu)}(g)^\ast\,\chi(g)=\frac{1}{|G|}\sum_{C}N_C\,\chi^{(\mu)}(C)^\ast\,\chi(C) \qquad (2)

where: nμn_\mu — how many times irrep μ\mu appears in Γ\Gamma; ∣G∣|G| — the group order; χ(g)\chi(g) — character of the reducible representation; χ(μ)(g)\chi^{(\mu)}(g) — character of irrep μ\mu; NCN_C — number of elements in class CC. Reading it: dot the reducible character into each irreducible one, weight by class size, divide by ∣G∣|G|; the answer is always a non-negative integer — a built-in check. Selection rules are the second pillar and follow from the same orthogonality. A matrix element cannot depend on the choice of frame, so it must equal its own group average, and that average projects onto the trivial (totally symmetric) representation Γ1\Gamma_1.

⟨f∣O^∣i⟩=0unlessΓ1⊂Γf∗⊗ΓO⊗Γi(3)\langle f|\hat{O}|i\rangle=0\quad\text{unless}\quad \Gamma_1\subset\Gamma_f^\ast\otimes\Gamma_O\otimes\Gamma_i \qquad (3)

where: Γ1\Gamma_1 — the trivial representation; Γi,Γf\Gamma_i,\Gamma_f — irreps of the initial and final states; ΓO\Gamma_O — irrep under which the operator transforms; ⊗\otimes — direct product of representations. Reading it: form the triple product, decompose it with equation (2), and check whether n1≥1n_1\ge 1; if the invariant is absent the amplitude is exactly zero. To use equation (3) you need the symmetry of the operator. The electric-dipole operator is proportional to the position vector, so it transforms as the coordinates (x,y,z)(x,y,z); the Raman (polarizability) operator transforms as the quadratic forms.

ΓOE1=Γ(x,y,z),ΓORaman=Γ(x2,y2,z2,xy,yz,zx)(4)\Gamma_O^{\text{E1}}=\Gamma(x,y,z),\qquad \Gamma_O^{\text{Raman}}=\Gamma(x^2,y^2,z^2,xy,yz,zx) \qquad (4)

where: E1 — electric-dipole process, O^∝r\hat{O}\propto\mathbf{r}; Raman — polarizability process, O^\hat{O} quadratic in coordinates. Reading it: a mode or transition is infrared active if its irrep appears among (x,y,z)(x,y,z) in the character table's basis column, and Raman active if it appears among the quadratic forms. For molecular vibrations the reducible representation is the 3N3N-dimensional space of Cartesian displacements. Only atoms left fixed by an operation contribute to its character, and each contributes the trace of the 3×33\times3 orthogonal matrix acting on its (x,y,z)(x,y,z).

χ3N(g)=Nunmoved(g) (±1+2cos⁡θg)(5)\chi_{3N}(g)=N_{\text{unmoved}}(g)\,\big(\pm 1+2\cos\theta_g\big) \qquad (5)

where: Nunmoved(g)N_{\text{unmoved}}(g) — atoms left in place by gg; θg\theta_g — its rotation angle; the sign is ++ for proper operations (E,CnE,C_n) and −- for improper ones (σ,Sn,i\sigma,S_n,i). Reading it: a proper rotation contributes 1+2cos⁡θ1+2\cos\theta per fixed atom, an improper one −1+2cos⁡θ-1+2\cos\theta; multiply by the count of fixed atoms and you have χ3N\chi_{3N} class by class. The physical vibrations are what is left after removing the three rigid translations and three rigid rotations, whose irreps are read straight off the character table's basis columns.

Γvib=Γ3N⊖Γtrans⊖Γrot(6)\Gamma_{\text{vib}}=\Gamma_{3N}\ominus\Gamma_{\text{trans}}\ominus\Gamma_{\text{rot}} \qquad (6)

where: Γ3N\Gamma_{3N} — full displacement representation; Γtrans\Gamma_{\text{trans}} — irreps spanning (x,y,z)(x,y,z); Γrot\Gamma_{\text{rot}} — irreps spanning (Rx,Ry,Rz)(R_x,R_y,R_z); ⊖\ominus — subtraction of multiplicities. Reading it: what remains has dimension 3N−63N-6 (or 3N−53N-5 for a linear molecule), the correct count of normal modes. The discrete story carries over unchanged to the continuous rotation group SO(3), whose irreps are the angular-momentum multiplets D(ℓ)D^{(\ell)}. A rotation by θ\theta is diagonal in the ∣ℓm⟩|\ell m\rangle basis, and summing the diagonal phases gives a compact closed form.

χ(ℓ)(θ)=∑m=−ℓℓeimθ=sin⁡ ⁣[(ℓ+12)θ]sin⁡(θ/2)(7)\chi^{(\ell)}(\theta)=\sum_{m=-\ell}^{\ell}e^{im\theta}=\frac{\sin\!\big[(\ell+\tfrac12)\theta\big]}{\sin(\theta/2)} \qquad (7)

where: ℓ\ell — angular-momentum quantum number, with dim⁡D(ℓ)=2ℓ+1\dim D^{(\ell)}=2\ell+1; θ\theta — rotation angle; mm — magnetic quantum number. Reading it: evaluating this at the rotation angles of a finite point group is exactly what subduces an atomic multiplet to a crystal site. Finally, combining two angular momenta is the decomposition of a tensor product of SO(3) irreps — the Clebsch–Gordan series, the continuous-group version of equation (2).

D(j1)⊗D(j2)=⨁J=∣j1−j2∣j1+j2D(J)(8)D^{(j_1)}\otimes D^{(j_2)}=\bigoplus_{J=|j_1-j_2|}^{j_1+j_2}D^{(J)} \qquad (8)

where: j1,j2j_1,j_2 — the angular momenta combined; JJ — total angular momentum, stepping by 1. Reading it: the dimensions must balance, (2j1+1)(2j2+1)=∑J(2J+1)(2j_1+1)(2j_2+1)=\sum_J(2J+1), another automatic check.

Derivation ASelection rules from the vanishing-integral theorem

The goal is to prove, from symmetry alone, that a matrix element I=⟨f∣O^∣i⟩=∫ψf∗ O^ ψi dτI=\langle f|\hat{O}|i\rangle=\int\psi_f^\ast\,\hat{O}\,\psi_i\,\mathrm{d}\tau is forced to vanish unless the integrand contains the trivial representation, then to specialize to the electric dipole and recover the atomic rule Δℓ=±1\Delta\ell=\pm1. The whole argument is that an integral over all space is a number, and a number cannot change under a relabelling of coordinates by a symmetry of the system, so II must equal its own average over the group — a projection onto the invariant part of the integrand.

Step 1. Because the volume measure dτ\mathrm{d}\tau is invariant under every U(g)U(g), the integral is unchanged if we act with gg on all three factors and then average over the group: I=1∣G∣∑g∫[U(g)ψf]∗ [U(g)O^U(g)−1] [U(g)ψi] dτI=\frac{1}{|G|}\sum_g\int[U(g)\psi_f]^\ast\,[U(g)\hat{O}U(g)^{-1}]\,[U(g)\psi_i]\,\mathrm{d}\tau. Why: this rewrites the single integral as the group average of itself, which is the trick that exposes what part of the integrand can survive.

Step 2. The product ψf∗ O^ ψi\psi_f^\ast\,\hat{O}\,\psi_i is an object built from three pieces living in Γf∗\Gamma_f^\ast, ΓO\Gamma_O and Γi\Gamma_i, so under GG it transforms in the tensor product Γf∗⊗ΓO⊗Γi\Gamma_f^\ast\otimes\Gamma_O\otimes\Gamma_i. Why: products of representation-carrying quantities transform in the product representation, which is the only fact about tensor products we need.

Step 3. Averaging a representation over the group is the projector onto its trivial-irrep content, because 1∣G∣∑gχ(μ)(g)=δμ,1\frac{1}{|G|}\sum_g\chi^{(\mu)}(g)=\delta_{\mu,1} — every non-trivial irrep sums to zero by orthogonality with the constant character. Why: the group average annihilates everything that is not invariant, leaving only the totally symmetric part standing.

Step 4. Therefore I∝n1(Γf∗⊗ΓO⊗Γi)I\propto n_1\big(\Gamma_f^\ast\otimes\Gamma_O\otimes\Gamma_i\big): if the trivial representation does not appear in the triple product, the projection is zero and I=0I=0 identically, not merely small. Why: this is the selection rule in its most general form, valid for any group and any operator.

Step 5. Specialize to the electric dipole, O^=r\hat{O}=\mathbf{r}, which transforms as D(1)D^{(1)} for a free atom. Then D(1)⊗D(ℓ)=D(ℓ−1)⊕D(ℓ)⊕D(ℓ+1)D^{(1)}\otimes D^{(\ell)}=D^{(\ell-1)}\oplus D^{(\ell)}\oplus D^{(\ell+1)}, so a non-zero amplitude to a state D(ℓ′)D^{(\ell')} requires ℓ′∈{ℓ−1,ℓ,ℓ+1}\ell'\in\{\ell-1,\ell,\ell+1\}; parity (r\mathbf{r} is odd) then kills Δℓ=0\Delta\ell=0, leaving Δℓ=±1\Delta\ell=\pm1. Why: the abstract rule collapses to the textbook atomic result the instant you feed it the rotation group and the dipole operator. *The abstract rule specializes cleanly: for atoms it returns Δℓ=±1, Δm=0,±1\Delta\ell=\pm1,\ \Delta m=0,\pm1, and for molecules it becomes the point-group rule that a mode is IR active if and only if its irrep lies in Γ(x,y,z)\Gamma(x,y,z).*

Go deeper · Derivation B · 2 minNormal modes and IR/Raman activity of water
Derivation BNormal modes and IR/Raman activity of water

The goal is to count and label every vibration of H₂O and decide which are infrared and which Raman active, using nothing but the C2vC_{2v} character table. Water has N=3N=3 atoms and point group C2v={E,C2,σv(xz),σv′(yz)}C_{2v}=\{E,C_2,\sigma_v(xz),\sigma_v'(yz)\} with ∣G∣=4|G|=4; we place the molecule in the xzxz plane with the C2C_2 axis along zz, the standard convention that makes the antisymmetric stretch come out as B1B_1. The plan is to build the 3N=93N=9 displacement character from unmoved-atom counts, decompose it with the magic formula, and strip away the translations and rotations.

Step 1. Count atoms left in place by each operation: EE fixes all three; C2C_2 (the zz axis) fixes only the oxygen; the molecular plane σv(xz)\sigma_v(xz) fixes all three; the perpendicular plane σv′(yz)\sigma_v'(yz) fixes only the oxygen. So Nunmoved=(3,1,3,1)N_{\text{unmoved}}=(3,1,3,1). Why: only atoms that stay put contribute to the character, so this count is the whole geometric input.

Step 2. Multiply by the per-atom contribution ±1+2cos⁡θ\pm1+2\cos\theta: for EE it is 1+2cos⁡0=31+2\cos0=3, for C2C_2 it is 1+2cos⁡π=−11+2\cos\pi=-1, and for each reflection it is −1+2cos⁡0=1-1+2\cos0=1. Hence χ3N=(3⋅3, 1⋅(−1), 3⋅1, 1⋅1)=(9,−1,3,1)\chi_{3N}=(3\cdot3,\ 1\cdot(-1),\ 3\cdot1,\ 1\cdot1)=(9,-1,3,1). Why: this turns the geometry of Step 1 into the reducible character that equation (2) will digest.

Step 3. Apply the magic formula against the four C2vC_{2v} irreps A1(1,1,1,1)A_1(1,1,1,1), A2(1,1,−1,−1)A_2(1,1,-1,-1), B1(1,−1,1,−1)B_1(1,-1,1,-1), B2(1,−1,−1,1)B_2(1,-1,-1,1), each class of size 1: nA1=14(9−1+3+1)=3n_{A_1}=\tfrac14(9-1+3+1)=3, nA2=14(9−1−3−1)=1n_{A_2}=\tfrac14(9-1-3-1)=1, nB1=14(9+1+3−1)=3n_{B_1}=\tfrac14(9+1+3-1)=3, nB2=14(9+1−3+1)=2n_{B_2}=\tfrac14(9+1-3+1)=2, giving Γ3N=3A1⊕A2⊕3B1⊕2B2\Gamma_{3N}=3A_1\oplus A_2\oplus 3B_1\oplus 2B_2. Why: the nine displacement coordinates are now sorted into symmetry species, and 3+1+3+2=93+1+3+2=9 confirms the dimension.

Step 4. Remove the rigid motions read from the table: Γtrans=(z,x,y)=A1⊕B1⊕B2\Gamma_{\text{trans}}=(z,x,y)=A_1\oplus B_1\oplus B_2 and Γrot=(Rz,Ry,Rx)=A2⊕B1⊕B2\Gamma_{\text{rot}}=(R_z,R_y,R_x)=A_2\oplus B_1\oplus B_2. Subtracting, Γvib=(3−1)A1⊕(1−1)A2⊕(3−1−1)B1⊕(2−1−1)B2=2A1⊕B1\Gamma_{\text{vib}}=(3-1)A_1\oplus(1-1)A_2\oplus(3-1-1)B_1\oplus(2-1-1)B_2=2A_1\oplus B_1. Why: the three translations and three rotations carry no vibrational information, and removing them leaves exactly 3N−6=33N-6=3 genuine modes.

Step 5. Read off the activities: A1A_1 carries zz (IR) and x2,y2,z2x^2,y^2,z^2 (Raman); B1B_1 carries xx (IR) and xzxz (Raman). Both species are therefore simultaneously infrared and Raman active — the symmetric stretch and bend (A1A_1) and the antisymmetric stretch (B1B_1) all appear in both spectra. Why: with no inversion centre nothing forbids mutual activity, and this is exactly what is measured for water. *As a check: χ3N=(9,−1,3,1)\chi_{3N}=(9,-1,3,1) decomposes to 3A1⊕A2⊕3B1⊕2B23A_1\oplus A_2\oplus 3B_1\oplus 2B_2, and after removing translations and rotations Γvib=2A1⊕B1\Gamma_{\text{vib}}=2A_1\oplus B_1 — three fundamentals, all IR and Raman active, in agreement with the observed bands at 3657, 1595 and 3756 cm⁻¹.*

Go deeper · Derivation C · 1 minSplitting of a d level in an octahedral crystal field
Derivation CSplitting of a d level in an octahedral crystal field

The goal is to compute how a five-fold degenerate free-ion dd level splits when the ion sits at the centre of an octahedron, by subducing the SO(3) irrep D(2)D^{(2)} to the finite octahedral rotation group OO. This is the founding calculation of ligand-field theory: it turns the continuous rotational symmetry of a free atom into the discrete cubic symmetry of a crystal site and reads the level pattern off a character table. The group OO has order 24 and classes E, 8C3, 3C2, 6C4, 6C2′E,\,8C_3,\,3C_2,\,6C_4,\,6C_2' with rotation angles 0, 2π/3, π, π/2, π0,\,2\pi/3,\,\pi,\,\pi/2,\,\pi.

Step 1. Evaluate the SO(3) character χ(2)(θ)=sin⁡(5θ/2)/sin⁡(θ/2)\chi^{(2)}(\theta)=\sin(5\theta/2)/\sin(\theta/2) at each class angle. At θ=0\theta=0 it is the dimension, 55; at 2π/32\pi/3, sin⁡(5π/3)/sin⁡(π/3)=−1\sin(5\pi/3)/\sin(\pi/3)=-1; at π\pi, sin⁡(5π/2)/sin⁡(π/2)=1\sin(5\pi/2)/\sin(\pi/2)=1; at π/2\pi/2, sin⁡(5π/4)/sin⁡(π/4)=−1\sin(5\pi/4)/\sin(\pi/4)=-1. So over (E,8C3,3C2,6C4,6C2′)(E,8C_3,3C_2,6C_4,6C_2') the characters are χ(2)=(5,−1,1,−1,1)\chi^{(2)}=(5,-1,1,-1,1). Why: these five numbers are the reducible character of the dd level regarded as a representation of the smaller group OO.

Step 2. Decompose with the magic formula, using class sizes (1,8,3,6,6)(1,8,3,6,6) and ∣G∣=24|G|=24. Against E:(2,−1,2,0,0)E:(2,-1,2,0,0) one gets nE=124(1⋅2⋅5+8⋅(−1)(−1)+3⋅2⋅1+0+0)=2424=1n_E=\tfrac{1}{24}(1\cdot2\cdot5+8\cdot(-1)(-1)+3\cdot2\cdot1+0+0)=\tfrac{24}{24}=1. Why: the doublet EE appears once, and the class-weighted inner product gives the integer directly.

Step 3. Against T2:(3,0,−1,−1,1)T_2:(3,0,-1,-1,1) one gets nT2=124(1⋅3⋅5+0+3⋅(−1)⋅1+6⋅(−1)(−1)+6⋅1⋅1)=124(15−3+6+6)=1n_{T_2}=\tfrac{1}{24}(1\cdot3\cdot5+0+3\cdot(-1)\cdot1+6\cdot(-1)(-1)+6\cdot1\cdot1)=\tfrac{1}{24}(15-3+6+6)=1; the same arithmetic against A1,A2,T1A_1,A_2,T_1 returns zero. Why: only EE and T2T_2 survive, so the five-dimensional level splits into a doublet and a triplet with nothing left over.

Step 4. Assemble the result, D(2)↓O=E⊕T2D^{(2)}\big\downarrow_{O}=E\oplus T_2, with dimensions 2+3=52+3=5. In the full group OhO_h (adding inversion) both keep the even label from ℓ=2\ell=2, giving eg⊕t2ge_g\oplus t_{2g}. Why: the barycentre-preserving split of the dd orbitals into the doublet {dz2,dx2−y2}\{d_{z^2},d_{x^2-y^2}\} and the triplet {dxy,dyz,dzx}\{d_{xy},d_{yz},d_{zx}\} is the celebrated 10 Dq10\,Dq ligand-field splitting. *As a check: χ(2)=(5,−1,1,−1,1)\chi^{(2)}=(5,-1,1,-1,1) over the classes of OO decomposes to E⊕T2E\oplus T_2; in the full group OhO_h the labels become ege_g (upper, +6Dq+6Dq) and t2gt_{2g} (lower, −4Dq-4Dq), separated by the measured 10Dq10Dq.*

Limiting cases

How the general formalism specializes at the extremes of symmetry:

CaseResultMeaning
Full rotational symmetryD(ℓ), dim⁡=2ℓ+1D^{(\ell)},\ \dim=2\ell+1Free atom: SO(3) irreps give the familiar (2ℓ+1)(2\ell+1)-fold degeneracy of orbital angular momentum.
Cubic crystal field (OhO_h)D(2)↓Eg⊕T2gD^{(2)}\downarrow E_g\oplus T_{2g}A dd level (ℓ=2\ell=2, five-fold) splits into a doublet EgE_g and a triplet T2gT_{2g} — the ligand-field splitting 10Dq10Dq.
No symmetryG={e}G=\{e\}Only the trivial one-dimensional irrep exists: every level is non-degenerate and every transition is symmetry-allowed.
Centrosymmetric moleculeIR active ⇔\Leftrightarrow Raman inactiveAn inversion centre imposes the mutual exclusion rule — no fundamental is both IR and Raman active.

Worked examples

1 · Vibrations of ammonia (fully worked, C₃ᵥ)

Ammonia NH₃ is a trigonal pyramid with N=4N=4 atoms and point group C3v={E,2C3,3σv}C_{3v}=\{E,2C_3,3\sigma_v\}, ∣G∣=6|G|=6. Find the number, symmetry species, and IR/Raman activity of its normal modes.

Step 1. Count unmoved atoms. EE fixes all four; each C3C_3 (about the principal axis through N) fixes only the nitrogen; each mirror plane σv\sigma_v contains N and one hydrogen, fixing two atoms. So Nunmoved=(4,1,2)N_{\text{unmoved}}=(4,1,2) over (E,2C3,3σv)(E,2C_3,3\sigma_v).

Step 2. Per-atom contributions: EE gives 1+2cos⁡0=31+2\cos0=3; C3C_3 gives 1+2cos⁡(2π/3)=1−1=01+2\cos(2\pi/3)=1-1=0; σv\sigma_v gives −1+2cos⁡0=1-1+2\cos0=1. Multiply: χ3N=(4⋅3, 1⋅0, 2⋅1)=(12,0,2)\chi_{3N}=(4\cdot3,\ 1\cdot0,\ 2\cdot1)=\mathbf{(12,0,2)}.

Step 3. Decompose with the magic formula, class sizes (1,2,3)(1,2,3), against A1(1,1,1)A_1(1,1,1), A2(1,1,−1)A_2(1,1,-1), E(2,−1,0)E(2,-1,0): nA1=16(12+0+3⋅2)=3n_{A_1}=\tfrac16(12+0+3\cdot2)=3, nA2=16(12+0−3⋅2)=1n_{A_2}=\tfrac16(12+0-3\cdot2)=1, nE=16(2⋅12+0+0)=4n_{E}=\tfrac16(2\cdot12+0+0)=4. Thus Γ3N=3A1⊕A2⊕4E\Gamma_{3N}=3A_1\oplus A_2\oplus 4E, and the dimension checks: 3+1+8=123+1+8=12.

Step 4. Remove translations Γtrans=(z)+(x,y)=A1⊕E\Gamma_{\text{trans}}=(z)+(x,y)=A_1\oplus E and rotations Γrot=(Rz)+(Rx,Ry)=A2⊕E\Gamma_{\text{rot}}=(R_z)+(R_x,R_y)=A_2\oplus E. Then Γvib=(3−1)A1⊕(1−1)A2⊕(4−1−1)E=2A1⊕2E\Gamma_{\text{vib}}=(3-1)A_1\oplus(1-1)A_2\oplus(4-1-1)E=\mathbf{2A_1\oplus 2E}, of dimension 2+4=6=3N−62+4=6=3N-6.

Step 5. Read activities: A1A_1 carries zz and z2,x2+y2z^2,x^2+y^2 — IR and Raman active; EE carries (x,y)(x,y) and (x2−y2,xy),(xz,yz)(x^2-y^2,xy),(xz,yz) — IR and Raman active. So ammonia has four distinct fundamental frequencies (2A₁ + 2E), every one of them both infrared and Raman active — the symmetric stretch and the umbrella bend (A1A_1), and the degenerate stretch and degenerate bend (EE), observed near 3337, 950, 3444 and 1627 cm⁻¹.

2 · Mutual exclusion in carbon dioxide (linear, D∞ₕ)

CO₂ is linear and centrosymmetric, N=3N=3. Show that it has four vibrational modes and that no fundamental is both IR and Raman active.

Step 1. A linear molecule has 3N−53N-5 vibrations, so CO₂ has 9−5=49-5=\mathbf{4} modes. Analysing the displacement representation of the D∞hD_{\infty h} group sorts them into a symmetric stretch Σg+\Sigma_g^+, an antisymmetric stretch Σu+\Sigma_u^+, and a doubly degenerate bend Πu\Pi_u — a total of 1+1+2=41+1+2=4, as required.

Step 2. Classify by parity. The subscript gg (gerade, even under inversion) or uu (ungerade, odd) is decisive: quadratic forms are always gg, and the coordinates (x,y,z)(x,y,z) are always uu. Hence Raman activity lives on gg modes and IR activity on uu modes.

Step 3. Apply it. The symmetric stretch Σg+\Sigma_g^+ is Raman active only (it appears among x2+y2,z2x^2+y^2,z^2) and IR silent — the two C=O bonds lengthen in phase, so the dipole never changes. The antisymmetric stretch Σu+\Sigma_u^+ (2349 cm⁻¹) and the bend Πu\Pi_u (667 cm⁻¹) are IR active only. No mode is both: this is the mutual exclusion rule, and observing a band in only one of the two spectra is a direct fingerprint of a centre of symmetry.

Physics of the portrait

Group theory did not walk into physics on its own; it was carried there, against considerable resistance, by a handful of people who saw that symmetry was not decoration but computation. The clearest single figure is Eugene Wigner, whose 1931 book made representation theory a working tool of the quantum physicist.

Origins · 1 minEugene Wigner and how this was discovered
EW
Eugene Wigner
1902–1995 · group theory and symmetry in quantum mechanics

Eugene Paul Wigner, a Hungarian-born chemical engineer turned theoretical physicist, was the person who most decisively made group theory indispensable to physics — at a time when many dismissed it as the "Gruppenpest," the group-theory plague. Drawn into the new quantum mechanics in the late 1920s, Wigner realized that the degeneracies and selection rules of atomic spectra were not accidents of the Schrödinger equation but consequences of rotational symmetry, expressible entirely through the representation theory of the rotation group. His 1927 papers introduced the representation-theoretic analysis of atomic terms, and the resulting Wigner–Eckart theorem — that every matrix element of an irreducible tensor operator factorizes into a geometric Clebsch–Gordan coefficient times a single reduced matrix element — remains the most-used consequence of Schur's lemma in all of spectroscopy. His 1931 monograph Gruppentheorie und ihre Anwendung auf die Quantenmechanik der Atomspektren systematized the field and taught a generation to compute with characters. Wigner also proved that symmetries must be represented by unitary or antiunitary operators, pioneered time-reversal, and laid the groundwork for classifying relativistic particles by the irreducible representations of the Poincaré group. He shared the 1963 Nobel Prize in Physics "through the discovery and application of fundamental symmetry principles" — the citation that made "symmetry principle" a household phrase in physics.

History
1927
— Wigner applies representation theory to the structure of atomic spectra
1928
— Weyl's Gruppentheorie und Quantenmechanik brings the machinery to a broad audience
1929
— Bethe computes crystal-field term splittings in his Termaufspaltung in Kristallen
1931
— Wigner's monograph systematizes group theory for the quantum mechanics of atomic spectra
1963
— Wigner shares the Nobel Prize for fundamental symmetry principles

Going further

Texts. Tinkham, Group Theory and Quantum Mechanics (Dover) — the classic physicist's treatment of characters and selection rules · Hamermesh, Group Theory and Its Application to Physical Problems (Dover) · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed., ch. 16–17 · Stone & Goldbart, Mathematics for Physics (CUP 2009), ch. on group representations and their applications · Hassani, Mathematical Physics, ch. on group representation theory · Cotton, Chemical Applications of Group Theory, 3rd ed. (for the spectroscopy). The key numerical results are: (1) the H₂O displacement character χ3N=(9,−1,3,1)\chi_{3N}=(9,-1,3,1) in the xzxz-plane convention decomposes to 3A1⊕A2⊕3B1⊕2B23A_1\oplus A_2\oplus 3B_1\oplus 2B_2, giving Γvib=2A1⊕B1\Gamma_{\text{vib}}=2A_1\oplus B_1; (2) the NH₃ character (12,0,2)(12,0,2) decomposes to 3A1⊕A2⊕4E3A_1\oplus A_2\oplus 4E and Γvib=2A1⊕2E\Gamma_{\text{vib}}=2A_1\oplus 2E; (3) the SO(3) character χ(2)=(5,−1,1,−1,1)\chi^{(2)}=(5,-1,1,-1,1) over the classes of OO gives D(2)→E⊕T2D^{(2)}\to E\oplus T_2; (4) the Clebsch–Gordan dimensions satisfy (2j1+1)(2j2+1)=∑J(2J+1)(2j_1+1)(2j_2+1)=\sum_J(2J+1). Note the convention subtlety: with the molecule in the yzyz plane the same physics labels the antisymmetric stretch B2B_2 — the mode count and activities are convention-independent, only the Mulliken label flips.

Prerequisites. Group Representations · Reducible and Irreducible Representations · Schur's Lemmas · Characters and Orthogonality Theorems · Character Tables and Decomposition of Representations · Tensor Products and Clebsch–Gordan Decomposition. Leads to. The Rotation Group SO(3) · Lie Algebras and the Exponential Map · Tensor Products and Clebsch–Gordan Decomposition.

Ask the tutor: Tell me a molecule or an ion in a crystal field and I will build its character, decompose it, and read off the degeneracies, normal modes, and IR/Raman selection rules step by step.

Why does a symmetry of the Hamiltonian create degeneracy?

If H^\hat{H} commutes with every U(g)U(g), and ψ\psi is an eigenstate with energy EE, then U(g)ψU(g)\psi is also an eigenstate with the same EE, because H^(U(g)ψ)=U(g)(H^ψ)=E U(g)ψ\hat{H}(U(g)\psi)=U(g)(\hat{H}\psi)=E\,U(g)\psi. The whole set {U(g)ψ}\{U(g)\psi\} therefore lies in one energy level and forms an invariant subspace — a representation of GG. Generically that representation is irreducible (if it were reducible, an extra symmetry would be hiding), so the degeneracy of the level equals the dimension of an irrep. This is why the allowed degeneracies of a system are dictated in advance by its symmetry group: for a point group they are the irrep dimensions listed in the character table, and nothing else can occur without an accident.

What is the "magic formula" really doing?

It is an orthogonal projection. The irreducible characters form an orthonormal basis for class functions under ⟨χ′,χ⟩=1∣G∣∑gχ′(g)∗χ(g)\langle\chi',\chi\rangle=\frac{1}{|G|}\sum_g\chi'(g)^\ast\chi(g). Any reducible character is an integer combination χ=∑μnμχ(μ)\chi=\sum_\mu n_\mu\chi^{(\mu)}, so the inner product with χ(ν)\chi^{(\nu)} picks out its coefficient nνn_\nu by orthogonality — exactly as projecting a vector onto a basis direction returns its component. Summing class-by-class with the class sizes NCN_C as weights is a bookkeeping speed-up, since the character is constant on each class. That every nμn_\mu comes out a non-negative integer is a strong check: a fractional or negative answer means an arithmetic slip.

How do I know which operator symmetry to use for a transition?

Match the operator to the physical process and read its transformation from the character table's basis columns. Electric-dipole (E1) transitions use O^∝r\hat{O}\propto\mathbf{r}, transforming as (x,y,z)(x,y,z) — so a mode or state is IR active when its irrep sits beside xx, yy, or zz. Raman scattering probes the polarizability, quadratic in the field, so its operator transforms as x2,y2,z2,xy,yz,zxx^2,y^2,z^2,xy,yz,zx, and Raman activity is read from the quadratic-form column. Magnetic-dipole transitions use the rotations (Rx,Ry,Rz)(R_x,R_y,R_z), and electric-quadrupole transitions the quadratic forms again. The table lists all of these precisely so you can look up the operator's irrep and feed it into the triple-product selection rule.

Why place water in the xz plane rather than the yz plane?

Purely to match the most common textbook label. The physics is identical either way: H₂O has three modes, two of species A1A_1 and one of the BB type, all IR and Raman active. But the name of that BB mode depends on which in-plane axis you call xx. In the xzxz-plane convention the antisymmetric stretch comes out B1B_1; in the yzyz-plane convention the identical mode is called B2B_2, because B1B_1 and B2B_2 are distinguished only by their behaviour under the two vertical mirrors, and swapping the planes swaps the labels. Mulliken's conventions do not fully fix the choice for C2vC_{2v}, so both appear in the literature. Always state your axis choice; never treat the B1B_1-versus-B2B_2 label as physical.

Can group theory tell me how strong an allowed transition is?

No — and this is its single most important limitation. Group theory is a binary oracle: it certifies only that an amplitude is either forced to zero or not forced to zero. When a transition is allowed, the actual magnitude depends on the overlap of the real wavefunctions and the dynamics, which no character table contains. A symmetry-allowed line can still be extremely weak, and a formally forbidden line can appear once the assumptions relax — through vibronic coupling, higher multipoles, or spin–orbit mixing. The value of the method is in the zeros it guarantees and the degeneracies it enforces; for intensities you hand the surviving matrix elements to perturbation theory or a direct computation.

How does the same framework handle both molecules and atoms?

Through the identical toolkit — characters, orthogonality, decomposition — applied to different groups. For a molecule the group is a finite point group with a handful of tabulated irreps, and you decompose the 3N3N displacement representation. For a free atom the group is the continuous rotation group SO(3), whose irreps are the multiplets D(ℓ)D^{(\ell)} of dimension 2ℓ+12\ell+1, with the character formula χ(ℓ)(θ)=sin⁡[(ℓ+12)θ]/sin⁡(θ/2)\chi^{(\ell)}(\theta)=\sin[(\ell+\tfrac12)\theta]/\sin(\theta/2) replacing a finite table. Crystal-field theory is the bridge: take the atomic D(ℓ)D^{(\ell)}, evaluate its SO(3) character only at the rotation angles present in the crystal's point group, and decompose with the same magic formula — the continuous problem reduced to a finite one.

Before this

Where this leads