Theme

The thread: A basis is not a thing

Hybrid orbitals, localised bonds and canonical molecular orbitals are related by transformations that leave every observable unchanged. Arguing about which is real is arguing about a coordinate system.

101 essays carry this theme. The first 24 are shown with their figures; the rest are listed without, so that the page does not grow without limit.

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. When the molecule does not stop

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.

methane — Td. The molecule with the point group found from its coordinates rather than looked up: every candidate operation was applied and kept when it permuted the atoms among themselves. What is taught wrongly

Hybridisation does not explain

Methane's photoelectron spectrum has two bands, not one. Four equivalent sp³ bonding orbitals cannot produce that, and the resolution is that hybridisation was never a claim about what a measurement would find.

sulfur hexafluoride — Oh. The molecule with the point group found from its coordinates rather than looked up: every candidate operation was applied and kept when it permuted the atoms among themselves. Beyond the octet

Hypervalency without d orbitals

Sulfur hexafluoride is not d²sp³ hybridised. The d orbitals are far too high in energy to contribute meaningfully, the bonding is three-centre four-electron, and the textbook account has been known to be wrong for fifty years.

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. What a spectrum settles

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.

benzene — D6h. The molecule with the point group found from its coordinates rather than looked up: every candidate operation was applied and kept when it permuted the atoms among themselves. Beyond the octet

Delocalisation

Benzene does not alternate between two structures. It has one structure, and the two Kekulé forms are basis functions in a description of it — which is a different and much less exciting claim than the one usually made.

The same d shell in four fields. The five d energies in octahedral, tetrahedral, cubic, square planar coordinations, computed and drawn in units of the octahedral splitting. The dashed line is the barycentre, which the five levels sum to whatever the arrangement. What the shape is for

Two models, one ratio

A tetrahedron splits a d shell by four ninths of what an octahedron does. Two models that share nothing but the ligand directions — an integral over a point-charge potential and a rotated diagonal matrix — both produce that number to eight decimal places, and neither was told it.

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. Beyond the octet

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.

sp3 hybrids. The directions the hybrids point, with the angle between them computed from the coefficients rather than quoted. The set is an orthogonal transformation of the atomic orbitals, so it describes the same space in different coordinates. Bonding models

Hybrids are a basis

An sp³ hybrid set is an orthogonal matrix applied to the atomic orbitals. Rotating a basis changes no observable, so asking whether the electrons are really in hybrids is asking which coordinate system nature prefers.

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. When the molecule does not stop

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.

Splitting goes with overlap. For each pair, the atomic levels on the outside and the combinations they form in the middle, with the splitting drawn in proportion to the computed overlap. A pair that symmetry forbids does not split at all, because its overlap is exactly zero. Bonding models

Molecular orbital and valence bond

Two frameworks, taught as rivals, describing the same molecules. One starts from delocalised orbitals and localises; the other starts from localised bonds and delocalises. Pushed far enough they meet.

Which rings close a shell. Each ring is filled with its own number of pi electrons and asked whether the highest occupied shell came out full. Of the rings drawn here, C6 and C10 close — at 6 and 10 electrons — which is Hückel's 4n+2, produced here rather than recalled. Beyond the octet

Where two-centre bonding stops

A bond between two atoms is a special case, not the general one. Rings, clusters and metals are held together by orbitals spread over many centres, and the arithmetic that describes them is the arithmetic already used for benzene.

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. What the shape is for

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.

sp3 hybrids. The directions the hybrids point, with the angle between them computed from the coefficients rather than quoted. The set is an orthogonal transformation of the atomic orbitals, so it describes the same space in different coordinates. Where the atoms go

Bent's rule, against the substituent series

More electronegative substituents get more p character, so the angle between them shrinks. The rule predicts a trend, the trend is observed, and the quantity it depends on turns out not to be one quantity.

Hückel levels of cyclobutadiene. 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. What is taught wrongly

Delocalisation is stabilising, and other things that are false in general

Spreading electrons over more centres is supposed to lower their energy. Cyclobutadiene is delocalised in exactly the same sense as benzene, gains exactly nothing by it, and the arithmetic says so before any experiment does.

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. Bonding models

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.

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. Beyond the octet

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.

