Concept

Reduced mass — where it appears

The effective mass of a two-body vibration, m₁m₂/(m₁+m₂). It fixes the frequency for a given force constant, which is why isotopic substitution moves every frequency by an amount arithmetic can predict exactly.

Named by 19 essays across 3 fields — each of them below, with the objects they name alongside it.

H₂O and its D isotopologue. Every frequency of H₂O joined to the frequency the same force field gives when every H is replaced by D, with the ratio on each join. The force constants were not refitted and could not be: they do not depend on mass. The product of all the ratios is fixed by the masses and the moments of inertia alone, and is checked against that identity while this figure is drawn.

The isotope shift is arithmetic

Replace hydrogen with deuterium and every frequency drops. The usual rule says by a factor of the square root of two — and of water's three modes, not one of them does that. What is exact is a different identity, and the force constants cancel out of it.

spectra · Normal mode

A structure out of a spectrum. Two rotational constants and two bond lengths, three times over, from three pairs of carbonyl sulfide isotopologues. Above them, the same inversion run on moments computed from a known structure, which returns it to twelve figures. The measured pairs disagree with one another by 7.4 milliangstrom, which is the difference between a ground-state average and an equilibrium geometry.

A bond length out of a spectrum

A linear triatomic has two bond lengths and one moment of inertia, so one measurement cannot determine it. Substituting an isotope gives a second measurement on the same structure, and two equations in two unknowns have a solution — which comes out differently depending on which isotope is used.

spectra · Rotation

6 fitted force fields. Every valence force field fitted here, ordered by the size of its bond stretching constant, with the stretching frequencies of the molecule beside it. The two orders are not the same, which is the whole of what separates a force constant from a frequency. The last two columns say how many constants were fitted to how many observed frequencies, and a field with as many of the first as the molecule has distinct frequencies fits exactly and reports nothing.

The frequency is not the bond strength

Sulfur dioxide's S–O force constant is larger than water's O–H constant, and its stretching bands sit at a third of the frequency. A vibrational frequency carries a mass as well as a force, and the two cannot be separated by looking at a spectrum.

wrong · Normal mode

The six motions that are not modes

Three N minus six is quoted in every textbook and the six are almost never computed. Diagonalising a mass-weighted Hessian gives six eigenvalues at arithmetic noise, and projecting their vectors onto the three translations and three rotations written down from the geometry alone accounts for every one of them.

spectra · Normal mode

A microwave constant predicted from an infrared one. For four diatomics: the rotational constant and the stretching frequency, both measured, and the centrifugal distortion constant predicted from them by 4B³/ω² — then the constant a microwave spectroscopist fits to the line positions. The prediction and the fit agree within a few per cent across three orders of magnitude in the quantity, and nothing connects the two measurements except the assumption that the bond stretching under rotation is the same bond that vibrates.

The rotor that stretches

A rigid rotor's lines are evenly spaced, and a real molecule's are not — it pulls itself apart as it spins. How much is not a fitting parameter: it follows from the stretching frequency by one relation, and the prediction agrees with the measured constant to a few per cent across four molecules spanning three orders of magnitude.

spectra · Rotation

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 $water, drawn with its nuclei the size they are. The molecule with a disc round each nucleus whose radius is the computed root-mean-square displacement of that nucleus in the vibrational ground state. The hydrogens' discs are a substantial fraction of the bond length, and this is at no temperature at all.$

The atoms are not at the points

Every structure in this collection is a set of points, and every bond angle it argues about is a property of that set. Computed from the fitted force fields it already has, water's hydrogens are 0.094 Å from where they are drawn and its bond angle has a spread of 8.9 degrees — larger than the difference between 104.5 and the tetrahedral value that half the essays here are about. This is at no temperature at all.

shape · Approximation

H³⁵Cl: the well, its states and their averages. The Morse potential built from H³⁵Cl's measured vibrational constants, with the lowest four states drawn at their computed energies and the average separation of each marked. Every average lies to the right of the minimum, because the well is not symmetric — and they move outward as the state rises.

The bond length that depends on the isotope

