Concept

Eigenvalue — where it appears

A number λ for which an operator returns λ times a vector it acts on. Every orbital energy, vibrational frequency and rotational level computed from a matrix is one, and the eigenvector beside it says what is moving.

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

Chains of 2, 4, 8, 16, 40: the levels crowd, the edges do not move. Every level of a chain of 2, 4, 8, 16, 40 sites, drawn at its computed energy on one axis. Adding sites adds levels inside a fixed interval rather than widening it, which is what makes the large system a band rather than a wider molecule.

A solid is a molecule that did not stop

Diagonalise a chain of two atoms, then four, then forty. Nothing new happens at any point, and by forty the levels are a band. The passage from molecule to solid is not a change of subject; it is the same matrix at a different size, and every step of it can be watched.

solids · Bands in a solid
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
⟨x²⟩ is the average number of neighbours. For each structure, the mean coordination counted off the edge list beside the mean of x² measured off the eigenvalues. They are equal by an identity about graphs, not by a limit — the full band width beside them obeys no such rule.

The width of a band is a count of neighbours

The mean of the squared level energies equals the average coordination — exactly, for any structure, with no limit taken and no periodicity assumed. It is the one statement in this field that is arithmetic rather than physics, and the usual textbook formula for band width is a special case of something weaker.

solids · Bands in a solid
The density of states of a chain of 2000. The 2000 levels of a linear chain, binned into 34 intervals across the band, with the closed-form density drawn through them. The density piles up at both edges because that is where the level spacing turns over, and nothing periodic was assumed to get it.

A density of states is not a spectrum

A molecule's spectrum is a list of positions and a solid's is a shape, and the shape is not the density of states. Between the two sits everything the count leaves out — which transitions are allowed, how strongly, and from where to where.

solids · Bands in a solid
Hückel levels of hexatriene. The orbital energies of the pi system, computed as the eigenvalues of the molecule's adjacency matrix. Degenerate levels share a line. The electrons are placed by aufbau with Hund's rule, so a half-filled degenerate shell shows as two unpaired spins.

Conjugation, and its limits

Every extra double bond in a chain lowers the gap between the highest occupied and lowest empty orbital, which is why long conjugated molecules are coloured. The trend has a limit, and the limit is not where the arithmetic says it should be.

beyond · Delocalisation
Chains of 4, 16, 64, 160: the levels crowd, the edges do not move. Every level of a chain of 4, 16, 64, 160 sites, drawn at its computed energy on one axis. Adding sites adds levels inside a fixed interval rather than widening it, which is what makes the large system a band rather than a wider molecule.

A band with no structure in it

Everything in this field is computed from a finite matrix with no periodicity assumed, which is a real method and a real limitation. It produces a band and cannot produce a band structure — and the difference between those two words is worth an essay, because it is the boundary of what a finite matrix can honestly say.

solids · Bands in a solid
Hückel levels of benzene. The orbital energies of the pi system, computed as the eigenvalues of the molecule's adjacency matrix. Degenerate levels share a line. The electrons are placed by aufbau with Hund's rule, so a half-filled degenerate shell shows as two unpaired spins.

Hückel theory and what it gets right

A conjugated system's orbital energies are the eigenvalues of a matrix of ones and zeroes. Nothing about carbon enters, no geometry enters, and the results that survive are exactly the ones that depend on neither.

bonding · Delocalisation
Rings of 6, 10, 20, 60 and the band at 2000. The Hückel levels of rings of 6, 10, 20, 60 atoms, all of them inside the same interval from −2 to +2, beside the density of states of a ring of 2000. The histogram is the computed levels; the line through it is the closed-form density, which diverges at both band edges.

The band limit

Every level of a ring of n atoms lies between −2β and +2β, however large n gets. The levels do not spread out as the molecule grows; they crowd into a fixed interval — and that crowding, computed, is a band with its density of states diverging at both edges.

beyond · Multicentre
How long a chain has to be before its ends stop mattering. The difference in energy per site between a ring and a chain of the same length, against that length. It falls as one over the length, which is what it means for the difference to be an end effect, and the size at which it drops below a thousandth of a β is printed.

Where a molecule stops being one

