Concept

Isotope substitution — where it appears

Replacing an atom by a heavier isotope, which changes every mass-dependent quantity and leaves the potential energy surface untouched. It is the standard way to locate an atom in a rotational structure determination and to assign a vibrational mode, since the shifts follow from masses alone.

Named by 10 essays across one field — each of them below, with the objects they name alongside it.

H₂O with 1→D: what each mode is made of. H₂O with 1→D. Each row is one distinct frequency and each column one internal coordinate; the bar is the share of the motion in that coordinate. A mode whose largest share reaches nine tenths is a motion of one bond or one angle and is named for it. 3 of 3 here are.

When a mode becomes a bond stretch

Water's two stretching modes are each exactly half in one O–H bond and half in the other, which is why neither of them belongs to a bond. Change one hydrogen to deuterium and the same force field at the same geometry gives two modes that are 99.5 and 99.7 per cent in a single bond each. Nothing about the bonding changed; a mass did.

spectra · Normal mode
Three parameters, two numbers, and a curve of answers. seven structures of formaldehyde, every one of which reproduces the measured rotational constants A and B exactly. The C=O length runs from 1 to 1.3 ångström, the C–H length from 1.56 down to 0.95, and the HCH angle from 74.58 to 164.47 degrees. The third constant is not a third number: for a planar molecule it is fixed by the other two, and it comes out at 1.14 for every member.

Three numbers is not a structure

Formaldehyde's rotational spectrum gives three constants, of which a planar molecule's are only two independent numbers, and its structure has three parameters. Seven structures are computed here that reproduce A and B to the last digit the solver carries: the C=O length runs from 1.000 to 1.300 ångström, the C–H length from 1.557 down to 0.952, and the HCH angle from 74.6 degrees to 164.5.

spectra · Rotation
Every sign, lost. Formaldehyde in its own principal axes. Open circles are the atoms where they are; filled ones are where Kraitchman's equations put them, from the change in the three moments when each atom in turn is made heavier. The two agree to 7.6e-8 ångström — the equations are an identity for a rigid structure — but they return the square of each coordinate, so the two hydrogens at b = ±0.9348 both come back at +0.9348 and land on the same point.

The coordinate an isotope reports

Kraitchman's equations return an atom's position from the change in the moments when that atom alone is made heavier, and for a rigid structure they are an identity — formaldehyde's four atoms come back to a part in ten million. What they return is the square of each coordinate, so both hydrogens at b = ±0.9348 come back at +0.9348; every out-of-plane coordinate comes back imaginary at a moment error of one part in a hundred thousand; and the famous error cancellation, measured at a factor of thirteen, still leaves the answer two and a half times worse than a direct fit.

spectra · Rotation
The factor of thirteen was generous. An invented error model against the one computed from a force field. Its invented mismatch of five per cent made the correlation between the parent and the substituted species worth a factor of 12; the computed mismatch of 35 per cent makes it worth 2.5. The substitution structure went from being 2.5 times worse than a direct fit to 26, and at the computed size of the correction its worst coordinate is out by 10.7 per cent.

The correction that was invented

A standard error model puts the zero-point error in a rotational constant at a few tenths of a per cent, shared between the three moments by invented weights, with the parent and its deuterated form differing by five. Computed from a force field it is 1.88 per cent, one of the three shares is negative, and the mismatch is 35 — so the cancellation the substitution method rests on is worth a factor of 2.5 and not thirteen.

spectra · Rotation
The coefficient each principal axis needs. The dimensionless coefficient the molecule-averaged expression requires, evaluated on each principal axis separately rather than on the mean of the three. Across the twelve axes the positive ones span a factor of 2.33, against the 1.36 the molecule-averaged version spans — and one of them is negative, which no positive constant can be.

The axis that goes the other way

One dimensionless coefficient turned a zero-point correction into an expression a spectroscopist could evaluate, and the three principal axes were averaged over to get it. Split by axis it gets worse, not better — the coefficients span 2.3 where the molecule-averaged ones span 1.36 — and water's smallest moment does not grow at all. It shrinks.

