When the molecule does not stop

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.

Worth reading first: Two bands, if the chain is short enough · A defect is a level in the gap.

An equal binary alloy has a gap that depends on how the two components are arranged, even though its mean and its variance cannot. Three arrangements of one composition: the ordered one’s gap is exactly 2δ2\delta, the segregated one’s exactly 2δ42\delta - 4, and the random one’s has no closed form and falls as the chain grows — from 0.185 at 512 sites to 0.088 at 32,768 — because a longer chain holds a longer run of like sites, and a longer run is a narrower sub-band.

That was an argument about where a band edge is. It said nothing about what is at the edge, and the two questions have different answers: a defect level, which is a state living on one site and its neighbours, and a band edge is ordinarily a state living everywhere.

The states at this edge are the first kind, and what they are is unusually clean.

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.
Fig. 1 The forty levels nearest the gap centre of five random alloys, against the length of the run of like sites each sits on. The line is a closed form with nothing fitted in it. Two of the forty are somewhere else, and they are the interesting ones.

What a run is, taken on its own

A run of LL like sites is a piece of pure material, and a piece of pure material with ends is a chain of LL sites. Its levels are the ones computed in the band limit: 2cos(kπ/(L+1))2\cos(k\pi/(L+1)), and its states are

ψk(j)sin ⁣(kπjL+1),\psi_k(j) \propto \sin\!\left(\frac{k \pi j}{L+1}\right),

which is the particle in a box, arriving here because a run really is one.

The quantity that says how much of a chain a state occupies is the participation ratio, (ψ2)2/ψ4(\sum \psi^2)^2 / \sum \psi^4 — one for a state on a single site, LL for a state spread evenly over LL. For a sine profile the two sums are elementary:

sin2=L+12,sin4=3(L+1)8,\sum \sin^2 = \frac{L+1}{2}, \qquad \sum \sin^4 = \frac{3(L+1)}{8},

so the participation ratio is 23(L+1)\tfrac{2}{3}(L+1), with no parameter in it at all. A box state of a run of thirteen sites occupies 9.33 of them, and that is a prediction rather than a fit.

What the alloy’s gap states are

Take a 400-site chain, half its sites raised by δ=4\delta = 4 and half lowered, arranged at random, and look at the eight levels nearest the middle of the gap. For each, find where its weight is centred and how long the run of like sites is at that place.

energy occupies run 23(L+1)\tfrac{2}{3}(L+1) ratio weight in the run
2.0522 9.16 13 9.33 0.982 100.0%
−2.1278 5.85 8 6.00 0.975 99.9%
2.1627 5.19 7 5.33 0.973 99.8%
2.2055 9.20 13 9.33 0.986 99.9%
−2.2128 8.30 6 4.67 1.779 65.1%
−2.2153 8.27 6 4.67 1.771 65.1%
−2.2921 3.88 5 4.00 0.971 99.6%
2.2921 3.88 5 4.00 0.971 99.6%

Six of the eight sit on the closed form to within three per cent, keep essentially all of their weight inside their own run, and are box states of a box the shuffle happened to make.

The measurement makes no reference to the disorder at all. Nothing here fits a localisation length, and nothing needs to know how the sites were arranged: the run is read off the chain, the closed form follows from its length, and the state matches.

Across five chains that is 38 of 40 gap states, with a worst departure of 3.1 per cent.

One state on one run, and one shared between two. The weight of a gap state on each site, drawn over the 46 sites around where it lives, with the sites of the raised kind shaded. The first occupies 9.16 sites on a run of 13 — the closed form says 9.33 — and keeps 100.0 per cent of its weight inside it. The second occupies 8.30 sites on a run of 6, which is 1.78 times the closed form, and keeps only 65.1 per cent inside: it is shared with a second run nearby, and it arrives with a partner split from it by a thousandth.
Fig. 2 Two of them drawn on the sites they live on, with the raised sites shaded. The first is a half-sine on a run of thirteen and stops where the run stops. The second is on a run of six and is nothing like a half-sine.

The two that are not

The pair at −2.2128 and −2.2153 occupies 8.3 sites on a run of six, which is 1.78 times what a box of that size holds, and keeps only 65 per cent of its weight inside.

