Concept

Correlation energy — where it appears

The exact energy less the best single-determinant one, which is a difference rather than a property. Its magnitude grows without limit with the repulsion, because a mean field cannot stop paying for coincidences a spread-out density forces.

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

Two is the ordinary answer, and one is a different kind of correlation. The local exponent of the correlation energy in the repulsion, against the repulsion, for three systems at half filling. The two closed-shell systems tend to two as the repulsion vanishes, which is ordinary perturbation theory. The ring of four tends to one, at repulsions fifty times smaller than the hopping.

A third kind of correlation

The two-site model gave an exact identity — the power of the correlation energy in the repulsion equals the occupation of the bonding natural orbital — and asked whether anything like it survives with more orbitals. It does not, and the way it fails is better than the identity was: a ring of four gives a power of one where every closed-shell system gives two, at repulsions fifty times weaker than the hopping, because its reference was never a single state.

beyond · Correlation
The error, against the repulsion it is an error about. The energy a mean field misses, for two electron counts on 4 sites, against the strength of the repulsion. Each is a straight line at large repulsion and the dashed line through it is not a fit: its slope is the count of coincidences a uniform density forces, computed from the electron number and the site number alone.

A mean field cannot get out of the way

The energy a mean field misses grows without limit as the repulsion rises — 27.82 at U = 32 for four electrons on four sites, against an exact energy of −0.2946, so the error is ninety-four times the answer. The rate it grows at is not an energy at all: it is n²/4N, a count of the coincidences a spread-out density cannot avoid.

beyond · Correlation
One exact state, two correlation energies. The energy each of two mean fields misses, against the repulsion, for one system whose exact energy is a single smooth curve. Below U = 2 the unrestricted search returns the restricted answer and the two definitions agree to the last bit. Above it they part: at U = 32 the restricted reference reports -27.82 and the unrestricted one -0.11, a factor of 259.43, and one is growing while the other falls.

The reference decides the correlation

The correlation energy of one exact state, measured against two references that are both called Hartree–Fock, is −27.82 and −0.107 at the same repulsion — a factor of 259, on a system whose exact energy is a single smooth curve. Below the instability at U = 2 the two agree to the last bit; above it one grows without limit while the other falls, and the reference that reports almost no correlation has ⟨S²⟩ = 1.99 where a singlet is zero.

beyond · Correlation
Nearly all of the error cancels, and the answer gets worse. For each repulsion: the error a spin-paired mean field makes in the total energy of one four-site system and of two two-site ones with the same number of electrons, and the error left in the difference between them. The cancellation improves from 83 to 97 per cent along the axis. The residue as a share of the quantity being computed goes the other way, from 1 to 423 per cent, because the reaction energy shrinks faster than what survives.

Two wrong numbers and a right difference

A mean field gets the total energy of a four-site system wrong by 12.11 and of two two-site systems wrong by 12.49, and 96.9 per cent of that error cancels out of the difference between them. The residue is 0.38 — and the reaction energy it is a residue of is 0.09, so the cancellation improves and the answer gets worse at the same time. Change the pair being compared to a singlet and a triplet and nothing cancels at all: the sign goes.

wrong · Approximation
A correction that stops belonging to the system it is added to. The exact ground state of a four-site Hubbard ring, the unrestricted mean field's, and the composite: the mean field plus the correlation correction computed on the symmetric molecule. At ε = 0 the two systems are the same one and the composite is exact. As the sites are made unlike, the transferred correction stops being the right one — the exact correlation energy shrinks towards nothing while the transferred number does not — and the last rows are the recipe adding a correction almost as large as the error it is meant to remove.

The correction that was computed somewhere else

Every composite method rests on one assumption: that an expensive correction computed on a small case can be added to a cheap calculation on a large one. Tested on four sites where both answers are exact, it removes 99.74 per cent of the cheap method's error near the reference and −1450 per cent of it further away — at which point the recipe is fifteen times as wrong as the calculation it was improving.

wrong · Approximation
Where the electrons are, without subtracting anything. The opposite-spin pair distribution of a half-filled ring of 6 at six repulsions, by separation, each divided by what uncorrelated electrons of the same density would give. At no repulsion it is one everywhere; at a repulsion of 16 the chance of finding two electrons on one site is 0.0430 of that, and what is missing has turned up next door. Nothing here is a difference between two calculations.

Where the electrons are, without subtracting anything

A correlation energy is an exact energy less a mean-field one, so the answer depends on which mean field was picked — by a factor of 259. The pair distribution says the same thing with no subtraction in it, and two systems matched to the same correlation energy turn out to have electrons in visibly different places.

beyond · Correlation
A cheap number that predicts an expensive failure. Thirty-two systems. Along the bottom, how far the mean field's own symmetry breaking has moved between the reference and the target — a quantity available before any exact calculation. Up the side, how wrong the transferred correction turns out to be. They rank together at 0.902, and the open marks are the systems whose broken solution has collapsed entirely, which is where the diagnostic stops being a scale and becomes a warning.

The warning a cheap calculation gives

