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 15 essays across 2 fields — 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
No substitution that keeps the symmetry can tell the two fields apart. Ten isotopologues of boron trifluoride, real and hypothetical, each computed under the fitted force field and under the second field a boron-10 substitution cannot exclude, with the largest difference between their frequencies on a logarithmic axis. The six that keep all three fluorines alike give the same spectrum under both fields to a millionth of a wavenumber, the precision of the calculation. The four that replace one or two fluorines differ by four to seven tenths of a wavenumber.

The third mass that cannot help

Boron trifluoride's E′ block has three force constants and two frequencies, so its exact fits form a closed curve, and a boron-10 substitution cuts that curve to two fields. The obvious remedy is a third isotopologue, and it was named as the first data that could pick one of the two. It is not, if the third substitution keeps the molecule's symmetry: every such isotopologue gives the same spectrum under both fields to a millionth of a wavenumber, because its sum rule is one plane written in two reciprocal masses and its product rule is Teller and Redlich's. Only a substitution that breaks the symmetry adds an equation — and it separates the two fields by half a wavenumber.

spectra · Normal mode
Ammonia's E verdicts across the fields its spectrum allows. The concentration of each distortion on one occurrence of ammonia's E species — the share of its E projection on the favoured vibration — for every field in the walk, with the stretch–bend constant held and the rest refitted to both isotopologues. The band is where the fit cannot tell the fields apart and the dotted line is the published fit. Across the band the angle under the B row reads 0.502 to 0.579, the angle at least cost 0.840 to 0.923, the bond at least cost 0.890 to 0.957.

The worst case was the loosest constant

A symmetry label's verdict on a distortion is read off mode shapes, and mode shapes are set by a fitted force field. Walked along the direction its own frequencies see least, ammonia's field moves the verdict on its bond angle from 0.840 to 0.923 without the fit getting measurably worse — and boron trifluoride, the census's worst case at 0.573, has a whole curve of fields that reproduce its spectrum exactly, along which the same verdict runs from a coin to 0.862.

symmetry · Point group
Carbon monoxide's fourth constant picks the fit that turns over. Along the Morse curves with a linearly bent range that hold carbon monoxide's measured ωₑ and ωₑxₑ, the ωₑyₑ read from the lowest four levels, against the bend λ. The two members that also have the measured αₑ are marked: the one at λ = −0.221 dissociates and predicts 0.0182 cm⁻¹; the one at λ = −0.453 turns over at 2.10 bond lengths and predicts 0.0110. The measured value, read from the published level expansion the same way, is 0.0110 — the second, to 0.0 per cent.

The fourth constant chose a curve, not a depth

Read from its lowest four levels the way every potential here is read, carbon monoxide's measured ωₑyₑ is 0.01101 cm⁻¹ — and of the two bent-range Morse curves that fit its first three constants, the one that predicts 0.01100 is the one that never dissociates. Let the range bend and stay bounded and every molecule, the hydrides included, is fitted to all four constants by a curve that does dissociate. The depths those curves carry are 6 to 31 per cent too deep, further from the measurement than the plain Morse curve's on all four, and hydrogen chloride and deuterium chloride — one potential — are fitted by two different ones.

spectra · Rotation
Ammonia's two Coriolis constants along the family of fields the spectrum cannot tell apart. ζ₃ and ζ₄ of NH₃ on every member of ammonia's family of force fields, each refitted to the frequencies of NH₃ and ND₃ with its stretch–bend constant held, against that constant. The shaded band is where the fit is no worse than its own noise, -1.000 to -0.634. Across it ζ₄ runs from -0.807 to -0.641 and ζ₃ from 0.418 to 0.584, while their sum stays at -0.222 on every member.

Two Coriolis constants are one number

Ammonia's frequencies leave its force field free along a band, and across that band a symmetry label's verdict on the bond angle changes. The measurement proposed to settle it was the Coriolis constant of each E vibration, which reads mode shapes directly — and it would settle it, sixteen times more sharply than the frequencies do. But the two constants are not two measurements. Their sum is fixed by the geometry and the masses, ζ₃ + ζ₄ = I∥/2I⊥ − 1 (the moments of inertia about the axis and across it), on every field there is, so the spectrum's Coriolis constants carry exactly one number about the field between them.

symmetry · Point group
Around boron trifluoride's loop of exact fits, ζ₃ takes every value — each one twice. The Coriolis constant ζ₃ of ¹¹BF₃'s E′ stretch on every field of the closed loop that reproduces its frequencies exactly, against the B–F stretch constant. It runs over its whole possible range, −1 to +1. The level of the fitted field, 0.791, meets the loop at exactly two points — B–F constants 7.425 and 3.966 — which are the fitted field and the partner field a ¹⁰B spectrum cannot tell from it.

The Coriolis constant was an isotope shift

Boron trifluoride's E′ block has a closed loop of force fields that fit its spectrum exactly, a ¹⁰B spectrum cuts the loop to two, and no third isotopologue that keeps the symmetry can choose between them. The Coriolis constant depends on mode shapes rather than frequencies, so it looked like the measurement that could. Around the loop it takes every value from −1 to +1 — and under the two fields it is the same to nine decimals, because for a planar top it is exactly the boron's share of the motion, which is what a boron isotope shift measures.

spectra · Normal mode

Named alongside it

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

Force constantModel limitNormal modeConventionRotational constantDegeneracyLeast-squaresMoment of inertiaUnderdeterminationValence force fieldZero-point energyIrreducible representations

All concepts