Orbitals at the 90 per cent contour. Several orbitals drawn at the same enclosed fraction and, unless stated otherwise, at the same scale — so the sizes on the page are the sizes. Each contour was solved for separately by integrating that orbital's own density. Contours drawn: 2pz at 90% of its density, |ψ| = 9.48e-3; 2px at 90% of its density, |ψ| = 9.49e-3; 2py at 90% of its density, |ψ| = 9.49e-3. Orbitals

Complex harmonics against real ones

The p orbitals every chemist draws are not eigenfunctions of anything. They are real combinations of the complex solutions, chosen because they point along axes — and the choice is invisible until a magnetic field makes it matter.

One metal orbital, two ligands competing for it. Metal–ligand bond orders in a three-orbital model as the left-hand ligand's interaction is turned up. Its own bond order rises and the bond order to the ligand opposite falls, from 0.62 at equal strengths to 0.42 at the strongest. Nothing else in the model can carry the effect: switch the second bond off and it vanishes exactly. What the shape is for

The trans influence is an overlap argument

Two ligands on opposite sides of a metal both bond through the same metal orbital, and there is only one of it. Strengthen one and the bond order to the other falls — computed exactly on three levels, and measured as a bond length that grows by a tenth of an ångström.

Hückel levels of three-centre four-electron. 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. Beyond the octet

Three-centre bonding, computed

Three orbitals in a line and four electrons: a bonding level, a level with exactly zero amplitude on the central atom, and an empty antibonding one. The middle atom never exceeds an octet, and the ligands carry the charge — which is why every molecule that needs this arrangement has electronegative ligands.

π bond orders in naphthalene. Each bond drawn with a thickness set by its computed pi bond order, and the number printed beside it. The orders come from the eigenvectors and cost nothing extra once the matrix is diagonalised. Bonding models

Bond order from the eigenvectors

Having diagonalised the matrix, the coefficients are already there. Two sums over them give every bond's order and every atom's charge, and naphthalene's three kinds of bond come out in the order the measurements find them.

The Td character table. The irreducible representations of the molecule's point group. The class headings carry the number of operations found by generating the group from the coordinates, and those counts were produced before this table was opened. What symmetry decides

Degeneracy is a group theorem

How many orbitals can share an energy is decided before any energy is computed. The dimensions of a group's irreducible representations are the only degeneracies it permits, and a molecule with no representation larger than one cannot have a degenerate level at all.

Oh descending to D4h. Every irreducible representation of Oh restricted to D4h and reduced there — A1g, A2g, Eg, T1g, T2g, A1u, A2u, Eu, T1u, T2u. A representation that arrives in more than one piece is a degeneracy the lower symmetry cannot hold, so a level carrying it must split when the molecule distorts. What symmetry decides

Descent in symmetry

Lower a molecule's symmetry and its labels stop being available. Which of them survive, which split, and into what, is decided by restricting characters to the operations that are left — arithmetic, not a table to be looked up.

pyrrole against cyclopentadienyl anion. The Hückel levels of pyrrole beside those of cyclopentadienyl anion, 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. Beyond the octet

Six electrons in a ring that is not all carbon

Pyrrole and the cyclopentadienyl anion hold the same six π electrons in the same five-membered ring. The anion spreads them perfectly evenly; pyrrole's nitrogen keeps 1.720 of them and furan's oxygen keeps 1.791 — the more electronegative atom donates less, and the bond orders to it fall with it.

The s character a bond angle requires. The fraction of the one s orbital each of two equivalent hybrids must carry in order to meet at a given angle, from pure p at 90 degrees to half s at 180. The marks are measured bond angles — H₂O, NH₃, H₂S, PH₃ — and the column beside them is the budget each implies: what one bond hybrid takes, and what is left for the lone pairs to share. The relation is orthogonality, not a model of bonding, and nothing here is fitted. Where the atoms go

The angle does not fix the hybridisation

Two equivalent hybrids are orthogonal only at one relation between their s character and the angle between them. Run it backwards on measured bond angles and water's bonds come out sp³·⁹⁹, hydrogen sulfide's sp²⁷, and phosphine's lone pair takes eighty-four per cent of the one s orbital.