There is no size at which a molecule becomes a solid, and the useful question is a different one — how large must it be before a given property has stopped changing? The answers differ by a factor of several hundred between one property and the next, and every one of them is a measurement.

solids · Bands in a solid
A half-filled ring's cheapest excitation goes to zero. The energy of the smallest available excitation of a half-filled ring, against the number of atoms, on log axes. Every doubling at least halves it, so in the limit there is no smallest excitation — which is what a metal is, before any band picture is drawn.

What a metal actually is

Not shiny, not a good conductor, not an element on the left of the table. A metal is a system with excitations of arbitrarily small energy, and that definition can be checked on a sequence of finite rings without drawing a band diagram or mentioning conduction at all.

solids · Metal
A chain of 100 is more stable alternating. The electronic energy gained by alternating the bonds of a half-filled chain, the elastic cost of doing it, and their sum, all per site and all against the alternation. The sum has its best value away from zero, so the evenly spaced chain is not the stable structure.

A chain cannot stay even

Diagonalise a half-filled chain, alternate its bonds slightly, and diagonalise again. The electrons gain more than the springs lose, and they do so for every spring constant whatever — because the gain is steeper than a parabola near the origin and a logarithm beats any constant.

solids · Peierls distortion
cyclobutadiene: what alternation costs and gains. The π energy of cyclobutadiene as its bonds are alternated by β(1 ± δ), the elastic cost of alternating them at a stated stiffness, and the sum. The π curve falls linearly 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 vibration that lowers the symmetry

A molecule in a degenerate electronic state distorts until the degeneracy is gone. Which distortion it needs is a direct product; whether it wins is a race between a π energy falling linearly and a σ frame resisting quadratically, and both powers are measured here.

beyond · Aromaticity
The gap against size: uniform against δ = 0.15. The HOMO–LUMO gap of a half-filled chain plotted against the number of sites, on log axes, for a uniform chain and for one whose bonds alternate. The uniform sequence falls without limit; the alternating one settles at four times the alternation.

The gap is not the band width

Two numbers describe a band and they answer different questions. The width is set by how many neighbours an atom has; the gap is set by how unequal they are. A structure can have a wide band and no gap, a narrow band and a large one, and changing one leaves the other alone.

solids · Peierls distortion
Two orbitals, 4 electrons, S = 0 and S = 0.25. Two interacting orbitals with 4 electrons in them, drawn twice: once with the overlap set to zero and once with it kept at 0.25. Dropping the overlap makes the two shifts equal, which is the picture usually taught; keeping it makes the upper level rise by more than the lower falls, which is why four electrons in two orbitals is a repulsion.

The antibonding level goes up more

The two-level diagram every course draws is symmetric, and the symmetry is an artefact of setting the overlap to zero. Keep it, and the upper level rises further than the lower one falls — which is why helium has no molecule and why closed shells push each other apart.

bonding · Overlap
pyridine against benzene. The Hückel levels of pyridine beside those of benzene, which is the same graph with one diagonal entry and the bonds touching it changed. The parent's levels are symmetric about α because its matrix has nothing on the diagonal; the substituted system's are not, and the asymmetry is the size of the fitted parameter rather than a result.

Hückel with a heteroatom

A hydrocarbon's Hückel matrix contains no fitted number at all. Put one nitrogen in the ring and two arrive at once, the levels stop being symmetric about α, and the delocalisation energy stops being quotable — all from changing three entries in a matrix.

bonding · Delocalisation
Levels of a ring closed with a half turn in it. The orbital energies of a ring with one resonance integral reversed in sign, which is what half a turn in the ribbon of p orbitals does to it. The levels come in degenerate pairs from the bottom up rather than singly, so the count that closes a shell is 4n rather than 4n + 2. Degenerate levels share a line. The electrons are placed by aufbau with Hund's rule, so a half-filled degenerate shell shows as two unpaired spins.

The ring with a twist in it

Reverse the sign of one resonance integral in a ring and its levels stop being one nodeless orbital above degenerate pairs and become degenerate pairs all the way up. The count that closes a shell changes from 4n+2 to 4n, and every aromatic ring size and antiaromatic ring size exchange places.