Hydrogen chloride and deuterium chloride have the same potential energy curve, and their measured equilibrium lengths agree to three hundredths of a milliångström. Their average bond lengths differ by 4.34 mÅ — a hundred times more — because a lighter atom explores more of a well that is not symmetric.

wrong · Approximation

What boron trifluoride's bands are strong in, and what they move. Every infrared-active mode of boron trifluoride, with its band strength and the root-mean-square displacement of its atoms in the zero point, each scaled to its own largest. The two do not order the modes the same way — the rank correlation between them is 0.2 — and the strongest band belongs to the mode at 719 cm⁻¹, in which 89.93 per cent of the motion is the lightest atom's. Every mode moves the same weighted amount of mass, exactly, so that is not what separates them either.

The mode that moves least radiates most

Boron trifluoride's strongest infrared band is the one in which the fluorines barely move: ninety per cent of the motion belongs to the boron, which is a fifth of the molecule's mass. The mode that moves the most mass is nine and a half times weaker. Across five molecules the rank correlation between band strength and how far the atoms actually go runs from +1 to −0.66, and every normal mode carries exactly the same weighted motion by construction.

shape · Dipole

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

NH₃: four states in a well the molecule does not sit at the bottom of. The umbrella coordinate of NH₃ — the signed distance of the N atom from the plane of its three H atoms — with the quartic well that has its minima at the measured 0.3816 ångström and its barrier at the quoted 2020 wavenumbers. The lowest 4 states are drawn at their computed energies. The lowest sits 587 wavenumbers above the bottom, which is 29.1 per cent of the way up the barrier, so the state is far from the harmonic bottom that a drawing of a pyramid implies.

A barrier is not what a splitting measures

Ammonia's inversion barrier is quoted everywhere as 2020 wavenumbers. Put that number into the simplest double well its own measured geometry allows and the ground-state splitting comes out at 1.3508 against a measured 0.7935, and the excited one at 68.37 against 35.81. Both are too large because a splitting is an area under a barrier and a height is only one of its two dimensions.

shape · Inversion

One per cent on the barrier is 3.6 per cent on the splitting. The ground inversion splitting of NH₃'s quartic well against the barrier height, both logarithmic, with the geometry and the reduced mass held at their measured values. The curve is visibly bent: its local slope is -3.56 at the published barrier and steepens either side, so a power law is a tangent to it rather than a description of it. The measured 0.7935 wavenumbers is reached at 2330, which is 15.3 per cent above the quoted 2020 — so a splitting wrong by a factor of 1.70 is a barrier wrong by a sixth. The same derivative read the other way is what makes a barrier quoted to ten per cent useless for predicting a splitting.

The exponent that runs both ways

How hard does a splitting depend on a barrier? Locally, as the power −3.5628 — and the local slope runs from −2.53 to −6.15 across the same sweep, so there is no power law. What is exact is stranger: rescaling the equation forces the mass exponent to be one below the barrier's and the geometry exponent to be twice the mass's, so the model's three sensitivities are one number and the arithmetic reproduces both identities to six decimals.

shape · Inversion

The mass is worth a factor of 1.6, and the other two 1e+4 and 9e+4. Ammonia's umbrella well, with each of phosphine's three differences substituted into it one at a time and then all together. The reduced mass is 11 per cent larger and costs a factor of 1.61. The pyramid is 2.01 times taller and costs 9.6e+3; the barrier is 6.1 times higher and costs 9.4e+4. Phosphine's own splitting is below what the arithmetic resolves, so it is drawn at that bound.

It was never the mass

Phosphine does not invert, and the reason given is that phosphorus is heavier than nitrogen. Three things about phosphine differ from ammonia. Substituting each into ammonia's own well one at a time, the reduced mass costs a factor of 1.61, the pyramid height costs 9,600 and the barrier 94,000 — and the mass is the smallest of the three by four orders of magnitude.

shape · Inversion

Three constructions of one number, spanning a factor of 5.6. The reduced mass of NH₃'s umbrella coordinate under each construction, against position along the coordinate. Two of the three are constants and the third is not: if the bonds are held at their measured length, the ligands must slide outward as the apex descends, and their radial motion adds to the mass. It runs from 2.4866 at the plane to 2.9874 at the pyramid — 20 per cent — and the coordinate itself stops existing at one bond length, which is where the curve ends.

