The collection

Every essay — page 15

Page 15 of 36, continuing through the fields in the same order.

Orbitals Where the atoms go Bonding models What symmetry decides Beyond the octet What a spectrum settles When the molecule does not stop What the shape is for What is taught wrongly Series Named objects Orbitals Refutations Search

Orbitals

A one-electron wavefunction, drawn as a contour surface at a level somebody chose. Nodes, signs, and what the picture is a picture of.

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.

8 figures

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.

7 figures

A ninety per cent surface with a neighbour beside it. The contour enclosing 90 per cent of a one-electron ion's density at an effective charge of 1.6, drawn with no field as a circle and in a field of 0.05 atomic units as the closed curve. The surface moves out by 47.3 millibohr on the side the field pulls the density towards and in by the same amount on the far side — 2.84 per cent of its own radius. What the sphere encloses does not change to first order; only where the surface is does.

The surface a neighbour moves

An ion in a crystal sits in the field of the ion next to it, and that field moves its contour. The displacement has a closed form, it is checked against the polarisability it implies, and it turns out to be almost perfectly anti-correlated with the discrepancy it was proposed to explain — the pairs that need the most correction get the least.

7 figures

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.

6 figures $The two halves change places. What fraction of the counterpoise correction belongs to the lighter of two unlike atoms, against their separation. At a bonding distance it is 1.26 per cent — essentially the whole correction is the heavier atom's — and by 9.0 bohr it is 67. The two change places at 6.29 bohr. A symmetric pair's share is exactly a half everywhere, which is what makes one frozen number a complete description there and nowhere else.$

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.

6 figures

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.

7 figures

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.

6 figures

The trimer's correction, and the sum of its pairs. The counterpoise correction of a three-fragment system computed directly — each fragment's energy alone less its energy in the whole trimer's basis — against the sum of the three pairwise corrections, on a logarithmic axis. The sum is the larger everywhere the difference is above the solver's noise: 30.4 per cent at 1.6 bohr and nothing by six.

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.

6 figures

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.

6 figures $The correction collapses and the error made assembling it grows. At a separation of 2 bohr, two quantities against the number of Gaussians a centre. The trimer's own counterpoise correction falls by a factor of 2717 from one function to six — a bigger basis has less to borrow. The fraction by which summing the pairwise corrections overshoots it rises from 5.3 per cent to 38.0. Improving the calculation makes the assembly proportionally worse.$

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.

6 figures $Both quantities, and the line they have to be read against. At the closest separation, the counterpoise correction itself and the error a pairwise assembly of it makes, against the number of Gaussians a centre. Both fall — the correction by a factor of 3631, the error by 793 — and the fraction rises by exactly the ratio of those two. The dashed line is a kilocalorie a mole. The only basis where the correction is above it and the error below it is two.$

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.

6 figures

The residue is four times below the model's own error. The quantity whose sign is wanted, beside the accuracy of the numbers it is a difference of. The measured residues are 0.031 and 0.056 ångström; the model gets a single separation right to 0.242 on average. A fourfold alternating difference of quantities known that badly cannot resolve something that small, and that arithmetic was available before any of this was computed.

The residue is below its own noise

Two interaction terms, both negative, leave the sign a coin toss — and a larger block would settle it. There is a larger block: a model that needs no measured separations supplies twenty-one. It cannot settle anything, because a fourfold alternating difference of distances known to a quarter of an ångström cannot resolve three hundredths of one.

6 figures