What is taught wrongly

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.

Bonding models

The localisation transformation, demonstrated

Four equivalent bonds or one a₁ and three t₂ — two descriptions of methane's bonding electrons that disagree about everything except the electron density, which they agree about to the last bit a double can hold.

What symmetry decides

An infinite group, worked in a finite one

A linear molecule has infinitely many symmetry operations, and every formula in character theory divides by the number of them. The standard device is to work in a finite subgroup — and it is worth computing what that trade costs rather than putting it in a footnote.

Beyond the octet

What one pair can hold together

Put a single electron pair into a ring of any size and it supplies a total bond order of exactly two and a π energy of exactly 4β — three atoms, eight atoms or six hundred. Spreading a pair over more centres divides the bonding among them; it neither creates nor destroys any.

When the molecule does not stop

A defect is a level in the gap

Change one site's energy in a chain of a hundred and sixty and a level leaves the band, carrying a state that lives on a handful of atoms. In one dimension it happens for any change however small — the threshold measured on chains of forty, eighty and a hundred and sixty halves with every doubling, so there is no threshold at all.

Bonding models

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.

Beyond the octet

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.

What is taught wrongly

What a photoelectron spectrum measures

The bands of a photoelectron spectrum are routinely read off as orbital energies. They are ionisation energies, which is a different quantity — and the identification rests on two errors of about an electronvolt each that happen to have opposite signs.

Orbitals

A filled shell has no shape

Sum the angular densities of a complete p shell and the answer is 3/4π in every direction, to sixteen decimal places. A filled d shell gives 5/4π. The lobes are in the decomposition and not in the density, and nothing that measures a closed-shell atom can see them.

What is taught wrongly

A half-filled band is not always a metal

Every band picture rests on an approximation that a whole class of materials refuses — each electron moving in an average field, never seeing another one individually. Where the repulsion is strong enough, a half-filled band describes an insulator, and no amount of care with the band fixes it.

What symmetry decides

Every group a molecule can fall to

A distortion takes a molecule's symmetry away, and what is left is not a free choice — it has to be a group. Closing subsets of methane's twenty-four operations returns thirty of them, which is exactly the number the symmetric group on four letters has.

Where the atoms go

Hybrids that were never orthogonal

Two sp³ hybrids are orthogonal at 109.47° and nowhere else. Water's are drawn at 104.5° and overlap by 0.062; cyclopropane's are drawn at 60° and overlap by 0.625, which is not a small correction to a basis but a description that has stopped being one.

Bonding models

Water's lone pairs are not a pair

Every course draws two equivalent lone pairs on water, pointing away from the hydrogens like a pair of ears. Its photoelectron spectrum shows the two bands they would produce at 12.6 and 14.7 electronvolts, two point one apart, in different symmetry species.

Orbitals

The atom does not bring its own orbital

Build a one-electron diatomic from two hydrogen 1s functions and it comes out 25 per cent too long and 37 per cent too weakly bound. Let the molecule choose how large those functions are and the bond length is right to three figures, at an exponent of 1.238 — the orbital contracts by a quarter when the bond forms.

When the molecule does not stop

What holds a solid together

Four kinds of interaction between two units of matter, each computed from a stated model and put on one logarithmic scale. The ordering is not the one usually taught — an ion pair at contact beats a shared pair, which is a fact about the comparison rather than about the numbers.

What symmetry decides

Why a character table stops where it stops

A character table is square, and both of its dimensions are forced. The number of representations equals the number of classes, the squares of their dimensions sum to the order of the group, and between them there is no room for another row.

What a spectrum settles

Fewer bands than electrons

Methane has eight valence electrons and two photoelectron bands. Ammonia has eight and three; water has eight and four. The count is not of electrons, not of bonds and not of orbitals — it is of the symmetry species the occupied orbitals fall into, and it falls as the symmetry rises.

Beyond the octet

Hypervalency is about the ligands

