Concept

Least-squares — where it appears

Fitting parameters by minimising the sum of squared residuals. The parameters that come out depend on which quantity is fitted and over what range, and a fit that reproduces its own data well may extrapolate to something quite wrong.

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

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
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
Two humps and a dip, which is not what a trend looks like. The measured enthalpy of hydration of the first transition series, in kilojoules per mole, against a straight line fitted through it. The measurements do not fall on the line: they rise and dip at manganese, rise and dip again at zinc. Both dips are at configurations with no ligand field stabilisation — d⁵ high spin and d¹⁰ — and the line alone accounts for only 73 per cent of the variation.

The double hump and what removes it

The hydration enthalpies of the first transition series do not lie on a line — they rise, dip at manganese, rise and dip again at zinc. Subtract the ligand field stabilisation computed from the same model that describes their spectra and what is left is a line, with the one fitted parameter landing inside the range a spectrum measures.

applied · Ligand field
The electrons repel each other less inside the complex. Seven chromium(III) complexes. Δ is the first band; B is solved from the second in closed form; β is B against the free ion's 918 cm⁻¹, which is measured on the gaseous ion. Every β is below one — the electrons in a complex repel each other less than the same electrons in the free ion, because they have more room. The two orderings are different: fluoride splits least and reduces the repulsion least, cyanide does both most, and the middle of the two series is not the same middle.

The electrons repel less inside the complex

Two measured bands determine two parameters in closed form, and one of them is the repulsion between the d electrons — which comes out below the free ion's value for every complex, by between two and forty-seven per cent. The ligands that split most are not the ligands that reduce the repulsion most, so a complex is characterised by two numbers rather than one.

applied · Ligand field
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 far benzene is from its own group, against how far it has been pushed. The measure — the distance to the nearest structure with the ideal group — against the size of the distortion, for two ways of distorting the same molecule. Each curve is smooth and quadratic; the markers on it are where the point-group search changes its verdict, which happens at one step and says nothing about the steps either side of it.

How much symmetry is left

A point group is a verdict and every real structure fails it. Two distortions of benzene that are two per cent apart in how far they sit from D6h need tolerances 1.88 times apart before the search will call either of them D6h — and one tolerance, applied to six molecules, admits amounts of asymmetry differing by a factor of fifty-six.

symmetry · Point group
Flat rings: the angles, the torsions and what is measured. For each ring from three to eight, held flat: its interior angle, how far that is from tetrahedral, the angle strain that follows, the torsional strain of having every bond eclipsed, and the measured strain energy. The five-ring is the row the essay is about — its angles are almost ideal and it is strained.

The strain that is not in the angles

Cyclopentane's flat bond angles are 108°, a degree and a half from tetrahedral, and its angle strain computed from a standard bending constant is 0.4 kJ mol⁻¹. Its measured strain is twenty-six. The missing sixty kilojoules are torsional — one ethane barrier for every bond in the ring, which no account built on bond angles mentions.

shape · Strain
The same sample, fitted over four temperature ranges. A pair coupled at -50 cm⁻¹, its susceptibility computed exactly, fitted to a Curie–Weiss law over four ranges. The moment and the Weiss temperature the fit reports both depend on which range was used, and the quality of the fit does not warn about it.

The moment a fit invents

One coupled pair of spins, its susceptibility computed exactly, fitted to a Curie–Weiss law over four temperature ranges. The moments reported are 2.471, 2.535, 2.566 and 3.590 Bohr magnetons, and the Weiss temperatures −51, −78, −78 and −292 K — from one sample, measured perfectly, with three of the four fits agreeing with their own data to better than a part in three hundred.

applied · Magnetism
seven sets of parameters, one spectrum. Benzene's π levels from seven sets of the three Hückel parameters, each fitted to reproduce the two measured ionisation energies exactly. The two occupied levels sit at the same energy in every column, because that is what was fitted. The empty level moves from -3.15 to 4.69 electronvolts across the family, and the resonance integral from -3.05 to -7.66.

One spectrum, a line of models

Hückel theory has three parameters and benzene's photoelectron spectrum supplies two numbers, so the fit has a curve of solutions rather than a point. Along it the resonance integral runs from −3.05 to −7.66 electronvolts, the empty level moves by eight, the terminal-to-central bond order ratio in butadiene goes from 2.000 to 2.671 — and the delocalisation energy is 6.100 electronvolts in every member.