They arrive as a pair, split by 0.0025. That is the shape of the answer: the state is shared between two runs close enough to each other that a state on one reaches the other, and two coupled boxes give a symmetric and an antisymmetric combination split by twice the coupling. It is the two-impurity problem, arriving unbidden because a random arrangement eventually puts two runs near one another.

So the answer to how localised are these states is: as localised as their run, until there are enough runs that they start to find each other, at which point they stop being defect levels and start being a defect band. On a chain of 400 that has happened to one pair out of eight; on a chain of 32,768 with the same composition it will have happened to a great many more, and the sub-band edge computed from the longest run alone will have acquired structure the longest run does not describe.

Two impurities at 2 to 16 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.
Fig. 3 The same physics on purpose rather than by accident: two impurities in a chain, and the splitting between the symmetric and antisymmetric combinations against how far apart they are. A pair split by a thousandth is a pair about nine sites apart.

Why a run and not a cluster

There is a question hiding in the phrase the run the state sits on, and it is worth answering before the closed form is trusted anywhere.

A run of like sites in a one-dimensional chain is unambiguous: walk left until the site energy changes, walk right until it changes, and the stretch between is the run. That is a definition with no parameter in it. It is also a definition that only exists in one dimension — in two, the sites of one kind form a connected region of some shape, and the box it makes is no longer a length but an outline whose levels depend on its whole geometry.

So the closed form here is not a general fact about disordered alloys. It is a general fact about one-dimensional ones, and it holds because the only property a one-dimensional region has is its length.

What survives in higher dimensions is the argument rather than the formula: a connected region of one component is a piece of that material with a boundary, its levels are the levels of that region, and a level near the middle of the gap belongs to a large region. The gap edge is still set by the largest region and the states there still live on it. What is lost is the ability to write down where the level is without solving the region, which is exactly what makes the one-dimensional case worth having: it is the case in which the whole account is arithmetic.

There is one more thing the definition buys. Because a run is found from the site energies alone, the two halves of every row in the table above come from different objects — one from a list of site energies, the other from an eigenvector — and neither computation can see the other. A stabilisation measured from somewhere is this collection’s standing warning about quantities computed against references chosen after the fact; there is no such freedom here, because the reference is an integer read off the chain.

The contrast does nothing, which was not expected

A gap state is a bound state, and a bound state in a deeper well ought to be more tightly bound. Turning the contrast up from δ=2.5\delta = 2.5 to δ=8\delta = 8 takes the gap from 1.11 to 12.10 — a factor of eleven — and the weight the state keeps inside its own run goes from 91.3 per cent to 91.0.

It does not move. The prediction was wrong and the reason is worth having.

The barrier confining the state is not the depth of the well. It is the ends of the run, where the next site is of the other kind, and a site of the other kind is off-resonance by 2δ2\delta — which is already large compared with the hopping at the smallest contrast that opens a gap at all. Making it larger changes a barrier that was already effectively complete.

What the contrast does change, slightly, is the participation ratio: 9.04 sites at δ=2.5\delta = 2.5 and 9.25 at δ=8\delta = 8, creeping up towards the box value of 9.33 from below. A state in a finite well is slightly compressed relative to the infinite one, and a deeper well compresses it less. The whole effect is two per cent.

A deeper well does not hold the state any better. How much of a gap state's weight stays inside the run that carries it, against the contrast between the two components. It sits at 91.3 per cent at δ = 2.5 and 91.0 at δ = 8, across a gap that has grown from 1.11 to 12.10. The confinement is set by the run's ends, where the neighbours are of the other kind, and that barrier is already complete at the smallest contrast that opens a gap at all.
Fig. 4 The weight a gap state keeps inside its own run, against the contrast between the components. A gap eleven times deeper, and a flat line.

The control, which had to be the ordered arrangement

A participation ratio of nine on a chain of four hundred is only evidence of localisation if the same measurement returns something large for a state that is not localised — and the obvious control does not work.

Every state of a disordered one-dimensional chain is localised. Taking a level from the middle of a sub-band of the same random alloy and measuring it returns a small number too, which proves nothing about the gap state in particular.