beyond · Aromaticity
The radial set of a 6-vertex cage. The energies of the 6 radial orbitals of a deltahedron — one per vertex, pointing at the centre — as the eigenvalues of the cage's own adjacency matrix. The top level is nodeless and is the only one that is, which is Perron's theorem about a connected graph rather than an observation. Degenerate levels share a line. The electrons are placed by aufbau with Hund's rule, so a half-filled degenerate shell shows as two unpaired spins.

A cage needs one pair more than it has corners

Every closed borane holds n+1 skeletal electron pairs for n vertices, and the extra one is a theorem about connected graphs rather than an observation about boron. A cage's radial orbitals have exactly one nodeless combination, always, whatever its shape.

beyond · Multicentre
Two orbitals, 2 electrons, S = 0 and S = 0.2. Two interacting orbitals with 2 electrons in them, drawn twice: once with the overlap set to zero and once with it kept at 0.2. Dropping the overlap makes the two shifts equal, which is the picture usually taught; keeping it makes the upper level rise by more than the lower falls, which is why four electrons in two orbitals is a repulsion.

Overlap is not interaction

Two orbitals interact by the square of their coupling divided by the distance between them in energy. A coupling half again as large, with a gap four times worse, buys a third less stabilisation — so the pair that overlaps best is often not the pair that bonds best.

applied · Overlap
A 6-ring at 111°: the twist-boat. A closed ring of 6 equal bonds meeting at 111°, drawn from the coordinates the closure conditions produce. Its torsions are 16.9°, -63.6°, 44.2° and repeat; its puckering amplitude is 0.508 bond lengths at a phase of 344°, with q₃ exactly zero, so it lies on the equator. Turning the ring changes none of those numbers, which is what makes them a description of the shape rather than of the view.

The ring that cannot hold still

Cyclohexane's chair is rigid and its boat is not, and that is a statement about the rank of a matrix rather than about strain. Hold every bond length and every bond angle fixed and count what is left: the chair has nothing, and the boat sits on a continuous loop of shapes with the same bonds and the same angles.

shape · Strain
The wrong shape, fitted as well as it can be. The exact hydrogen 1s orbital and the best sums of one, two, three and six Gaussians, each with its exponents optimised for the energy. Three of them already reproduce the exact function to 99.94 per cent by overlap, which is why the method works at all — and the two places it goes wrong, at the nucleus and far out, are exactly where the other faces of this figure look.

A Gaussian is the wrong shape

Sixty years of molecular calculation are built on functions that get the two ends of an orbital wrong. A Gaussian has no cusp at the nucleus and dies too fast far away, and no number of them fixes either — while three of them already reproduce hydrogen's 1s to better than 99.9 per cent by overlap, and that is why the method works.

orbitals · Basis
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
The cost of a ratio decided somewhere else. Four bases containing exactly the same six primitives, differing only in how many of the linear coefficients the calculation may choose. The horizontal axis is the effective nuclear charge, which is this one-electron problem's only knob for a different environment; the contraction was fitted at one. At 1.238 — the exponent H₂⁺ chooses when a bond forms — the fully contracted basis is 0.0283 hartree above what the same six functions could give, and one freed coefficient removes most of it.

A contraction is a decision made once

Every published basis set freezes its primitive functions into fixed combinations, on an isolated atom, before any molecule is in sight. Freeing one coefficient recovers three quarters of what that costs — and freeing it at the other end of the basis recovers one per cent.

orbitals · Basis
Two impurities at 2 to 14 sites apart. The splitting between the two levels a pair of impurities of strength -2β pulls out of a chain of 61, against how far apart they are, on a logarithmic scale. It falls by a constant factor per site of separation, and that factor is the decay of the isolated bound state computed from its energy alone. The two levels close on the single impurity's level as the pair separates.

Two defects, and the level between them

One impurity pulls a level out of a band at exactly −√(h²+4). Two of them pull out a pair, split by an amount that falls by a factor of 0.41421 for every site of separation — which is √2 − 1, predicted from the isolated level's energy alone and measured to six figures.

solids · Defect
How far the ground state is from being one determinant. The occupations of the two natural orbitals of the Hubbard dimer against the repulsion. At zero they are two and nothing, which is a single determinant exactly. As the repulsion grows they converge on one and one, which is a state no single determinant has — the failure is in the description rather than in the number. The entropy of the occupations rises to ln 2, one bit: the two determinants of the singlet.

