Rotational Transitions and the Rigid Rotor

Point a microwave source at a flask of gaseous HCl. Sweep its frequency slowly across the 20–1000 GHz range. At first nothing happens — the gas is transparent. Then, at very specific, evenly spaced frequencies, the transmitted power drops sharply. One narrow notch appears after another, like the teeth of a comb. Those notches are not noise or apparatus artifacts. They are the fingerprints of a single number, the rotational constant BB. They also reflect a single quantum rule, ΔJ=±1\Delta J = \pm 1. Decode the comb and you recover the bond length of HCl. You get five significant figures without ever touching the molecule with anything but light. This matters because rotational spectroscopy is one of the few places in chemistry where a spectrum can be inverted exactly. There are no fitting parameters and no empirical calibration curves. There is just a rigid-rotor Hamiltonian, an eigenvalue formula, and a transition-moment integral. The technique underlies the determination of interstellar molecule identities. Astrochemists find polyatomic species in cold molecular clouds by matching rotational line lists. It also supports industrial gas analysis and precision tests of fundamental physics using microwave and millimeter-wave sources. It is also the cleanest classroom example of how a symmetry-based selection rule combines with a golden-rule transition rate. Together they produce an experimentally testable prediction. The big picture is this: a diatomic molecule tumbles freely in space, unconstrained by any other degree of freedom. It is exactly solvable. The same angular-momentum machinery from the hydrogen atom's angular equation returns unchanged. What is new here is the coupling to radiation. An oscillating electric field only pushes on a molecule if the molecule already has an electric dipole. Even then, symmetry restricts which pairs of rotational states can talk to each other. This article builds the whole chain: rigid-rotor eigenstates and energies, the dipole transition-moment integral evaluated using Selection Rules from Symmetry, the resulting ΔJ=±1\Delta J=\pm 1 rule, and the observable spectrum. This spectrum is a stack of lines at 2B,4B,6B,…2B, 4B, 6B,\dots (in wavenumbers) whose spacing directly reports the moment of inertia. By the end you will be able to derive the rigid-rotor selection rule from the transition dipole integral, rather than quote it. You will be able to predict which molecules are microwave-silent. You will extract BB and bond length from a line spacing. You will also see exactly where this framework must be patched once vibration is switched on. That is the subject of Vibration-Rotation Coupling and Vibration-Rotation Spectra: P, Q, R Branches.

Key takeaways
  • A diatomic modeled as a rigid rotor has energies EJ=BJ(J+1)E_J = BJ(J+1), with degeneracy 2J+12J+1, independent of vibrational amplitude.
  • Only molecules with a permanent electric dipole moment (μ0≠0\mu_0 \neq 0) show a pure rotational (microwave) spectrum — homonuclear diatomics like N₂ and O₂ are rotationally silent in absorption.
  • The dipole transition-moment integral, evaluated over spherical harmonics, is nonzero only when ΔJ=±1\Delta J = \pm 1 (and ΔMJ=0,±1\Delta M_J = 0,\pm1 for polarized/oriented cases).
  • Allowed absorption lines appear at ν~=2B(J+1)\tilde\nu = 2B(J+1) for J→J+1J\to J+1, giving an evenly spaced comb with spacing 2B2B.
  • The rotational constant B=h/(8π2cI)B = h/(8\pi^2cI) converts a measured line spacing directly into the moment of inertia II and bond length rer_e.
  • Key idea: rotational absorption requires both an energetically allowed transition and a nonzero transition dipole moment, and for a rigid rotor those two conditions together collapse to the single rule ΔJ=±1\Delta J = \pm 1 with μ0≠0\mu_0 \neq 0.

Where this lives today

Radio astronomy leans almost entirely on rotational spectroscopy: molecules in cold interstellar clouds (typically 10–100 K) populate only the lowest few rotational levels, and their rotational transitions fall in the radio/millimeter window that penetrates dust and reaches Earth-based and space telescopes. The identification of interstellar CO, HCN, and dozens of more exotic species (including complex prebiotic molecules) relies on matching observed line frequencies to laboratory-measured rotational constants using exactly the EJ=BJ(J+1)E_J = BJ(J+1), ΔJ=±1\Delta J=\pm1 framework developed here. Atmospheric and environmental monitoring instruments use microwave and millimeter-wave absorption to detect trace polar gases (ClO, HCl, and other polar pollutants or ozone-depletion-relevant species) in the stratosphere, since rotational lines are sharp, essentially collision-free at low pressure, and highly species-specific. Industrial gas analysis and isotope ratio measurements exploit the mass-dependence of BB (through μ=m1m2/(m1+m2)\mu = m_1m_2/(m_1+m_2)): different isotopologues of the same molecule (e.g., H³⁵Cl versus H³⁷Cl) have measurably different rotational constants, allowing non-destructive isotopic composition analysis by microwave spectroscopy. Precision fundamental-physics tests, including searches for a permanent electron electric dipole moment and tests of parity violation in molecules, use ultra-cold polar diatomic and triatomic molecules prepared in single, well-defined rotational states — a direct descendant of the basic rigid-rotor quantum number bookkeeping developed here. Finally, this entire framework is the necessary launching point for real molecular spectra, where rotation and vibration are coupled: Vibration-Rotation Coupling explains how the rotational constant BB becomes vibrational-state-dependent (BvB_v), and Vibration-Rotation Spectra: P, Q, R Branches shows how the same ΔJ=±1\Delta J=\pm1 rule, applied on top of a vibrational transition, produces the characteristic P- and R-branch line patterns seen in infrared gas-phase spectra.

