Field

How it is known

A crystal's structure is read from the pattern it scatters. The reciprocal lattice, systematic absences, and what an experiment can and cannot see.
A lattice and its reciprocal. The reciprocal lattice is where a crystal scatters. A long axis in the crystal gives closely spaced spots and a short one gives widely spaced spots, so the diffraction pattern is the structure turned inside out.

The reciprocal lattice

Nobody has seen a space group. Crystals are read from where they scatter, and where they scatter is a second lattice in which long has become short and short has become long.

What pg scatters. The diffraction pattern computed from the atom positions alone. Where a glide plane is present, alternate reflections along a row cancel exactly — and those missing spots are how the glide is identified in an experiment, since the glide itself is never seen.

Systematic absences

The most informative part of a diffraction pattern is the part that is not there. A glide plane cancels alternate reflections along a row, exactly, and those missing spots are how a symmetry nobody can see is identified.

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.

168 reflections, 26 rings, one axis. Every reflection of a square lattice out to indices of 6, drawn as a point of the reciprocal lattice with the rings of equal length it lies on. A single crystal gives each point its own place on a detector; a powder is every orientation at once, so what is measured is the ring rather than the point, and the 168 reflections arrive as 26 lines on the single axis at the right, with height the number sharing each. The rings are found by sorting the lengths and grouping equal ones, and every point is checked to sit on the ring it was given. 1 of the rings are marked: they hold reflections that no operation of the lattice's point group relates, so the collapse has added together two genuinely independent measurements and no experiment can take them apart again.

What a powder pattern loses

Grind a crystal up and every orientation is present at once, so a two-dimensional pattern of spots collapses onto a single axis. Reflections that had their own places arrive together, and some of the coincidences are exact and have nothing to do with symmetry.

What p3 scatters, and what the scattering shows. The structure on the left has point group 3, of order 3. The intensities it scatters, on the right, have point group 6, of order 6 — more symmetric than the thing that produced them. Reversing the sign of both indices conjugates every term in the sum and leaves the modulus alone, so a diffraction pattern always acquires a centre of symmetry, and in the plane a centre is a half turn. Both numbers are measured: the left from the operations, the right by testing each candidate against the computed intensities.

The symmetry diffraction adds

A diffraction pattern is always more symmetric than the crystal that made it. The extra symmetry is not a mistake in the experiment and no care removes it — it is a property of what a detector records, and it collapses the seventeen groups onto six.

p4m, displaced and then measured. Every atom of a p4m pattern moved by up to 1.2 per cent of a cell edge, and the resulting positions examined at 26 tolerances. At zero tolerance only the identity survives, so the structure has no exact symmetry whatever. Between 0.021 and 0.094 the count sits at 8, which is the group that was displaced. Above that it climbs to 15, accepting operations no version of this pattern has. The climb is not even steady: at 4 of the 25 steps the count falls as the tolerance is loosened, because operations accepted separately at one threshold merge into one at the next. The correct answer is a step on a staircase and nothing in the coordinates says which step.

Near-symmetry, and the tolerance that is not here

Every claim on this site is decided by integer arithmetic, so no threshold is ever chosen. Measured coordinates do not arrive that way, and the moment a tolerance is introduced the answer stops being a fact about the structure and becomes a fact about the threshold.

The h0l layer of P2₁/c. The h0l reflections of P2₁/c out to 4 in each index, with each spot decided by summing the structure factor over the group's operations: 44 survive and 36 vanish identically, whatever the atoms are. The pattern of holes is the condition h0l: l even, read back off the spots rather than imposed on them.

The reflections that are not there

A screw axis and a glide plane leave no mark on the intensity of any reflection. What they do is delete some, exactly, for every possible arrangement of atoms — and the pattern of deletions is computed here from the sum a crystallographer writes down, rather than read from a table.

A structure, and the vectors between its atoms. On the left, 4 atoms in a cell. On the right, every one of the 16 vectors between them, each drawn from a common origin: 13 distinct positions, with the 4-fold peak at the origin being each atom paired with itself. That right-hand picture is what a Patterson map shows, and it is the thing a diffraction experiment gives without phases. It has more peaks than the structure has atoms — n² against n — which is why interpreting one is hard, and why it is always symmetric about its centre.

The map that needs no phases

