The Rigid Rotor Model

Point a radio telescope at a cold molecular cloud and the loudest thing in the millimetre sky is a single sharp line from carbon monoxide at 115 GHz. That line, and the whole comb beside it, is the fingerprint of a molecule turning end over end — and the theory that predicts where the rungs fall is the rigid rotor. It borrows two results wholesale: from the Born–Oppenheimer approximation, the licence to freeze the electrons and treat the nuclei as a mechanical object; from the quantum theory of angular momentum, the eigenvalue L^2YJM=ℏ2J(J+1) YJM\hat L^2 Y_J^{M}=\hbar^2 J(J+1)\,Y_J^{M}, the single fact that turns a spinning dumbbell into a ladder of discrete energies. The move is deliberately brutal: hold the bond rigidly at its equilibrium length ReR_e, switch off vibration, and let only the two orientation angles survive. What remains is a mass point gliding on a sphere — one of the handful of exactly solvable quantum problems, its wavefunctions the spherical harmonics that shape atomic orbitals. The payoff is EJ=BJ(J+1)E_J=BJ(J+1), with the rotational constant BB fixed entirely by the moment of inertia I=μRe2I=\mu R_e^2. Because BB depends only on geometry, the logic runs backwards: measure the line spacing and you have read the bond length to four or five significant figures — the most precise ruler in molecular physics. The (2J+1)(2J+1) degeneracy sets how thermal population spreads across the ladder, sculpting the envelope of every band. Everything real — centrifugal stretching, vibration–rotation coupling — is a correction bolted onto this rigid skeleton.

Key takeaways
  • The rigid rotor freezes a diatomic at its equilibrium bond length ReR_e, discarding vibration so only the two angles (θ,ϕ)(\theta,\phi) remain — a point mass on a sphere.
  • The angular Schrödinger equation gives EJ=ℏ22IJ(J+1)E_J=\dfrac{\hbar^2}{2I}J(J+1) with J=0,1,2,…J=0,1,2,\dots and I=μRe2I=\mu R_e^2; the eigenfunctions are the spherical harmonics YJMY_J^{M}.
  • Each level is (2J+1)(2J+1)-fold degenerate in MJM_J, so levels spread apart as 2B(J+1)2B(J+1) while the pure-rotation lines stay equally spaced by 2B2B.
  • Because I=μRe2I=\mu R_e^2, one spacing yields the bond length: CO's 115.271 GHz line (B~≈1.92 cm−1\tilde B\approx1.92\ \mathrm{cm^{-1}}) gives Re≈1.13 A˚R_e\approx1.13\ \text{Å}.
  • The thermal envelope (2J+1)e−BJ(J+1)/kBT(2J+1)e^{-BJ(J+1)/k_BT} peaks at a nonzero Jmax⁡≈kBT/2B−12J_{\max}\approx\sqrt{k_BT/2B}-\tfrac12 — for CO at room temperature near J≈7J\approx7, not J=0J=0.
  • Key idea: freezing the bond length turns a diatomic into a particle on a sphere whose angular momentum is quantized, giving EJ=ℏ22IJ(J+1)E_J=\dfrac{\hbar^2}{2I}J(J+1) set entirely by the moment of inertia I=μRe2I=\mu R_e^2.
The core insight

Freezing the bond length reduces a diatomic to a particle on a sphere; its angular momentum is quantized as L^2→ℏ2J(J+1)\hat L^2\to\hbar^2 J(J+1), giving EJ=ℏ22IJ(J+1)E_J=\dfrac{\hbar^2}{2I}J(J+1) fixed entirely by the moment of inertia I=μRe2I=\mu R_e^2 — so the spacing of rotational lines directly reads out the bond length.

Where this lives today

The rigid rotor is quietly everywhere a molecule is counted by its rotation. In radio and millimetre astronomy it is foundational: the 115.271 GHz115.271\ \mathrm{GHz} CO J=1→0J{=}1\to0 line worked above is the single most important tracer of molecular gas in galaxies, and instruments like ALMA map star-forming clouds by imaging the CO ladder rung by rung, reading gas mass, density, and temperature from the level energies and thermal envelope — every such measurement resting on EJ=BJ(J+1)E_J=BJ(J+1) and its distortion correction. Closer to home, laboratory microwave and terahertz spectrometers run the same physics in reverse — fixing molecular geometry to a thousandth of an ångström, resolving isotopologues, and identifying trace species by their unique comb of lines (the "fingerprinting" used to detect organics in space and monitor atmospheric gases). Rotational Raman scattering extends the reach to homonuclear N2\mathrm{N_2} and O2\mathrm{O_2}, letting lidar sound atmospheric temperature. In each case the rigid rotor supplies the scaffold; distortion, vibration–rotation coupling, and nuclear-spin statistics are the refinements on top.