So the control is the ordered arrangement of the same composition: the same site energies, the same δ\delta, in alternating order rather than shuffled. Its band-edge state occupies 133.7 sites of 400 — a third of the chain, and growing with the chain, which is what an extended state does.

That contrast is the check, and it is the one that makes the participation ratio a measurement rather than a formula: the arrangement decides, and it decides by a factor of fifteen.

One composition, three arrangements, three shapes — at a contrast of 4. The density of states of a chain of 300 sites, half of them raised by 4 and half lowered by the same, arranged three ways: alternating, at random, and segregated into two blocks. All three have the same mean and the same variance to nine decimal places — 17.993333 — because both are traces and a trace does not know which site is where. The shapes are not alike, and the ordered arrangement is the only one of the three with a gap at this contrast.
Fig. 5 The three arrangements, at the contrast used here. The ordered one has a clean gap; the random one has a gap with things in it, and this essay is about what the things are.

What this says about a large alloy

The random alloy’s gap, computed from the longest run, leaves open what a very large alloy is best described as. The answer supported here is:

A split band with a defect band in the gap. The two sub-bands are the two components’ bands, narrowed — and narrowed rather than shifted, which is the distinction a band’s width being a count of neighbours makes exactly: a run of six has fewer neighbours available at its ends than a bulk site does. Between them are box states of runs — one level per run per box quantum number, at energies set by the run’s length through the closed form — and as the chain grows the runs become numerous enough that neighbouring ones couple and the levels smear into a band of their own.

That is not a description anybody draws, and the reason is that it is three descriptions at once: a band structure, a set of defect levels, and a coupling between them. Each of the three is something computed separately, and the alloy is the case where all three are present in one system with no parameter distinguishing them.

Why the gap closes as slowly as it does

The box picture explains more than the states; it explains the rate at which the gap shrinks, and the rate is the distinctive part of the whole phenomenon.

A run of LL like sites carries levels at 2cos(kπ/(L+1))2\cos\bigl(k\pi/(L+1)\bigr), and the one nearest the gap is the run’s own band edge, at 2cos(π/(L+1))2\cos\bigl(\pi/(L+1)\bigr). Expanding that gives

2π2(L+1)2+,2 - \frac{\pi^2}{(L+1)^2} + \dots,

so a run’s edge approaches the infinite-run limit as the inverse square of its length. The gap of the whole chain is set by the longest run it happens to contain, so the gap approaches its limiting value at that rate too.

Now ask how the longest run grows. In a random chain at equal composition, a specified run of LL like sites occurs with probability 2L2^{-L} per starting position, so the longest one in a chain of NN sites is near log2N\log_2 N — a length that grows with the logarithm of the chain rather than with any power of it.

Putting the two together: the gap’s excess over its limit falls as 1/(log2N)21/(\log_2 N)^2. That is an extraordinarily slow decay, and it is the reason the measured gap behaves the way it does. Between 512 sites and 32,768 the chain grew by a factor of 64 and the longest run only from about nine sites to about fifteen, so the predicted shrinkage of the gap’s excess is somewhere between 2.5 and 3.1 depending on how the run length is counted. The measured factor is 2.10.

Two points cannot separate that from a weak power law, and the agreement in magnitude is the most the argument can claim on this evidence. What the argument does supply, and no fit could, is a reason for the gap to be so reluctant to close: doubling the chain adds one site to its longest run, and one site to a run of fifteen is a small change to a quantity that already enters as an inverse square.

It also says the limit is approached and never reached. A chain of 102310^{23} sites has a longest run near seventy-six, and its gap is still above the segregated arrangement’s by a term of order π2/762\pi^2/76^2. A random alloy in the thermodynamic limit has the segregated arrangement’s gap, and any real sample has a slightly larger one, by an amount that depends on how big the sample is — which is the sort of size dependence that is easy to mistake for a preparation effect.

What is quoted, and what is computed

Nothing is quoted. There is no material here: a chain of 400 sites, half of them raised by a stated amount, arranged by a seeded shuffle. Every number is computed from that.

