Concept

Direct methods — where it appears

Estimating phases from amplitudes alone, using the fact that a crystal's density is non-negative and made of atoms. The constraints are probabilistic, they weaken as the square root of the number of atoms, and using them is a search rather than a calculation.

Named by 6 essays across one field — each of them below, with the objects they name alongside it.

The phase problem. The same structure rebuilt from its diffraction three ways: with the true amplitudes and phases, with the phases scrambled, and with the amplitudes discarded but the phases kept. The atoms survive the loss of the amplitudes and do not survive the loss of the phases, which is what an experiment throws away.

The phase problem

A detector records how much light arrives and not when it arrives, so half of every diffraction measurement is thrown away before it is written down. The half that is lost turns out to be the half that carries the structure.

diffraction · Diffraction
⟨cos Φ⟩ against κ, 379 triplets. The mean cosine of the triplet, binned by the concentration κ = 2|E₁E₂E₃|/√N, for the 379 triplets of a structure of 24 atoms whose reflections all exceed |E| = 1.2. The curve is Cochran's I₁(κ)/I₀(κ), computed from the distribution and not fitted to anything; the points are measured, with the number of triplets in each bin printed above. They agree to 0.1 root-mean-square. The measured points sit slightly above the curve throughout, which is the finite structure showing: Cochran's derivation assumes atoms placed at random and there are only 24 of them.

Three phases that do not move when the origin does

A phase is a property of the description, not of the crystal: shift the origin and every one of them changes. A sum of three phases whose indices add to zero does not change, because the shifts cancel. That sum is the smallest thing about a structure that a diffraction experiment could in principle know, and it is not distributed at random.

diffraction · Direct methods
the most consistent answer has 100 per cent of the signs. 60 runs of the sign procedure from 60 different random starts, each plotted at its self-consistency — a figure computed without any knowledge of the answer — against the fraction of its signs that are in fact right. Throwing out the 1 run that reached the uniform solution — every sign the same, perfectly consistent and physically a single peak — the highest consistency belongs to a run with 100 per cent of the signs right. The ranking works here, and the reason it works is that the cell is small. Nothing in the plot's horizontal axis knows the answer, which is the only reason a procedure of this kind is a procedure at all.

The formula that has the answer already

The tangent formula rebuilds each phase from all the others, and the true phase set is very nearly a fixed point of it — hand it the answer and it hands the answer back. Start it anywhere else and it does not arrive. Having a fixed point and finding it are different problems, and the second is where the subject spent twenty years.

diffraction · Direct methods
A map from amplitudes alone, 0.59 grid steps out. The density after 150 cycles of flipping, with the atoms that produced the data drawn as rings — moved into the origin and the handedness the solution chose, because a phase set does not fix either and comparing without allowing for them measures the arbitrariness of the description. Every peak of the map is an atom and every atom has a peak. Nothing about the arrangement went into the calculation: the input was a list of amplitudes and a random set of phases.

The solver that knows no symmetry

Compute a map from amplitudes and random phases, reverse the sign of everything below a small threshold, transform back and keep the phases. Repeat. The structure appears — and so does its space group, which was never supplied.

diffraction · Direct methods
The sign turns when the cross terms go weak. Every quartet among the strong reflections of a small structure, sorted by the mean of its three cross terms, with the mean cosine of the phase sum in each bin. Where the cross terms are strong the quartet behaves like a triplet and the cosine is near one. Where they are weakest it is -0.698 — negative — and 93 per cent of those quartets have a cosine below nought. The four reflections of the quartet are equally strong in every bin; what changes is three reflections that are not in the sum at all.

The relation that can say no

Every phase relation before this one pushes a sum towards zero, so none of them can contradict another — a phase set satisfying all of them badly is still satisfying them in the same direction. A quartet can be estimated at π instead, and it is when its three cross terms are weak, so the information arrives from the reflections nobody would have thought worth measuring.

diffraction · Direct methods
An equality, with a known factor in it. Sayre's identity says the structure factor is proportional to the convolution of the structure factors with themselves — an equality, not a probability — with a factor that depends on the atoms and the resolution and on nothing else. Computed for a structure of Gaussian atoms, the ratio of the two sides falls with resolution exactly as predicted: the logarithm of the ratio is a straight line in the squared index, its slope is minus half the width in the atomic factor, and the ratio is very nearly real and positive, so the identity relates the phases and not only the amplitudes.

The relation that is an equality

Every relation of this kind so far is a probability — right nine times in ten, useless applied once and decisive applied ten thousand times. One is not. For a structure of equal, resolved atoms the squared density has peaks in the same places, so every structure factor is exactly a convolution of all the others, with a factor that depends only on the atom. It is exact, and on its own it is useless.

diffraction · Direct methods

Named alongside it

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

Phase problemStructure factorTripletFigure of meritNormalised structure factorStatisticsStructure invariantElectron densityRefinementCase analysisCharge flippingFourier synthesis

All concepts