The three-centre four-electron bond is offered as the reason sulfur hexafluoride needs no d orbitals. Read it forwards instead of backwards and it is a requirement rather than a permission — the arrangement puts half an electron onto each ligand before any electronegativity difference is applied, which is why the hypervalent compounds are fluorides and SH₆ is not a compound.

When the molecule does not stop

The lattice sum that depends on the order of adding

An ionic solid's binding is the sum of every pair of charges in it, and the series does not converge absolutely — rearranged, it gives a different answer. That is a genuine mathematical difficulty rather than a technicality, and it is the clearest example of something a real-space, neighbour-by-neighbour method cannot compute at all.

Bonding models

Where molecular orbital theory dissociates

The molecular orbital description of a two-electron bond puts both electrons on the same atom half the time — at every bond length, including infinite. The exact answer falls from a half to 0.0039 as the atoms separate, and the point where the two standard models are equally wrong is exactly U = 4t.

Orbitals

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.

Where the atoms go

The lone pair is not the missing term

Ammonia's dipole is 1.47 debye and nitrogen trifluoride's is 0.235, although the N–F bonds are far more polar than the N–H ones. The bond sums explain the reversal exactly and point in opposite directions — and the lone pair that is supposed to make up the difference has to be worth 0.58 debye in one molecule and at least 1.48 in the other.

What symmetry decides

The projector is unique, the basis is not

A reduction says how many times each representation appears. It cannot say which combinations of orbitals they are — that takes a projection operator, and for a degenerate representation the operator returns a different pair of orbitals depending on which function it is handed first. Both pairs are correct, they span the same space, and every energy computed from either is identical.

What a spectrum settles

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.

Orbitals

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.

What is taught wrongly

A ranking is not a difference

Four electronegativity scales agree about order to a Spearman coefficient of 0.99 and disagree about size by a factor of sixty-two. Put each on its own range and the O–H bond is 39 per cent of the way across on three scales and 4.5 on the fourth.

Bonding models

Two pictures, one plane

Molecular orbital theory and valence bond theory are taught as rival descriptions of a two-electron bond. In a model small enough to solve exactly they are two vectors in a two-dimensional space, the exact answer lies in the plane they span at every repulsion, and it is neither of them at any repulsion but two.

Bonding models

A method that is not additive

Put a molecule next to a copy of itself, far enough away that they do not interact, and the exact energy doubles exactly. A truncated calculation does not — because the truncation forbids both halves being excited at once, which is something the pair can do and neither half can.

Beyond the octet

Four centres, and the pair that will not localise

The occupied orbitals of a molecule can be mixed freely without changing anything observable, and the freedom is usually spent on making them as local as possible. For the methyllithium tetramer the answer is four centres — one carbon and the three lithiums of the face it caps — and no mixing of the four pairs reduces it.

Beyond the octet

Hypervalency does not stop at three centres

The three-centre four-electron bond is written up everywhere as an arrangement peculiar to hypervalent molecules. It is the first member of a family — five centres and six electrons, seven and eight — and the family predicts alternating bond strengths that the polyiodide crystal structures have.

Orbitals

The orbital in momentum space

Every orbital has a second picture as complete as the first and almost never drawn. Nothing is added by taking it — it is the same function in the other variable — but the uncertainty product falls out of it, and the functions quantum chemistry is built from turn out to be the only ones that attain the bound.

Orbitals

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.

Bonding models

Two bent bonds, or a σ and a π

A carbon–carbon double bond is drawn two ways and the two look like rival claims about what is there. They are one occupied space in two bases, related by a rotation of exactly forty-five degrees: the density is identical to the last bit a double holds, the bent pair are sp⁵ hybrids at 50.77° to the axis, and their charge sits 0.235 Å off the plane of the molecule where neither canonical orbital's does.

Orbitals

The measurement a basis was not fitted to

Six Gaussians reproduce hydrogen's energy to eleven parts in a hundred thousand and its Compton profile to three parts in a thousand — twenty-five times worse, on a quantity an X-ray scattering experiment measures directly. The gap between the two errors widens as the basis is improved, because the energy is the one property a variational fit is best at.

Bonding models

A weight that depends on how it is weighed