The physical picture

Picture a dumbbell: two nuclei of masses m1m_1 and m2m_2 joined by a stiff, massless rod of length ReR_e. Classically it tumbles freely about its centre of mass; quantum mechanically its orientation is a wavefunction smeared over a sphere. Removing the centre-of-mass motion collapses the two-body problem to one particle of reduced mass μ=m1m2/(m1+m2)\mu=m_1m_2/(m_1+m_2) pinned to a sphere of radius ReR_e, whose allowed wavefunctions are precisely the spherical harmonics YJM(θ,ϕ)Y_J^{M}(\theta,\phi) — the same angular functions that dress the hydrogen atom's orbitals. In this sense the rigid rotor is the hydrogen atom's angular problem with the radial part amputated, and the whole apparatus of L^2\hat L^2 quantization comes along for free. The consequence is a ladder whose rungs are not evenly spaced. Since EJ∝J(J+1)E_J\propto J(J+1), the gap from JJ to J+1J+1 is 2B(J+1)2B(J+1), widening as you climb. Yet — the point everyone trips over — the pure-rotation spectrum is a comb of equally spaced lines, because each line sits at a gap 2B(J+1)2B(J+1) and the next at 2B(J+2)2B(J+2), differing by a fixed 2B2B. Heavy, long-bonded molecules have large II, small BB, a densely packed ladder; light molecules like H2\mathrm{H_2} have tiny II, enormous BB, rungs flung far apart — which is why H2\mathrm{H_2}'s first transition lies deep in the infrared while CO rotates in the microwave.

Figure
The diatomic frozen at bond length R_e tumbles as a rigid dumbbell of reduced mass μ, and its energy ladder E_J = B̃·J(J+1) climbs with ever-widening rungs — for CO with B̃ ≈ 1.92 cm⁻¹ the levels J = 0,1,2,3,4,5 sit at 0, 3.85, 11.54, 23.07, 38.45, 57.68 cm⁻¹, spaced by 2B̃, 4B̃, 6B̃, ...
Particle on a sphere → an unevenly spaced ladder E_J = B̃ J(J+1) rigid rotor (reduced mass μ) m₁ m₂ centre of mass R_e I = μ R_e² μ = m₁m₂/(m₁+m₂) E J=0 0 J=1 3.85 J=2 11.54 J=3 23.07 J=4 38.45 J=5 57.68 cm⁻¹ 2B̃ 4B̃ 6B̃ gaps grow as 2B̃(J+1); each level is (2J+1)-fold degenerate

The (2J+1)(2J+1) degeneracy counts the ways the angular-momentum vector can tip relative to a lab axis — the projections MJ=−J,…,+JM_J=-J,\dots,+J. With no field these share one energy, but they act as a statistical weight: the population of level JJ goes as (2J+1)e−BJ(J+1)/kBT(2J+1)e^{-BJ(J+1)/k_BT}, a rising degeneracy against a falling Boltzmann factor. That product peaks at a nonzero Jmax⁡J_{\max}, so a warm gas does not sit in its ground state — it spreads across dozens of levels, and its brightest line comes from well up the ladder. This is what turns a rotational band into a thermometer.

Figure
Pure-rotation lines fall at equal intervals of 2B̃ ≈ 3.85 cm⁻¹, but their strengths trace the thermal envelope P_J ∝ (2J+1)·e^(−B̃J(J+1)/kT); for CO at T = 300 K this crests at J ≈ 7 rather than at J = 0, because the rising 2J+1 degeneracy fights the falling Boltzmann factor
Thermal envelope of the CO rotational lines at 300 K — peak at J ≈ 7, not J = 0 rotational quantum number J (lines equally spaced by 2B̃ = 3.85 cm⁻¹) relative population P_J J_max = 7 (2J+1)·e^(−B̃J(J+1)/kT) 0 6 8 12 17

Finally, the geometry that makes this a measuring instrument: because I=μRe2I=\mu R_e^2 and nothing else, a single number — the spacing 2B2B — pins down ReR_e to a precision no other simple model can match. A caricature so crude it forbids the molecule to vibrate still delivers bond lengths good to a thousandth of an ångström.

The core equations