Two kinds of correlation, and only one is small

Correlation energy is defined as a subtraction, and the definition hides that the thing subtracted is not one thing. In the two-site model the local power of the correlation energy in the interaction falls from two to one — and at every interaction strength it equals, exactly, the occupation of the bonding natural orbital.

beyond · Correlation
The characters of C6, and the levels they are. The 6 representations of the ring's rotation group, drawn as points on the unit circle at 2πk/6. Each level is twice the horizontal coordinate: a representation and its complex conjugate have the same real part, so they are degenerate, and the levels pair up automatically. Only k = 0 — and k = n/2 when n is even — lands on the real axis, so one level is unpaired at the bottom and the shell closes at 4n + 2. Nothing here has been diagonalised.

The ring's levels are its group's characters

The Hückel energies of a cyclic system are twice the real part of the characters of its own rotation group, so the level pattern — one level, then pairs, then one more if the ring is even — is a group theorem rather than a calculation. Hückel's 4n+2 rule is a statement about the representations of a cyclic group and about nothing else.

symmetry · Aromaticity
The cliff, and the slope leading to it. The smallest eigenvalue of the overlap matrix, and what an extra function is worth, as that function is brought towards one already in the basis. Both fall together, and the energy stops improving long before the matrix stops being invertible.

A function that is already there

Adding a function to a basis set can only lower the energy, so a bigger basis is a better one. What that leaves out is that the same act makes the functions less independent: an optimised basis's smallest overlap eigenvalue halves with every function added, and a function placed on top of one already present buys less than a millionth of what a well-placed one buys while driving that eigenvalue to 4×10⁻¹⁰ — past which the calculation is refused outright.

orbitals · Basis
The same difference, at five repulsions. The charge transferred to the more electronegative of two atoms against the difference in their orbital energies, at five strengths of the repulsion between the two electrons. Only the topmost curve is the two-level result; every other one moves far less charge at the same difference.

A difference does not make a transfer

Every electronegativity scale reports one thing: how far apart two atoms are in their appetite for electrons. Put that difference into a model with repulsion in it and the charge it actually moves is not determined at all — the same difference of one moves 0.447 of an electron with no repulsion and 0.0019 with a strong one, a factor of two hundred and forty-two, and nothing on any scale distinguishes the two cases.

bonding · Electronegativity
The two figures, per electron. For each ring, the delocalisation energy against isolated double bonds and against the same ring with one bond deleted, divided by the number of pi electrons. The upper bar of each pair is the usual reference; the lower one is the cycle's own contribution, and the two are nowhere proportional.

A stabilisation is measured from somewhere

Benzene's delocalisation energy is 2β against three isolated double bonds and 1.0121β against the same six carbons with one bond deleted. Cyclobutadiene's is exactly zero against the first reference and −0.4721β against the second, so one of the two references records the most famously destabilised ring in the subject as neutral.

beyond · Delocalisation
Four fillings, four periods. four fillings of a ring of 120, and for each of them what a distortion of every available period is worth. Every period pays the same elastic cost, so the bars compare what the electrons give back and nothing else. The winner is one over the filling in every case.

The distortion the filling chooses

A half-filled chain of equal bonds is unstable and alternates — long, short, long, short. That is the case everyone is shown, and it is one case. Fill the chain a third of the way instead and the alternation is worthless: what wins is a pattern that repeats every three bonds, and the period is one over the filling at every filling tried.

solids · Peierls distortion
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
water: every level up to J = 4, from a matrix. The rotational levels of water at κ = -0.4322, each J diagonalised in the symmetric-top basis. A symmetric top would show one level per K with everything above K = 0 doubly degenerate; here every degeneracy is split, and the size of each splitting is what the third constant is measured from.

The top that reports all three

A symmetric top hides one of its two rotational constants in every line of its spectrum. Break the symmetry and the hiding stops: for water, twenty-three of the twenty-five levels up to J = 4 move when A is changed, and the two that do not are the ground state and the one at B + C. There is no formula for any of them.

