CSIR-NET Rotational and Vibrational Spectroscopy — Bond Lengths and Force Constants
Structure elucidation from IR, NMR and mass spectra is a pattern-recognition skill. This article is about the other half of spectroscopy — the quantitative half, where a microwave spectrum gives you a bond length to three decimal places and an infrared band gives you a force constant in newtons per metre. These questions are pure physical chemistry, the formulae are few, and the marks are reliable once you have practised the unit handling. Every number below is computed step by step.
Rotation: the rigid rotor
- J — rotational quantum number, 0, 1, 2, … ; each level is (2J+1)-fold degenerate.
- B — rotational constant. Expressed in cm−1 when c is in cm s−1; in Hz if c is omitted.
- I — moment of inertia, kg m2. μ — reduced mass, kg (convert from u with 1 u = 1.66054 × 10−27 kg).
- Selection rule ΔJ = ±1, and the molecule must have a permanent dipole moment. So HCl and CO give microwave spectra; N2, O2 and CO2 do not.
The transition J → J+1 appears at ν̃ = 2B(J+1), so the pure rotational spectrum is a ladder of lines separated by 2B. That even spacing is the fingerprint of a rigid rotor and the fastest way to extract B from a printed spectrum.
Worked example 1 — bond length of 12C16O from B = 1.9313 cm−1.
Step 1 — reduced mass.
μ = (12.000 × 15.9949) ÷ (12.000 + 15.9949) = 191.939 ÷ 27.9949 = 6.8562 u.
In kg: 6.8562 × 1.66054 × 10−27 = 1.13850 × 10−26 kg.
Step 2 — moment of inertia. Rearranging B = h/(8π2cI):
8π2 = 78.9568.
Denominator = 78.9568 × (2.99792 × 1010 cm s−1) × 1.9313 =
4.5715 × 1012.
I = (6.62607 × 10−34) ÷ (4.5715 × 1012) =
1.4494 × 10−46 kg m2.
Step 3 — bond length.
r = √(I/μ) = √[(1.4494 × 10−46) ÷ (1.13850 × 10−26)] =
√(1.2731 × 10−20) = 1.1283 × 10−10 m =
112.8 pm.
Cross-check. The accepted equilibrium bond length of carbon monoxide is about 112.8 pm, so the arithmetic and the unit conversions are both sound. Note that c had to be in cm s−1 because B was in cm−1 — mixing that up is the single most common way to be out by a factor of 100.
Line spacing: 2B = 3.86 cm−1, which is where in the far infrared or microwave you would look for these lines.
Two refinements the paper likes
Centrifugal distortion. A real bond stretches as the molecule spins faster, so I grows and the levels fall slightly below the rigid-rotor prediction:
The observable consequence is that the line spacing slowly decreases at high J rather than staying exactly 2B. A question that says "the spacing is not quite constant" is asking for this term.
Which line is strongest? The intensity of a rotational line depends on the population of the starting level, which is a product of a rising degeneracy (2J+1) and a falling Boltzmann factor. The maximum occurs at
Worked example 2 — most populated level of CO at 300 K.
kT = (1.380649 × 10−23) × 300 = 4.1419 × 10−21 J.
hcB = (6.62607 × 10−34)(2.99792 × 1010)(1.9313) =
3.8364 × 10−23 J, so 2hcB = 7.6729 × 10−23 J.
Ratio = 4.1419 × 10−21 ÷ 7.6729 × 10−23 = 53.98.
Jmax = √53.98 − 0.5 = 7.35 − 0.5 = 6.85, so J = 7 is the most
populated level.
This is why a rotational band peaks in the middle rather than at J = 0: the ground level has the largest Boltzmann factor but the smallest degeneracy.
Vibration: the harmonic oscillator
- v — vibrational quantum number; note the zero-point energy ½hν, which does not vanish at v = 0.
- k — force constant, N m−1; a measure of bond stiffness, not bond strength, though the two usually run together.
- Selection rule Δv = ±1 for a harmonic oscillator, and the vibration must change the dipole moment for it to be infrared active.
Worked example 3 — force constant of H35Cl from ωe = 2990 cm−1.
Step 1 — reduced mass.
μ = (1.00783 × 34.96885) ÷ (1.00783 + 34.96885) = 35.2427 ÷ 35.9767 = 0.97960 u
= 0.97960 × 1.66054 × 10−27 = 1.62666 × 10−27 kg.
Step 2 — frequency in s−1.
ν = c ν̃ = (2.99792 × 1010) × 2990 = 8.96378 × 1013 s−1.
Step 3 — force constant.
k = 4π2 ν2 μ = 39.4784 × (8.96378 × 1013)2 ×
(1.62666 × 10−27)
= 39.4784 × (8.0349 × 1027) × (1.62666 × 10−27)
= 317.21 × 1.62666 = 516 N m−1.
That agrees with the standard quoted value for the H–Cl bond, which is the check worth making before you move on.
Worked example 4 — the isotope shift, D35Cl. The force constant is a property of the electronic structure and is essentially unchanged by isotopic substitution, so only μ changes:
μ(DCl) = (2.01410 × 34.96885) ÷ (2.01410 + 34.96885) = 70.4308 ÷ 36.9830 = 1.90441 u.
ν̃(DCl)/ν̃(HCl) = √(μHCl/μDCl) = √(0.97960 ÷ 1.90441) = √0.51439 =
0.71721.
ν̃(DCl) = 2990 × 0.71721 = 2145 cm−1.
A near-halving of the wavenumber on going from H to D is the standard qualitative result, and it is why C–H stretches are deliberately shifted out of the way by deuteration in mechanistic work.
Anharmonicity — where overtones come from
A real potential is not a parabola: it flattens out towards dissociation. The Morse approximation gives
Fundamental (0→1): ν̃ = ωe − 2ωexe
First overtone (0→2): ν̃ = 2ωe − 6ωexe
Two consequences follow immediately. Levels get closer together as v rises, and Δv = ±2, ±3 transitions become weakly allowed — which is why overtones exist at all, and why the first overtone appears at slightly less than twice the fundamental rather than exactly twice.
Worked example 5 — HCl with ωe = 2990.9 cm−1 and ωexe = 52.8 cm−1.
Fundamental = 2990.9 − 2(52.8) = 2990.9 − 105.6 = 2885.3 cm−1,
matching the observed HCl fundamental near 2886 cm−1.
First overtone = 2(2990.9) − 6(52.8) = 5981.8 − 316.8 =
5665.0 cm−1 — clearly less than 2 × 2885.3 = 5770.6.
Dissociation energy, and an honest caveat. The Morse model gives De = ωe2 ÷ (4 ωexe) = (2990.9)2 ÷ (4 × 52.8) = 8 945 483 ÷ 211.2 = 42 356 cm−1. Using 1 cm−1 = 11.9627 J mol−1, that is about 507 kJ mol−1. The experimentally accepted dissociation energy of HCl is appreciably lower than this. The discrepancy is not an arithmetic error — the Morse formula systematically overestimates De for real molecules, because the true potential deviates from the Morse shape at large separation. Quote the method with the number, and say that it is an upper estimate.
Rovibrational spectra: P, Q and R branches
In the gas phase, a vibrational transition is accompanied by a change in J, which splits the band into branches.
| Branch | ΔJ | Position relative to the band centre | When it appears |
|---|---|---|---|
| P | −1 | Lower wavenumber | Always, for a diatomic |
| Q | 0 | At the band centre | Only when the molecule has vibrational angular momentum — e.g. a bending mode; absent in a diatomic stretch |
| R | +1 | Higher wavenumber | Always, for a diatomic |
The lines within each branch are spaced by roughly 2B, and there is a gap of about 4B across the missing centre where the Q branch would be. Being asked to explain the "missing central line" of the HCl gas-phase infrared band is a standard question, and the answer is simply that ΔJ = 0 is forbidden for that transition.
Raman, and how many modes there are
Raman scattering depends on a change in polarisability, not dipole moment, so it sees vibrations that the infrared cannot. Selection rules are ΔJ = 0, ±2 for rotation (giving lines spaced by 4B) and Δv = ±1 for vibration.
CO2 is the textbook illustration. It is linear with three atoms, so it has 3N − 5 = 4 vibrational modes: a symmetric stretch (Raman active, infrared inactive, because the dipole stays zero), an antisymmetric stretch (infrared active) and a doubly degenerate bend (infrared active). Water, being bent, has 3N − 6 = 3 modes and all three are active in both techniques, because it has no centre of symmetry. Benzene, with 12 atoms and non-linear, has 3(12) − 6 = 30 vibrational modes.
Stokes lines (the molecule gains energy) are more intense than anti-Stokes lines (it loses energy) simply because the lower level is more populated at ordinary temperatures — a Boltzmann argument, not a selection rule.
Mistakes that cost marks
- Using c in m s−1 with B in cm−1. Match the units first; this error changes the answer by a factor of 100 and is by far the commonest.
- Using molar masses instead of atomic masses in μ. Divide by NA, or work in u and convert with 1.66054 × 10−27 kg per u.
- Forgetting the dipole requirement. Homonuclear diatomics have no pure rotational or infrared vibrational spectrum — but they are Raman active.
- Assuming the strongest rotational line is J = 0. Degeneracy pushes the maximum out to Jmax.
- Treating the first overtone as exactly twice the fundamental. That is only true for a perfect harmonic oscillator, in which case overtones would not appear at all.
- Assuming a force constant changes on isotopic substitution. It does not, to a very good approximation — only μ does.
- Using 3N − 6 for a linear molecule. Linear molecules have only two rotational degrees of freedom, so it is 3N − 5.
- Expecting a Q branch in a diatomic stretch. It is absent, and the gap is the evidence.
Where this appears in the paper
| Sub-topic | Typical question form |
|---|---|
| Rigid rotor | Bond length from B, or B from a line spacing |
| Populations | Most intense rotational line; effect of temperature |
| Harmonic oscillator | Force constant from a wavenumber, or the reverse |
| Isotope effects | Predict ν̃ for the deuterated or 13C-substituted species |
| Anharmonicity | Fundamental and overtone positions; dissociation energy estimate |
| Band structure | Explain P, Q, R branches and the missing central line |
| Activity rules | Which modes are infrared or Raman active; mutual exclusion |
| Mode counting | 3N − 6 versus 3N − 5 for a stated molecule |
Use this as a checklist of sub-topics rather than a forecast. How many questions appear and how they are weighted is set by the official notification for the session you are sitting, and that document is the only reliable source.
Every calculation above needs the same four numbers. Planck's constant, the speed of light, the Boltzmann constant and the atomic mass unit appear in all five worked examples, and one mistyped exponent ruins the answer. The Scientific Constants tool lists them with their exact accepted values and units, so you can copy them correctly instead of trusting memory.
Open the Scientific Constants Tool →Preparing for CSIR-NET, GATE, IIT-JAM or CUET-PG? ABC Chemistry runs dedicated competitive-exam batches at its coaching centre and online for students across India — details at abcchemistry.in.