Concept

Force constant — where it appears

The second derivative of the energy with respect to a displacement, which sets a vibrational frequency. Usually it is fitted to observed frequencies rather than computed, and the fit is what limits what the frequencies can then be used to prove.

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

H₂O: 3 distinct modes. The displacement of every atom in 3 normal modes of H₂O, drawn from the eigenvectors of the mass-weighted Hessian. Under each is the internal coordinate with the largest share of the motion and how large that share is; where no coordinate holds nine tenths, the mode is not a motion of one bond or one angle and is labelled so.

Normal modes are not bond stretches

Water has two stretching frequencies and two O–H bonds, and it is almost irresistible to pair them off. Computed, each mode is exactly half in one bond and half in the other, and neither frequency belongs to a bond at all.

spectra · Normal mode
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
¹¹BF₃: what each mode is made of. ¹¹BF₃. 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. 1 of 4 here are.

Group frequencies, and where they stop

A carbonyl band sits near 1,700 wavenumbers in every ketone anybody looks at, and that regularity is real. Computed for a set of small molecules, four of twenty-one distinct frequencies belong to a single internal coordinate — and the four are exactly the coordinates symmetry leaves alone.

spectra · Normal mode
Donation strengthens the C–O bond, back-donation weakens it. Two interactions computed as two-level problems: 0.39 of an electron donated out of the ligand's σ orbital, worth 0.06 on the C–O bond order, and 0.5 donated back into π*, worth -0.4. Beside them, six measured stretching frequencies: five isoelectronic species differing only in charge, and free CO.

Back-bonding is two interactions

A carbon monoxide molecule bound to a metal donates from an orbital that is slightly antibonding and accepts into one that is strongly antibonding, so the two halves of the bonding move its stretching frequency in opposite directions. Five isoelectronic complexes differing only in charge settle which wins — and one of them stretches above free CO.

applied · Electron count
Every force field the spectrum permits. 5 force fields for H₂O, each one refitted with rr:O:H,H held at the value on the axis, and every one of them reproducing all of the molecule's own frequencies. The heavy line is the error each makes on the isotopologue it was never fitted to. The spectrum cannot choose between these fields; the isotopologue can.

The force field is not in the spectrum

Water has three vibrational frequencies and its force field has four constants, so the spectrum cannot determine the field. A whole family of quite different fields reproduces all three frequencies exactly, and only a second isotopologue can tell them apart.

wrong · Normal mode
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
methane, substituted 5 ways. CH₄, CH₃D, CH₂D₂, CHD₃, CD₄ — the same molecule with the same force constants, differing only in which nuclei are heavy. Every column is a consequence of which of the parent's operations survive the mass pattern: the surviving set is a subgroup, it is named from its own conjugacy classes, and the vibrations are reduced in it. The allowed and computed columns come from group theory and from the mass-weighted Hessian respectively, and they agree in every row.

A spectrum that changes when only a mass does

Methane has four distinct vibrational frequencies and two infrared bands. Replace two of its hydrogens with deuterium and it has nine and eight — with every force constant identical, every nucleus where it was, and the potential energy surface unchanged.

spectra · Spectrum
benzene: what alternation costs and gains. The π energy of benzene as its bonds are alternated by β(1 ± δ), the elastic cost of alternating them at a stated stiffness, and the sum. The π curve falls quadratically in δ, which is measured here and is the whole difference between a molecule that distorts and one that does not. The stiffness is a stated parameter, so where the minimum falls is not a claim; which power of δ each curve goes as is.

The hexagon is the frame's doing

Benzene's delocalisation energy is quoted as 2β and read as the reason its bonds are equal. Let the bonds alternate and the π energy goes down, not up — at every ring size, for 4n+2 as much as for 4n. The π electrons are not what keeps the hexagon regular; the σ frame is, and the margin between them is uncomfortably thin.

beyond · Delocalisation
H₂O: 3 distinct modes. The displacement of every atom in 3 normal modes of H₂O, drawn from the eigenvectors of the mass-weighted Hessian. Under each is the internal coordinate with the largest share of the motion and how large that share is; where no coordinate holds nine tenths, the mode is not a motion of one bond or one angle and is labelled so.

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
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
boron trifluoride at 1454 cm⁻¹, shared out four ways. One mode of boron trifluoride, with each internal coordinate's share of it computed by four conventions. Each convention is a bar in every group; the groups are the coordinates. A number quoted for this band without saying which bar it is has not said much.

How much of a band is a bond stretch

Boron trifluoride's 1454 cm⁻¹ band is 36.7 per cent B–F stretch, or 49.9, or 97.9, depending on which of four standard ways of sharing a mode out among internal coordinates is used. All four are defensible, all four sum to one, and the spread between them is sixty-one percentage points on one band of one molecule.

spectra · Normal mode
How many coordinates, and how many motions. For each of six molecules, one square per internal coordinate the valence set carries — bond stretches, angle bends and an out-of-plane wag where there is one — with a rule drawn at the number of vibrational degrees of freedom. three of them have more coordinates than motions, and which ones is decided by shape rather than by size: ammonia's three angles are independent and boron trifluoride's are not, and the only difference is that one is flat.

More coordinates than motions

Methane has ten internal coordinates and nine ways to vibrate, boron trifluoride seven and six, formaldehyde seven and six. The excess is not bookkeeping: it is a combination of coordinates that describes no displacement of any atom, and adding twenty-five units of force constant along it moves every frequency by five parts in a hundred million.