The closed form 23(L+1)\tfrac{2}{3}(L+1) is derived from two elementary sums and checked against the participation ratio of an actual eigenvector, which is the only sense in which it is verified — the two are computed by different routes and must agree.

The runs are found by walking the chain from the state’s centroid, which is a property of the site energies alone and knows nothing about any eigenvector. The state and the run are measured independently, which is what makes their agreement evidence.

What this cannot say

It is one dimension. Localisation in one dimension is special: every state is localised at any disorder, which is why the control had to be the ordered arrangement rather than a state from the middle of a band. In three dimensions weak disorder leaves states extended and the whole question changes shape.

There are no electrons in it. These are one-electron levels, and which of them is occupied has not been asked. A gap state that is filled and a gap state that is empty are very different things for a material, and what a filled band is and is not is where that question lives.

And the runs are a caricature of an alloy. A real solid solution has short-range order — a tendency for like or unlike neighbours — so its runs are not the runs of a fair shuffle. That changes the distribution of run lengths and therefore the distribution of gap levels, and it does not change the closed form for any one of them.

The state on a defect at site 31 of 60. The amplitude of one eigenvector at each site of a 60-site chain whose site 31 has its energy raised by 1.6β. The two colours are the two signs. A state belonging to the whole chain would reach across it; this one does not.
Fig. 6 The simplest case: one site changed, one level pulled into the gap, and a state on that site and its neighbours. Everything above is that picture with the single site replaced by a run the arrangement made.

What was checked

Most of the gap states are box states of their run, stated as a majority of a census over five chains rather than as an observation about one.

And for those, the closed form needs no fitting — a worst departure of 3.1 per cent between a participation ratio computed from an eigenvector and one computed from an integer.

Some exceed it, which is checked as well: a census in which every state matched would mean the shuffle had been too kind to produce the shared pair, and the shared pair is half the finding.

The ordered arrangement’s band-edge state occupies more than eight times as many sites and more than a twentieth of the whole chain — the control, without which a small number would be arithmetic rather than evidence.

And the weight inside the run does not depend on the contrast, checked as a bound on the variation — two points across a factor of three in δ\delta — because the prediction that it would rise was wrong and a check written to the prediction would have failed.

The gap a long enough chain does not have. How far the random alloy's gap sits above the segregated alloy's closed form of 2δ − 4, at a contrast of 4, against the length of the chain it was measured on. It falls as the chain grows — from 0.1849 at 512 sites to 0.0876 at 32768 — and the open marks are what a run of the longest stretch of like sites the chain happens to contain would put at the band edge. The mechanism is that a longer chain holds a longer run, and a longer run is a narrower sub-band.
Fig. 7 The instrument that decides whether a separation in a finite spectrum is a gap at all: the same alloy at four lengths, asking whether the split survives as the chain grows or falls away as a level spacing. It survives, and that is what entitles everything above to call it a gap — a single computed spectrum could not have told the two cases apart.

Still open: the gap’s density of states in closed form

The obvious open question is the distribution rather than the individuals. A run of length LL appears in a random chain with a known probability, and each run contributes levels at energies the closed form gives — so the density of states inside the gap is computable in closed form from the run-length distribution, with no diagonalisation anywhere. Comparing that against the levels an actual diagonalisation puts in the gap would be a strong test of the whole picture, and it would say what fraction of the gap states the box description accounts for on a chain far too large to solve.

The nearer question is the coupling that made the shared pair. Two runs share a state when they are close enough, and close enough is a distance already measured for two impurities — the splitting against separation is drawn above. Working out at what chain length the typical gap between neighbouring long runs falls below that distance would give the size at which the defect levels become a defect band, as a number rather than as an expectation, and it is the size at which the longest-run account stops being the whole story.

What links here

Computed from the collection rather than written here: the essays that point at this one.

Reads more easily once this is understood

Essays that name this one as worth reading first.

Shares its objects with

Essays naming at least two of the same things, that neither author linked.

Named objects

A dashed tag is an object no other essay names yet.

Bands in a solidBand gapClosed formDefect levelDegeneracyDisorderEigenvalueLocalisationModel limitOne-electron modelsParticipation ratioProbability densityThermodynamic limitTight-binding models