A diffraction experiment measures intensities and loses phases, so the electron density cannot be computed from it. One map can be: the transform of the intensities, whose peaks are not atoms but the vectors between them — every ordered pair, brought to a common origin.

What each group extinguishes. The extinction conditions of 9 space groups, each derived by summing the structure factor over that group's own operations and reading the surviving rule off the result: P1 — nothing; P1̅ — nothing; P2 — nothing; Pm — nothing; P2/m — nothing; P222 — nothing; Pmm2 — nothing; P4 — nothing; P23 — nothing. 9 of the 9 extinguish nothing at all, and diffraction alone cannot distinguish those from each other.

Where the experiment runs out

Absences narrow the space group down and often not to one. Two groups can extinguish exactly the same reflections and scatter with exactly the same symmetry, and telling them apart needs something the diffraction pattern does not contain.

P2₁2₁2₁: the sections its symmetry forces. The Patterson cell of P2₁2₁2₁ with the sections marked. Each operation (M, t) sends an atom at x to Mx + t, so the vector between them is (I − M)x − t; where I − M is singular that vector cannot leave a plane, and the plane's equation comes from the left null space in integers. This group has 3 such operations, giving the sections w = 0.5, u = 0.5, v = 0.5. A heavy atom's vector to its own image is somewhere on one of them, which is what made structure solution possible before computers: a plane can be searched by eye and a volume cannot.

Where symmetry stacks the vectors

A Patterson map of a real structure is a blur with thousands of overlapping peaks. A screw axis rescues it: the vectors between symmetry-related atoms cannot leave a plane, so the search for a heavy atom is a search of a section rather than of a volume — and which plane it is falls out of the operation's matrix in integers.

Friedel's law, as an equality rather than a resemblance. Each pair of bars is a reflection and its opposite for a structure of three atoms in no particular arrangement. They are the same height, and not approximately: with real scattering factors, negating the indices conjugates the structure factor, and conjugation does not change a modulus. The test behind this figure requires the largest difference over 40 pairs to be below 10⁻⁹ and it comes back exactly zero. This is why a diffraction pattern is centrosymmetric whatever the crystal is, and why the thirty-two classes collapse to eleven before a structure is even proposed.

The law that hides handedness

With real scattering factors, negating the indices conjugates the structure factor and leaves the intensity exactly alone — so every diffraction pattern is centrosymmetric whatever the crystal is. The escape is an imaginary component that the negation does not touch, and it is how the handedness of a molecule is measured.

Every vector between every pair of atoms. The Patterson map of a four-atom structure: the transform of the intensities with every phase set to zero, so it is computable from a measurement and nothing else. Its peaks are not atoms but the vectors between them, and the ringed one is the strongest that is not the origin — taken here as an interatomic vector exactly as a crystallographer takes the vector between two heavy atoms. That it really is one of the structure's own vectors is checked rather than assumed: it lands within 0.0031 of a cell of a difference of two positions, and it sits 0.67 of a cell from the origin, well outside the skirt of the tall peak there. A candidate taken too close to the origin is the same peak seen again and the whole method fails quietly.

Solving from the vector set

A Patterson map contains a copy of the structure laid over every atom in turn. Shift it by one interatomic vector, take the pointwise minimum with itself, and the copies that fail to coincide are cut away — leaving the structure, together with its inverse, from a measurement that carries no phases at all.

N(z): the fraction of reflections weaker than z. The cumulative distribution of normalised intensities, measured on two structures built from the same atoms — one with an inversion centre, one without — and drawn against the two closed forms, 1 − e^(−z) without a centre and erf(√(z/2)) with one. The curves are furthest apart at small z, which is the useful end: a centrosymmetric structure has far more nearly-absent reflections, because its structure factor is a single real number that can pass through zero rather than a complex one that rarely does.

Whether there is a centre is a statistic

Everything else on this site is decidable: a pattern has a symmetry or it does not, and the detector settles it in integers. Whether a structure has an inversion centre is not like that. No single reflection carries the answer — the distribution of all of them does.

Two candidates, and the sign that chooses. The phase of reflection (2, 3), recovered from three measurements and no model. The circle is every complex number of the measured amplitude; the isomorphous difference fixes the cosine of the angle between the unknown phase and the heavy atom's, leaving the two candidates marked; the anomalous difference fixes the sine, which picks one. The recovered phase agrees with the true one to fifteen decimal places, and the true phase was never used in the calculation.