The mass nobody chose

Every one-dimensional treatment of ammonia's inversion needs a mass, and the measurement does not supply one. Three constructions are defensible and they give 1.35075, 0.94420 and 0.0000055 wavenumbers. The honest one is not a constant at all, and it moves the answer thirty per cent towards the measurement — which means the usual choice is the wrong one.

shape · Inversion

The orderings in use move the splitting by 0.47 per cent between them. The change in NH₃'s ground inversion splitting under each ordering of the kinetic operator, relative to BenDaniel–Duke, with the bond-conserving mass throughout. The bars are exact solves and the ticks are first-order perturbation theory. The five span 0.469 per cent, from −0.407 to 0.060; the change from a constant mass to the bond-conserving one, in the same well and box, is 43.1 per cent, 92 times as large.

An ordering worth half a per cent

A mass that varies along a coordinate has no unique quantum kinetic energy, and the choice among the Hermitian orderings in use was the one thing left that could undo a forty-three per cent correction to ammonia's splitting. It cannot. The five orderings anybody uses span 0.47 per cent between them, a ninety-second of the correction, and the family only reaches the measurement at exponents three times larger than any of them.

shape · Inversion

Every prediction of the isotope ratio overshoots, and the mass decides nothing. The ratio of NH₃'s ground inversion splitting to ND₃'s, predicted six ways, against the measured 14.94. At the published barrier the usual mass gives 15.77 and the bond-conserving mass 17.79. With a quartic fitted to NH₃'s splitting they give 19.26 and 19.00; with a well whose shape is fitted to both of NH₃'s lines, 17.62 and 18.48. The bond-conserving mass is nearer the measurement in one of the three pairs and further in two, and the difference within any pair is smaller than the distance of either from the measurement.

Deuterium cannot tell the masses apart

A reduced mass built by holding ammonia's bonds rigid predicts a deuterated molecule differently from any constant mass, and ND₃'s splitting is measured. Run as a test, it cannot choose. Every well and every mass needs a barrier for ND₃ several per cent lower than for NH₃, every prediction of the isotope ratio from a well fitted to NH₃ overshoots by eighteen to twenty-nine per cent, and the two masses differ by less than either misses — in opposite directions in the two wells.

shape · Inversion

The term two of the four molecules have and two do not. The radial coefficient of the bond-conserving reduced mass for the four isotopologues: the sum of the ligand masses, and what is left after the asymmetric correction. Three ligands at a hundred and twenty degrees on a circle whose radius changes as the apex descends move their own horizontal centre of mass outward — unless the three masses are equal. A frame that does not translate has to subtract that motion, and the amount is half the sum of the squared mass differences over the total mass. It is zero at both ends of the series and the same number in the middle.

The two that are not on the line

Ammonia and its fully deuterated twin are two points, and two points cannot show a curve. Putting the partly deuterated molecules between them needs a term neither symmetric one has — three ligands of unequal mass move their own centre of mass sideways as the apex descends — and it moves the prediction by half a per cent, which is what a whole change of mass construction was worth.

shape · Inversion

The covariant operator is BenDaniel–Duke plus this. The difference between the Laplace–Beltrami operator — the one a one-dimensional manifold with metric μ(x) distinguishes, carried across to the flat measure by the unitary map ψ ↦ μ^(¼)ψ — and the BenDaniel–Duke ordering, divided by the function it was applied to, at forty-one positions inside the molecule's own range. The curve drawn through the marks is the two-function fit every ordering is a combination of, and it passes through them to a part in ten million.

The ordering a manifold picks

A position-dependent mass leaves the kinetic energy with no unique quantum form, and an earlier sweep of the five orderings in use found half a per cent between them. A one-dimensional reduction is a one-dimensional manifold, a manifold has a distinguished Laplacian, and carrying it to the flat measure lands on exactly one of those five — not the one with no extra potential, and not the one anybody reaches for.

shape · Inversion

Named alongside it

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

Model limit Isotopologue Inversion splitting Convention Double well Tunnelling Zero-point energy Born–Oppenheimer separation Harmonic approximation Normal mode Bond length Force constant

All concepts