The central result is the energy ladder. It follows from writing the rotational energy as L^2/2I\hat L^2/2I and inserting the angular-momentum eigenvalue; the derivation below fills in the steps, but the formula is the thing to memorize — every downstream result hangs on it.

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

where: EJE_J — energy of level JJ; ℏ\hbar — reduced Planck constant; II — moment of inertia about the centre of mass; JJ — rotational quantum number, a non-negative integer. Physical reading: the energy climbs as J(J+1)J(J+1), so J=0J=0 sits at exactly zero — no rotational zero-point energy — and the rungs fan apart as you go up. The moment of inertia is where geometry enters. Removing the centre-of-mass motion recasts the two nuclei as one particle of reduced mass μ\mu at fixed radius ReR_e, with inertia μRe2\mu R_e^2. This is the only place the masses and bond length appear — why measuring II measures ReR_e.

I=μRe2,μ=m1m2m1+m2(2)I=\mu R_e^2,\qquad \mu=\frac{m_1 m_2}{m_1+m_2} \qquad (2)

where: μ\mu — reduced mass; ReR_e — fixed (equilibrium) internuclear separation; m1,m2m_1,m_2 — nuclear masses. Physical reading: a light-atom pair (small μ\mu) or a short bond gives a small II and thus a large, widely spaced spectrum. Spectroscopists rarely quote BB in joules. The energy constant B=ℏ2/2IB=\hbar^2/2I is converted to a wavenumber B~\tilde B (in cm−1\mathrm{cm^{-1}}) through E=hcB~E=hc\tilde B, the form catalogued for every molecule. Both are called "the rotational constant".

B=ℏ22I  (energy),B~=Bhc=ℏ4πc I=h8π2c I  (cm−1)(3)B=\frac{\hbar^2}{2I}\ \ (\text{energy}),\qquad \tilde B=\frac{B}{hc}=\frac{\hbar}{4\pi c\,I}=\frac{h}{8\pi^2 c\,I}\ \ (\mathrm{cm^{-1}}) \qquad (3)

where: BB — rotational constant in energy units; B~\tilde B — in wavenumbers; hh — Planck constant; cc — speed of light. Physical reading: B~∝1/I\tilde B\propto1/I, so a big, heavy molecule shows a small B~\tilde B and a finely spaced comb. Because EJ∝J(J+1)E_J\propto J(J+1), adjacent levels are not equally spaced. Subtracting two neighbours gives a gap that grows linearly with JJ — the algebraic fact behind the "levels spread, lines don't" confusion.

EJ+1−EJ=2B(J+1)(4)E_{J+1}-E_J=2B(J+1) \qquad (4)

where: J+1J+1 — the upper member of the pair. Physical reading: the gaps run 2B2B, 4B4B, 6B6B; but since the ΔJ=+1\Delta J=+1 absorption lines land on these gaps, consecutive lines differ by a constant 2B2B. Each level carries a hidden multiplicity: the energy is blind to the orientation quantum number MJM_J, which runs over 2J+12J+1 values, so every level is (2J+1)(2J+1)-fold degenerate — the statistical weight that shapes band intensities.

gJ=2J+1,MJ=−J,−J+1,…,+J(5)g_J=2J+1,\qquad M_J=-J,-J+1,\dots,+J \qquad (5)

where: gJg_J — number of degenerate sublevels; MJM_J — projection of angular momentum on the lab axis. Physical reading: higher levels are "wider" targets for thermal population, which is why the envelope climbs before the Boltzmann factor drags it back down. Underneath the ladder is a differential equation. With the bond frozen the Hamiltonian is purely angular, and the Schrödinger equation is the eigenvalue problem for L^2\hat L^2 on a sphere — solved by the spherical harmonics.

−ℏ22I[1sin⁡θ∂∂θ ⁣(sin⁡θ∂∂θ)+1sin⁡2θ∂2∂ϕ2]ψ=Eψ,ψJ,M=YJM(θ,ϕ)(6)-\frac{\hbar^2}{2I}\left[\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\!\left(\sin\theta\frac{\partial}{\partial\theta}\right)+\frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\right]\psi=E\psi,\qquad \psi_{J,M}=Y_J^{M}(\theta,\phi) \qquad (6)

where: θ,ϕ\theta,\phi — orientation angles of the bond axis; ψ\psi — rotational wavefunction; YJMY_J^{M} — spherical harmonic of degree JJ, order MM. Physical reading: the bracket is exactly −L^2/ℏ2-\hat L^2/\hbar^2, so the equation says "energy = L^2/2I\hat L^2/2I", and its eigenvalue ℏ2J(J+1)\hbar^2 J(J+1) hands back equation (1). The last relation governs what a spectrum looks like. Weighting each level by its degeneracy and Boltzmann factor gives the relative population; its derivative set to zero locates the most populated level.