bonding · Models
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
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
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 coupling a fit reports, and the coupling the sample has. Exact susceptibilities of Heisenberg chains of two, four, six and eight spins, every one of them coupled at -50 cm⁻¹, each fitted with the two-spin expression over 80–600 K. The two-spin sample returns its own coupling exactly; every longer chain returns one too large, by more the longer it is, up to 20.9 per cent. Every one of those fits has an R² above 0.99, so nothing in the fit reports that anything is wrong.

The model is what is fitted

Fit a pair of coupled spins with the two-spin expression and it hands back the coupling exactly, from any temperature range. Fit a chain of eight with the same expression and it hands back −66.7 where the sample has −50, with a residual of 0.998 and a g factor of 1.973 — three numbers of which only the last says anything is wrong, and it is the one nobody looks at.

applied · Magnetism
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
A correction computed at 3 bohr and used everywhere. Three binding curves for H₂⁺ in 2 Gaussians a centre: uncorrected, properly counterpoise corrected at every separation, and corrected once at 3 bohr with that value subtracted throughout. The frozen curve is the uncorrected one shifted down by a constant, so its minimum sits at 2.2270 bohr — exactly where the uncorrected minimum is, and 4.2 millibohr from where the full correction puts it. The depth moves and the structure does not.

A correction computed at one length

The counterpoise correction is expensive, so it is evaluated once at a reference geometry and subtracted across a whole potential surface. A constant does not move a minimum — so a frozen correction returns the uncorrected bond length exactly, at every reference geometry and in every basis, and everything the correction does to a structure is the part that has just been thrown away.

orbitals · Basis
Two molecules, one eigenvalue, and 1.27 eV between them. Six alternant hydrocarbons placed by the Hückel eigenvalue of their highest occupied level — computed by diagonalising each molecule's own adjacency matrix — against the measured first π ionisation energy. Ethene and benzene share an eigenvalue of exactly 1 and their measurements differ by 1.27 eV; butadiene and naphthalene share 0.618 and differ by 0.94. A model that reads only the eigenvalue is a function of it, so it must give each pair one answer, and the two vertical pairs are the whole of its error.

A parameter that never finds a value

Show a two-parameter model six measurements instead of two and it stops being underdetermined and starts being wrong. Adding the third parameter improves the fit by one part in eighty, moves the resonance integral by a factor of four, and never finds a best value at all — because nine tenths of the error is a term the model does not have.

bonding · Models
Four bonds, or one a₁ and three t₂. The same four occupied orbitals written in two bases. On the left each orbital sits on one bond; on the right one is shared over all four hydrogens and three follow the Cartesian directions. The transformation between them is orthogonal, so the density is unchanged.

The answer a search is most likely to give

If a cage has ten equally good localised descriptions, why does the literature agree about its picture? The hoped-for answer was that one basin is very large. Counting four hundred and eighty starting points says it is not: two independent searches agree one time in ten, and the description they most often return is not the best one.

beyond · Multicentre
What the second channel buys, and what it cannot. Five measured splittings against three models. Adding the π overlap takes the error from 1747 to 1402 cm⁻¹, and telling the model which ligand is an acceptor takes it to 961 — so the fact about occupation is worth more than twice the integral. The two halides are the pair that fixes which model is which, and no single one of the three gets both them and cyanide right.

The integral that cannot count electrons

Adding the π channel to a ligand-field splitting means one more overlap integral over the same two orbitals at the same distance. It removes a fifth of the error. Telling the model which ligand is a π acceptor — one word per ligand, quoted rather than computed — removes forty-five per cent, because an overlap cannot know whether the orbital it reaches is full or empty.

applied · Ligand field
A control that ranked better than the mechanism. Rank correlations against the additivity shortfall, over eight ion pairs. The overlap of the two closed shells ranks at 0.8571 — but the cation's formal charge, which cannot be a mechanism, ranks at 0.9524, so the set is confounded: its eight pairs split four and four by charge and everything else rises with it. Held fixed within a charge group the overlap still ranks at 0.80 — and so does the softness, at -1.00. Four pairs cannot separate two candidates.

A control that outranked the mechanism

Ruling polarisation out left one candidate, and the closed-shell overlap ranks at 0.857 against the additivity shortfall — which looked like the answer until the control was read. The cation's formal charge, which cannot be a mechanism, ranks at 0.9524. Eight pairs split four and four by charge cannot separate anything, and within a charge group two candidates both rank perfectly.