A correlation correction computed on one system and carried to another works until it does not, and nothing in the scheme says in advance which. The mean field's own symmetry breaking says: it collapses at a definite field, and the transfer fails where it goes. Across thirty-two systems the two rank together at 0.90.

wrong · Approximation
The same hole, priced three ways. The correlation hole of a ring of 6 at a repulsion of 8, weighted by three interactions. With an on-site interaction the answer is 100 per cent at separation zero — as an identity, since the interaction is zero everywhere else. With one that reaches a neighbour, the enhancement at separation one costs rather than pays, and gives back 44.0 per cent of the on-site saving; with a Coulomb tail, 46.9. Everything beyond one neighbour is worth under a twentieth of the on-site term.

Half of it is given back at one bond

The correlation hole of a ring removes 1.2653 pairs from zero separation and puts 1.1126 of them one site away. With an on-site repulsion that second number costs nothing, because the interaction is zero there — but with anything that reaches a neighbour it costs 44 per cent of what the hole saved, and with a Coulomb tail 46.9.

beyond · Correlation
A ring of 6 as the neighbour repulsion is turned up. At an on-site repulsion of 8, three quantities against the nearest-neighbour repulsion: the alternating structure factor, the double occupancy, and the nearest-neighbour opposite-spin pair distribution. The rise is steepest at V = 4.5, which is 0.563 times the on-site repulsion. Far past it the ring is charge ordered — nearly every electron paired on alternate sites, which is what a double occupancy approaching a half means.

The give-back that turned into a saving

Take a wavefunction optimised for an on-site repulsion, weight its correlation hole with an interaction that reaches one neighbour, and the enhancement at one bond gives back forty-four per cent of the on-site saving. Putting that neighbour term into the Hamiltonian and solving exactly does not shrink the give-back. It reverses its sign.

beyond · Correlation
Where the broken solution appears. The mean field's spin polarisation against the on-site repulsion, at half filling and no site-energy modulation, for three systems. Two of them are symmetric below a threshold and polarised above it — 1.672 for a chain of four and 2.355 for a ring of six. The third is polarised at every repulsion tested, because its half-filled shell is degenerate and the symmetric solution is unstable however small the repulsion is.

The half of the square a ring of four cannot show

There is a warning that says in advance whether a transferred correction will hold: the mean field's own symmetry breaking, which collapses at a definite site-energy modulation and takes the transfer with it. The other axis of the same square has a threshold too — on the other side — and the ring of four it was all measured on is the one system with no threshold to find.

wrong · Approximation
One sign change is a root and the other is not. The solved give-back against the neighbour repulsion on a ring of 6 at U = 8. It crosses zero at V = 3.895, where the price of everything beyond contact really is nothing, and changes sign again at V = 4.159, where the quantity it is a fraction of has vanished instead. The curve is broken at the second because it is an asymptote and not a crossing; four of the 12 points fall outside the band drawn here.

A sign change is not always a zero

The solved give-back changes sign somewhere between a neighbour repulsion of two and one of four, and a bisection looks like the way to find the value where the structure beyond contact contributes exactly nothing. It changes sign twice. One crossing is that value; at the other the quantity the fraction is a fraction of has vanished instead, and a bisection reports the two in identical words.

beyond · Correlation
The second number is the first one, rearranged. The composite's error against the change in the correlation energy, at every point on both axes of the square. They lie on the diagonal because they are the same quantity: the composite is the target's mean field plus the reference's correlation energy, so its error is the reference's correlation energy minus the target's. The largest departure across 20 points is 2.2e-16, which is the arithmetic's own precision and not a measurement.

The second number is the error, rearranged

A cheap diagnostic for a composite method leaves a scatter it cannot explain, and the number that ought to close it is the change in the correlation energy, already computed at every point, so the test is arithmetic rather than a calculation. It is arithmetic, and the arithmetic is the answer. The composite's error is that change with a sign on it.

wrong · Approximation
Five candidates, five failures, five different places. Every point of the square, with each cheap diagnostic's failing pair joined by a line. The five tested candidates fail on five different pairs involving 10 different points — no line shares an end with another. Had they all failed on one corner the lines would have converged on it, and the honest conclusion would have been that composites are safe away from that corner.

Five failures in five different places

Seven quantities a mean field produces for nothing have been tried as diagnostics, and none of them is usable. The question left is whether that is one finding or seven — whether the same awkward corner of the square breaks every candidate, or each is broken somewhere else. Each is broken somewhere else. Five candidates, five failing pairs, ten systems, and not one of them appearing twice.

wrong · Approximation
Eighty systems, and only sixteen of them on an arm. The composite's error at every combination of repulsion and staggered site energy for a four-site ring, with the correction transferred from one reference — the point at the top left, where it is exact by construction. The area of each mark is the error. Every system the diagnostics have previously been tested on lies along the top edge or the middle column; everything else is new, and it is most of the square.

The arms were the kindest part of the square