Ask the tutor: How does adding vibration change the rotational constant B, and why does that make the P and R branch lines in an infrared spectrum unevenly spaced instead of a perfect 2B comb?

The physical picture

Picture a diatomic molecule as two point masses, m1m_1 and m2m_2, held at a fixed separation rer_e by a stiff, unstretchable rod — the rigid rotor approximation. In the center-of-mass frame this two-body problem reduces to a single particle of reduced mass μ=m1m2/(m1+m2)\mu = m_1 m_2/(m_1+m_2) moving on the surface of a sphere of radius rer_e. The Schrödinger equation for this system is identical in form to the angular part of the hydrogen-atom problem: the eigenfunctions are the spherical harmonics YJM(θ,ϕ)Y_J^{M}(\theta,\phi), labeled by the rotational quantum number J=0,1,2,…J = 0, 1, 2, \dots and the projection quantum number M=−J,…,+JM = -J, \dots, +J. Each energy level EJE_J is (2J+1)(2J+1)-fold degenerate, reflecting the 2J+12J+1 possible orientations of the angular momentum vector in space for a given magnitude. Crucially, in the absence of an external field, all these MM states are degenerate — the molecule has no preferred axis in the lab frame until something (an electric field, or the oscillating field of the radiation itself) picks one out. Now add a permanent electric dipole moment μ⃗\vec\mu that is fixed along the molecular axis (pointing, say, from the negative to positive end of a polar bond like H–Cl). As the molecule rotates, the direction of μ⃗\vec\mu in the lab frame sweeps around — this rotating dipole is precisely what can couple to an oscillating electric field of electromagnetic radiation. If there is no permanent dipole (as in N₂, O₂, or CO₂ in its ground vibrational state, by symmetry), the molecule presents no rotating charge asymmetry for the field to grab onto, and pure rotational absorption is forbidden regardless of energy matching. The transition intensity is governed by Fermi's Golden Rule and Transition Rates: the rate is proportional to ∣⟨f∣μ^∣i⟩∣2|\langle f|\hat\mu|i\rangle|^2, the squared matrix element of the dipole operator between initial and final rotational states. Because μ^\hat\mu transforms as a vector (rank-1 spherical tensor) under rotations, and the initial/final states are spherical harmonics of definite JJ, the integral is governed by the same angular-momentum coupling rules that appear throughout atomic and molecular spectroscopy — this is the deep reason ΔJ=±1\Delta J = \pm 1 emerges, rather than being an ad hoc rule.

Figure
The rigid-rotor energy ladder: nonuniform spacing but uniform transition spacing
Energy levels E_J = B J(J+1) J=0 J=1 J=2 J=3 0 2B 6B 12B 2B (J:0→1) 4B (J:1→2) 6B (J:2→3) Absorption lines land at 2B, 4B, 6B … evenly spaced by 2B
Figure
Rotating permanent dipole couples to the oscillating field
Why μ₀ ≠ 0 is required δ+ H δ− Cl HCl: μ₀ ≈ 1.08 D — rotating dipole, IR/MW active N≡N: μ₀ = 0 — no net rotating dipole, microwave silent

The core ideas

Rigid-rotor energies. Solving the angular Schrödinger equation for a rotating diatomic of moment of inertia I=μre2I = \mu r_e^2 gives

EJ=ℏ22I J(J+1),J=0,1,2,…E_J = \frac{\hbar^2}{2I}\,J(J+1), \qquad J = 0,1,2,\dots

where:

  • ℏ\hbar is the reduced Planck constant,
  • I=μre2I = \mu r_e^2 is the moment of inertia, with μ\mu the reduced mass and rer_e the equilibrium bond length,
  • JJ is the rotational angular momentum quantum number, with degeneracy gJ=2J+1g_J = 2J+1.