The ionic character of a two-electron bond is quoted as a percentage. For one wavefunction at hydrogen's bond length, three conventions in the literature give 18.73, 34.74 and 5.88 per cent — a factor of six — and on a wavefunction with no ionic structure in it at all, one of them still reports a quarter.

What symmetry decides

How far, and along which coordinate

A continuous symmetry measure returns one number: how far a structure is from a shape. Projecting the same displacement onto the twelve symmetry species of benzene's group turns it into a list that sums to 1.000000000000, says which coordinates the structure left along, and shows that the totally symmetric part — 3.4 per cent of this one — moves every atom by 0.0200 ångström and lowers the symmetry measure by nothing at all.

Orbitals

The property that gets worse

One Gaussian fitted to hydrogen gets the mean radius exactly right — 1.500000, against an exact 1.5. Two Gaussians get it wrong by 1.4 per cent, which is three hundred and thirty-two thousand times further out, while the energy improves fivefold. The variational principle bounds one number and says nothing whatever about any other, and the sequence of errors in everything else need not even be monotone.

Bonding models

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.

Beyond the octet

One scale, from two centres to a cage

How many centres a pair of electrons holds together is two different numbers, and they separate exactly where the bonding is most deficient: the methyllithium tetramer's pairs sit on four atoms each and have a participation number of 2.54. Run the same measurement up the scale and the twelve-vertex borane refuses it — its localisation has at least ten maxima differing by six parts in a thousand, so the number of centres is not an output for it at all.

Where the atoms go

Three bent bonds, and the same hybrid

A triple bond localises into three bent bonds at 101.537 degrees to one another, each an sp⁵ hybrid with exactly one sixth s character. A double bond's two bent components are sp⁵ hybrids at 101.537 degrees. The two are the same hybrid, built out of different frameworks and in different numbers, and the reason is that a third of a half is a half of a third.

Bonding models

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.

Beyond the octet

How many descriptions a cage has

A localisation is a maximisation, and running it once reports the maximum it reached rather than the maximum there is. Run to exhaustion on a twelve-vertex borane it finds ten answers and then three batches of twelve starts in a row that find nothing new — and none of the ten is another one seen from a different side, because a symmetry of the cage cannot change a multiset of participation numbers and all ten multisets differ.

Orbitals

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

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.

What the shape is for

An integer nobody measured

The oxidation state of chromium in the hexacarbonyl is zero. Its charge, computed from the same wavefunction, is anywhere between −3.04 and −0.98 depending on how the shared electrons are divided — and the integer sits outside that whole range. The electron count, meanwhile, is eighteen at every point on it.

Beyond the octet

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.

Bonding models

The hybrids that point outside the bonds

Coulson's relation says two equivalent hybrids sharing an s orbital are orthogonal only between ninety and a hundred and eighty degrees. Cyclopropane's carbons make sixty, so its ring hybrids cannot point at the atoms they bond to — and the same relation says by how much they miss: 22.75 degrees each, falling to 0.03 in cyclohexane.

Orbitals

A correction that is two functions

A symmetric pair's counterpoise correction splits exactly in half, at every separation, to the last digit — which is why one frozen number describes it. Give the two atoms different charges and the split runs from 0.03 per cent to 67, changing places at 6.29 bohr: the correction a single number was standing in for is two functions of different shapes.

Bonding models

The angle that does not have to be searched for

A carbonyl's two bent components have no symmetry making them equivalent, so the mixing that best localises them looks like something to search for and their s characters look like two different numbers. Neither happens. The localisation functional is a quadratic with no linear term, so its maximum is at forty-five degrees or at the ends — bent bonds or canonical ones, decided by one inequality, with nothing in between.

Beyond the octet

The square that wastes an orbital

The orphan count that prices hypervalency was treated as a property of a molecule's composition — ligands plus lone pairs minus four. Run on twenty-six arrangements of the same ten molecules it is right for twenty-three and wrong for three, and all three are flat. A planar arrangement gives a main-group centre three usable orbitals rather than four, so the count is a property of the shape.

Orbitals

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.

Beyond the octet

Where the count stops being an effort