spectra · Rotation
Which numbers survive a change of frame, and which are coordinates. Every quantity this field quotes, against the four things that can be changed without changing the molecule: the zero of energy, the unit, the reference state a stabilisation is measured from, and what a "per" quantity is divided by. A filled mark is a quantity that moves. three of the 9 survive all four, and every one of them is a property of the eigenvectors rather than of the energies. The first two changes are exact symmetries of the model, so a quantity that moves under either is a coordinate and not a quantity at all.

Which numbers carry a frame

A Hückel calculation is written in two numbers nobody computes and quoted against reference states nobody measures, so every quantity it prints is a quantity in a frame. Nine of them, tested against four changes of frame: three survive all four, and the one that survives both exact symmetries and looks safest — a dimensionless ratio between two molecules — is the most fragile of all, because against one reference its denominator is exactly zero.

bonding · Delocalisation
How anisotropic a structure has to be. The fourth moment of a layered band, divided by the square of its second so that the width drops out, against the coupling between layers. The curve is the closed form that follows from adding cumulants — 3 − 3(2 + λ⁴)/(2(2 + λ²)²) — with nothing fitted; the points are eigenvalues. At λ = 1 it is the three-dimensional value 2.5 and at λ = 0 the two-dimensional 2.25, and it is halfway between them at λ = 0.47 — so the question of when a band is effectively two-dimensional has a number rather than an opinion.

Two bands, and the shape of each

A band's shape has a closed form when the structure it belongs to factorises. Put a second band beside it and the form survives exactly while the two do not talk, and fails as the square of the coupling once they do. Make one direction weaker instead and the form survives everywhere — 3 − 3(2 + λ⁴)/(2(2 + λ²)²), agreeing with eigenvalues to nine decimal places — so how anisotropic a structure has to be before its band is two-dimensional is a number: λ = 0.46518.

solids · Bands in a solid
Two errors, opposite signs, four orders of magnitude apart. A finite basis makes H₂⁺'s binding too large by letting each atom borrow the other's functions, and too small by describing the molecule incompletely. Both are computed here against the exact binding of 0.102634 hartree. The second is thousands of times the first at every basis size, and it is the first that counterpoise removes — so the corrected number is further from the true one than the uncorrected at every row of this table.

The basis the other atom lent

Two atoms in a molecule are described in each other's functions and the separated atoms are not, so the molecule is treated better than the pieces and the binding comes out too large. That is the basis set superposition error, it is removed by a standard correction, and for H₂⁺ in four Gaussians a centre it is six tenths of a microhartree against an incompleteness error of twelve millihartree — a factor of eighteen thousand the other way.

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
The current does not divide equally between equal rings. The current each ring of an acene carries under a uniform field, ring by ring, for four acenes. Naphthalene's two rings are equal by symmetry; anthracene's middle ring carries 1.180 times what its outer ones do, and tetracene's inner rings 1.222 times. Every ring has the same area and the same six carbons, and the response is a matrix rather than a set of parallel loops.

The current does not divide

A fused ring system's response was computed from the areas of its rings, and the obvious next question was whether the current divides between them the way it divides between two resistors. Giving each ring its own flux and taking the second derivatives says no: the response is a matrix, its off-diagonal entries are nearly half its diagonal ones, and anthracene's middle ring carries 1.18 times what its outer rings do.

symmetry · Aromaticity
The shared bond's response falls and reverses, at every coupling. How far the relaxation moves the most interior cross bond of an acene, against the number of rings, at three couplings. Every series falls monotonically and every one changes sign — at λ = 0.2 between 3 and 4 rings, at λ = 0.4 between 3 and 4 rings, at λ = 0.6 between 4 and 5 rings. Where it crosses is a property of the coupling; that it crosses is not. The natural reading — that the effect grows with the number of rings feeding a bond — is refused by every one of these curves.

An anomaly that is not the first of a series

Naphthalene's shared bond behaves differently when the geometry is allowed to answer back, and the obvious reading is that two rings feeding one bond is what does it. Run it up the acenes and the effect falls at every step and changes sign: +0.02036 for two rings, +0.00678 for three, −0.00029 for four, −0.00603 for five. Two rings feeding a bond is not the beginning of anything.