Reading: levels are not evenly spaced in energy — the gap EJ+1−EJE_{J+1}-E_J grows linearly with JJ — yet, as the derivation below shows, the transitions that are actually allowed produce an evenly spaced spectrum. It is conventional to define the rotational constant in wavenumber units,

B≡h8π2cI(cm−1),E~J=BJ(J+1) (cm−1)B \equiv \frac{h}{8\pi^2 c I} \quad \text{(cm}^{-1}\text{)}, \qquad \tilde E_J = B J(J+1)\ \text{(cm}^{-1}\text{)}

where:

  • hh is Planck's constant, cc the speed of light,
  • E~J=EJ/(hc)\tilde E_J = E_J/(hc) expresses the energy as a spectroscopic term value in cm⁻¹,
  • BB packages all the molecule-specific structural information (μ\mu, rer_e) into a single measurable number.

Reading: BB is large for light, short molecules (H₂, HF) and small for heavy, long ones (I₂, ICl) — rotational spacings directly encode molecular size and mass distribution.

Transition dipole moment and selection rule. The rate of absorption between rotational states ∣J,M⟩|J,M\rangle and ∣J′,M′⟩|J',M'\rangle is governed, via Fermi's Golden Rule, by the square of

μ⃗fi=⟨J′M′∣μ⃗^∣JM⟩=μ0∫YJ′M′∗(θ,ϕ) n^(θ,ϕ) YJM(θ,ϕ) dΩ\vec\mu_{fi} = \langle J' M' | \hat{\vec\mu} | J M \rangle = \mu_0 \int Y_{J'}^{M'*}(\theta,\phi)\, \hat n(\theta,\phi)\, Y_J^{M}(\theta,\phi)\, d\Omega

where:

  • μ0\mu_0 is the magnitude of the permanent dipole moment along the molecular axis,
  • n^(θ,ϕ)\hat n(\theta,\phi) is the unit vector along the bond axis, whose Cartesian components are proportional to the ℓ=1\ell=1 spherical harmonics,
  • dΩ=sin⁡θ dθ dϕd\Omega = \sin\theta\, d\theta\, d\phi is the solid-angle element,
  • the integral is purely angular because μ0\mu_0 (the radial/bond-length dependence) is factored out as a rigid-rotor constant.

Reading: the whole selection-rule problem reduces to a three-spherical-harmonic overlap integral — exactly the type of integral evaluated using Selection Rules from Symmetry and angular-momentum coupling (the Wigner–Eckart theorem in miniature).

Derivation

Derivation AThe Δ J = ±1 selection rule from the transition dipole integral

Goal. Show that ⟨J′M′∣μ⃗^∣JM⟩≠0\langle J'M'|\hat{\vec\mu}|JM\rangle \neq 0 requires μ0≠0\mu_0\neq 0 and J′=J±1J' = J\pm 1 (with M′=M,M±1M' = M, M\pm1 depending on polarization), starting only from the rigid-rotor eigenfunctions and the vector nature of μ^\hat\mu.

Step 1 — Write the dipole operator in terms of spherical harmonics. The permanent dipole points along the molecular axis, so its lab-frame Cartesian components are

μz=μ0cos⁡θ,μx±iμy=μ0sin⁡θ e±iϕ.\mu_z = \mu_0 \cos\theta, \qquad \mu_x \pm i\mu_y = \mu_0 \sin\theta\, e^{\pm i\phi}.

These are, up to normalization constants, exactly the ℓ=1\ell=1 spherical harmonics:

cos⁡θ∝Y10(θ,ϕ),sin⁡θ e±iϕ∝Y1±1(θ,ϕ).\cos\theta \propto Y_1^{0}(\theta,\phi), \qquad \sin\theta\, e^{\pm i\phi} \propto Y_1^{\pm1}(\theta,\phi).

So μ⃗^\hat{\vec\mu} behaves, angularly, exactly like a rank-1 (vector) spherical tensor operator — the same object that appears in the dipole selection rules for atomic transitions.

Step 2 — Reduce the matrix element to a three-YℓmY_\ell^m integral. The transition moment for, say, the zz-component becomes

μfi(z)=μ0∫YJ′M′∗ Y10 YJM dΩ.\mu_{fi}^{(z)} = \mu_0 \int Y_{J'}^{M'*}\, Y_1^{0}\, Y_J^{M}\, d\Omega.

Integrals of a product of three spherical harmonics are standard objects (they reduce to Clebsch–Gordan/Wigner 3-j coefficients) and vanish unless the three angular momenta satisfy a triangle condition: J1,J2,J3J_1, J_2, J_3 must be able to add like vectors, here 11, JJ, J′J'. Because the third harmonic has ℓ=1\ell=1, the only nonzero results occur when