PJ∝(2J+1) e−BJ(J+1)/kBT,Jmax⁡≈kBT2B−12(7)P_J\propto (2J+1)\,e^{-BJ(J+1)/k_BT},\qquad J_{\max}\approx\sqrt{\frac{k_BT}{2B}}-\tfrac12 \qquad (7)

where: PJP_J — relative population of level JJ; kBk_B — Boltzmann constant; TT — absolute temperature; Jmax⁡J_{\max} — most populated level. Physical reading: Jmax⁡∝TJ_{\max}\propto\sqrt{T}, so a hotter gas populates higher rungs and its band crests farther from J=0J=0.

Derivation AFrom the two-body rotor to the J(J+1) ladder

The goal is to descend from a real diatomic — two nuclei bound in a potential well — to the eigenvalue formula EJ=ℏ2J(J+1)/2IE_J=\hbar^2 J(J+1)/2I, making every idealization explicit. We separate the centre-of-mass motion to collapse the two-body problem to one particle, freeze the bond so the radial coordinate drops out, promote the classical energy to an operator, and borrow the one result from angular-momentum theory that does the quantizing. The point is that the rigid rotor imports its quantization rather than generating it: the spectrum is the L^2\hat L^2 spectrum in disguise.

Step 1. Split the nuclear coordinates into the centre-of-mass position and the relative vector r=r2−r1\mathbf r=\mathbf r_2-\mathbf r_1. The kinetic energy separates, and the relative motion is that of one particle of reduced mass μ=m1m2/(m1+m2)\mu=m_1m_2/(m_1+m_2) in the bond potential. Why: the centre-of-mass drifts freely and carries no spectroscopic information, so discarding it is exact — two bodies become one.

Step 2. Impose rigidity: fix ∣r∣=Re|\mathbf r|=R_e. The radial kinetic energy and V(Re)V(R_e) become constants (absorb V(Re)V(R_e) into the energy zero), leaving only motion on a sphere, with E=L2/2IE=L^2/2I and I=μRe2I=\mu R_e^2. Why: freezing the bond is the defining idealization — it trades vibration for an exactly solvable problem, valid because rotational quanta are far smaller than vibrational ones.

Step 3. Promote to quantum mechanics: L2→L^2L^2\to\hat L^2, so H^=L^2/2I\hat H=\hat L^2/2I. The eigenfunctions of L^2\hat L^2 on the sphere are the spherical harmonics YJMY_J^{M}, with J=0,1,2,…J=0,1,2,\dots and M=−J,…,+JM=-J,\dots,+J. Why: with no potential and no radial part, the eigenfunctions are fixed entirely by the sphere's geometry — a solved problem inherited from angular-momentum theory.

Step 4. Insert L^2YJM=ℏ2J(J+1) YJM\hat L^2 Y_J^{M}=\hbar^2 J(J+1)\,Y_J^{M} to read off EJ=ℏ2J(J+1)/2IE_J=\hbar^2 J(J+1)/2I. Since EJE_J is independent of MM, each level is (2J+1)(2J+1)-fold degenerate. Why: this eigenvalue is where all the ladder's physics lives — the J(J+1)J(J+1) shape, the zero-point-free ground state, and the degeneracy fall out at once. Cross-check: at J=0J=0, E0=0E_0=0 and the first transition costs 2B2B, matching 2B(J+1)2B(J{+}1) at J=0J=0.

Go deeper · Derivation B · 1 minReading the bond length off the rotational constant
Derivation BReading the bond length off the rotational constant

The goal is to turn a spectrum into a geometry: from the equally spaced comb of a pure-rotation spectrum, extract ReR_e to several significant figures. We locate the line positions from the ΔJ=+1\Delta J=+1 rule, identify the constant spacing with 2B~2\tilde B, express B~\tilde B through the moment of inertia, and invert I=μRe2I=\mu R_e^2 using isotope masses known to far higher precision than the spectrum. This is the workhorse calculation of microwave spectroscopy — why a rigid caricature yields real molecular structure.

Step 1. With the electric-dipole rule ΔJ=+1\Delta J=+1, absorption from level JJ lands at ν~J→J+1=2B~(J+1)\tilde\nu_{J\to J+1}=2\tilde B(J+1), so the lines sit at 2B~, 4B~, 6B~,…2\tilde B,\,4\tilde B,\,6\tilde B,\dots with constant spacing 2B~2\tilde B. Why: the observable is a spacing, not an absolute energy, so reading it off the comb needs no knowledge of the zero point.