bonding · Delocalisation
However hard the π channel is driven, the counted orbital stays the metal's. The metal's share of the filled T₂g orbital against the π coupling, with the eighteen-electron count drawn beside it. The share falls from 1 to 0.5467 across a coupling range of 80,000 cm⁻¹ and approaches a half from above without reaching it: the lower eigenvector of a two-level problem always carries more of the lower basis function, whatever the coupling. The count is eighteen at every point.

The count that cannot be broken by strength

Back-donation puts electrons into orbitals that are not the metal's, and the eighteen-electron rule counts the metal's nine. Turning the π channel up as far as it will go never breaks it: the counted orbital's metal share falls from 100 per cent to 54.67 and approaches a half from above without reaching it. What does flip it is not strength but order.

applied · Electron count
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 state in the gap is a particle in a box the alloy happened to make. The participation ratio of the 40 levels nearest the gap centre, over 5 chains of 400 sites, against the length of the run of like sites each one sits on. The line is 2(L+1)/3, the participation ratio of the ground state of an isolated chain of L sites, with nothing fitted. 38 of 40 lie on it to within 3.1 per cent. The ones above it are states shared between two runs close enough to talk, which is a defect band beginning.

A particle in a box the alloy made

A random alloy's gap is set by the longest run of like sites. What sits at the edge of that gap turns out to be the simplest state in quantum mechanics: a particle in a box of L sites, occupying 2(L+1)/3 of them, to within three per cent and with nothing fitted — except for the two states in forty that found a second run to share.

solids · Defect
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
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
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
The pair's regime is not a property of the pair. How much the pair's bonding responds to its own overlap — the ratio of what it is bonded by at an overlap of 0.4 to what it is bonded by at 0.1 — with and without a third orbital coupled to both. Alone it is 11.83, which is the regime in which bonding tracks overlap. With a third orbital present it falls to 1.46, 0.92, 0.74 — and two of those are below one, meaning a fourfold increase in the overlap between the two atoms buys them less bonding rather than more. The coupling comes from the overlap by the Wolfsberg–Helmholz rule with K = 1.75, which is fitted rather than derived. Every stabilisation here inherits that; the shape of the curve against separation does not, because K is a constant.

A regime that belongs to the neighbours

Two orbitals at the same energy are bonded in proportion to their overlap and two far apart are barely bonded at all — two regimes, and the natural question is whether the regime is a property of the pair. It is not. Put a third orbital beside them and the pair's response to its own overlap falls from twelvefold to less than one: more overlap buys less bonding.

wrong · Overlap
Two events, not one — and both below the n = 2 shell's. The field at which each coupled pair's shift stops being quadratic and starts being linear, for the n = 3 shell and for the n = 2 shell. The n = 3 shell has two, a factor of 3.66 apart, and both are far below n = 2's single one — so the linear effect returns is two events in a shell with a d, and it happens at a field thirty times weaker than in a shell without one.

Two events where there was one

A broken symmetry returns at a field where the coupling matches the gap it has to overcome, and in a shell with two levels that field is a single number. The next shell up has three, two gaps and two dipoles — so there are two crossovers, a factor of 3.66 apart, both of them below the single one of the shell below, and the exponent takes more than a decade of field to travel between them.

symmetry · Representation
Two bands, square below, triangular above. The 128 levels of a structure carrying two orbitals on every one of its 64 sites, where the lower orbitals hop on a square net and the upper ones on a triangular net, separated by 12 and coupled at 0.06. The two bands have different widths and different shapes before anything mixes, which is the situation a single net cannot construct. Every level here is also reproduced by a two-by-two problem, one per mode, to 3e-13.

The constant that belonged to one net

Divide one second-order departure by another and the coupling cancels, leaving one constant times the gap — which reads as arithmetic. It was arithmetic about a square net. On a triangular one — the same graph for both bands, nothing else changed — the quotient drifts by a hundred and twenty-eight per cent.

solids · Bands in a solid
Every crossover a shell has. The field at which each coupled pair's linear behaviour returns — the gap between the two levels divided by twice the dipole joining them — on a logarithmic axis. The n = 3 shell has 3 distinct fields rather than four, because its m = ±1 half is a two-level ladder with one coupled pair. All three are below the n = 2 shell's single one.

Three events, and a ratio of two dipoles

The m = 0 half of a shell has two crossovers because it has three levels and two coupled pairs. The other half has two levels and one, so the whole shell has three distinct fields rather than four — and two of the three share a gap exactly, which makes the ratio between them a ratio of two dipoles, 2/√3.