orbitals · Contour
How much of a curve each extra parameter has left to work with. The singular values of the design matrix for a susceptibility curve, for two, three and four parameters fitted to the same data, on a logarithmic scale. With four they run 12.411, 2.026, 0.149, 0.025 — a span of 500 — so one per cent data fix the first two to under 3 per cent and the last to 37. Each value is what is left of the measurement after the directions above it have taken their share, so a short bar is not a hard parameter but an absent one.

How many parameters a curve is worth

A susceptibility curve routinely carries four fitted parameters and the question of whether it can support them is never asked. It has an arithmetic answer: the four directions the fit sees span a factor of five hundred, so one per cent data fix the first two to under three per cent and the last to thirty-seven — and forty points reaching two kelvin are worth more than sixteen thousand starting at twenty.

applied · Magnetism
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
The slope goes to a half, and a window fit stops short of it. The local slope of the alternation against the reduced temperature, between each neighbouring pair of points, on a ring of 40 at K = 1.6. It rises monotonically from 0.4115 to 0.5053 as the transition is approached, crossing a half at about a part in a thousand of the reduced temperature. The fitted 0.44 is the average of the left-hand end of this curve; the exponent is one half, which is what a free energy analytic in one order parameter is obliged to give.

The exponent was the window's

A fit over the last decade before a distortion vanishes gives an exponent of 0.44, and running it on larger rings should say whether the number belongs to the transition or to a forty-site ring. It belongs to neither. The local slope runs to 0.5020 as the transition is approached, and 0.44 is what a fit over that particular decade returns — on every ring size and every stiffness, because the whole curve is one curve.

wrong · Metal
Two curves that cross, and two that do not. Every pair of parameters that reproduces one measured number, for three numbers, with the true system marked at Δ = 12 and t⊥ = 0.06. The lower band's shape is a function of the ratio of the two, so its curve is a straight line through the origin; the excess gap is a function of the coupling squared over the separation, so its curve bends. The two cross at one point. The upper band's shape draws a line almost on top of the first, because it is a function of the same ratio — a second measurement lying along the first fixes nothing the first had not already fixed.

Two ways of being second order

A band's shape and a band's gap are both second order in the coupling that mixes two bands, which sounds like a reason to measure only one of them. They are second order in different ways — one goes as the square of the ratio and the other as the square over the separation — and that single difference of one power is what turns a curve of possible answers into a point.

solids · Bands in a solid
What each of these predictors cannot see. Four predictors built elsewhere in this collection, each audited by its ties: the share of the variation in what was measured that the predictor demonstrably cannot account for, because two systems it assigns the same number to were measured to differ by that much. the spin-only moment leaves 23 per cent; the VSEPR angle leaves 92 per cent; the highest occupied Hückel eigenvalue leaves 41 per cent. The control has no ties at all and the instrument returns nothing for it, which is what it must do. No fitting anywhere: a tie is a claim a model of that form cannot escape.

Two systems a model cannot tell apart

Two molecules that share a Hückel eigenvalue but were measured to differ bound a whole family of models at 0.456 eV. That is a reusable instrument, and chemistry has plenty of predictors of exactly the same shape. Turned on themselves: VSEPR cannot account for 92 per cent of the variation in the four angles it predicts, and no fitting is involved anywhere.

bonding · Models
A length that keeps growing, and one that stops. The fitted decay length of the ring-current response, against the number of rings. The bare acene's runs 1.397, 1.669, 1.896, 2.104, 2.302 — up by a factor of 1.65 and still climbing — while its own gap falls from 0.590 to 0.1102. With a gap held open the same measurement gives 0.678, 0.642, 0.636, 0.636, 0.639, which has stopped moving by the third molecule. There is a magnetic reach, and an acene is too nearly gapless to have one.

A reach that has no length

A ring current's response to a neighbouring ring falls with distance, which invites asking for the length. Every acene computed gives a longer one — 1.397, 1.669, 1.896, 2.104, 2.302 rings — because the gap that would set the length is closing at the same time. Give the same molecule a gap that stays open and the number settles at 0.636 by the third one and does not move.

symmetry · Aromaticity
Eight wells, and where each one puts its pair. The total energy of each pair against separation, with the measured distance marked on every curve. The wells are deep and their minima are in the right region — tenths of an ångström from the measurements — which is what makes the comparison worth making. What they are not is closer to the measurements than the sum of two tabulated radii, and that is the result.