Seven cheap quantities have been tried as diagnostics for when a transferred correlation correction will fail, and every one of them was tested on systems lying along one of two arms — a repulsion swept at zero site energy, or a site energy swept at one repulsion. Filling in the square makes every one of them worse, by up to a factor of twelve, so the diagonal is not a region where they might work but the region where they fail hardest. And a grid asks a question a cross cannot: whether a candidate even sorts the systems the way the error does. One of them does not.

wrong · Approximation
A ratio is not a score until the gap it was taken at is beside it. Each candidate placed by how far its nearest comparison on the other axis is — the median gap in diagnostic value — and by the worst error ratio it produces. The bottom left is where a working diagnostic sits: close comparisons and ratios near one, which is where the control sits. The two candidates the banded rule could not score at all sit far to the right, so they were blank because they separate the two families rather than because they are good.

A blank is not a pass

Two of the seven cheap diagnostics for a transferred correlation correction could not be scored at all, because a scoring rule that only compares systems whose diagnostic values agree to fifteen per cent finds no such pair for them. The census reported them as blanks, which is correct and reads as the best entry in the table. Dropping the band and reporting the gap beside the ratio scores both — and the two turn out to be unscoreable for opposite reasons, one because it separates the two families as a monotone function of something that already failed, and one because it is unrelated to anything.

wrong · Approximation
The conventional start misses the lowest mean field on one band, and nowhere else. A four-site ring at every combination of repulsion and staggered site energy on a 12 by 12 grid. A small dot is a system where the conventional broken-symmetry start reaches the lowest mean-field solution; a circle is one where it lands above it, with area proportional to how far. There are 9, every one with a site energy one or two below the repulsion, and the largest miss is 1.124 at a repulsion of six and a site energy of four.

The worst system in the square was the solver's

Every score in the transferability square rests on an unrestricted mean field, and every one was computed from a single conventional start. On a band of nine systems, each with a site energy one or two below the repulsion, that start lands above the lowest solution the method has — by 1.12 at a repulsion of six and a site energy of four, which was the largest error in the whole square. Recomputed on the lowest field, the factor of forty-six that first showed the polarisation to be no rule falls to 2.61 and passes; the grid still fails every candidate, by a factor of fifteen rather than forty-eight; and the typical comparison, which no worst case could show, is far better than chance for every candidate but one.

wrong · Approximation
Two contributions of opposite sign, and where they cancel. The two terms of the give-back's numerator on a six-site ring, with the interaction reaching two bonds at half the strength of one. The first term is what the correlation hole returns at one bond and the second is what it costs at two — and they have opposite signs, because penalising pairs at two bonds depletes them where a one-bond interaction enhanced them. The numerator vanishes where the two are equal and opposite, which is not where either is zero.

A cancellation between two separations

The give-back's numerator has one term in it, so its root is the repulsion at which the pair count at one bond returns to its free value. Let the interaction reach two bonds and the numerator has two terms with opposite signs, and its root becomes a cancellation between an enhancement at one separation and a depletion at the next — at a repulsion where neither is at its free value. Which side of the old zero the new root falls on reverses between a reach of a third and a half, because whether two bonds are enhanced or depleted is the interaction's own doing.

bonding · Correlation
Where the correction comes from moves the worst case more than the start did. The polarisation difference's worst error ratio on the square's grid, with the correlation correction carried from eight references: the published one, three beside the band with no site energy, and the four systems where the two mean-field starts disagree. On the lowest mean field it runs from 9.7 to 213.2; the start alone moved the published reference's from 48.0 to 14.7. Rings are the conventional start; the dashed line is the census threshold of three.

The worst case measured the reference

Every score on the transferability square carries one correlation correction, fitted at a repulsion of eight and no site energy. Carried instead from seven other systems — three beside the band where the mean field's start misses, four on it — the polarisation's worst error ratio on the lowest mean field runs from 9.7 to 213, where fixing the start had moved it only from 48 to 14.7. Its median stays between 1.21 and 1.55. A worst ratio divides by the smallest comparable error it can find, and the reference decides where that is. And a reference fitted to the wrong solution passes the census at 2.05 while every system's error is above a half.

wrong · Approximation
Three terms in the numerator, and one place they cancel. The three contributions to the give-back's numerator on a six-site ring whose interaction falls as one over the distance round the ring — a half at two bonds, a third at the antipode — with their sum. The sum crosses zero once, at V = 5.9198, where the one-bond and antipodal terms are negative and the two-bond term positive: the two-bond enhancement is paying for depletions on either side of it.

A root that moves when nothing else does

Give the ring of six an interaction that falls off as one over the distance and the give-back's numerator has three terms. It still has one root, and at that root every separation's pair change has the opposite sign from the two-bond case. But the root has a weakness the earlier essays never tested: lower the interaction by the same amount at every separation and not one state changes, while the root moves — by as much as the reach of the interaction did, and out of existence altogether if the zero is put at one bond.

bonding · Correlation

Named alongside it

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

Exact diagonalisationHubbard modelReference stateElectron correlationModel limitOn-site repulsionSymmetry breakingHartree–FockComposite methodMany-electron wavefunctionsMean-field approximationTransferability

All concepts