symmetry · Representation
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
One geometry, four electron counts, four answers. The bond order between two orbitals with exactly no overlap and no resonance integral, as the third orbital's energy is swept, at every count the trio can hold. With none it is identically zero. With two it is positive and rises past one. With four it is negative and reaches -0.954. With six it is a horizontal line — the third orbital's energy stops mattering entirely.

A filled shell is not an empty statement

A bond order of −0.954 between two orbitals with no overlap and no resonance integral invites the prediction that at six electrons — every level occupied, the sum over a complete set — it would be exactly zero. It is exactly one seventh, and the reason is that a complete set in a non-orthogonal basis sums to the inverse of the overlap matrix, which has entries where the overlap has none.

wrong · Overlap
What happens to the level that was exact. The three levels of the trio as one of the two outer orbitals is raised. At zero detuning the middle one sits at -13.6 exactly — it is the antisymmetric combination, and nothing of its symmetry exists for it to mix with. The moment the two are made inequivalent that statement is gone: the level leaves linearly, and the other two barely move by comparison.

A symmetry holds or it does not

One level of a three-orbital trio sits at the free-atom energy exactly, at every third-orbital energy, because the antisymmetric combination of the pair has nothing of its own symmetry to mix with. Detuning one of the two by a twentieth of an electron volt moves it by half of that — first order, immediately, with no protected regime at all.

wrong · Overlap
The A–C overlap changes by a fifth and the exact level does not move. The three levels of the trio as the overlap between A and C is raised from 0.25 by up to 0.2, with B's overlap to C held and both site energies at -13.6 eV. The two outer orbitals stop being equivalent at the first step. The lowest level falls by 0.80 eV and the highest rises by 5.24, and the middle one stays at -13.6 eV — its largest departure over the whole sweep is 2.7×10⁻¹⁴ eV, which is rounding.

A level no symmetry was protecting

A three-orbital trio keeps one level at the free-atom energy exactly, and the reason given was that its two outer orbitals are equivalent. Make them inequivalent by changing one overlap rather than one energy and the level does not move at all — not to first order, not to any order, at any energy of the third orbital. It was never the symmetry. It is allyl's non-bonding orbital, held by a count.

bonding · Overlap
The verdicts hold and the number does not. For each predictor, what survives thirteen strictly increasing reparameterisations of its own scale. The count of exact ties, the share of the variation those ties leave unexplained, and the count of discordant pairs are identical under every one of them — to the last bit, because each asks only about the order of the predicted values and an increasing map preserves order. The slope floor asks for a ratio of differences, and an increasing map does not preserve differences.

A floor on models written in one scale

A tie asks whether two predicted values are equal and a discordance asks whether two differences have the same sign. Both are questions about order, and a strictly increasing change of scale preserves order — so both are exactly invariant under thirteen reparameterisations of all three predictors. The slope floor asks for a ratio of differences, and it moves by factors of nineteen, ten and eleven thousand.

bonding · Models
The unpaired electron's level stays put while its spin moves to the far end. Above, the trio's three levels as the overlap of A with C is raised from 0.25 to 0.45 while B's stays at 0.25. The lowest falls from -16.26 to -17.07 eV and the highest rises from -8.02 to -2.78 eV; the middle one, which holds the radical's unpaired electron, stays at −13.6 eV. Below, the spin on A and on B over the same change: from a half each to 0.236 on A and 0.764 on B, with the ratio of squared overlaps drawn as open circles on top.

The spin the count does not hold

Three orbitals in a row keep one level at the free-atom energy however the overlap of one end is changed, and at three electrons that level holds the radical's unpaired electron. Its energy does not move. Its spin does: from half on each end to 0.236 and 0.764 as one overlap goes from 0.25 to 0.45, exactly the squared ratio of the two overlaps, at every energy of the middle orbital. A coupling between the ends moves the level and cannot move the spin. The energy and the spin are answering to different things.

bonding · Overlap

Named alongside it

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

Model limitDegeneracyHückel theoryClosed formConventionBands in a solidApproximationReference stateOverlap integralTight-binding modelsConjugationDelocalisation

All concepts