spectra · Rotation
The frame H₂O → HDO turns. H₂O → HDO: the parent's principal axes and the daughter's, drawn on the same nuclei. Two of the three turn by 21.12° and the third does not move at all, because it is the normal to a plane no substitution can tilt. A per-axis comparison between these two species is comparing moments about lines this far apart, which is a rotation of the frame rather than a correction to a number.

Two moments about two different lines

Asked axis by axis, the substitution method's near-cancellation gives numbers as large as 163 per cent. The arithmetic is the smaller half of the answer. A principal axis is an eigenvector of a tensor built from the masses, so one deuterium turns water's frame by 21.12° — and the two moments being compared are not moments about the same line.

spectra · Rotation
The two terms in a vibrationally averaged rotational constant. A rotational constant averages 1/r², not r², so the moment it reports carries +2⟨Δr⟩/rₑ from the anharmonicity and −3⟨Δr²⟩/rₑ² from the harmonic spread. The two have opposite signs in every molecule here, and the anharmonic one — which is exactly zero in any symmetric well and therefore absent from every harmonic force field — is larger by a factor of 1.94 to 2.58.

The term a harmonic field cannot produce

The usual zero-point correction to a moment of inertia comes from a harmonic force field, which contains the mean square displacement and nothing else. A rotational constant does not average that. It averages one over r squared, whose leading correction is the mean displacement — zero in any symmetric well — and which enters with the opposite sign and about twice the size.

spectra · Rotation
A Morse curve's αₑ falls 4 to 15 per cent short of the measurement. The vibration–rotation constant αₑ of each diatomic, averaged over the states of a Morse curve built from its measured ωₑ, ωₑxₑ and rₑ, as a fraction of the tabulated value. H³⁵Cl: 0.27747 against 0.3072 cm⁻¹, ×0.903; D³⁵Cl: 0.10240 against 0.1133 cm⁻¹, ×0.904; ¹²C¹⁶O: 0.01674 against 0.0175 cm⁻¹, ×0.957; H¹⁹F: 0.68164 against 0.798 cm⁻¹, ×0.854. Every one is short, and HCl and DCl — one potential with two masses on it — are short by the same fraction.

The cubic a Morse curve guesses

The two terms in a vibrationally averaged rotational constant were computed on Morse curves built from measured constants, and the vibration–rotation constant αₑ is the measurement that tests them. In all four molecules the Morse curve's αₑ is short, by four to fifteen per cent. The averaging is not the error — it matches the closed form to four parts in ten thousand. The curve's cubic is, and the measurement asks for more of exactly the term a harmonic field cannot produce.

spectra · Rotation
One isotopologue converges and two run away. Methane's force field refitted with the bend–bend constant restored, to the light molecule alone and to both isotopologues. The first converges to an ordinary field with a C–H constant of 5.42. The second walks along the direction no frequency can see until two constants reach minus ninety-seven thousand, equal to four figures and opposite in effect. Adding data made it worse.

Adding data made it worse

A projection leaves a residual spread across four starting points, and an earlier essay attributed it to the eigenvalue problem rather than to the coordinate set — an attribution rather than a measurement. Refitting with two isotopologues was the proposed test. Run, it produces a runaway to minus ninety-seven thousand that one isotopologue alone does not, and the reason is that the direction it walks along is flat at every mass.

spectra · Normal mode
Every exact fit lies on one closed curve. The one-parameter family of force fields that reproduce all six of boron trifluoride's frequencies exactly, drawn in the B–F stretch constant and the stretch–bend coupling. It is a closed curve: the stretch constant runs from 3.94 to 12.11 millidyne per ångström and the coupling from -0.36 to 6.71. The four fits from four starting points sit together in one small patch of it.

The residual was a loop

Projecting away boron trifluoride's flat direction left four fits disagreeing by a tenth in the B–F stretch constant, and the disagreement was put down to the search or to the arithmetic. It is neither. Two frequencies cannot fix three constants in one symmetry block, so the exact fits form a closed curve along which the stretch constant runs from 3.9 to 12.1, and the four fits are four points on a short arc of it.

spectra · Normal mode

Named alongside it

The objects these essays reach for when they reach for this one.

ConventionModel limitMoment of inertiaRotational constantZero-point energyLeast-squaresDegeneracyForce constantIsotopologueNormal modeUnderdeterminationApproximation

All concepts