Step 2. Convert the energy constant B=ℏ2/2IB=\hbar^2/2I to wavenumbers via E=hcB~E=hc\tilde B: B~=ℏ/4πcI=h/8π2cI\tilde B=\hbar/4\pi c I=h/8\pi^2 c I, so a measured B~\tilde B fixes II immediately. Why: II is the only molecular quantity in B~\tilde B — a clean one-to-one map from spectrum to inertia, no free parameters.

Step 3. Insert I=μRe2I=\mu R_e^2 and solve: Re=ℏ/(4πc μ B~)R_e=\sqrt{\hbar/(4\pi c\,\mu\,\tilde B)}, with μ\mu from tabulated isotope masses. Why: because μ\mu is known to many digits and B~\tilde B to parts in 10610^6, the geometric output is limited only by the rigid-rotor idealization, not the arithmetic. Dimensional check: ℏ/(c μ B~)\hbar/(c\,\mu\,\tilde B) has units (J s)/(m s−1⋅kg⋅m−1)=m2\mathrm{(J\,s)/(m\,s^{-1}\cdot kg\cdot m^{-1})}=\mathrm{m^2}, so its root is a length.

Go deeper · Derivation C · 1 minWhy the band peaks away from J = 0
Derivation CWhy the band peaks away from J = 0

The goal is to explain quantitatively why a warm gas glows brightest not from its rotational ground state but from a level several rungs up, and to derive Jmax⁡≈kBT/2B−12J_{\max}\approx\sqrt{k_BT/2B}-\tfrac12. The tension is between the degeneracy 2J+12J+1 that rewards high JJ and the Boltzmann factor e−BJ(J+1)/kBTe^{-BJ(J+1)/k_BT} that punishes it. We treat JJ as continuous, differentiate the product, and locate the maximum — statistical mechanics riding on the rigid-rotor levels.

Step 1. Write PJ∝(2J+1) e−BJ(J+1)/kBTP_J\propto(2J+1)\,e^{-BJ(J+1)/k_BT} and treat JJ as continuous, valid when many levels are populated. Why: with Jmax⁡J_{\max} of order ten, smoothing the discrete ladder introduces only a fractional error ∼1/Jmax⁡\sim1/J_{\max}, fine for an estimate.

Step 2. Set dPJ/dJ=0dP_J/dJ=0. With J(J+1)≈(J+12)2J(J+1)\approx(J+\tfrac12)^2, the logarithmic derivative gives (2J+1)2=2kBT/B(2J+1)^2=2k_BT/B. Why: the zero of the derivative is where the rising degeneracy and falling exponential balance — the crest of the envelope.

Step 3. Solve: 2Jmax⁡+1=2kBT/B2J_{\max}+1=\sqrt{2k_BT/B}, hence Jmax⁡=kBT/2B−12J_{\max}=\sqrt{k_BT/2B}-\tfrac12. Why: the T\sqrt T scaling is the fingerprint of the competition — it makes band envelopes usable as gas thermometers over orders of magnitude in temperature. For CO (B~≈1.92 cm−1\tilde B\approx1.92\ \mathrm{cm^{-1}}) at 300 K300\ \mathrm{K} (kBT≈208 cm−1k_BT\approx208\ \mathrm{cm^{-1}}), Jmax⁡≈208/3.85−0.5≈6.9J_{\max}\approx\sqrt{208/3.85}-0.5\approx6.9 — the band crests near J=7J=7, as the figure shows.

Limiting cases

How the rigid rotor behaves in the physically important regimes:

CaseResultMeaning
Ground stateE0=0E_0=0No rotational zero-point energy — unlike the harmonic oscillator, a non-rotating state Y00=constY_0^0=\text{const} is allowed, with L=0L=0 and completely uncertain orientation.
Low-JJ, cold limitE1−E0=2BE_1-E_0=2BOnly the lowest rungs are populated; the spectrum is a few well-separated lines, and kBT≳2Bk_BT\gtrsim2B is needed to populate J=1J=1 at all.
High-JJ, classical limitEJ→ℏ2J22I=L22IE_J\to\dfrac{\hbar^2 J^2}{2I}=\dfrac{L^2}{2I}For large JJ, J(J+1)→J2J(J+1)\to J^2 and the ladder becomes a quasi-continuum recovering the classical energy L2/2IL^2/2I — the Bohr correspondence limit.
Non-rigid correctionEJ=BJ(J+1)−D J2(J+1)2E_J=BJ(J+1)-D\,J^2(J+1)^2Centrifugal stretching lengthens the bond at high JJ, lowering levels through the distortion constant D>0D>0; the rigid rotor is the D→0D\to0 limit.
Homonuclear moleculeno permanent dipoleN2,O2,H2\mathrm{N_2},\mathrm{O_2},\mathrm{H_2} still obey EJ=BJ(J+1)E_J=BJ(J+1), but have no pure-rotation absorption — their levels show up only in rotational Raman scattering.