A size a confound cannot supply

A rank correlation of 0.857 was beaten by a control that cannot be a mechanism, so a size is the next thing to ask for: does a closed-shell repulsion of the computed magnitude displace two ions by the tenths of an ångström the additive radii are wrong by. It does not. The balance of a Madelung attraction against six computed repulsions predicts six separations to 0.242 ångström where adding two tabulated radii predicts them to 0.183, and the displacement it produces ranks at 0.14 against the shortfall it was proposed to explain.

orbitals · Contour
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
The four things a susceptibility curve measures, in order. The four directions in the parameter space, best fixed first, each written as the product of powers it is. The Jacobian is logarithmic, so a direction is a set of exponents and a combination is a product — which is why the answer can be printed. The best-determined is g · J^-0.33, fixed to 0.08 per cent by a curve measured to one per cent; the worst is tip · rho^-0.40, fixed to 40. Neither is one parameter's own axis.

The product a curve measures

A susceptibility curve's fourth parameter is undetermined, and the question is which combination the free direction actually is. It is a product of powers, because the Jacobian is logarithmic — and at the usual window it is the temperature-independent term divided by the 0.40 power of the monomer fraction, fixed to forty per cent, while the product one place up is fixed to 6.7. A paper could print that instead of four numbers.

applied · Magnetism
A fourth data point, and 3 negative electronegativities. Each element's electronegativity from three points on its energy curve, the cubic coefficient a fourth point adds, and what the fourth point leaves. The cubic coefficient is one sixth of a second difference of the ionisation series, so it is largest where that series has a kink — and the alkali metals, whose second electron comes out of a closed shell, are pushed to Li -7.91, Na -3.42, K -1.49 eV. The last column is where the two roots of the fixed-point equation collide, beyond which the atom has no solution at all.

Where a closed form stops being one

What happens when the quadratic energy is not enough? A cubic makes the equalisation condition a quadratic with two roots, and something has to choose between them. The choice is easy and the finding is somewhere else — a cubic fitted through the dication gives lithium an electronegativity of −7.907 eV, a capacity of 0.073 of an electron, and a molecule of two alkali metals no solution at all.

wrong · Electronegativity
Two lengths off the same chains, and only one of them has an exponent. The length an end's influence reaches into a chain of 320, against the gap the bulk has opened, over ten elastic constants and a factor of twenty in the gap. Fitted over the first twelve bonds — as a short-window fit does — the exponent is -0.476. Taken from the local decay rate extrapolated to a bond infinitely far from the end, it is -1.029, and every neighbouring pair of points gives between -1.06 and -0.91. The argument says −1. The two lines are the same ten profiles read two ways.

A decay that keeps slowing down

A healing length fitted over the first twelve bonds of each relaxed chain goes as the gap to the power −0.60, against an argument that says −1. A fitted exponent can belong to its window, and this one does — but not because the window is too short for the chain. The profile is not an exponential at all, so there is no single length for an exponent to be of.

solids · Peierls distortion
One half is a curve and the other is not. The two halves of a counterpoise correction for an unequal pair, against separation, on a logarithmic axis. The heavier centre's falls smoothly over two decades; the lighter centre's scatters over more than one decade between neighbouring points. It is not a rough function — it is a difference of two energies of order a hartree whose difference is a millionth, and the solver does not have seven figures to spare.

The half that cannot be computed

How many points does each half of a counterpoise correction need to interpolate? The natural expectation is two different numbers. The answer is that the question is not yet askable: the lighter centre's half is a difference of two energies agreeing to six figures, its second differences are seven per cent of its own value, and no interpolation of it means anything. The third thing worth checking — the symmetric-pair check — works perfectly.

orbitals · Basis
Two instruments, and they do not agree. Each of this collection's four predictors under both tests. The exact-tie instrument reports the share of the variation a predictor demonstrably cannot account for, and it ranks VSEPR worst and has nothing at all to say about the ring-strain control. The near-tie instrument reports how much steeper any model must be somewhere than the set is overall, and it ranks the control worst and cannot speak about VSEPR, whose predictor is a label rather than a number. Neither instrument is the general one, and a predictor that passes one has not been audited.

The pair that is not a tie