One experiment gives the cosine, the other gives the sine

Friedel's law holding exactly is what makes the phase unreachable. Its breaking is what hands it back: an isomorphous difference fixes the cosine of the phase and leaves two candidates, and the anomalous difference fixes the sine, which chooses.

Everything measurable at λ = 1.54 Å, and it is a finite set. A section through the reciprocal lattice of a cubic cell of 10 Å, with the limiting sphere drawn. Bragg's law reaches a reflection only if its spacing is at least half the wavelength, so the measurable reflections are the lattice points inside a sphere of radius 2/λ and the points outside it are not merely unmeasured — no experiment at this wavelength can reach them. In this section 516 points lie inside; in the full sphere there are 9,092.

How many reflections there are

The reciprocal lattice is infinite and a measurement is not. A wavelength cuts a sphere out of it, the number of points inside is the sphere's volume times the cell's, and symmetry then removes all but a fraction — though never quite the fraction a division would give.

Three atoms, four resolutions. A one-dimensional Fourier synthesis of the same three atoms, cut at four different resolutions. Nothing is approximate except the edge: every amplitude and every phase used is exact, and the only information withheld is the reflections outside the sphere. The peaks broaden as the cut-off comes in, and beside every peak sits a negative ripple that the coarsest map cannot distinguish from a real absence of density. Both effects are the transform of the sphere rather than anything about the structure.

As sharp as the sphere is wide

A map made from a truncated sum is not a blurred picture of the structure. It is the structure convolved with the transform of the sphere — so peaks acquire a width proportional to the resolution, and a negative ripple of twenty-two per cent that no improvement in the data ever reduces.

Two structures on 8 sites with the same vectors. Two arrangements of 4 atoms on a ring of 8 positions. They are not the same arrangement — no rotation of the ring and no reflection carries one onto the other — and every interatomic vector occurs the same number of times in both. The bars below are the shared vector counts, which is the Patterson function of each: the tall one at the origin is the atom count and carries no information, and everything else is what a diffraction experiment measures. Their diffraction patterns are identical in every intensity, so no measurement of intensities, at any resolution, distinguishes them.

Two structures, one Patterson

Eight arrangements of four atoms on a ring of eight sites, and only seven distinct sets of interatomic vectors between them. Two of the arrangements are genuinely different and no measurement of intensities can tell them apart — at any resolution, for ever.

Where a homometric pair comes from. A set that factors as a sumset gives its own partner. If every point of A is a sum b + c with b in B and c in C, and every sum arises once, then reversing C produces a different set with the same vectors — because reversing a factor and reversing its conjugate cancel in the product that the vector set is. Both factors must be asymmetric, which is the constraint that decides where the construction can be used: a two-point set is its own reflection up to a translation, so the smallest useful factorisation is three points by three points, and the smallest structure it builds has nine atoms.

Where the pairs come from

A structure whose atoms are the sums of two smaller sets has a partner: reverse one factor and the interatomic vectors do not notice. The construction is Patterson's own, it explains why homometry exists, and the smallest structure it can build has nine atoms for a reason worth following.

pg scatters as pmm. The orbit of a motif under pg, and the set of all ordered differences between its points brought to a common origin. The vector set was handed to the detector with no indication of where it came from, and came back as pmm: 31 peaks from 6 atoms, on the same lattice.

Seventeen groups, seven vector sets

A map of interatomic vectors is more symmetric than the structure it came from, twice over: it always acquires a centre, and it loses every translation part. So a glide becomes a mirror, seventeen plane groups collapse onto seven — and the collapse is verified by handing the vectors to a detector that has never heard of Patterson.

Observations per unknown, against resolution. Unique reflections divided by refinable parameters, for a triclinic cell, no angle a right angle, with three coordinates and six displacement parameters for every atom and one non-hydrogen atom per 18 ų, in a molecular crystal. The scale is logarithmic because the fall is a cube: 16.3 at 0.8 Å and 0.30 at 3 Å. The line at one is where a determination stops being over-determined, and it is crossed at about 2.0 Å.

The unknowns against the observations

A structure determination is a fit of some number of parameters to some number of measurements, and both counts can be worked out before any data exist. The ratio turns out not to depend on how large the crystal's cell is, or on how symmetric it is — only on the resolution, and on that as a cube.