Worked examples

1 · The bond length of CO from its radio line (confidence builder)

CO's lowest pure-rotation transition J=0→1J{=}0\to1 is observed at 115.271 GHz115.271\ \mathrm{GHz} — the brightest millimetre line in the interstellar medium. Find CO's moment of inertia and bond length, using μCO=6.856 u\mu_{\mathrm{CO}}=6.856\ \mathrm{u} (12C16O^{12}\mathrm{C}^{16}\mathrm{O}), 1 u=1.6605×10−27 kg1\ \mathrm{u}=1.6605\times10^{-27}\ \mathrm{kg}.

Step 1. The J=0→1J{=}0\to1 line lies at 2B~2\tilde B, so B~=ν/(2c)=115.271×109/(2×2.998×1010)=1.9225 cm−1\tilde B=\nu/(2c)=115.271\times10^9/(2\times2.998\times10^{10})=\mathbf{1.9225\ cm^{-1}} (equivalently B=hν/2=3.82×10−23 J≈0.24 meVB=h\nu/2=3.82\times10^{-23}\ \mathrm{J}\approx0.24\ \mathrm{meV}, squarely microwave).

Step 2. Get the inertia from I=h/(8π2c B~)I=h/(8\pi^2 c\,\tilde B). With B~=192.25 m−1\tilde B=192.25\ \mathrm{m^{-1}}: I=6.626×10−34/(78.96×2.998×108×192.25)=1.456×10−46 kg m2I=6.626\times10^{-34}/(78.96\times2.998\times10^8\times192.25)=\mathbf{1.456\times10^{-46}\ kg\ m^2}.

Step 3. Invert I=μRe2I=\mu R_e^2. With μ=6.856×1.6605×10−27=1.1385×10−26 kg\mu=6.856\times1.6605\times10^{-27}=1.1385\times10^{-26}\ \mathrm{kg}: Re=I/μ=1.279×10−20=1.131×10−10 mR_e=\sqrt{I/\mu}=\sqrt{1.279\times10^{-20}}=1.131\times10^{-10}\ \mathrm{m}. R_e ≈ 1.131 Å, within a percent of the accepted 1.128 A˚1.128\ \text{Å} — the gap is the vibrational-averaging correction the rigid rotor ignores. One radio line has just measured a chemical bond.

2 · Where CO's band peaks, and which line that is

For the same CO (B~=1.9225 cm−1\tilde B=1.9225\ \mathrm{cm^{-1}}) at T=300 KT=300\ \mathrm{K}, find the most populated level and the wavenumber of the line it emits. Take kBT=208.5 cm−1k_BT=208.5\ \mathrm{cm^{-1}}.

Step 1. With everything in wavenumbers, Jmax⁡=kBT/2B−12=208.5/(2×1.9225)−0.5=54.2−0.5=6.9J_{\max}=\sqrt{k_BT/2B}-\tfrac12=\sqrt{208.5/(2\times1.9225)}-0.5=\sqrt{54.2}-0.5=6.9. Rounding, J_max = 7 — the band crests seven rungs above the ground state, not at J=0J=0.

Step 2. The line from the crest is J=7→8J{=}7\to8, at ν~=2B~(J+1)=2×1.9225×8=30.76 cm−1\tilde\nu=2\tilde B(J+1)=2\times1.9225\times8=\mathbf{30.76\ cm^{-1}} (about 922 GHz922\ \mathrm{GHz}). It is strongest not because it is highest but because (2J+1)e−BJ(J+1)/kBT(2J+1)e^{-BJ(J+1)/k_BT} maximizes there.

Step 3. Sanity check: at J=0J=0 the weight is 11, while at J=7J=7 it is 15×e−1.9225⋅56/208.5=15×e−0.516=8.915\times e^{-1.9225\cdot56/208.5}=15\times e^{-0.516}=8.9 times larger. Rising degeneracy wins until the Boltzmann factor overwhelms it near J=7J=7. The band peaks at J = 7, near 30.8 cm⁻¹, about nine times the ground-level weight.

