On a high, dry plateau, dozens of radio antennas point at a dark cloud between the stars. Among the frequencies they collect is one at 115.271GHz: carbon monoxide molecules, at a few kelvin, turning from their first rotational level to the lowest. From that single number a chemist obtains the length of the C–O bond to a fraction of a picometre; from the strengths of the next lines, the temperature of the cloud. The rotations and vibrations of molecules are the subject of this chapter: their levels, from the rotor and the oscillator of Chapter 1; the selection rules, from Chapter 5; and what spectra measure — bond lengths, force constants, dissociation energies.
You already know
Chapter 1 gave the rigid rotor, EJ=J(J+1)ℏ2/2I with degeneracy2J+1, and the harmonic oscillator, Ev=(v+21)hν, for a diatomic of reduced massμ. Chapter 5 stated the infrared and Raman selection rules. The Year 1 volume measured infrared spectra in wavenumbers and defined the polarisability of a molecule.
Antennas of a radio interferometer on a high plateau: they record, among many others, the rotational lines of carbon monoxide from interstellar clouds. Photograph ESO/B. Tafreshi (twanight.org), CC BY 4.0.
6.1 Rotational spectra
Definition 6.1(Rotational constant, isotopologue)
The rotational constant of a linear molecule is B=h/(8π2cI), in cm−1 (or B=h/8π2I in hertz), so that its rotational levels are F(J)=BJ(J+1) in wavenumbers. Molecules that differ only in the isotopes of their atoms are isotopologues (X12X2212CX16X2216O, X13X2213CX16X2216O); in the Born–Oppenheimer approximation they share the same bond lengths and force constants.
Definition 6.2(Gross and specific selection rules)
A gross selection rule states which property a molecule must have to show a kind of spectrum at all; a specific selection rule states which changes of quantum number are allowed.
Proposition 6.3(Pure rotational spectra)
A molecule has a pure rotational (microwave) spectrum only if it has a permanent dipole moment; for a linear molecule the specific rule is ΔJ=±1. Absorption lines lie at
ν~=F(J+1)−F(J)=2B(J+1),J=0,1,2,…,
equally spaced by 2B.
Proof. The interaction with light is through the dipole moment: a molecule with no permanent dipole has no oscillating dipole when it rotates. The rule ΔJ=±1 comes from the vanishing of ∫YJ′M′∗cosθYJM unless J′=J±1 (the dipole component along z transforms like Y1,0), admitted. Then B(J+1)(J+2)−BJ(J+1)=2B(J+1). ∎
Method 6.4(Bond length from a rotational line)
From the line J+1←J at frequency ν, deduce B=ν/2(J+1) (in hertz).
Compute the moment of inertia I=h/8π2B.
With the reduced mass from isotopic masses, r=I/μ.
Example 6.5(Carbon monoxide)
The line J=1←0 of X12X2212CX16X2216O is at 115271.20MHz: B=57635.6MHz=1.92252cm−1, I=1.4560×10−46kgm2. With μ=12×15.99491/27.99491=6.85621u, r0=113.09pm. This r0 averages the bond over its zero-point vibration and is slightly longer than the equilibrium length re (Section 6.2).
Definition 6.6(Centrifugal distortion constant)
A rotating bond stretches; the levels become F(J)=BJ(J+1)−DJ2(J+1)2, where D is the centrifugal distortion constant, and the lines 2B(J+1)−4D(J+1)3 close up slowly at high J.
Proposition 6.7(The strongest line)
If line intensities follow the population of the lower level, (2J+1)e−hcBJ(J+1)/kT, the most populated level is near
Jmax≈2hcBkT−21.
Proof. Treat J as continuous and differentiate: dJd[(2J+1)e−aJ(J+1)]=[2−a(2J+1)2]e−aJ(J+1) with a=hcB/kT; it vanishes at (2J+1)2=2/a. ∎
The rotational absorption spectrum of carbon monoxide, lines spaced by 2B=3.845cm−1, intensities taken as the population of the lower level, at 30K (blue) and 300K (red, shifted slightly for visibility), each temperature scaled to its strongest line. At 30K the strongest line is J=3←2; at 300K, J=8←7.
6.2 Vibrations of real bonds
A real bond does not obey Hooke’s law far from equilibrium: it softens when stretched and breaks.
Definition 6.8(Morse potential)
The Morse potential is V(r)=De[1−e−β(r−re)]2: harmonic near re (force constant2β2De), and tending to the dissociation limit De at large r.
Proposition 6.9(Levels of the Morse oscillator)
In wavenumbers the vibrational levels of a Morse oscillator are
a ladder whose rungs close up and end at vmax=⌊ω~e/2ω~exe−21⌋.
Proof.Admitted at this level.∎
The solution of the Morse equation is treated in more advanced courses; the levels are checked numerically in this chapter’s figure data.
Definition 6.10(Anharmonicity constant, fundamental, overtone, hot band)
ω~exe is the anharmonicity constant. The fundamental transition is v=1←0; an overtone is v=n←0 with n≥2, weak because the gross rule Δv=±1 of the harmonic oscillator is only slightly broken; a hot band starts from an excited level, v=2←1, and grows with temperature.
The spectroscopic dissociation energy of a molecule is the depth of its potential, De, measured from the minimum, or D0=De−G(0), measured from the lowest level — per molecule at 0K, unlike the bond dissociation enthalpy of the Year 2 volume, a molar enthalpy at 298K.
Corollary 6.12(De and D0 of a Morse oscillator)
De=ω~e2/4ω~exe and D0=De−(21ω~e−41ω~exe).
Proof. The first is the relation of the proposition; the second is G(0). ∎
Proposition 6.13(Birge–Sponer extrapolation)
The gaps ΔG(v+21)=G(v+1)−G(v)=ω~e−2ω~exe(v+1) decrease linearly with v; summing them until they vanish gives D0, the area under the plot of ΔG against v.
Proof. Subtract the two expressions of G. D0 is the energy from v=0 to the limit, the sum of all the gaps. ∎
Example 6.14(Hydrogen chloride)
For HX35X2235Cl, ω~e=2990.95cm−1 and ω~exe=52.82cm−1: the fundamental is at ω~e−2ω~exe=2885.3cm−1 and the first overtone at 2ω~e−6ω~exe=5665.0cm−1, slightly less than twice the fundamental. The Morse model predicts De=42342cm−1 and D0=40860cm−1; the true D0, from the enthalpies of formation of H, Cl and HCl at 0K, is 35760cm−1 (4.43eV). The real potential flattens out faster than a Morse curve fitted at the bottom: Birge–Sponer extrapolations from the lowest levels overestimate dissociation energies.
Left: the Morse potential of HCl built from its spectroscopic constants (solid), the harmonic parabola with the same curvature (dashed), and every third vibrational level, drawn between its turning points; the levels close up towards De. Right: the Birge–Sponer plot of the gaps; the shaded area is D0.
6.3 Rovibrational bands
A vibrational transition of a gas is accompanied by rotational changes, ΔJ=±1 for a diatomic in a 1Σ state.
Definition 6.15(P, Q and R branches, band origin)
In a rovibrational band, the lines with ΔJ=−1 form the P branch, those with ΔJ=+1 the R branch, and those with ΔJ=0, when allowed, the Q branch. The band originν~0 is the vibrational energy difference without rotation.
Theorem 6.16(Combination differences)
With rotational constantsB0 and B1 in the lower and upper vibrational levels, R(J)=ν~0+B1(J+1)(J+2)−B0J(J+1) and P(J)=ν~0+B1(J−1)J−B0J(J+1), and
R(J−1)−P(J+1)=4B0(J+21),R(J)−P(J)=4B1(J+21).
Proof. The first two follow from E=G(v)+BvJ(J+1) for both levels. R(J−1) and P(J+1) end on the same upper level J, from lower levels J−1 and J+1: their difference is B0[(J+1)(J+2)−(J−1)J]=B0(4J+2). Likewise R(J) and P(J) start from the same lower level and end on J+1 and J−1. ∎
Method 6.17(Analysing a rovibrational band)
Locate the gap: for a 1Σ diatomic there is no Q branch, and ν~0 lies in the gap, between R(0) and P(1).
Number the lines outwards from the gap: R(0),R(1),… and P(1),P(2),….
Use the combination differences to get B0 and B1 separately, then αe=B0−B1 and Be=B0+21αe.
The fundamental band of gaseous HCl at 300K, computed from its spectroscopic constants: P and R branches either side of the gap at the band origin. Each line is a doublet, HX35X2235Cl and the weaker, slightly lower HX37X2237Cl (natural abundances 76 % and 24 %). The R lines close up, the P lines spread out, because B1<B0.
Example 6.18(The constants of HCl from its band)
With R(0)=2905.58cm−1 and P(2)=2842.94cm−1, one finds 6B0=62.64cm−1, so B0=10.440cm−1. With R(1)=2925.23 and P(1)=2864.43cm−1, B1=10.133cm−1. So αe=0.307cm−1 and Be=10.593cm−1, the values the band was built from: the method recovers them exactly.
6.4 Raman spectroscopy
Definition 6.19(Raman and Rayleigh scattering)
When light of frequency ν0 crosses a sample, most of the scattered light keeps the frequency ν0: Rayleigh scattering. A small fraction is shifted by the rotational or vibrational frequencies of the molecules: Raman scattering. Lines at lower frequency are Stokes lines (the molecule gained energy), at higher frequency anti-Stokes lines (it lost energy).
Proposition 6.20(Classical origin of the Raman effect)
If the polarisability of a molecule vibrating at νvib varies as α=α0+α1cos(2πνvibt), the dipole induced by a field E0cos(2πν0t) oscillates at ν0 and at ν0±νvib. A vibration is Raman active only if it changes the polarisability.
Proof.μ=αE=α0E0cos(2πν0t)+21α1E0[cos2π(ν0+νvib)t+cos2π(ν0−νvib)t], by cosacosb=21[cos(a+b)+cos(a−b)]. Each oscillating dipole radiates at its own frequency; without α1 only the Rayleigh term remains. ∎
Proposition 6.21(Rotational Raman lines)
A linear molecule with an anisotropic polarisability shows rotational Raman lines with ΔJ=±2, displaced from the exciting line by
∣Δν~∣=F(J+2)−F(J)=4B(J+23)=6B,10B,14B,…
— the first line at 6B, then a spacing of 4B.
Proof. The polarisability ellipsoid of a rotating linear molecule returns to the same orientation twice per turn, so it modulates the scattering at twice the rotation frequency: ΔJ=±2 (the quantum rule, admitted). Then B(J+2)(J+3)−BJ(J+1)=B(4J+6). ∎
Example 6.22(Nitrogen, seen by Raman)
NX2 has no dipole and no infrared or microwave spectrum, but its polarisability is anisotropic: its rotational Raman lines start at 6B0=11.94cm−1 from the exciting line and are spaced by 7.96cm−1. Their intensities alternate 2:1 between even and odd J: the two X14X2214N nuclei are identical, and nuclear-spin statistics (Chapter 11) give even-J levels twice the weight of odd ones.
Left: the rotational Raman spectrum of nitrogen computed at 300K: Stokes lines at positive shifts (right), anti-Stokes lines at negative shifts (left, slightly weaker), with the 2:1 alternation of intensities; the grey line at zero marks the much stronger Rayleigh line. Right: energy scheme of Rayleigh and Stokes scattering through a virtual level (not a real state of the molecule).
6.5 Polyatomic molecules
Definition 6.23(Types of rotors)
A molecule with three equal principal moments of inertia is a spherical top (CHX4, SFX6); with two equal, a symmetric top (NHX3, CHClX3, CX6HX6); with all three different, an asymmetric top (HX2O, most molecules); linear molecules have one moment zero.
A spherical top has no dipole and no rotational spectrum; a symmetric top has levels BJ(J+1)+(A−B)K2, where K counts the angular momentum about its axis, and lines that still fall at 2B(J+1) since ΔK=0; an asymmetric top has an irregular spectrum, fitted by computer. The vibrations of polyatomic molecules are their normal modes, labelled and selected by Chapter 5: the infrared spectrum shows the modes that transform like x, y, z, each with its own rotational structure.
In the lab— A gas cell in an infrared spectrometer
Gaseous HCl is admitted into a 10cm cell with windows transparent in the infrared, inside a Fourier-transform spectrometer. At 1cm−1 resolution the P and R branches appear as single lines; at 0.25cm−1 each splits into its HX35X2235Cl/HX37X2237Cl doublet, 2cm−1 apart. Hydrogen chloride is toxic and corrosive: the cell is filled on a vacuum line in a fume hood.
Safety
Hydrogen chloride gas: a gas under pressure, corrosive to skin, eyes and the respiratory tract, toxic if inhaled. Handled only in closed apparatus or a fume hood.
History— Raman, 1928
In 1928 C. V. Raman and K. S. Krishnan observed, in sunlight filtered through a violet filter and scattered by liquids, a faint light of different colour that a crossed green filter could isolate. Raman received the 1930 Nobel Prize in Physics; laser sources turned his effect into a routine analytical technique fifty years later.
6.6 Exercises
Exercise 6.1★
From B0=1.9225cm−1 for X12X2212CX16X2216O, compute the moment of inertia and the bond length r0.
Solution
Solution of Exercise 6.1.
I=h/(8π2cB)=6.62607×10−34/(8π2×2.99792×1010×1.9225)=1.4561×10−46kgm2. μ=6.85621u=1.13850×10−26kg, so r0=I/μ=113.09pm.
Exercise 6.2★
With B0=10.440cm−1, give the positions of the first four lines of the rotational spectrum of HX35X2235Cl and the J of the strongest line at 300K.
Solution
Solution of Exercise 6.2.
2B0(J+1): 20.88, 41.76, 62.64, 83.52cm−1. Jmax=kT/2hcB−21=208.5/20.88−0.5=2.7: the level J=3 is the most populated, and the strongest line is 4←3 at 83.5cm−1.
Exercise 6.3★
Which of HX2, HCl, NX2, COX2, CHX4, SFX6 and HX2O show a pure rotational absorption spectrum? A rotational Raman spectrum?
Solution
Solution of Exercise 6.3.
Microwave absorption needs a permanent dipole: HCl and HX2O only. Rotational Raman needs an anisotropic polarisability: HX2, HCl, NX2, COX2, HX2O; not the spherical topsCHX4 and SFX6.
Exercise 6.4★
What fraction of HCl molecules is in v=1 at 300 and 1000K (fundamental 2885.3cm−1)? When does the hot band become visible?
Solution
Solution of Exercise 6.4.
N1/N0=e−hcν~/kT: at 300K, e−13.84=9.8×10−7; at 1000K, e−4.15=0.016. The hot band becomes noticeable (a per cent of the fundamental) only near 1000K.
Exercise 6.5★★
The fundamental and first overtone of HX35X2235Cl lie at 2885.3 and 5665.0cm−1. Deduce ω~e and ω~exe.
Solution
Solution of Exercise 6.5.
ω~e−2ω~exe=2885.3 and 2ω~e−6ω~exe=5665.0: subtracting twice the first from the second, −2ω~exe=−105.6, so ω~exe=52.8cm−1 and ω~e=2990.9cm−1.
Exercise 6.6★★
With the constants of the previous exercise, compute De and D0 of the Morse model, and compare with the true D0=35760cm−1. Convert both D0 into kJ/mol.
Solution
Solution of Exercise 6.6.
De=2990.952/(4×52.82)=42341cm−1; G(0)=1495.47−13.20=1482.3cm−1; D0=40859cm−1=488.8kJ/mol, against the true 35760cm−1=427.8kJ/mol: the Morse model overestimates by 14 %.
Exercise 6.7★★
Four lines of the HX35X2235Cl fundamental are R(0)=2905.58, R(1)=2925.23, P(1)=2864.43 and P(2)=2842.94cm−1. Find B0, B1, αe and the band origin.
For NX2, B0=1.9896cm−1. Give the shifts of the first three Stokes rotational Raman lines, and explain the alternation of their intensities.
Solution
Solution of Exercise 6.9.
Shifts 4B0(J+23): 11.94, 19.90, 27.85cm−1 (from J=0, 1, 2). The two X14X2214N nuclei are identical bosons of spin 1: the even-J levels have six nuclear-spin states, the odd-J three, so lines from even J are twice as strong as their neighbours.
Exercise 6.10★★★
The first three rotational lines of X12X2212CX16X2216O are at 115271.2018, 230538.0000 and 345795.9899MHz. Show that a rigid rotor does not fit them, and determine B and the centrifugal distortion constantD.
Solution
Solution of Exercise 6.10.
A rigid rotor would put the lines at ν1, 2ν1, 3ν1: 230542.40 and 345813.61MHz, above the measured ones by 4.4 and 17.6MHz. With νJ+1←J=2B(J+1)−4D(J+1)3: ν1=2B−4D, ν3=6B−108D, so 3ν1−ν3=96D: D=0.1835MHz and B=(ν1+4D)/2=57635.97MHz. Check: ν2=4B−32D=230538.0MHz.
Exercise 6.11★★★
Show that the Birge–Sponer sum of the gaps ΔG(v+21) of a Morse oscillator from v=0 to the last bound level is close to De−G(0), and evaluate the sum for HCl.
Solution
Solution of Exercise 6.11.
ΔG(v+21)=ω~e−2ω~exe(v+1) vanishes near v=ω~e/2ω~exe−1; the sum of an arithmetic progression of n terms from ω~e−2ω~exe down to almost 0 is about nω~e/2≈ω~e2/4ω~exe−ω~e/2≈De−G(0). For HCl, the 28 gaps (v=0 to 27) add up to 40857cm−1, against De−G(0)=40859cm−1.
Exercise 6.12★★★
The X16X2216O nucleus has spin zero, so the levels of X16X2216OX12X2212CX16X2216O with odd J do not exist. What is the spacing of its rotational Raman lines, and the shift of the first one, in terms of B? Why does COX2 have no microwave spectrum?
Solution
Solution of Exercise 6.12.
With only even J, the transitions J→J+2 start at J=0,2,4,…: shifts 6B,14B,22B,…, spacing 8B. COX2 has a centre of symmetry and no permanent dipole: no microwave spectrum.
6.7 Problem: Listening to a Cold Cloud
Problem 6.1
Weekend problem — the bond length of carbon monoxide from its first rotational line, the temperature of an interstellar cloud, the carbon-13 isotopologue, and the infrared band
Data: X12X2212CX16X2216O lines J=1←0 at 115.2712GHz and 2←1 at 230.5380GHz; X13X2213CX16X2216O line 1←0 at 110.2014GHz; masses X12X2212C 12 (exactly), X13X2213C 13.00335, X16X2216O15.99491u; natural abundance of X13X2213C 1.1 %; for X12X2212CX16X2216O, ω~e=2169.81cm−1, ω~exe=13.288cm−1, and the equilibrium length re=112.83pm. h=6.62607×10−34Js, u=1.66054×10−27kg, c=2.99792×1010cm/s, hc/k=1.43878cmK.
Part I — The bond.
Compute B0 in hertz and in cm−1.
Compute the reduced mass of X12X2212CX16X2216O in kilograms.
Compute the moment of inertia.
Deduce r0.
Compare with re. Why is r0 longer?
Predict the 2←1 line for a rigid rotor. Comment on the measured value.
Part II — The temperature of the cloud.
Give the energies, in cm−1 and in kelvin (E/k), of the levels J=0, 1, 2.
Write the ratio of the populations N2/N1 at temperature T.
The observations give N2/N1=0.85 (data of the problem). Find T.
Which level is the most populated at that temperature?
Why can such a cloud not be seen in the infrared vibrational band?
Why is the 1←0 line the most used tracer of cold gas?
Which molecule, X12X2212CX16X2216O or X13X2213CX16X2216O, has the lower band origin?
State the result: the bond length r0 of carbon monoxide obtained from one radio line.
Solution
Solution of Problem 6.1.
1.B0=ν/2=57.6356GHz; B0/c=1.92252cm−1. 2.μ=12×15.99491/27.99491=6.85621u=1.13850×10−26kg. 3.I=h/8π2B0=1.4560×10−46kgm2. 4.r0=I/μ=113.09pm. 5.r0 exceeds re=112.83pm by 0.26pm: B0 averages 1/r2 over the zero-point vibration in an anharmonic well, which spends more time at long distances. 6.2ν1←0=230.5424GHz; the measured line is 4.4MHz lower: centrifugal distortion. 7.0, 3.845 and 11.535cm−1, i.e. E/k=0, 5.53 and 16.60K. 8.N2/N1=35exp[−(16.60−5.53)K/T]. 9.ln(35/0.85)=0.673=11.06/T: T=16K. 10.Jmax=16.4/(2×1.43878×1.9225)−0.5=1.2: J=1. 11. A vibrational quantum corresponds to hcν~/k≈3080K: at 16K no molecule is vibrationally excited, so the cloud emits no vibrational light. 12. It needs only 5.5K of excitation, is excited even in the coldest gas, and CO is abundant and has a dipole moment. 13.μ=13.00335×15.99491/28.99826=7.17241u. 14.B∝1/μ: 115.2712×6.85621/7.17241=110.1893GHz. 15. The prediction is 12MHz (0.011 %) below the measured line: B0=Be−21αe, and the vibration–rotation term αe scales with a different power of μ; the equilibrium structure is the same, the zero-point average is not. 16.98.9/1.1=90. 17. The X12X2212CX16X2216O line is saturated (optically thick): it no longer grows with the amount of gas. 18. The rare X13X2213CX16X2216O line stays proportional to the amount of gas, the X12X2212CX16X2216O line gives the temperature: together they measure both. 19.ν~0=2169.81−2×13.288=2143.24cm−1, 4.666µm. 20.2ω~e−6ω~exe=4259.89cm−1. 21.G(0)=21ω~e−41ω~exe=1081.58cm−1. 22.R(0)−P(1)=2(B0+B1)≈4B≈7.7cm−1. 23.X13X2213CX16X2216O: heavier reduced mass, lower frequency. 24.r0(CO)=113.1pm, from a single radio frequency.