-2 0 -3 is forbidden and is reached 64 ways. A layer of the reciprocal lattice of Cc: pale spots are reflections the group extinguishes, solid ones are allowed. The path shows a detour — the beam diffracts once at -3 -3 -3 and again at 1 3 0, and the two together send it in exactly the direction a single reflection at -2 0 -3 would. That reflection is forbidden, so intensity arrives where the symmetry said none could. There are 64 such routes to this one spot inside this window alone.

The absence that fills itself in

A systematic absence is the strongest evidence this subject has: a whole zone of reflections cancelling exactly, for reasons of symmetry rather than of arithmetic accident. The exactness belongs to a model — that the beam scatters once. A beam that has already been diffracted can be diffracted again, and the two events together land where the group said nothing could.

pmg: 80 of 80 restricted. The reflections of pmg inside a window of ±4, with the ones whose phase symmetry restricts to two values picked out. Every solid spot has a structure factor that must be real up to a fixed rotation — a sign, in effect — whatever the atoms turn out to be, and the pale ones have a phase symmetry says nothing about. The 4 palest spots carry two incompatible restrictions at once, which leaves them nothing to be but zero. Which spots these are was computed from the operations, before any structure existed.

The zones that behave as if there were a centre

The phase problem is usually stated as though symmetry had nothing to say about phases. It is true of most reflections and false of some, and which is decidable from the group alone: where an operation carries a reflection onto its own negative, the phase is confined to two values half a turn apart, computed from that operation's translation.

⟨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.

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.

B = 3.45 against 3.4, K = 0.37 against 0.37. The mean intensity of each resolution shell of a cell of 1194 reflections, divided by Σf² computed from the cell's content alone, and logged. The points fall on a line whose slope gives B = 3.45 against the 3.4 put in, and whose intercept gives a scale of 0.37 against 0.37 — both recovered before a single atom has been placed. The atoms here are independent, so the line is straight at every resolution; the shells below the cutoff are marked in the second colour.

The average that knows the atoms and not where they are

Square a structure factor and average it over a shell of reflections at one resolution. The cross terms — every one of which carries a fact about the arrangement — cancel, and what is left is a sum over the *content* of the cell with no position in it anywhere. A scale and a temperature factor come out of that before a single atom has been placed.

One crystal, two rows: one streaks and one does not. Two rows of reflections from the same faulted crystal, at a fault rate of 0.05, each drawn against the same row from a perfect one. The upper row has h − k not divisible by three, so each layer contributes a different cube root of unity and the sequence of layers enters the sum: the sharp peaks collapse into a streak. The lower row has h − k divisible by three, the phase factor is one, every layer scatters in step, and the peaks are exactly as sharp as in the perfect crystal. The sorting is an integer condition — a reflection either can see the stacking or cannot, decided by h − k modulo three — which is the most direct evidence there is that the disorder is in the stacking and not in the layers. The profiles are averaged over 16 independently faulted crystals, because one crystal gives speckle rather than a diffuse profile.

The streaks a faulted stack makes

Close packing settles two directions and leaves the third to chance. A crystal that chooses wrongly now and then has a lattice in the plane of its layers and none across them — and its diffraction pattern says so, with some rows of spots as sharp as ever and others smeared into streaks, sorted by an integer condition.

What each Laue class buys, in measurements per reflection. The eleven Laue classes, with how many distinct reflections a block of indices holds under each and how many times a data set measures the average one. The redundancy is always below the order of the class and the gap is the special reflections. This is the number an experiment is planned around: repeated measurements of what symmetry says must agree are the only estimate of precision that does not come from a model, so a triclinic crystal has to be turned through far more of the sphere than a cubic one to be measured as well.

Every reflection, several times over

A diffraction experiment does not measure each reflection once. Symmetry relates a reflection to the others of its orbit, and those are the same reflection seen from another direction — so a hundred thousand measurements may contain twelve thousand reflections, each observed eight times.