Physics of the portrait

The rigid rotor is at once a homework exercise and a precision instrument — and the person who did most to turn its ladder into the language of chemistry read molecules line by line.

Origins · 1 minGerhard Herzberg and how this was discovered
GH
Gerhard Herzberg
1904–1999 · molecular spectroscopy and the structure of diatomics

Gerhard Herzberg was, by common consent, the twentieth century's master of molecular spectra. Forced to emigrate from Germany in 1935, he built at Saskatchewan and later Canada's National Research Council a school of spectroscopy whose output — above all the three-volume Molecular Spectra and Molecular Structure — became the definitive reference for rotational and vibrational structure. His contribution to the rigid rotor was not the equation but its exhaustive application: he and his co-workers tabulated the constants BB and DD, and hence the bond lengths, of a vast catalogue of molecules, showing that the J(J+1)J(J+1) ladder — corrected for stretching and vibration — described real molecules to spectroscopic accuracy. He identified the spectra of fleeting radicals like CH and CH2_2, extending the bookkeeping to molecules no bottle could hold. His 1971 Nobel Prize in Chemistry recognized his work on "the electronic structure and geometry of molecules, particularly free radicals" — geometry read, at its foundation, from the line spacings this article derives. His books taught generations that a rotational spectrum is a ruler, and the rigid rotor is where it is calibrated.

History
1912
— Bjerrum applies quantum ideas to molecular rotation, predicting that infrared bands carry a discrete rotational fine structure
1926
— Schrödinger's wave mechanics solves the rigid rotor exactly, its eigenfunctions the spherical harmonics
1927
— The J(J+1) energy ladder replaces the old quantum theory's flawed rotational term
1927
— Dennison explains hydrogen's heat capacity through the ortho/para nuclear-spin statistics of its rotational levels
1946
— Postwar microwave spectroscopy resolves pure-rotation lines, turning B into bond lengths good to five figures

Going further

Texts. Foot, Atomic Physics (molecular rotation) · Demtröder, Laser Spectroscopy (rotational structure of diatomics) · Herzberg, Spectra of Diatomic Molecules (the classic reference) · Atkins & Friedman, Molecular Quantum Mechanics (the rigid-rotor solution) · Saleh & Teich, Fundamentals of Photonics.

Provenance. The energy formula EJ=ℏ2J(J+1)/2IE_J=\hbar^2 J(J+1)/2I, the inertia I=μRe2I=\mu R_e^2, the degeneracy 2J+12J+1, the spacing 2B(J+1)2B(J+1), and the peak Jmax⁡≈kBT/2B−12J_{\max}\approx\sqrt{k_BT/2B}-\tfrac12 are standard results verifiable against Foot, Demtröder, Herzberg, and Atkins & Friedman. The eigenvalue L^2YJM=ℏ2J(J+1)YJM\hat L^2 Y_J^{M}=\hbar^2 J(J+1)Y_J^{M} is recalled from the quantum theory of angular momentum, and the Born–Oppenheimer separation from aq-optics-born-oppenheimer-approximation. The worked numbers were recomputed here: the CO J=1→0J{=}1\to0 line at 115.271 GHz115.271\ \mathrm{GHz} gives B~=1.9225 cm−1\tilde B=1.9225\ \mathrm{cm^{-1}}, I=1.456×10−46 kg m2I=1.456\times10^{-46}\ \mathrm{kg\,m^2}, Re=1.131 A˚R_e=1.131\ \text{Å} (accepted 1.128 A˚1.128\ \text{Å}); Jmax⁡=6.9J_{\max}=6.9 at 300 K300\ \mathrm{K} with kBT=208.5 cm−1k_BT=208.5\ \mathrm{cm^{-1}}. A verifier should reconfirm: (i) the L^2\hat L^2 eigenvalue; (ii) B~=h/8π2cI\tilde B=h/8\pi^2 cI; (iii) the dimensional check on ReR_e; (iv) Jmax⁡J_{\max} at 300 K300\ \mathrm{K}.

Builds on. The Born–Oppenheimer Approximation · The Chemical Bond and Molecular Potential Energy Curves · Angular Momentum in Quantum Mechanics. Leads to. Rotational Spectra and the Rotational Constant · Vibration–Rotation Spectra and P, Q, R Branches · Raman Scattering and the Rotational Raman Effect · Centrifugal Distortion (the non-rigid rotor).

