Skip to content

Molecular Symmetry

Molecular symmetry turns a geometric observation into exact statements about a specified molecular model. Once the equilibrium geometry, isotopic composition, Hamiltonian, and coordinate convention are fixed, symmetry can:

  • divide electronic, vibrational, and rotational problems into independent blocks;
  • label states and orbitals without relying on a particular visualization;
  • identify degeneracies protected by the model;
  • prove that selected matrix elements vanish;
  • organize normal modes and spectroscopic activity;
  • expose when a numerical calculation has broken an intended symmetry.

The useful question is therefore not merely “what shape is the molecule?” It is:

Which transformations leave the relevant Hamiltonian unchanged, how do the states and operators transform, and what follows from those transformation laws?

A molecular point group is the practical fixed-geometry answer. It consists of rotations, reflections, inversion, and improper rotations that map the nuclear framework onto itself while leaving at least one point fixed. Point groups are powerful, but they are not the whole symmetry of a freely rotating, vibrating, permuting molecule.

This page owns the practical molecular workflow:

geometry↓point group↓representation↓irreducible labels↓physical consequences.\begin{gathered} \text{geometry} \\ \downarrow \\ \text{point group} \\ \downarrow \\ \text{representation} \\ \downarrow \\ \text{irreducible labels} \\ \downarrow \\ \text{physical consequences}. \end{gathered}

It develops point-group identification, character reduction, molecular selection rules, normal-mode classification, and symmetry-adapted orbitals.

Related topics have separate canonical homes:

The emphasis here is not memorizing tables. It is learning to reconstruct and audit the reasoning.

After removing overall translation, describe the nuclear framework by body-fixed positions

R=(R1,…,RN).\mathbf R = \left( \mathbf R_1,\ldots,\mathbf R_N \right).

An orthogonal transformation g∈O(3)g\in O(3) is a point-group operation when it maps this framework onto itself, allowing exchange only among identical nuclear species. Schematically,

GR={g∈O(3):gR=PgR},G_{\mathbf R} = \left\{ g\in O(3): g\mathbf R = P_g\mathbf R \right\},

where PgP_g is a permutation that preserves nuclear species. For an isotopically specified rovibrational problem, operations that exchange different isotopes need not remain symmetries of the nuclear kinetic energy.

The point-group label therefore belongs to a particular geometry and model. It is commonly assigned at an equilibrium structure, a transition structure, or another stationary geometry. A generic displaced geometry usually has less symmetry than the stationary point even when the displacement is a normal coordinate classified by the higher-symmetry reference group.

Schoenflies notation is standard in molecular spectroscopy and quantum chemistry:

  • C1C_1 contains only the identity.
  • CsC_s has one mirror plane, and CiC_i has inversion.
  • CnC_n, CnvC_{nv}, and CnhC_{nh} have a principal nn-fold axis, with the indicated vertical or horizontal mirror planes.
  • DnD_n, DnhD_{nh}, and DndD_{nd} also have nn twofold axes perpendicular to the principal axis.
  • S2nS_{2n} groups are generated by an improper rotation.
  • TdT_d, OhO_h, and IhI_h describe ideal tetrahedral, octahedral, and icosahedral frameworks.
  • C∞vC_{\infty v} and D∞hD_{\infty h} are the linear-molecule groups without and with inversion, respectively.

Ideal examples include:

  • bent water, C2vC_{2v};
  • pyramidal ammonia, C3vC_{3v};
  • trigonal-planar boron trifluoride, D3hD_{3h};
  • methane, TdT_d;
  • benzene, D6hD_{6h};
  • hydrogen chloride, C∞vC_{\infty v};
  • carbon dioxide, D∞hD_{\infty h}.

These labels assume idealized isolated structures. Substitution, isotopic labeling, an external field, a crystal environment, or a distorted geometry can reduce the symmetry.

For a finite equilibrium geometry:

  1. Check whether the molecule is linear.
  2. Check for one of the high-symmetry polyhedral groups.
  3. Find the highest-order proper rotation axis and call it the principal axis.
  4. Look for nn twofold axes perpendicular to that principal axis.
  5. Test for a horizontal mirror plane, vertical or dihedral mirror planes, an inversion center, and improper axes.
  6. Verify closure by listing every operation, not only every visible symmetry element.

The last step matters. A symmetry element is a geometric object such as an axis or plane. A symmetry operation is a transformation such as C3C_3 or C32C_3^2. One threefold axis supplies two nonidentity rotations.

Point groups are stabilizers, not molecular trajectories

Section titled “Point groups are stabilizers, not molecular trajectories”

The point group of Re\mathbf R_e contains operations that leave the reference framework invariant as a set. It does not include arbitrary internal motions that carry the molecule through other geometries. Internal rotation, inversion tunneling, pseudorotation, and exchange of identical nuclei through feasible paths may require a molecular symmetry or permutation–inversion group.

That distinction is essential in high-resolution spectroscopy. A rigid equilibrium point group is often the correct first tool, but it is not always the exact symmetry group of the full rovibronic Hamiltonian.

The standard spatial operations are:

E:identity,Cnk:rotation through 2πk/n,σ:reflection in a plane,i:r↦−r,Snk:rotation followed by reflection.\begin{aligned} E &: \text{identity}, \\ C_n^k &: \text{rotation through } 2\pi k/n, \\ \sigma &: \text{reflection in a plane}, \\ i &: \mathbf r\mapsto-\mathbf r, \\ S_n^k &: \text{rotation followed by reflection}. \end{aligned}

For a chosen principal axis:

  • σh\sigma_h is perpendicular to the axis;
  • σv\sigma_v contains the axis;
  • σd\sigma_d contains the axis and bisects two perpendicular C2C_2 axes.

With the principal axis chosen as zz, one may write

Cn=Rz ⁣(2πn),Sn=σhCn.C_n = R_z\!\left(\frac{2\pi}{n}\right), \qquad S_n = \sigma_h C_n.

The order of factors is harmless here because the rotation about zz commutes with reflection in the plane perpendicular to zz. For more general products, operator order must be kept explicit.

At fixed nuclear geometry, a spatial operation acts on every electronic coordinate. For a spinless coordinate-space wavefunction,

(U^gψ)(r1,…,rNe)=ψ(g−1r1,…,g−1rNe).\begin{aligned} \left( \widehat U_g\psi \right) &\left( \mathbf r_1,\ldots,\mathbf r_{N_e} \right) \\ &= \psi \left( g^{-1}\mathbf r_1,\ldots, g^{-1}\mathbf r_{N_e} \right). \end{aligned}

with the nuclei simultaneously mapped onto symmetry-equivalent positions. If the clamped-nuclei electronic Hamiltonian respects the operation,

[H^e(R),U^g]=0for every g∈GR.\left[ \widehat H_e(\mathbf R), \widehat U_g \right] = 0 \qquad \text{for every }g\in G_{\mathbf R}.

The inverse in the wavefunction argument makes the operators compose in the same order as the geometric transformations.

Take water in the yzyz plane, with the C2C_2 axis along zz. The two vertical planes are then:

  • σv(xz)\sigma_v(xz), which exchanges the two hydrogens;
  • σv′(yz)\sigma_v'(yz), the molecular plane, which fixes all three nuclei.

The C2(z)C_2(z) rotation also exchanges the hydrogens. All four operations