The twin fraction, recovered from a moment and nothing else (12 atoms). A structure of 12 atoms twinned at each of 6 fractions, with the second moment of its intensity distribution measured and the fraction solved back out of it. The recovery is within a few hundredths as far as thirty per cent — 4 rows here — and 4 of the 6 fractions get a number at all. Beyond thirty per cent the relation flattens: the derivative of 2α(1−α) vanishes at a half, the two roots meet, and a small error in the moment becomes a large one in the fraction. Where the sampled moment falls below 1.5 the quadratic has no real root and the estimate refuses rather than clamping, which is why a nearly perfect twin is the hard case in practice rather than the easy one.

A twin hides in the statistics

A twinned crystal scatters as two orientations at once and the detector cannot separate them. What arrives is a sum of two intensities — and adding two independent quantities narrows a distribution, which is a signature no model of the structure is needed to read.

The diffuse intensity of an alloy with α₁ = -0.46. The diffuse part of the scattering across the wavevectors an 12 × 12 block can be asked about, one square per wavevector with darkness the intensity. The marked squares are where the average structure scatters — the sharp part, which is what a Bragg reflection is. The diffuse maximum here is at (0.50, 0.50), which is the zone boundary: the alloy is trying to alternate, and a crystal that succeeded would put a sharp reflection exactly there. Summed over every wavevector, the intensity is exactly one per site whatever the correlations are — order moves scattering about, it does not create it.

The average scatters sharply and the rest does not

A crystal whose lattice is perfect and whose occupation is not scatters in two parts: the average structure gives Bragg reflections, and the variance is spread over everything between them. The split is exact, the total is one unit per site whatever the disorder does, and an ensemble of n arrangements mislays exactly a fraction 1/n of it.

How much one site knows about another, by separation. The Warren–Cowley parameters: the average of the product of the occupations of two sites a given vector apart, over every pair in every configuration. The value at the origin is exactly one — a site always agrees with itself — and it falls away with distance, alternating in sign where the alloy prefers unlike neighbours. These numbers are the whole of what the diffuse scattering measures: its intensity at a wavevector is their Fourier transform, computed here separately and agreeing to the last bits of the arithmetic. Nothing about them requires the crystal to be ordered, and their falling away is what short-range order means.

The order a diffuse pattern measures

Where a diffuse maximum sits says what the crystal is trying to become, and its shape is the Fourier transform of how much each site knows about its neighbours. The correlations are a small array of numbers, the intensity is their transform, and neither route to the other loses anything.

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.

One number, and it is the fraction. Data simulated from crystals that are nought, a quarter, a half, three quarters and wholly inverted, each fitted for the single parameter. The fitted values sit on the diagonal to better than five parts in a hundred, which is what makes the parameter a measurement of composition rather than a test of a hypothesis: a crystal is allowed to be part one hand and part the other, and a value near a half is a real answer about the specimen rather than a failure of the determination.

How much of it is the other hand

A crystal of one enantiomer is a hypothesis, not an observation. What the diffraction actually measures is a fraction — how much of the specimen is the inverted structure — and the useful part of that measurement is the uncertainty on it.

Two populations, and neither of them empty. The reflections of a structure in which three quarters of the atoms are paired by a half-cell shift, sorted by the parity of h + k and each class scaled by its own mean. The two histograms have the same shape, which is the point: each class on its own is an ordinary acentric distribution. What differs is the scale — the odd class is a sixth of the even one on average — and no odd reflection is absent, so no extinction rule fires and nothing about the space group is affected.

A translation that is nearly there

Half a structure copied onto the other half by a half-cell shift, with nothing exact about it. No reflection vanishes, so no extinction rule fires — and the test for a centre of symmetry answers yes about a structure that has none.

66 reflections that no rotation reaches. Every reflection inside the limiting sphere of an orthorhombic cell, 7 × 11 × 13 Å at 1.4 ångström, plotted by its distance from the rotation axis against its height along it. The ones marked are those a rotation about that axis can never bring into diffracting position: turning the crystal moves a point on a circle at fixed height, so a point too close to the axis can never acquire the component along the beam that the Ewald condition demands. The blind region is a cusp about the axis, it is 4.3 per cent of the sphere here, and nothing but remounting the crystal removes it.

What one turn of the crystal reaches

Every reflection inside the limiting sphere is measurable by some orientation. A crystal on a spindle has one axis, and a region around it never reaches the Ewald sphere at all — however patiently the crystal is turned.