Ask the tutor: Want to see how the same BB that fixes CO's bond length reappears as the spacing of the P and R branches in its infrared vibration–rotation band?

Why does the rigid rotor have no zero-point energy when the harmonic oscillator does?

Because the two confinements differ in kind. The oscillator confines a coordinate to a well, and the uncertainty principle then forbids the particle from sitting still at the bottom, forcing a minimum 12ℏω\tfrac12\hbar\omega. The rotor confines nothing radially — the bond length is fixed by fiat — and its "position" is an orientation. The state Y00Y_0^0 is constant over the whole sphere: the axis is completely undetermined, consistent with zero angular momentum. No coordinate is being squeezed, so there is no zero-point penalty. Hence E0=0E_0=0 exactly, and the ground state is a genuinely non-rotating molecule.

If the energy levels spread apart, why is the spectrum a set of equally spaced lines?

Both are true at once. The levels obey EJ=BJ(J+1)E_J=BJ(J+1), so the gap from JJ to J+1J+1 genuinely grows: 2B2B, 4B4B, 6B6B. But a spectral line is a transition, and with ΔJ=+1\Delta J=+1 the line from level JJ sits at that gap 2B(J+1)2B(J+1) and the next at 2B(J+2)2B(J+2). Subtract: consecutive lines differ by a constant 2B2B. The levels fan out while the lines march in lockstep — and measuring that 2B2B spacing is how BB, and then the bond length, is extracted.

What exactly is the reduced mass doing here?

The reduced mass turns a two-body rotation into a one-body problem. Two nuclei orbiting their common centre of mass is awkward directly, but the relative coordinate obeys the same mathematics as a single particle of mass μ=m1m2/(m1+m2)\mu=m_1m_2/(m_1+m_2) at distance ReR_e from a fixed centre; its inertia μRe2\mu R_e^2 equals the true two-body inertia, so nothing is lost. For a lopsided molecule like HCl, μ\mu is close to the hydrogen mass because the heavy chlorine barely moves; for a symmetric one it is half the atomic mass. Getting μ\mu right — with the correct isotope masses — is what makes the extracted bond length trustworthy to five figures.

Why can't homonuclear molecules like N₂ show a pure-rotation spectrum?

Because pure-rotation absorption requires a permanent electric dipole for the radiation to grip. A homonuclear diatomic — N2\mathrm{N_2}, O2\mathrm{O_2}, H2\mathrm{H_2} — is perfectly symmetric, so its dipole is zero and stays zero as it rotates; with no oscillating dipole to drive, the ΔJ=+1\Delta J=+1 absorption lines do not appear. The rotational levels are still there, exactly as EJ=BJ(J+1)E_J=BJ(J+1) predicts — just invisible to direct absorption. Rotational Raman scattering, which couples through polarizability rather than dipole (with ΔJ=±2\Delta J=\pm2), is how air's main gases are measured in practice.

How good is the rigid approximation, and when does it break?

Surprisingly good at low JJ, predictably worse as it spins faster. The model assumes a bond that never stretches, but under rotation a real molecule feels a centrifugal pull that lengthens ReR_e, raising II and lowering the high-JJ levels below the rigid prediction. This enters as −D J2(J+1)2-D\,J^2(J+1)^2, with DD typically 10−410^{-4}–10−610^{-6} times BB — negligible for the first rungs, measurable once JJ reaches a few tens. A second correction: BB depends weakly on the vibrational state, because a vibrating bond has a larger average ReR_e, which is why the spectroscopic B0B_0 differs slightly from the equilibrium BeB_e. Both are small and well understood — refinements on a sound skeleton, not failures of it.

How does a rotational spectrum become a thermometer?

Through the population envelope. Each line's strength scales with the population of its lower level, (2J+1)e−BJ(J+1)/kBT(2J+1)e^{-BJ(J+1)/k_BT} — rising (degeneracy) then falling (Boltzmann), peaking at Jmax⁡≈kBT/2B−12J_{\max}\approx\sqrt{k_BT/2B}-\tfrac12. Because Jmax⁡∝TJ_{\max}\propto\sqrt T, the brightest line in a band directly reads temperature: colder gas peaks at lower JJ, hotter at higher JJ, and fitting the whole envelope pins TT down tightly. This is how astronomers derive the temperature of interstellar clouds and cometary comae from CO lines, and how lidar sounds atmospheric temperature from the rotational Raman envelope of N2\mathrm{N_2}.

Before this

Where this leads