{E, C2(z), σv(xz), σv′(yz)}\left\{ E,\, C_2(z),\, \sigma_v(xz),\, \sigma_v'(yz) \right\}

map the nuclear framework onto itself and form C2vC_{2v}.

Water molecule with the C2 axis and two vertical reflection planes labeled

Axis convention for the worked C2vC_{2v} example. The molecular plane is yzyz; C2(z)C_2(z) and σv(xz)\sigma_v(xz) exchange the hydrogens, while σv′(yz)\sigma_v'(yz) fixes each nucleus.

Axis conventions are not cosmetic. If the molecule is instead placed in the xzxz plane, labels called B1B_1 and B2B_2 in one convention can be exchanged in another. The displacement pattern and measurable polarization rule remain the same.

Ordinary point groups act on spatial coordinates. If spin–orbit coupling is important, spinors transform under a double group, and a 2π2\pi rotation cannot simply be represented as the identity on the spin state. Magnetic fields also require care because axial vectors and antiunitary time reversal enter the symmetry analysis.

An external electric field restricts the group to operations that preserve the field direction and polarity. A calculation performed in a field must therefore use the symmetry of “molecule plus field,” not the field-free molecular point group.

Let Γ\Gamma label an irreducible representation of GG. A basis of one copy of that irrep transforms as

U^g∣Γ,a,μ⟩=∑ν=1dΓDνμ(Γ)(g)∣Γ,a,ν⟩.\widehat U_g \lvert\Gamma,a,\mu\rangle = \sum_{\nu=1}^{d_\Gamma} D_{\nu\mu}^{(\Gamma)}(g) \lvert\Gamma,a,\nu\rangle.

Here:

  • dΓd_\Gamma is the irrep dimension;
  • μ\mu labels components inside one irrep;
  • aa distinguishes repeated copies of the same irrep.

If the Hamiltonian commutes with the group, its matrix is block diagonal in irreducible sectors:

H=⨁Γ(CmΓ⊗VΓ).\mathcal H = \bigoplus_\Gamma \left( \mathbb C^{m_\Gamma} \otimes V_\Gamma \right).

States belonging to inequivalent irreps cannot mix under a symmetry-preserving Hamiltonian. Repeated copies of the same irrep can mix, so a symmetry label alone need not identify a unique state.

Every component of a single irrep has the same energy when the Hamiltonian is group invariant. A multidimensional irrep therefore enforces degeneracy:

dΓ>1⟹at least dΓ componentsin one symmetry-protected multiplet.\begin{gathered} d_\Gamma>1 \quad\Longrightarrow\quad \text{at least }d_\Gamma\text{ components} \\ \text{in one symmetry-protected multiplet}. \end{gathered}

The converse is false. Distinct one-dimensional irreps can cross, and unrelated levels can be accidentally degenerate. Conversely, numerical use of an Abelian subgroup can print the components of a parent multidimensional irrep under different labels even though the exact parent symmetry keeps them degenerate.

For common molecular point groups:

  • AA and BB usually denote one-dimensional irreps;
  • EE denotes a two-dimensional irrep;
  • TT denotes a three-dimensional irrep;
  • subscripts distinguish irreps with different behavior under selected operations;
  • gg and uu denote even and odd parity under inversion;
  • primes and double primes denote even and odd behavior under a designated horizontal reflection.

The precise subscript convention depends on the group and chosen axes. Read the table header before interpreting a label. The symbol EE can denote either the identity operation or a two-dimensional irrep; context distinguishes them.

A chosen set of coordinates, atomic orbitals, or displacement vectors usually transforms reducibly. If

Γred≅⨁αnαΓα,\Gamma_{\mathrm{red}} \cong \bigoplus_\alpha n_\alpha\Gamma_\alpha,

then nαn_\alpha is the number of copies of irrep Γα\Gamma_\alpha. This decomposition is the bridge from geometry to useful basis functions.

The character of a representation matrix is its trace:

χ(Γ)(g)=Tr⁡D(Γ)(g).\chi^{(\Gamma)}(g) = \operatorname{Tr} D^{(\Gamma)}(g).

Characters are invariant under a basis change and are constant on conjugacy classes. A character table therefore lists:

  • conjugacy classes in columns;
  • irreducible representations in rows;
  • irreducible characters in the central block;
  • common coordinate, rotation, and quadratic functions as interpretive aids.

The functions printed beside an irrep are examples of objects with that transformation law. They are not an exhaustive basis and do not mean that the irrep “is” a coordinate.

For finite groups, irreducible characters obey

1∣G∣∑g∈Gχ(α)(g)∗χ(β)(g)=δαβ.\frac{1}{|G|} \sum_{g\in G} \chi^{(\alpha)}(g)^* \chi^{(\beta)}(g) = \delta_{\alpha\beta}.

Taking the inner product of a reducible character with an irreducible character gives the multiplicity:

nα=1∣G∣∑C∣C∣ χ(α)(C)∗χ(red)(C).n_\alpha = \frac{1}{|G|} \sum_C |C|\, \chi^{(\alpha)}(C)^* \chi^{(\mathrm{red})}(C).

The sum is over conjugacy classes CC, and ∣C∣|C| is the number of operations in the class. A valid reduction must give nonnegative integers and satisfy the dimension check

dim⁡Γred=∑αnαdα.\dim\Gamma_{\mathrm{red}} = \sum_\alpha n_\alpha d_\alpha.

Failure of either check usually means that an atom-fixed count, vector trace, class size, or axis convention is wrong.

The projector onto the full isotypic sector of irrep Γ\Gamma is

P^(Γ)=dΓ∣G∣∑g∈Gχ(Γ)(g)∗U^g.\widehat{\mathcal P}^{(\Gamma)} = \frac{d_\Gamma}{|G|} \sum_{g\in G} \chi^{(\Gamma)}(g)^* \widehat U_g.

Applying it to trial orbitals or internal coordinates produces symmetry-adapted combinations. This character projector selects all copies of Γ\Gamma. Resolving individual rows and repeated copies requires the fuller matrix-element projectors or an additional orthogonalization and diagonalization within the projected space.

Products of states or operators transform in tensor-product representations. Their characters multiply:

χ(Γ⊗Λ)(g)=χ(Γ)(g)χ(Λ)(g).\chi^{(\Gamma\otimes\Lambda)}(g) = \chi^{(\Gamma)}(g) \chi^{(\Lambda)}(g).

The resulting character can be reduced with the same formula. For one-dimensional real irreps, direct products are especially simple: multiply the signs operation by operation.

The compact C2vC_{2v} and C3vC_{3v} tables used below are collected at Character Tables.

Each of the three nuclei contributes three Cartesian displacement coordinates, so water begins with a nine-dimensional representation Γ3N\Gamma_{3N}. Its character under operation gg is

χ3N(g)=∑A: gA=ATr⁡Rg.\chi_{3N}(g) = \sum_{A:\,gA=A} \operatorname{Tr}R_g.

Only nuclei left at their own positions contribute to the trace. A nucleus moved to another site contributes zero to the diagonal trace, even when the destination contains an identical nucleus.

For the axis convention in the figure:

χ3N(E)=3(3)=9,χ3N ⁣(C2)=1(−1)=−1,χ3N ⁣(σv(xz))=1(1)=1,χ3N ⁣(σv′(yz))=3(1)=3.\begin{aligned} \chi_{3N}(E) &=3(3)=9, \\ \chi_{3N}\!\left(C_2\right) &=1(-1)=-1, \\ \chi_{3N}\!\left(\sigma_v(xz)\right) &=1(1)=1, \\ \chi_{3N}\!\left(\sigma_v'(yz)\right) &=3(1)=3. \end{aligned}

The factors 33, −1-1, and 11 are traces of the corresponding three-dimensional vector transformations. Reducing the character

χ3N=(9,−1,1,3)\chi_{3N} = \left( 9,-1,1,3 \right)

against the C2vC_{2v} table gives

Γ3N=3A1⊕A2⊕2B1⊕3B2.\Gamma_{3N} = 3A_1 \oplus A_2 \oplus 2B_1 \oplus 3B_2.

The dimensions add to 99, as required.

With zz along the C2C_2 axis and the molecule in the yzyz plane,

Γtrans=Γ(x,y,z)=A1⊕B1⊕B2,Γrot=Γ(Rx,Ry,Rz)=B2⊕B1⊕A2.\begin{aligned} \Gamma_{\mathrm{trans}} &= \Gamma(x,y,z) = A_1\oplus B_1\oplus B_2, \\ \Gamma_{\mathrm{rot}} &= \Gamma(R_x,R_y,R_z) = B_2\oplus B_1\oplus A_2. \end{aligned}

Therefore

Γvib=Γ3N−Γtrans−Γrot=2A1⊕B2.\begin{aligned} \Gamma_{\mathrm{vib}} &= \Gamma_{3N} - \Gamma_{\mathrm{trans}} - \Gamma_{\mathrm{rot}} \\ &= 2A_1 \oplus B_2. \end{aligned}

The two A1A_1 modes are conventionally the symmetric stretch and bend; the B2B_2 mode is the asymmetric stretch in this axis convention. All irreps are one dimensional, so C2vC_{2v} does not enforce a vibrational degeneracy.

The subtraction is subtraction of irrep multiplicities after decomposition. It is not subtraction of sets or individual displacement arrows.

For a transition operator O^\widehat O, consider

Mfi=⟨ψf∣O^∣ψi⟩.M_{fi} = \langle\psi_f| \widehat O |\psi_i\rangle.

The matrix element can be nonzero only if

Γf∗⊗ΓO⊗Γi⊃Γts,\Gamma_f^* \otimes \Gamma_O \otimes \Gamma_i \supset \Gamma_{\mathrm{ts}},

where Γts\Gamma_{\mathrm{ts}} is the totally symmetric irrep. Equivalently, the multiplicity of the totally symmetric irrep is

nts=1∣G∣∑g∈Gχf(g)∗χO(g)χi(g).n_{\mathrm{ts}} = \frac{1}{|G|} \sum_{g\in G} \chi_f(g)^* \chi_O(g) \chi_i(g).

If nts=0n_{\mathrm{ts}}=0, symmetry forces the matrix element to vanish within the stated model. If nts>0n_{\mathrm{ts}}>0, the transition is symmetry allowed, but its intensity can still be extremely small.

The electric-dipole components transform like xx, yy, and zz. In the water convention:

Γ(x)=B1,Γ(y)=B2,Γ(z)=A1.\begin{aligned} \Gamma(x)&=B_1, & \Gamma(y)&=B_2, \\ \Gamma(z)&=A_1. \end{aligned}

From an A1A_1 initial state, an xx-, yy-, or zz-polarized transition can reach B1B_1, B2B_2, or A1A_1, respectively, provided the other quantum numbers and approximations also permit the transition.

For a fundamental vibration from a totally symmetric vibrational ground state, the normal coordinate must transform like at least one dipole component. Water’s 2A1+B22A_1+B_2 fundamentals are therefore all infrared active in the electric-dipole and harmonic approximations.

The leading nonresonant Raman operator is the symmetric polarizability tensor. Its components transform like

x2, y2, z2, xy, xz, yz.x^2,\ y^2,\ z^2,\ xy,\ xz,\ yz.

A normal mode is first-order Raman active when its irrep occurs among those quadratic functions. This is a different condition from infrared activity.

For a centrosymmetric molecule, dipole components are ungerade and polarizability components are gerade. A harmonic fundamental from a gerade ground state therefore cannot be both electric-dipole infrared active and first-order Raman active. This mutual-exclusion rule depends on inversion symmetry and on the stated leading operators.

A nominally forbidden line can acquire intensity through:

  • vibronic or spin–orbit mixing;
  • Coriolis or anharmonic coupling;
  • isotopic or environmental symmetry breaking;
  • external fields;
  • magnetic-dipole or electric-quadrupole operators;
  • higher-order terms in the dipole or polarizability expansion;
  • intensity borrowing from a nearby allowed state.

The correct statement is always conditional: “forbidden by this symmetry and operator at this order.”

At a symmetry-preserving stationary geometry, the mass-weighted Hessian F\mathbf F commutes with the displacement representation:

D(g)†FD(g)=F.\mathbf D(g)^\dagger \mathbf F \mathbf D(g) = \mathbf F.

It can therefore be block diagonalized by irreducible representation. This gives three immediate consequences:

  1. modes from inequivalent irreps do not mix while the symmetry is preserved;
  2. multidimensional irreps enforce degenerate harmonic frequencies;
  3. transition-property derivatives can be screened by symmetry before an intensity calculation.

Modes of the same irrep can mix freely. If several modes are labeled A1A_1, the label does not uniquely mean “the symmetric stretch.” Frequencies, displacement vectors, and internal-coordinate projections are needed for a physical assignment.

For ammonia in C3vC_{3v}, an EE vibrational mode has two components. Any orthonormal basis spanning that two-dimensional eigenspace is valid. Numerical eigenvectors can rotate within the subspace as the method, basis, or geometry changes while the observable subspace and common harmonic frequency remain unchanged.

A small computed splitting of an expected EE pair can signal:

  • a geometry that is not exactly symmetry adapted;
  • loose optimization or integration thresholds;
  • use of a lower-symmetry computational subgroup;
  • a genuinely symmetry-breaking perturbation;
  • an electronically unstable high-symmetry reference.

The point-group label alone does not decide which explanation applies.

When a molecule distorts from group GG to subgroup HH, an irrep of GG may remain irreducible or split when restricted to HH:

Γ(G)↓H=⨁αnαΓα(H).\Gamma^{(G)} \downarrow H = \bigoplus_\alpha n_\alpha \Gamma_\alpha^{(H)}.

Correlation tables record this restriction. They are indispensable for tracking orbital and vibrational labels along reaction paths, symmetry-breaking coordinates, and Jahn–Teller distortions. Equal subgroup labels permit mixing; they do not guarantee that mixing is large.

The detailed normal-coordinate construction and intensity derivatives are developed in Normal Modes of Polyatomics. Raman Spectroscopy specializes the quadratic-function test to measured Raman tensors, depolarization ratios, and oriented samples.

Suppose the two hydrogen 1s1s functions in water are ∣HL⟩|H_L\rangle and ∣HR⟩|H_R\rangle. Their two-dimensional permutation representation has character

χH=(2,0,0,2),\chi_H = \left( 2,0,0,2 \right),

because C2(z)C_2(z) and σv(xz)\sigma_v(xz) exchange the functions, whereas EE and σv′(yz)\sigma_v'(yz) fix them. Character reduction gives

ΓH=A1⊕B2.\Gamma_H = A_1 \oplus B_2.

Normalized symmetry-adapted linear combinations are

∣HA1⟩=∣HL⟩+∣HR⟩2,∣HB2⟩=∣HL⟩−∣HR⟩2.\begin{aligned} |H_{A_1}\rangle &= \frac{ |H_L\rangle+|H_R\rangle }{\sqrt2}, \\ |H_{B_2}\rangle &= \frac{ |H_L\rangle-|H_R\rangle }{\sqrt2}. \end{aligned}

On oxygen, the 2s2s and 2pz2p_z orbitals transform as A1A_1, 2py2p_y as B2B_2, and 2px2p_x as B1B_1 in this convention. A totally symmetric one-electron Hamiltonian couples only matching symmetry sectors:

⟨ϕα∣h^∣ϕβ⟩=0if Γα≇Γβ.\langle\phi_\alpha| \widehat h |\phi_\beta\rangle = 0 \qquad \text{if } \Gamma_\alpha\not\cong\Gamma_\beta.

Thus the A1A_1 hydrogen combination can mix with oxygen A1A_1 functions, and the B2B_2 combination can mix with oxygen 2py2p_y. The B1B_1 oxygen 2px2p_x function has no partner in this minimal hydrogen σ\sigma basis.

Symmetry proves the zeros. It does not determine the nonzero coupling strengths, orbital ordering, occupation, or degree of bonding.

A molecular orbital has an irrep label, but a many-electron state is classified by the direct product of all occupied spin-orbital transformation laws, reduced with the required fermionic antisymmetry and spin coupling. The symmetry of a single frontier orbital is not automatically the symmetry of an excited state.

For open-shell systems or strong spin–orbit coupling, one may need:

  • spin-adapted configuration state functions;
  • double-group irreps;
  • total angular-momentum projections for linear molecules;
  • explicit treatment of near-degenerate electronic configurations.

Quantum-chemistry programs often exploit only an Abelian subgroup even when the geometry has a non-Abelian point group. This keeps all working irreps one-dimensional and simplifies integral storage. Components of a parent EE or TT irrep then receive separate subgroup labels.

When reporting a calculation, record:

  • the geometry and isotopologue;
  • the axis convention;
  • the full intended point group;
  • the subgroup actually used by the program;
  • whether symmetry was detected, imposed, or disabled;
  • how near-degenerate states or orbitals were matched across calculations.

Electronic Coulomb interactions depend on nuclear charges, not nuclear masses. Replacing one isotope can therefore leave the clamped-nuclei electronic potential symmetric at the same geometry while reducing the symmetry of the nuclear kinetic energy and the full rovibrational problem.

This is why an isotopic substitution can activate weak vibrational lines or split rovibrational structure without substantially changing the electronic potential-energy surface.

For tunneling between equivalent minima, internal rotation, umbrella inversion, or pseudorotation, the experimentally relevant states can extend over several geometries. A single equilibrium point group cannot represent every feasible permutation and inversion operation. Molecular symmetry groups and permutation–inversion groups provide the appropriate global classification.

A nonlinear molecule in a degenerate electronic state can be unstable to a symmetry-lowering distortion through the Jahn–Teller effect. The high-symmetry point group still organizes the parent electronic and vibrational spaces, but the distorted minima belong to subgroups.

At intersections of electronic potential-energy surfaces, point-group symmetry can permit or forbid couplings on special subspaces, while generic conical intersection topology depends on the local branching plane. See Conical Intersections for the canonical degeneracy geometry.

Near symmetry is often chemically informative, but “almost an irrep” is not an exact group-theory concept. If a perturbation is small, write

H^=H^0+λV^,\widehat H = \widehat H_0 + \lambda\widehat V,

where H^0\widehat H_0 has the higher symmetry. The labels of H^0\widehat H_0 then remain useful zeroth-order labels, and V^\widehat V controls mixing and splitting. Quantifying λV^\lambda\widehat V is more informative than assigning a higher point group by visual tolerance alone.

From coordinates to a defensible assignment

Section titled “From coordinates to a defensible assignment”
  1. Specify the object. State the geometry, isotopologue, charge, electronic state, and external environment.
  2. Remove arbitrary orientation. Choose and report body-fixed axes, usually guided by principal inertial axes and the highest-order symmetry axis.
  3. Find candidate operations. Apply each operation to the Cartesian coordinates and match atoms only within the same species or isotope class relevant to the model.
  4. Use a declared tolerance. Coordinate noise is not a symmetry principle. Compare assignments across tighter optimization and matching tolerances.
  5. Verify the group. Check products, inverses, and the full operation count.
  6. Construct the representation. Decide whether the basis consists of Cartesian displacements, orbitals, internal coordinates, or many-electron states.
  7. Compute characters and reduce. Check integer multiplicities and total dimension.
  8. Derive consequences. Block diagonalization, degeneracies, direct products, and selection rules follow only after the representations are fixed.
  9. Validate numerically. Inspect forbidden matrix elements, expected degeneracies, and stability under tighter thresholds or explicit symmetrization.

An automatic point-group detector solves a numerical matching problem. Its answer can change when:

  • an optimized coordinate differs by roundoff;
  • a shallow mode leaves the structure slightly distorted;
  • the tolerance merges genuinely distinct positions;
  • atom ordering or isotopic labels are inconsistent;
  • a high-symmetry stationary point is unstable.

Use the detector as evidence, not as an oracle. Report the tolerance or program default when the assignment affects degeneracies, state labels, or thermochemical symmetry numbers.

Bent, pyramidal, tetrahedral, and planar are geometric descriptions. A point group is the complete set of symmetry operations for a specified structure.

“Every visible axis gives one operation”

Section titled ““Every visible axis gives one operation””

An nn-fold axis contributes Cn,Cn2,…,Cnn−1C_n,C_n^2,\ldots,C_n^{n-1}, and some powers may belong to different conjugacy classes.

A character is the trace of a representation matrix. It is a sum of eigenvalues and can have magnitude larger than one.

“Different irreps can never have the same energy”

Section titled ““Different irreps can never have the same energy””

Symmetry prevents their mixing but does not prevent an accidental or parameter-tuned degeneracy.

Equal symmetry makes mixing permissible. The coupling matrix element and energy separation determine whether it is appreciable.

Symmetry only determines whether the leading matrix element is forced to zero. Radial overlap, property derivatives, populations, and dynamical factors set the intensity.

The line may become weakly allowed when the Hamiltonian, state mixing, or transition operator goes beyond the approximation used to derive the rule.

“The software label is coordinate independent”

Section titled ““The software label is coordinate independent””

Physical subspaces are coordinate independent; names such as B1B_1 and B2B_2 depend on axis conventions and subgroup choices.

“The equilibrium point group is the exact molecular symmetry”

Section titled ““The equilibrium point group is the exact molecular symmetry””

Overall rotation, nuclear permutation, large-amplitude motion, spin, and external fields can require a larger or smaller symmetry framework.

Before using a molecular symmetry label, ask:

  • What geometry and isotopologue does the label describe?
  • Which axis convention is being used?
  • Is the group exact for the Hamiltonian or only a useful reference symmetry?
  • What basis carries the representation?
  • Were characters calculated from fixed basis functions correctly?
  • Do reduction multiplicities and dimensions check?
  • Is a claimed degeneracy protected by a multidimensional irrep?
  • Does a selection rule specify the operator and approximation?
  • Was full symmetry or only a computational subgroup used?
  • Could a distortion, field, environment, or nonadiabatic coupling relax the conclusion?
  • A molecular point group is attached to a specified nuclear geometry and Hamiltonian.
  • Symmetry operations act on coordinates; representations describe their action on states, orbitals, or displacement spaces.
  • Characters provide a basis-independent route from reducible coordinates to irreducible sectors.
  • A symmetry-preserving Hamiltonian is block diagonal in inequivalent irreps.
  • Multidimensional irreps enforce degeneracy, but symmetry does not explain every degeneracy.
  • Selection rules are statements about a state–operator–state direct product, not universal declarations that a line is absent.
  • Normal modes and orbitals can be built from the same projection and reduction machinery.
  • Point groups are a controlled first approximation to molecular symmetry, not a substitute for permutation–inversion groups, double groups, or dynamical analysis when those are needed.

Exercise 1: Identify idealized point groups

Section titled “Exercise 1: Identify idealized point groups”

Assign the equilibrium point group of idealized HCl\mathrm{HCl}, CO2\mathrm{CO_2}, H2O\mathrm{H_2O}, NH3\mathrm{NH_3}, and CH4\mathrm{CH_4}. State the geometric feature that distinguishes each answer.

Solution
  • HCl\mathrm{HCl} is linear and heteronuclear, so it is C∞vC_{\infty v}.
  • CO2\mathrm{CO_2} is linear and centrosymmetric, so it is D∞hD_{\infty h}.
  • Bent H2O\mathrm{H_2O} has one C2C_2 axis and two vertical mirror planes, so it is C2vC_{2v}.
  • Pyramidal NH3\mathrm{NH_3} has one C3C_3 axis and three vertical mirror planes, so it is C3vC_{3v}.
  • Ideal tetrahedral CH4\mathrm{CH_4} has tetrahedral symmetry, TdT_d.

These are equilibrium-geometry labels. Isotopic substitution, distortion, or an external environment can reduce them.

Exercise 2: Reduce the water displacement representation

Section titled “Exercise 2: Reduce the water displacement representation”

Using the C2vC_{2v} character

χ3N=(9,−1,1,3),\chi_{3N} = \left( 9,-1,1,3 \right),

show that

Γ3N=3A1⊕A2⊕2B1⊕3B2.\Gamma_{3N} = 3A_1 \oplus A_2 \oplus 2B_1 \oplus 3B_2.

Then remove translations and rotations.

Solution

Each C2vC_{2v} class contains one operation. Taking the character inner product with each row of the canonical table gives

nA1=14(9−1+1+3)=3,nA2=14(9−1−1−3)=1,nB1=14(9+1+1−3)=2,nB2=14(9+1−1+3)=3.\begin{aligned} n_{A_1} &= \frac14 \left( 9-1+1+3 \right) =3, \\ n_{A_2} &= \frac14 \left( 9-1-1-3 \right) =1, \\ n_{B_1} &= \frac14 \left( 9+1+1-3 \right) =2, \\ n_{B_2} &= \frac14 \left( 9+1-1+3 \right) =3. \end{aligned}

With

Γtrans=A1⊕B1⊕B2,Γrot=A2⊕B1⊕B2,\begin{aligned} \Gamma_{\mathrm{trans}} &= A_1\oplus B_1\oplus B_2, \\ \Gamma_{\mathrm{rot}} &= A_2\oplus B_1\oplus B_2, \end{aligned}

the vibrational representation is

Γvib=2A1⊕B2.\Gamma_{\mathrm{vib}} = 2A_1\oplus B_2.

Its dimension is three, equal to 3N−63N-6 for nonlinear water.

Apply symmetry to the two water hydrogen 1s1s functions. Show that their representation is A1⊕B2A_1\oplus B_2 and construct normalized symmetry-adapted linear combinations.

Solution

The identity and molecular-plane reflection fix both functions, while C2(z)C_2(z) and σv(xz)\sigma_v(xz) exchange them. Thus

χH=(2,0,0,2),\chi_H = \left( 2,0,0,2 \right),

which reduces to A1⊕B2A_1\oplus B_2. The normalized combinations are

∣HA1⟩=∣HL⟩+∣HR⟩2,∣HB2⟩=∣HL⟩−∣HR⟩2.\begin{aligned} |H_{A_1}\rangle &= \frac{ |H_L\rangle+|H_R\rangle }{\sqrt2}, \\ |H_{B_2}\rangle &= \frac{ |H_L\rangle-|H_R\rangle }{\sqrt2}. \end{aligned}

Their relative signs describe transformation behavior. An overall sign change of either SALC has no physical effect.

Exercise 4: Polarization selection in water

Section titled “Exercise 4: Polarization selection in water”

Assume an A1A_1 initial state in C2vC_{2v}. Which final-state irreps can be reached by xx-, yy-, and zz-polarized electric-dipole transitions in the axis convention used on this page?

Solution

The coordinate components transform as

Γ(x)=B1,Γ(y)=B2,Γ(z)=A1.\begin{aligned} \Gamma(x)&=B_1, & \Gamma(y)&=B_2, \\ \Gamma(z)&=A_1. \end{aligned}

For one-dimensional real irreps, the product

Γf⊗Γ(μa)⊗A1\Gamma_f \otimes \Gamma(\mu_a) \otimes A_1

contains A1A_1 precisely when Γf=Γ(μa)\Gamma_f=\Gamma(\mu_a). Therefore

x:B1,y:B2,z:A1.x:B_1, \qquad y:B_2, \qquad z:A_1.

Other quantum numbers and the actual transition moment can impose additional restrictions or make an allowed line weak.

Exercise 5: Prove vibrational mutual exclusion

Section titled “Exercise 5: Prove vibrational mutual exclusion”

For a centrosymmetric molecule with a gerade vibrational ground state, show that a one-quantum harmonic fundamental cannot be both electric-dipole infrared active and first-order Raman active.

Solution

The electric dipole is ungerade. A nonzero infrared matrix element from a gerade ground state therefore requires an ungerade final vibrational state:

g⊗u⊗u=g.g\otimes u\otimes u = g.

The polarizability tensor is gerade. A nonzero first-order Raman matrix element from the same ground state requires a gerade final state:

g⊗g⊗g=g.g\otimes g\otimes g = g.

A normal mode of a centrosymmetric molecule is either gerade or ungerade, not both. Hence the same fundamental cannot satisfy both conditions within these approximations. Higher-order mechanisms or broken inversion symmetry can relax the observed rule.

Let ∣ϕα⟩|\phi_\alpha\rangle and ∣ϕβ⟩|\phi_\beta\rangle belong to inequivalent irreps of a finite point group, and let h^\widehat h be totally symmetric. Why must

⟨ϕα∣h^∣ϕβ⟩=0?\langle\phi_\alpha| \widehat h |\phi_\beta\rangle =0?
Solution

The matrix element transforms in

Γα∗⊗Γts⊗Γβ≅Γα∗⊗Γβ.\Gamma_\alpha^* \otimes \Gamma_{\mathrm{ts}} \otimes \Gamma_\beta \cong \Gamma_\alpha^* \otimes \Gamma_\beta.

By irreducible-character orthogonality, the totally symmetric irrep occurs in Γα∗⊗Γβ\Gamma_\alpha^*\otimes\Gamma_\beta only when the two irreps are equivalent. For inequivalent irreps, group averaging projects the matrix element to zero. This is the finite-group form of the block-diagonal statement behind Schur’s lemma.

Move the water molecule from the yzyz plane to the xzxz plane while keeping zz as the C2C_2 axis. What happens to the B1B_1 and B2B_2 labels? Does the infrared polarization physics change?

Solution

Interchanging the roles of xx and yy interchanges the two vertical reflection planes. Under the standard C2vC_{2v} convention, labels attached to xx and yy therefore exchange:

B1⟷B2.B_1 \longleftrightarrow B_2.

The asymmetric stretch receives the other printed label. Its displacement pattern and the statement that it couples to dipole polarization in the molecular plane perpendicular to the C2C_2 axis do not change. The relabeling is conventional; the observable relation is geometric.

A calculation on idealized C3vC_{3v} ammonia is run without imposed symmetry. Two modes expected to form an EE pair differ by a small frequency. Give a diagnostic sequence that distinguishes numerical symmetry breaking from a physical instability.

Solution

First inspect the optimized Cartesian geometry and quantify deviations from C3vC_{3v} under a stated atom-matching tolerance. Then:

  1. tighten geometry, integration-grid, self-consistent-field, and Hessian thresholds;
  2. symmetrize the geometry and recompute the analytic or finite-difference Hessian;
  3. project the two eigenvectors onto the expected EE displacement subspace;
  4. test whether the splitting shrinks with numerical convergence;
  5. inspect the electronic and vibrational stability of the high-symmetry stationary point;
  6. repeat with symmetry explicitly imposed and with it disabled.

A splitting that vanishes under convergence is numerical. A robust splitting accompanied by a lower-energy distorted stationary point indicates genuine symmetry lowering. The point-group label alone cannot settle the diagnosis.

  • C. J. H. Schutte, J. E. Bertie, P. R. Bunker, J. T. Hougen, I. M. Mills, J. K. G. Watson, and B. P. Winnewisser, “Notations and Conventions in Molecular Spectroscopy: Part 2. Symmetry Notation,” Pure and Applied Chemistry 69, 1641–1650 (1997), doi:10.1351/pac199769081641.
  • P. R. Bunker and P. Jensen, Molecular Symmetry and Spectroscopy, 2nd ed., NRC Research Press, 1998, ISBN 978-0-660-17519-5.
  • P. R. Bunker and P. Jensen, Fundamentals of Molecular Symmetry, IOP Publishing, 2005, doi:10.1201/9781315273334.
  • F. A. Cotton, Chemical Applications of Group Theory, 3rd ed., Wiley, 1990.
  • D. M. Bishop, Group Theory and Chemistry, Dover, 1993.
  • M. Tinkham, Group Theory and Quantum Mechanics, Dover, 2003.
  • E. B. Wilson Jr., J. C. Decius, and P. C. Cross, Molecular Vibrations: The Theory of Infrared and Raman Vibrational Spectra, McGraw–Hill, 1955.
  • G. Herzberg, Molecular Spectra and Molecular Structure II: Infrared and Raman Spectra of Polyatomic Molecules, Van Nostrand, 1945.
  • P. W. Atkins and R. S. Friedman, Molecular Quantum Mechanics, 5th ed., Oxford University Press, 2011.
  • International Union of Pure and Applied Chemistry, “Point Group,” Compendium of Chemical Terminology, 5th ed., doi:10.1351/goldbook.P04703.
  • NIST Computational Chemistry Comparison and Benchmark Database, species organized by molecular point group, Standard Reference Database 101.