Two plane structures with one Patterson. Two arrangements of 4 atoms on a 4 by 4 torus. No translation and no half-turn carries one onto the other, so they are different structures; every interatomic vector occurs the same number of times in both, so no measurement of intensities distinguishes them. The plane case is not the chain case with an extra index: a mirror in the plane sends a vector set to its mirror image rather than to itself, which is the one place the analogy with a cycle breaks.

Two structures on a torus, and one Patterson

Homometry was settled here on a ring of positions, which is a crystal in one dimension. Moving the same exhaustive search to a torus asks whether the coincidence is commoner or rarer when the vectors have a plane to land in — and the honest answer is that dimension is not what decides it.

How many reflections the centre test needs. The error rate of two tests for a centre of symmetry against the number of reflections used, measured on 60 centrosymmetric and 60 non-centrosymmetric structures at each point. The moment test — the one in every textbook, comparing ⟨|E|² − 1⟩ to its two theoretical values — reaches one error in twenty at 160 reflections and one in a hundred at 320. A likelihood ratio, which uses each reflection's own value instead of one average, reaches the same at 40 and 80. The gap is the price of summarising a distribution by its mean, and it is about a factor of four.

How many reflections it takes to know there is a centre

The test for a centre of symmetry compares one average of the intensities against two theoretical values a quarter apart. Whether that is a measurement depends on how many reflections went into the average, and the only honest way to find out is to run the test on structures whose answer is already known and count the mistakes.

One crossing or the other, and never both. Two events on a square patch of lattice: a path of occupied sites crossing from left to right, and a path of vacant sites crossing from top to bottom. On the triangular lattice exactly one of them happens in every configuration tested — the claim is combinatorial rather than statistical, so one counterexample would end it. On the square lattice both can fail at once, and do, in more than a quarter of the configurations. That difference is the whole of what follows.

The threshold a symmetry pins down

Occupy sites at random and somewhere the occupied ones first join up across the crystal. For almost every lattice that occupancy is known only to a few digits. For the triangular lattice it is exactly a half, and the reason is that on a lattice whose faces are all triangles an occupied path and a vacant path cannot slip past each other — a statement about one configuration at a time, with no probability in it.

Two kinds of atom, and more ambiguity rather than less. Exhaustive searches on rings of four sizes. The third column counts homometric groups when every atom is identical; the fourth counts them when each atom may be one of two kinds. The fourth is larger at every size, and at nine and ten sites the third is nothing at all — there is no pair of arrangements of four identical atoms that a diffraction experiment cannot separate, and there are six and four once the atoms may differ. Distinguishing the atoms adds information to the structure and adds ambiguity to the measurement.

When the atoms are not all the same

Every homometric pair found so far is a pair of point sets, where an atom is a point and counts once. Give the atoms different scattering powers and the ambiguity does not go away — it grows. On a ring of nine there is no pair of four identical atoms that diffraction cannot separate, and there are six once two kinds of atom are allowed.

The same structure, mapped from intensities and from differences. Left, the ordinary Patterson map of the structure: 14762 interatomic vectors, a continuous field of overlapping peaks, and the two vectors between the anomalous scatterers — circled — nowhere among its strongest. Right, the map from squared Bijvoet differences over the same reflections: two peaks after the origin, and they are those two vectors. The difference map is 7381 times smaller a problem to read.

A map of the atoms that break the law

Feed a Patterson synthesis the differences between the two halves of each Friedel pair instead of the intensities, and the map that comes back holds the vectors between the anomalous scatterers and nothing else. Two atoms among a hundred and twenty-two: 14,762 vectors become two.

Two symmetric structures a diffraction pattern cannot separate. Two arrangements of 6 atoms on a 6 × 6 torus, each invariant under the plane group p6m, drawn beside the Patterson they share. No translation and no inversion carries one onto the other, so they are different structures; every one of the thirty-six interatomic vector counts is the same, so every diffracted intensity is the same and no measurement at any resolution separates them. Of the 4 structures with this symmetry and this many atoms, there are only 3 Pattersons — so imposing the most symmetric of the seventeen plane groups has not removed the ambiguity.

Symmetry does not rescue a Patterson

Every homometric pair found so far sits on a bare ring with no operations imposed, and a real crystal sits in a space group. Impose one and the ambiguity does not go away: 12 of the 13 groups searched still have pairs, and at six atoms the hexagonal groups are indistinguishable two to three times as often as the general position.

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.

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.

All essays