A cage's localised descriptions were counted by a search that stopped when it stopped finding new ones, which makes the count a property of the stopping rule. Run to four thousand starts the count is fourteen and the last new description appears at start 124 — after which three thousand eight hundred and seventy-six starts add nothing. The four rarest are found six times in a thousand, which is what says nothing rarer is hiding.

What is taught wrongly

A capacity that is largest where there is none

A cubic through four charge states implies a largest amount of electron an atom can accept, and it is natural to ask whether the idea survives its own model. It does not. The three atoms whose anions are not bound — beryllium, magnesium and nitrogen — have the largest capacities in the set, two of them infinite; the three smallest capacities all belong to atoms whose anions are bound. And among the fourteen where the test cannot bite, the ordering is sensible.

Bonding models

A bond with nothing in the middle

Two head-on 2p functions have an overlap that changes sign at 5.03 bohr, and the picture of what that means is worth drawing. Below that separation the combination the model fills has a nodal plane through the midpoint of the bond, so two electrons in it put exactly nothing between the nuclei; above it the same model fills the other one. Which picture a bonding orbital has is decided by a separation.

Orbitals

The assembly that counts one share twice

A counterpoise correction is divided unequally between its two centres, and the first place that matters is a three-fragment system, where the pairwise corrections are added up and the assembly must double-count one share and undercount another. It does: the heavy centre is over-corrected by seventeen per cent and the light ones under-corrected by two and a half, and the two do not cancel.

Bonding models

The node that decided a picture

A hydrogenic calculation cannot say whether formaldehyde's localised description is two bent components or σ and π, because the integral that decides it came out at 0.00718 where it should be the largest in the molecule. Replacing the hydrogenic radial functions with nodeless Slater ones moves that overlap by a factor of ninety-seven, leaves the π overlap identical to the last bit, and flips the answer: the description is bent.

Beyond the octet

Three shapes from one search

A cage's fourteen localised descriptions fell into two groups with a factor of five between them and nothing in the gap, which made counting the big basins look like a stopping rule. Two more cases from the same family give a smooth tail over sixty-three descriptions and a single basin holding ninety-eight per cent. One search, one family, three shapes.

Orbitals

The correction that gets harder to assemble

Summed across a trimer, pairwise counterpoise corrections come to fifteen per cent more than the trimer's own. Three Gaussians a centre is a small basis, so the natural question was whether the fraction shrinks with a better one or stays put. It does neither. The correction falls by three orders of magnitude and the fraction grows fivefold.

Bonding models

An interior maximum a third orbital allows

Two orbitals related by a mirror plane give a Boys functional with no linear term, so it has no interior maximum and no angle is ever searched for. Add the oxygen lone pair and the search over a three-dimensional rotation group finds one — eight and a half degrees of lone pair mixed into two bent bonds, worth a further two and a half per cent that no pair of orbitals could reach.

Beyond the octet

Fifty descriptions of one molecule

A search whose largest basin takes ninety-eight per cent of its starts will report one description however long it is run, and nobody runs four thousand starts when the first fifty agree. Whether the rare ones are worse descriptions or merely rarer is one number per description, already computed and never looked at. On two of three cages they are not worse — they are the same answer, to parts per million.

Orbitals

One basis size where it is worth doing

The pairwise assembly's error grows with the basis — 5.25 per cent at one Gaussian a centre, 37.98 at six — while the correction itself falls by a factor of three thousand. Which of the two should a practitioner care about? Drawing a line at a kilocalorie a mole answers it: there is exactly one basis size at which the correction is worth computing and its pairwise assembly is accurate enough to use.

What is taught wrongly

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.

What the shape is for

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.

Beyond the octet

Counting was right except where it mattered

A cage whose localisation gives dozens of descriptions that are all the same answer raises a worry: if degeneracy is common across the family, counting descriptions is the wrong measure of ambiguity. Across forty-eight cage-and-filling pairs it is the right measure on eleven of the thirteen that have anything to count — and it fails on the one leaned on hardest.

Bonding models

The five figures were an identity