symmetry · Normal mode
The shallowest slope, and the steepest. For three molecules, the softest and stiffest vibrations, the species each belongs to, and what a unit distortion along each costs. The cost is not taken from the frequency: it is computed by resolving the distortion onto the normal coordinates in the mass-weighted metric and adding up ω²q². That it comes back as the frequency is the identity the whole comparison rests on — the species decomposition uses only the coordinates and the group, and the cost uses only the masses and the force constants, and the two have to agree on this case before they can be asked to disagree on any other.

The coordinate it was already soft along

A distorted structure can be resolved into the symmetry species of its ideal group, and that resolution cannot say which coordinate a molecule fell down by itself and which one something outside pushed it along. The force field answers that, and the two halves agree on the case where they must — a unit distortion along a normal mode costs exactly that mode's frequency, to a part in a million, computed from masses and force constants by one side and from coordinates and characters by the other.

symmetry · Point group
The 6-ring at 111.5°: the alternating form and what a search finds. The alternating ring — every atom displaced above or below the plane in turn — with its dihedral angles and its two strain terms, beside the lowest-torsion member a search constrained only by bond angles and closure returns. Both satisfy every geometric constraint; only one of them is staggered.

The explanation with the wrong sign

Two methyl groups on one carbon make a ring easier to close, by up to eleven thousand-fold, and the textbook reason is that they compress the ring's internal angle. Computed from a standard bending term, that compression helps a three-ring and a four-ring and hinders every ring from five up — predicting a slowing of 0.754-fold for the five-ring measured to speed up 250-fold. The account with the right sign is about rotations rather than angles, and it has a ceiling of 36.5 that the measurement is already above.

shape · Strain
The ceiling rises and the measurements fall, so they cross. The largest acceleration the rotamer account can produce, against the ring being closed, with the measured gem-dimethyl accelerations on the same axis. Closing a bigger ring means freezing more rotations, so the ceiling rises steeply; the measurements go the other way. The five-membered ring's 250-fold acceleration is above its own ceiling of 36.5 and the six-membered ring's tenfold one is far below its 121 — so the account is refused at one size and sufficient at the next.

A ceiling that rises where the measurements fall

The rotamer account of the gem-dimethyl effect has a largest possible acceleration, which looks like a limitation. It is a prediction: the ceiling is a closed form in the number of rotations a closure freezes, it rises steeply with ring size, and the measurements fall — so the account is refuted for the five-membered ring and more than sufficient for the six.

shape · Strain
10 of 55 force constants that no spectrum can see. The map from methane's 55 independent force constants to its Cartesian Hessian, as a spectrum: 45 directions the frequencies respond to and 10 they do not, out of 55. The null block is exactly the size the redundancy count predicts — 10 for 1 redundancy on 10 coordinates — and the two are computed by different routes, one a rank and one a closed form. A force field quoted to four figures is quoted along 45 directions that were measured and 10 that were chosen.

Ten directions no frequency can see

A redundant force field has one obvious flat direction. There are ten: methane has fifty-five independent force constants and a ten-dimensional subspace of them that no spectrum can touch. The literature's repair is to project — and the metric a vibrational analysis is naturally written in cannot define the projection at all, because a redundancy is a null vector of it.

symmetry · Normal mode
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
How much of each of methane's force constants the spectrum fixes. Each of methane's 55 force constants, grouped by kind, with the squared length of its projection onto the determined subspace. One means the constant is fixed on its own. All 4 stretching constants are; not one bending constant is, at 0.5357 each. The fractions sum to 45, which is the rank of the map from constants to the Hessian, and that sum is an arithmetic check rather than a result.

The forty-five that are fixed

Methane's fifty-five-dimensional space of force constants has ten directions no frequency can see, and the useful thing to report is the forty-five that are fixed. It is a table: every stretching constant is fixed on its own, no bending constant is, and projecting a fitted field onto the determined subspace takes four fits that span 0.88 mdyn per ångström down to four that span 0.0027.

spectra · Normal mode
One expression, four molecules, a factor of eleven. The computed zero-point correction to each molecule's moment of inertia against the expression A²·Σ(1/ν)/I, which has no fitted quantity in it. The dashed line is the mean dimensionless coefficient, 0.6343; the four points lie within 18 per cent of it, on corrections that span a factor of 11.4. The expression is evaluated from a moment of inertia and a list of wavenumbers, which is what a spectroscopist has before doing anything.

An expression for what was a warning

The zero-point correction to a moment of inertia comes out at 1.88 per cent where a few tenths had been assumed, and heavy molecules are safer. Written out, the correction is a mean curvature times the sum of reciprocal wavenumbers over the moment — one line, evaluated from things a spectroscopist has before starting. One coefficient serves four molecules whose corrections span a factor of eleven.

spectra · Rotation
The two molecules with a direction no spectrum can see. Methane's fifty-five independent force constants and boron trifluoride's twenty-eight, split into the combinations a spectrum determines and the ones it cannot touch. Both molecules have a redundant coordinate set — methane's six angles at a tetrahedral centre are five coordinates' worth, and boron trifluoride's three angles at a planar centre are two — so both have a flat space, and the smaller molecule's is the larger share of its field.

The second molecule with a blind spot

Methane's fifty-five force constants have ten directions no spectrum can touch, and the forty-five that are fixed were worth a table. Boron trifluoride is the only other molecule here with a redundant coordinate set, and it is not the same case: eight kinds of constant rather than five, a quarter of its field invisible rather than a fifth, and an out-of-plane coordinate the redundancy cannot reach.

spectra · Normal mode
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.

Normal modeVibrational modesValence force fieldHarmonic approximationModel limitLeast-squaresDegeneracyInternal coordinateConventionIsotopologueMode compositionUnderdetermination

All concepts