An exact tie bounds every model of a form with no fitting, and it is available only for predictors that read integers. A near tie bounds the model's derivative instead, and a pair the predictor orders the wrong way round refuses monotonicity outright. Run on four standard predictors, the new instrument ranks them in a different order from the old one — and its worst case is the control the old one could not see.

bonding · Models
What the gap alone would predict, and what was fitted. Each halide's π scale as a multiple of fluoride's: the value fitted to the spectrochemical series, against what the energy denominator alone gives with the metal orbital at the vacuum level — which is the weakest the denominator effect can be. It over-predicts at every ligand, and moving the metal level down makes it worse.

The gap that would have to be smaller

An angular overlap parameter is an overlap squared over an energy denominator, and the usual fit folds the denominator away. Put the measured ionisation energies back in and the denominator alone over-predicts the trend down the halide group at every metal level a donor permits — the smallest it can give is 1.67 against a fitted 1.43. So the overlap has to shrink down the group, which is the opposite of the usual expectation.

applied · Ligand field
6 rings fused two ways. Two catacondensed chains of 6 hexagons: the linear one, where every fusion continues the line, and the angular one, where the fusions alternate. They have the same formula and the same number of bonds and they are not the same graph. Ring centres are numbered; in the linear molecule two rings k steps apart have centres 1.7321k units apart and in the angular one they do not, which is the whole reason this pair can be asked the question.

Neither of the two separations

A linear acene cannot pose the question, because the number of fusions between two rings and the distance between their centres are the same variable there. Bending the molecule pulls them apart — and the response follows neither. Two pairs of rings the same distance apart differ by two thirds, and the larger one is at the greater distance.

symmetry · Aromaticity
Four ways of predicting the same six separations. How far each model's predicted separation is from the measured one, pair by pair. The best additive model misses by 0.0098 ångström on average, the tabulated radii by 0.0237, radii read off the ions' own densities by 0.1828, and the balance of a Madelung attraction against a computed repulsion by 0.2421.

The residue that is two numbers

Is a shortfall between a sum of radii and a measured separation a property of the pairs, or of the compromise a universal table makes? The eight separations are two complete two-by-two blocks, so what no assignment of radii can reproduce is not a residual at all — it is an alternating sum, computed by subtraction, and it comes to three hundredths of an ångström and six.

orbitals · Contour
One heteroatom, and how far the relaxation carries it. The change every ring's bonds undergo when one carbon of the end ring is given a site energy, on a chain of 12 fused hexagons with its gap held open. The upper curve measures each molecule's couplings from its own mean bond order, which is what this collection's relaxation has always done; it stops falling at 4.48e-5 and stays there. The lower curve measures both from the same mean and keeps falling to 7.89e-9. Nothing about the molecule differs between them.

The floor was in the bookkeeping

A heteroatom in one ring of a fused chain is a perturbation with a place, and the relaxation carries it along the molecule. Measured directly, the response stops falling after five rings and sits at four hundredths of a millionth for ever. That floor is not the molecule. It is the relaxation measuring each system's couplings from its own mean, and removing it recovers nine orders of magnitude.

bonding · Delocalisation
The end pair against the deepest pair, at each separation. For each separation, the response of the pair that touches an end divided by the response of the pair at the same separation sitting deepest in the molecule. A straight chain is below one at every separation and a zigzag is above one at every separation, so the end effect has opposite signs on the two shapes. The third chain — two straight arms meeting at one angular ring — is above two at three separations and below one at the fourth, which is a third behaviour and not an intermediate one.

An end effect with two signs

Neither of two separations accounts for the scatter in a fused ring system's response, and the natural guess is the end: pairs with more molecule outboard should behave differently from pairs at an edge. They do. In a straight chain an end pair responds a third less than an interior one, in a zigzag a quarter more, and in a chain of two straight arms meeting at one angular ring the anomaly is in the middle.

symmetry · Aromaticity
One turned fusion, moved along the chain. Chains of 9 rings differing in one integer: which fusion's direction is turned. Turned at the first or last fusion, that leaves one angular ring beside an end; anywhere between, it leaves two adjacent angular rings whose turns cancel. The widest ratio between two pairs at the same separation, against which fusion is turned. A straight chain gives 1.520 and every bent one gives more — from 3.052 to 4.600. The two ends of the curve are the one-ring members; every interior point is a two-ring step.

One integer, and everything it changes