Two lone pairs satisfy the two-orbital mixing criterion by a margin of 1.732128, against √3 = 1.732051 — agreement to five figures on a number assembled from three integrals over a numerical grid. It is exact. The criterion is the reciprocal of the p amplitude of the outward hybrid, and everything else in the molecule cancels out of it.

What the shape is for

A contraction that cannot reach three of them

The angular overlap model's own derivation gives a π/σ ratio that disagrees with the fitted parameters by up to sixfold, and the metal's contraction is the obvious candidate to account for it. The whole range Slater's rules allow moves the ratio by a factor of two. Two ligands need charges far outside it, one needs a charge past where the overlap rule can be trusted at all, and two are unreachable at any charge because a quotient of squared overlaps cannot be zero or negative.

Bonding models

The other end of the bracket is not a number

Every claim can be priced by the measurement error it would take to overturn it, assuming the two errors independent, and the correlated extreme ought to turn each price into a bracket. At perfect correlation between equally precise measurements the difference is exact and the price is infinite — and every price is multiplied by the same factor, so the ordering of prices cannot change.

What symmetry decides

The variation was the basis

Solved in the five functions a field in one plane couples, the tilted Stark problem has one avoided crossing at every tilt and a crossing field that moves by a third across ninety degrees. Solved in the whole nine-function shell the field does not move at all — the same number at every direction, to eleven decimal places — and the thirty-four per cent was the truncation.

What is taught wrongly

The worst of the six was the one we asked about

The capacity's correlation with the second ionisation energy explains 2.9 per cent, and the natural next step is the same pair test against every other candidate input. Every input fails it — but the second ionisation energy fails it by a factor of twenty-four more than the best, and the test itself had to be repaired first, because the version the question implied reports the capacity failing to be a function of itself.

What symmetry decides

Two levels cannot make a minimum

The one minimum between coupled levels in a Stark shell sits at 0.0196 of the s–p gap, and the nearest two-state estimate at 0.0192 — two per cent away, which reads as the estimates having been aimed at the right feature all along. They were not. Two coupled levels only ever separate, so no estimate can be where its own pair is closest. The minimum belongs to a third level, exists only while the d level sits within a quarter of the s–p gap, and meets the estimate by crossing it.

What the shape is for

The correction that moves three of them backwards

Every ligand radial function in the angular overlap sweeps was a neutral atom's, while three of the five donors carry a formal charge. Giving each one the charge it actually has moves three of the five computed ratios — and moves all three away from the fitted parameter, none towards it. The whole window each donor's own oxidation states allow sits above the value it was meant to reach.

What is taught wrongly

A size the fit was not made from

Every input tested against the chemical capacity so far has been inside the cubic that produced it, so the search was constrained to fail. A size is not: a hydrogenic orbital at the effective charge Slater's rules give an atom's valence shell has an exact radius, on one scale for every atom, and the cubic knows three energies and no length at all. It fails the test by sixteen times the floor.

Where the atoms go

The mass nobody chose

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

Bonding models

The criterion that has never said no

The two-orbital mixing criterion is a tautology on a carbonyl's lone pairs, and it is tempting to contrast that with the double bond, where the criterion seems to discriminate. Both halves are wrong. The double bond's margin is wider, ethene's is infinite because symmetry puts both centroids at the same point, and the criterion has never returned a negative on any case it is usually put to — though the same three orbitals supply two pairs it refuses outright.

Beyond the octet

A second criterion left a gap too

The localisation spreads are bimodal — an empty factor of two hundred around the degeneracy threshold — and the gap might belong to the criterion rather than to the cages. Boys localisation leaves a gap of seven decades on the same forty-eight pairs, so it belongs to the cages. But the two criteria disagree about six of them, in both directions, and the family's most ambiguous cage under one is exactly degenerate under the other.

Beyond the octet

The cage is on both sides

Two localisation criteria classify six cage-and-filling pairs differently, and a symmetry explanation for the six is the natural first guess. The icosahedron's graph has a hundred and twenty automorphisms, the most in the family, and supplies three of the six disagreements and nine of the agreements. What the six do have in common is sharper than a symmetry: in every one, one criterion's spread is not small but zero.

Where the atoms go

An ordering worth half a per cent

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

All themes