J′=J+1, J, or J−1.J' = J+1,\ J,\ \text{or } J-1.

Step 3 — Eliminate J′=JJ' = J by parity. Spherical harmonics have definite parity under inversion, YJM(π−θ,ϕ+π)=(−1)JYJM(θ,ϕ)Y_J^M(\pi-\theta,\phi+\pi) = (-1)^J Y_J^M(\theta,\phi). The operator cos⁡θ\cos\theta (or sin⁡θ e±iϕ\sin\theta\,e^{\pm i\phi}) is odd under inversion. The integral ∫YJ′M′∗(odd function)YJM dΩ\int Y_{J'}^{M'*}(\text{odd function}) Y_J^M\, d\Omega is an integral over all space of a function with overall parity (−1)J′×(−1)1×(−1)J=(−1)J+J′+1(-1)^{J'}\times(-1)^1\times(-1)^J = (-1)^{J+J'+1}. For the integral to survive, the integrand must be even, so

(−1)J+J′+1=+1 ⟹ J+J′ must be odd.(-1)^{J+J'+1} = +1 \ \Longrightarrow\ J+J' \text{ must be odd}.

This immediately rules out J′=JJ' = J (which would make J+J′=2JJ+J'=2J even) and confirms only J′=J±1J'=J\pm1 survive.

Step 4 — Confirm the surviving cases are actually nonzero (not just parity-allowed). Explicit evaluation of the standard integral

∫0π∫02πYJ+1M∗cos⁡θ YJM sin⁡θ dθ dϕ  =  (J+1)2−M2(2J+1)(2J+3)\int_0^\pi \int_0^{2\pi} Y_{J+1}^{M *}\cos\theta\, Y_J^{M}\, \sin\theta\, d\theta\, d\phi \;=\; \sqrt{\frac{(J+1)^2 - M^2}{(2J+1)(2J+3)}}

gives a nonzero numerical coefficient for every allowed (J,M)(J,M) pair (it is never accidentally zero), confirming J→J+1J\to J+1 (and by symmetry J→J−1J \to J-1) are genuinely allowed, not merely parity-permitted.

Step 5 — Reinstate the prefactor μ0\mu_0. Every step above carried μ0\mu_0 as a multiplicative constant pulled out front. If μ0=0\mu_0 = 0 — as for any homonuclear diatomic or any molecule whose bond-axis dipole vanishes by symmetry — the entire matrix element is identically zero for every J,M,J′,M′J,M,J',M', regardless of the angular integral. Hence the complete selection rule is

 μ0≠0andΔJ=±1 \boxed{\ \mu_0 \neq 0 \quad \text{and} \quad \Delta J = \pm 1\ }

with the accompanying (often ignored in unpolarized experiments) rule ΔM=0,±1\Delta M = 0, \pm 1 depending on the polarization of the exciting radiation relative to the lab-frame quantization axis.

Go deeper · Derivation B · 1 minLine positions and the spacing 2B
Derivation BLine positions and the spacing 2B

Goal. Convert the selection rule into predicted absorption frequencies and show the spectrum is an evenly spaced comb with spacing 2B2B, then show how a measured spacing yields the bond length.

Step 1 — Absorption from level JJ to J+1J+1. Absorption requires the molecule to gain energy, so only ΔJ=+1\Delta J = +1 contributes to an absorption spectrum (the −1-1 branch corresponds to stimulated emission from an already-populated higher level, relevant in emission/maser contexts but not to a simple absorption scan). The photon energy needed is

ΔE~=E~J+1−E~J=B(J+1)(J+2)−BJ(J+1).\Delta \tilde E = \tilde E_{J+1} - \tilde E_J = B(J+1)(J+2) - BJ(J+1).

Step 2 — Expand and simplify.

ΔE~=B[(J+1)(J+2)−J(J+1)]=B(J+1)[(J+2)−J]=2B(J+1).\Delta \tilde E = B\Big[(J+1)(J+2) - J(J+1)\Big] = B(J+1)\Big[(J+2)-J\Big] = 2B(J+1).

Step 3 — Tabulate the comb. For J=0,1,2,3,…J = 0,1,2,3,\dots (the quantum number of the lower state), the absorbed photon wavenumbers are

ν~J→J+1=2B, 4B, 6B, 8B,…\tilde\nu_{J\to J+1} = 2B,\ 4B,\ 6B,\ 8B,\dots

Consecutive lines are separated by exactly

Δ(ΔE~)=2B(J+2)−2B(J+1)=2B,\Delta(\Delta\tilde E) = 2B(J+2) - 2B(J+1) = 2B,

independent of JJ — a direct experimental signature that the molecule is (to a good approximation) behaving as a rigid rotor.

Step 4 — Invert to get structural information. Once BB is read off from the (constant) line spacing, the moment of inertia and bond length follow from

I=h8π2cB,re=Iμ.I = \frac{h}{8\pi^2 c B}, \qquad r_e = \sqrt{\frac{I}{\mu}}.

This is the entire logical chain that lets a microwave spectrometer report a bond length to five or six significant figures: measure a frequency spacing, divide by 22, invert two closed-form relations.

Worked examples

Example 1 — Foundational: energy levels and the first three lines of ¹²C¹⁶O. Carbon monoxide has a well known rotational constant B=1.9313 cm−1B = 1.9313\ \text{cm}^{-1} (a nice asymmetric, non-round number to keep every step honest). Compute the term values for J=0,1,2,3J=0,1,2,3 and the first three absorption line positions.

E~J=BJ(J+1)=1.9313 J(J+1) cm−1\tilde E_J = B J(J+1) = 1.9313\, J(J+1)\ \text{cm}^{-1}
  • J=0J=0: E~0=1.9313(0)(1)=0 cm−1\tilde E_0 = 1.9313(0)(1) = 0\ \text{cm}^{-1}
  • J=1J=1: E~1=1.9313(1)(2)=3.8626 cm−1\tilde E_1 = 1.9313(1)(2) = 3.8626\ \text{cm}^{-1}
  • J=2J=2: E~2=1.9313(2)(3)=11.588 cm−1\tilde E_2 = 1.9313(2)(3) = 11.588\ \text{cm}^{-1}
  • J=3J=3: E~3=1.9313(3)(4)=23.176 cm−1\tilde E_3 = 1.9313(3)(4) = 23.176\ \text{cm}^{-1}

Absorption lines, using ν~J→J+1=2B(J+1)\tilde\nu_{J\to J+1} = 2B(J+1):

ν~0→1=2(1.9313)(1)=3.8626 cm−1\tilde\nu_{0\to1} = 2(1.9313)(1) = 3.8626\ \text{cm}^{-1}
ν~1→2=2(1.9313)(2)=7.7252 cm−1\tilde\nu_{1\to2} = 2(1.9313)(2) = 7.7252\ \text{cm}^{-1}
ν~2→3=2(1.9313)(3)=11.5878 cm−1\tilde\nu_{2\to3} = 2(1.9313)(3) = 11.5878\ \text{cm}^{-1}

Check the spacing: 7.7252−3.8626=3.8626 cm−1=2B7.7252 - 3.8626 = 3.8626\ \text{cm}^{-1} = 2B, and 11.5878−7.7252=3.8626 cm−1=2B11.5878-7.7252 = 3.8626\ \text{cm}^{-1} = 2B — the constant spacing predicted in Step 3 of the second derivation holds exactly, as it must for a rigid rotor by construction. Limiting case check: as J→0J\to 0 the ground state has zero rotational energy (no zero-point rotational energy, unlike vibration) — this is a genuine feature of an angular-momentum ladder starting at J=0J=0, not an approximation artifact.

Example 2 — Applied: extracting the bond length of ¹H³⁵Cl from a measured microwave spectrum. Suppose a microwave absorption scan of gaseous HCl shows lines at (to five significant figures) 20.79420.794, 41.58841.588, and 62.382 cm−162.382\ \text{cm}^{-1} — an intentionally non-round data set as would come from a real recording. Extract BB, then II, then rer_e. Step 1 — Confirm even spacing and read off 2B2B:

41.588−20.794=20.794 cm−1,62.382−41.588=20.794 cm−141.588 - 20.794 = 20.794\ \text{cm}^{-1}, \qquad 62.382 - 41.588 = 20.794\ \text{cm}^{-1}

So 2B=20.794 cm−12B = 20.794\ \text{cm}^{-1}, giving

B=10.397 cm−1.B = 10.397\ \text{cm}^{-1}.

Step 2 — Convert to moment of inertia using B=h/(8π2cI)B = h/(8\pi^2 c I), so I=h/(8π2cB)I = h/(8\pi^2 cB). Using h=6.6261×10−34 J sh = 6.6261\times10^{-34}\ \text{J s}, c=2.9979×1010 cm/sc = 2.9979\times10^{10}\ \text{cm/s}:

I=6.6261×10−348π2(2.9979×1010)(10.397) J s2/cmI = \frac{6.6261\times10^{-34}}{8\pi^2 (2.9979\times10^{10})(10.397)} \ \text{J s}^2/\text{cm}

Evaluate the denominator: 8π2=78.9578\pi^2 = 78.957; then 78.957×2.9979×1010=2.3671×101278.957 \times 2.9979\times10^{10} = 2.3671\times10^{12}; then ×10.397=2.4611×1013\times 10.397 = 2.4611\times10^{13}.

I=6.6261×10−342.4611×1013=2.6924×10−47 kg m2I = \frac{6.6261\times10^{-34}}{2.4611\times10^{13}} = 2.6924\times10^{-47}\ \text{kg m}^2

(the units work out to kg·m² once the cm⁻¹ and cm/s are consistently converted; carrying the conversion through gives this standard textbook value region for HCl). Step 3 — Compute reduced mass. For ¹H (1.007825 u) and ³⁵Cl (34.96885 u):

μ=(1.007825)(34.96885)1.007825+34.96885 u=35.243635.9767 u=0.97962 u\mu = \frac{(1.007825)(34.96885)}{1.007825+34.96885}\ \text{u} = \frac{35.2436}{35.9767}\ \text{u} = 0.97962\ \text{u}

Convert: μ=0.97962×1.66054×10−27 kg=1.6266×10−27 kg\mu = 0.97962 \times 1.66054\times10^{-27}\ \text{kg} = 1.6266\times10^{-27}\ \text{kg}. Step 4 — Solve for bond length:

re=Iμ=2.6924×10−471.6266×10−27=1.6553×10−20 m=1.2866×10−10 mr_e = \sqrt{\frac{I}{\mu}} = \sqrt{\frac{2.6924\times10^{-47}}{1.6266\times10^{-27}}} = \sqrt{1.6553\times10^{-20}}\ \text{m} = 1.2866\times10^{-10}\ \text{m}
re≈1.2866 A˚r_e \approx 1.2866\ \text{Å}

This lands close to the accepted HCl bond length near 1.271.27–1.291.29 Å, illustrating that a handful of microwave line positions — obtained without touching the sample chemically — pin down a bond length to picometer precision. Limiting case check: if instead the two nuclei had equal mass (m1=m2m_1=m_2), the molecule would be homonuclear, μ0=0\mu_0=0 by symmetry, and this entire spectrum would not exist — consistent with the requirement derived in Step 5 of the first derivation.

Example 3 — Exam/competition level: intensity alternation vanishes, but selection-rule logic still predicts the silent partner. A common trap question: "N₂ has a large bond dipole moment change during vibration according to some naive electronegativity argument — does it show a microwave rotational spectrum?" Justify the answer using the transition dipole integral, and then contrast with the isotopically substituted, chemically identical but not-quite-symmetric case of ¹⁴N¹⁵N. Reasoning: N₂ is homonuclear: by the D∞hD_{\infty h} point group, the molecule possesses a center of inversion, and a permanent electric dipole moment (a polar vector) cannot be invariant under inversion unless it is exactly zero. Formally, μ0\mu_0 transforms as Σu+\Sigma_u^+ symmetry, while the totally symmetric ground electronic state is Σg+\Sigma_g^+; the dipole operator has ungerade parity while the molecule's charge distribution has gerade symmetry, forcing μ0=0\mu_0 = 0 identically (see Selection Rules from Symmetry for the general group-theoretic statement). Step 5 of the first derivation shows that once μ0=0\mu_0=0, the entire transition-moment integral is zero for every J,J′J,J' pair — there is no accidental "leakage" pure rotational transition, no matter how large the individual bond polarizabilities are. N₂ has zero pure-rotational microwave absorption spectrum. Now consider 14N15N^{14}\text{N}^{15}\text{N}. Chemically the bonding is unchanged — both nuclei are still nitrogen, same electronic structure, same electron cloud, hence still essentially zero electronic dipole moment to extremely good approximation. Isotopic substitution breaks the mechanical symmetry (unequal nuclear masses shift the center of mass off the geometric bond midpoint) but not the electronic symmetry; the equilibrium electron distribution is still (almost) mirror-symmetric about the original bond center. In practice 14N15N^{14}\text{N}^{15}\text{N} has an extremely tiny nonzero dipole moment (of order 10−310^{-3} D, arising from small non-Born–Oppenheimer/vibrational effects), so an extraordinarily weak pure rotational spectrum is in principle observable with sensitive enough instrumentation — but the intensity is many orders of magnitude below a normal polar molecule like CO, illustrating that the point-group-based selection rule is a robust zeroth-order truth, with only tiny corrections surviving beyond it. Numerical check on the "even spacing" logic surviving isotopic substitution: the rotational constant itself still changes measurably (since μ\mu changes with isotopic mass), so if ¹⁴N¹⁵N did show lines, they would still be spaced by 2B′2B' with B′=h/8π2cI′B' = h/8\pi^2cI' using the new reduced mass — the derivation's Step 3 (line spacing = 2B2B, independent of JJ) is a universal rigid-rotor result, unaffected by whether the transition happens to be allowed or merely extremely weak.

Limiting cases, scope & common traps

Regime / caseWhat happensWhy
Homonuclear diatomic (N₂, O₂, H₂, Cl₂)No pure rotational absorption spectrumμ0=0\mu_0 = 0 by inversion symmetry — kills the entire transition moment (Step 5)
Heteronuclear diatomic (CO, HCl, HF)Full comb of lines at 2B,4B,6B,…2B,4B,6B,\dotsμ0≠0\mu_0\neq0 and ΔJ=±1\Delta J=\pm1 both satisfied
J→JJ\to J (elastic, no rotational change)Forbidden in pure rotational absorptionParity argument in Step 3: J+J′J+J' must be odd, so ΔJ=0\Delta J=0 is always forbidden
ΔJ=±2,±3,…\Delta J = \pm 2, \pm3,\dotsForbidden for electric-dipole rotational transitionsTriangle rule with ℓ=1\ell=1 operator only connects J′=J,J±1J' = J, J\pm1; parity removes J′=JJ'=J
Heavy, long molecule (I₂, ICl)Very small BB, closely spaced lines, requires far-infrared/THz sources rather than microwaveB∝1/IB\propto 1/I and I=μre2I=\mu r_e^2 grows with mass and bond length
Light, short molecule (HF, HD)Very large BB, widely spaced lines, easily resolved by ordinary microwave sourcesSmall II
J→∞J\to\infty (very high rotational excitation)Rigid-rotor formula increasingly fails; centrifugal distortion stretches the bondReal bonds are not infinitely stiff; the rigid-rod assumption breaks down
Vibrationally excited moleculeBB becomes vibrational-state-dependent, BvB_v, shifting all lines slightlyForeshadows Vibration-Rotation Coupling — the average ⟨1/r2⟩\langle 1/r^2\rangle changes with vibrational amplitude

The chemist behind it

The rigid-rotor treatment of molecular rotation and its use in interpreting microwave spectra is inseparable from the broader development of quantum theory of angular momentum, but the experimental and instrumental breakthrough that made rotational spectroscopy a practical chemical tool is largely credited to Charles Hard Townes and collaborators in the years following World War II, when radar-derived microwave technology was redirected toward molecular spectroscopy.

Origins · 1 minCharles Hard Townes and how this was discovered
CH
Charles Hard Townes
1915–2015 · Post-war microwave spectroscopy and the birth of the maser/laser

Townes trained as a physicist at Caltech and Duke, then worked at Bell Labs during WWII on radar systems — technology that generated exactly the coherent microwave sources needed to probe molecular rotational transitions directly, rather than inferring them indirectly from older far-infrared methods. After the war he moved to Columbia University, where he and his students used surplus radar klystrons to record precision microwave absorption spectra of gases including ammonia and simple diatomics, extracting rotational constants, bond lengths, and nuclear quadrupole coupling constants with unprecedented accuracy. His deep familiarity with molecular energy-level populations and stimulated emission, developed through this microwave spectroscopy work, led directly to the invention of the maser (1953–1954) and, with Arthur Schawlow, the theoretical proposal for the laser (1958). Townes shared the 1964 Nobel Prize in Physics for this work, and his microwave spectroscopy textbook with A. L. Schawlow remains a canonical reference for the field.

History
- 1926 — Schrödinger's wave mechanics supplies the angular momentum eigenfunctions (spherical harmonics) later applied to molecular rotation.
- 1930s — Far-infrared and radiofrequency studies establish rotational constants for simple diatomics, but resolution and sensitivity remain limited.
- 1946 — Postwar release of radar klystron technology enables Townes and coworkers at Columbia to build the first modern microwave spectrometers.
- 1949–1955 — Systematic microwave rotational spectroscopy of diatomics and light polyatomics establishes bond lengths, dipole moments, and nuclear quadrupole couplings as routine measurements.
- 1964 — Townes, Basov, and Prokhorov awarded the Nobel Prize in Physics for maser/laser principles rooted in this microwave spectroscopy tradition.
Why doesn't a homonuclear diatomic show any rotational spectrum at all, even weakly?

A homonuclear diatomic like N₂ or O₂ has a center of inversion symmetry, and its electronic charge distribution is exactly mirror-symmetric about the bond midpoint in the absence of any distorting influence. A permanent electric dipole moment is a vector quantity that must reverse sign under inversion; since the molecule's ground-state charge density is invariant under inversion, the only value the dipole moment can consistently take is zero. Because the entire transition dipole matrix element carries μ0\mu_0 as a multiplicative prefactor, setting μ0=0\mu_0=0 forces every rotational transition moment to vanish identically, not just approximately. This is a strict symmetry statement, not a matter of degree — it is why homonuclear diatomics are used as inert bath gases in microwave spectroscopy experiments without interfering with the spectrum of the polar analyte.

Do rotational Raman spectra follow the same ΔJ = ±1 rule?

No — rotational Raman scattering follows a different selection rule, ΔJ=0,±2\Delta J = 0, \pm 2, because the relevant coupling is to the molecular polarizability rather than the permanent dipole moment, and polarizability is a rank-2 tensor rather than a rank-1 vector. This means homonuclear diatomics such as N₂ and O₂, which are microwave-silent because μ0=0\mu_0=0, are perfectly Raman-active in rotation because their polarizability still changes as they tumble. Rotational Raman spectroscopy is therefore the standard technique for obtaining bond lengths of symmetric, nonpolar diatomics that pure rotational microwave absorption cannot see at all.

Why is there no J = 0 rotational ground-state energy analogous to vibrational zero-point energy?

Rotational and vibrational zero-point behavior arise from different physics. Vibrational zero-point energy exists because the harmonic oscillator Hamiltonian has an uncertainty-principle floor — you cannot simultaneously localize a particle at the potential minimum and give it zero momentum. Rotational angular momentum, however, is quantized starting at J=0J=0, which corresponds to a spherically symmetric, structureless wavefunction (Y00Y_0^0 is a constant) with genuinely zero angular momentum and zero rotational kinetic energy. There is no analogous uncertainty-principle floor forcing J≥1J\ge 1; a molecule can, and typically does at low temperature, occupy the true J=0J=0 ground rotational state with exactly zero rotational energy.

How does temperature affect which rotational lines appear most intensely?

Line intensities in a rotational spectrum are governed by the Boltzmann population of the initial state combined with the (2J+1)(2J+1) degeneracy of each level, giving a population factor proportional to (2J+1)e−BJ(J+1)hc/kBT(2J+1)e^{-BJ(J+1)hc/k_BT}. At very low JJ the degeneracy factor is small, while at very high JJ the Boltzmann exponential suppresses the population; the result is that intensities rise from J=0J=0, pass through a maximum at some intermediate Jmax⁡J_{\max} that depends on the ratio B/(kBT)B/(k_BT), and then fall off. This is exactly the same population-weighting logic that later produces the intensity envelope of the P and R branches once vibration is included, in Vibration-Rotation Spectra: P, Q, R Branches.

Is the rigid-rotor model exact, or does it break down in real spectra?

The rigid-rotor model is an idealization that assumes a perfectly inextensible bond; real molecules have bonds that stretch slightly under the centrifugal force of rotation, especially at high JJ, which slightly lowers the energy relative to the rigid-rotor prediction. This is corrected by adding a centrifugal-distortion term, −DJ2(J+1)2-DJ^2(J+1)^2, to the energy expression, where DD is a small positive constant much smaller than BB. For most rotational spectroscopy at low to moderate JJ, the rigid-rotor formula EJ=BJ(J+1)E_J=BJ(J+1) is accurate to a very good approximation, and the distortion correction is treated as a perturbative refinement rather than a fundamental change in the selection rule or overall spectral pattern.

Why do we need Fermi's Golden Rule at all if the selection rule already tells us which transitions are allowed?

The selection rule (ΔJ=±1\Delta J=\pm1, μ0≠0\mu_0\neq0) only tells you whether the transition dipole moment integral is zero or nonzero — it is a yes/no gate. Fermi's Golden Rule is what converts the nonzero transition dipole moment into an actual, quantitative transition rate, incorporating the density of radiation states, the strength of the incident field, and the specific numerical value of the (nonzero) matrix element computed in the derivation above. In short, the selection rule determines which lines can appear at all, while the golden-rule machinery from Fermi's Golden Rule and Transition Rates determines how strong each allowed line actually is.

Why does the moment of inertia formula use the reduced mass rather than the total mass of the molecule?

The two-body rotation problem (two point masses connected by a rigid rod, rotating about their common center of mass) can be exactly transformed into an equivalent one-body problem: a single particle of mass μ=m1m2/(m1+m2)\mu = m_1m_2/(m_1+m_2) orbiting at the full bond distance rer_e from a fixed center. This reduction is a standard result of classical and quantum two-body mechanics and is why the reduced mass, not the total mass, enters I=μre2I=\mu r_e^2 and therefore BB. Using the total mass instead would badly overestimate the moment of inertia for asymmetric diatomics like HCl, where the light hydrogen atom does almost all of the orbiting around a nearly stationary chlorine.

Before this

Where this leads