Eight molecules with the same rings, the same carbons and the same graph distance between every pair, differing in which fusion's direction is turned — one angular ring when the turned fusion is at an end, two adjacent ones anywhere else. The scatter within a separation class runs from three to four and a half times, against a straight chain's one and a half — and the two ends of one molecule disagree by up to a factor of four.

symmetry · Aromaticity
A heteroatom on ring 5, and the profile in both directions. The response of each ring to a heteroatom placed on ring 5 of 12, on a logarithmic scale. Two rings share the peak, because an interior ring's outermost carbon belongs to two rings at once, and the profile falls away at the same rate on both sides — which is the answer wanted: the reach is the molecule's and not the end's. Every point is a mean over the six bonds of its ring.

The reach is the molecule's

Every measurement of the relaxation's reach so far puts the heteroatom on the end ring, because that gives the longest run to measure a decay over — and the response of a molecule to a perturbation at its end is not the response to one in its middle. Moved inward, the decay is the same in both directions and the same as the end's. The amplitude is not: it halves, and splits across two rings.

bonding · Delocalisation
The two overlaps, squared, at the measured bond lengths. For each chromium(III) donor: the σ overlap squared, a metal 3d(z²) against the donor's p(z), and the π overlap squared, a 3d(xz) against its p(x). Everything is computed — the radial functions from Slater's rules, the separation from the measured bond length, the integral by quadrature. Chloride's π overlap is 3.5 times fluoride's, which is the opposite of what the overlap argument required of it.

The overlap the model is not proportional to

Without a computed π overlap, the natural argument reasons about one instead: the denominator over-predicts the halide trend, so the overlap must shrink down the group to cancel part of it. Computed, it grows — 3.5 times from fluoride to chloride. And the fitted parameter changes sign across the series, which no ratio of squared overlaps can do.

applied · Ligand field
The exponent's sign is the ring's parity. The power the temperature-independent term carries in the free combination, for an open chain and for rings of five to ten spins. Every odd ring is negative and every even ring and the chain is positive — so the free product is ρ·χ_TIP raised to a power whose SIGN changes, which is a different combination rather than a shifted one. Nothing here is a near miss: the closest pair on either side of zero are +0.23 and −0.29.

The sign a frustrated ring changes

There is a sharper question than whether a low-temperature feature buys back a fourth parameter: does it change which combination is free? It does, and by a sign. Every odd ring of spins leaves free the monomer fraction times a negative power of the temperature-independent term, and every even ring a positive power, with no case in between.

applied · Magnetism
Adding one measurement, and six of fifteen change verdict. Each atom's chemical capacity from a cubic fitted through four electron counts and from a quartic fitted through five, on a logarithmic axis, with an unbounded capacity drawn at the right-hand margin. Five atoms go from a finite capacity to an unbounded one and one goes the other way — and the five are exactly the five largest the cubic reported. The seven that stay finite keep their order and change their values.

Six of fifteen change verdict

The chemical capacity is finite or infinite according to the sign of a coefficient that is one sixth of a second difference of three measurements, and the pair test asked how much of its reported ordering is the quantity and how much is the fit's resolution. Adding a fifth point changes the verdict on six of fifteen atoms — and the five it unbounds are the five it had largest.

wrong · Electronegativity
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
The one-bend family was one angular ring at its ends and two everywhere else. Filled circles: the worst spread within a separation class for the eight nine-ring chains built by bending one fusion, placed at the angular rings each actually has. The chains bent at the first and last fusions have one angular ring, next to an end; the six bent in between have two adjacent angular rings whose turns cancel. The lower line is a single angular ring at each interior position, running 3.05, 3.66, 3.46, 3.03, 3.46, 3.66, 3.05; the upper line is two adjacent ones, running 4.30, 4.60, 4.06, 4.06, 4.60, 4.30. Every published point lies on one line or the other.

The scatter counts angular rings

Eight nine-ring chains built by turning one fusion looked like one bend moved along a molecule. They were two families: one angular ring when the turned fusion is at an end, two adjacent angular rings everywhere else. Built from stated angular rings instead, twenty-seven chains show the ring-current scatter ignores which way a ring turns, does not grow with how many turn, and is mostly explained by one count per pair — each angular ring between two rings multiplies their response by 0.60, each one under them by 0.65.

bonding · Aromaticity

Named alongside it

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

Model limitConventionUnderdeterminationApproximationReference stateClosed formHückel theoryForce constantDegeneracyDelocalisationNormal modeOverlap integral

All concepts