Crystallography & Symmetry
patterns derived from their groups
Theme

The theme: From the diffraction back — page 3

Nobody has seen a space group. They are inferred from where a crystal scatters and, just as informatively, from where it does not.
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. How it is known

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. How it is known

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. How it is known

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.

 342 unlabelled spots, cell volume 52. A bag of 342 reflection positions with no indices on them, collected out to a bound of 3 on each index. Their pairwise differences generate the reciprocal lattice; a basis of that is taken by integer elimination and then reduced, and the reduced basis is printed. Its determinant is 52, which is the volume of the cell the reflections were computed from — so the cell has been recovered from positions alone, with no intensity used anywhere. Symmetry at work

A cell from a bag of spots

A single-crystal experiment returns a list of directions with no labels on them. Recovering the cell is recovering the lattice those directions generate, and the whole of it is take differences, reduce, read the answer. What no quantity of data settles is whether the lattice found is the true one or a sublattice of it.

 A peak that grows, and not fast enough. The strongest peak of three chains, divided by the square of the number of letters, as each chain is lengthened. A Bragg reflection is a sum of terms in phase, so its intensity grows as the square of the count and this number settles: the Fibonacci chain and the period-doubling chain both do, at exponents of about two. The Thue–Morse chain does neither — its strongest peak grows, so it is not diffuse scattering, and it grows more slowly than the square, so it is not a Bragg peak. The fitted exponents are printed beside each curve and no threshold enters the comparison. Order without repetition

Neither a peak nor a bump

A chain whose strongest reflection grows as the length to the power one and a half. A Bragg peak grows as the square and a diffuse bump grows as the length itself, so this is neither — and the essay that ruled out the first possibility could only say so by quoting a theorem.

 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. How it is known

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.

 p4m: freezing Γ3 leaves pmg. The same crystal three times. On the left, a pattern with the full symmetry of p4m. In the middle, the displacement each atom takes under the order parameter — the arrows are the mode, drawn in the first colour for one set of atoms and the second for the other where there are two. On the right, the atoms moved by a small multiple of those displacements. The group of the right-hand pattern is pmg, of index 4 in the parent, and it was found by the detector from the point set alone. The prediction — which operations carry the displacement field to itself — is made separately, from the mode and not from the points, and the two lists of operations are identical. Into space

The cell a zone-boundary mode doubles

An order parameter that alternates from cell to cell keeps only half the translations, so the frozen structure has a cell twice as large and reflections that were never there before. The phases at such a wavevector are ±1, so the whole computation stays in exact integers.

 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. How it is known

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. How it is known

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. How it is known

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 it is known

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. How it is known

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. How it is known

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. How it is known

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 it is known

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.

 The figure of merit a supercell always beats. One line list, indexed on the true cell and on five multiples of it. The mean discrepancy is not merely similar down the column, it is identical to every digit: the supercell's grid of allowed Q values contains the true cell's grid exactly, so each line lands on precisely the same place and misses by precisely the same amount. Any figure of merit built on agreement alone therefore returns one number for the whole family, and cannot prefer the true cell. What falls is the last column, and the only thing in it that the fit does not already contain is the count of lines the cell says should have been seen. Symmetry at work

The figure of merit a supercell always beats

Indexing a powder pattern returns a ranked list rather than an answer, and the ranking needs a number. The obvious number — how well the cell accounts for the lines — is exactly the number a supercell cannot lose on, because the supercell's grid contains the true cell's grid and the discrepancies are identical to every digit. What has to be paid for is the lines nobody saw.

 The alias accounts for every line and predicts more. The observed lines above, and below them the grid of a supercell that explains all of them. The full ticks are the observed lines, which the alias reproduces exactly; the faint ones are lines the alias predicts and nobody saw. That second set is the only thing that separates the two cells, and it is why an indexing criterion has to charge for unobserved lines rather than measure agreement. Symmetry at work

Every alias is a supercell

A cell that explains every line of a powder pattern is not a near miss and not a coincidence: its reciprocal grid contains the true one, which means its own cell is a superlattice of the true cell. So the ambiguity of indexing is the arithmetic of superlattices, and it can be counted — two cells with one unknown, sixteen with two, sixty-two with three, all of them accounting for the same twenty lines exactly.

 The cross, from the selection rule alone. The layer lines of a helix with the first maximum of each marked on both sides. Nothing here is a picture of a photograph: each mark is at the radius where the Bessel function of the lowest order the selection rule permits on that layer line first peaks, and that radius is proportional to the order. The order rises by one per layer line until the middle of the repeat, so the maxima lie on two straight lines through the origin — the X — and the larger marks are the layer lines that reach the axis. The classification

What a thread scatters

A helix with ten subunits in a turn is not a screw axis a crystal may have, and nothing about its diffraction pattern is lawless. The pattern lies on layer lines, and on each one only certain angular orders may contribute — a selection rule as hard as any extinction condition. The lowest permitted order rises by one per layer line, a Bessel function of order n does nothing until its argument is about n, and the maxima therefore lie on two straight lines through the origin.

 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. How it is known

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.

 A patch of hats scatters a pattern that repeats. The diffracted intensity of the 1217 points of a patch of 183 tiles laid out as the hat, over 2 by 2 cells of the kite grid's reciprocal lattice, whose edges are the faint lines. Every local maximum above a hundredth of the central peak is a disc with area proportional to its intensity; 12 reach the central peak's full height. 72 maxima are drawn. Adding a reciprocal lattice vector to the scattering vector changes the intensity by at most 1.1e-15 of the central peak, so each cell holds the same pattern. Order without repetition

How much of the hat is a crystal

Put a scatterer on every corner of a patch of hats and the diffraction pattern repeats exactly, because every corner sits on a lattice. Inside each repeat the strongest reflections are those of an ordinary crystal with partly filled sites, and by Parseval's identity they carry sixty-three per cent of what the pattern holds. The aperiodicity the hat is famous for lives in the remaining third, in reflections a hundred times weaker.

 Square ice scatters a pinch at the origin. The intensity scattered by the horizontal arrows of square ice, averaged over 100 configurations on a 32 by 32 torus, over the whole Brillouin zone with the origin at the centre; darker is more intense. Along the horizontal axis through the origin the intensity falls to zero — 2.9e-32 at the smallest wavevector — while along the vertical axis it stays near 1.75, so the two meet at the origin in a pinch. Order without repetition

The ice rule is a conservation law

Two arrows in and two out at every vertex is a statement that nothing flows in or out anywhere. That makes one half of the arrow field vanish identically, in every arrangement and not merely on average, and what is left scatters with a pinch at the origin: an intensity that approaches different values from different directions. Break the rule now and then and the pinch acquires a width Debye and Hückel predicted for a salt solution.

 The average is the site's orbit, with occupancies. A molecule at a site of symmetry mmm keeping a subgroup of order two takes four orientations, and the average over them is the site group's orbit of each of the molecule's atoms, every image at one over the length of its own orbit. Atoms in general positions give eight images at an eighth each, and give the same eight whichever subgroup the molecule keeps. Atoms on a locus the model keeps give a shorter orbit at a higher occupancy, drawn larger and darker, and those are the only atoms that differ between models. The total scattering is the same for every model, so all of them agree exactly at zero scattering angle. What a lattice forbids

The molecule size that hides a disorder

Two disorder models with the same occupancy leave averaged structures that differ only in a handful of partial atoms. The difference is 20% in structure factors for a ten-atom molecule and 3% for a sixty-atom one — so the data choose between the models for a small molecule and stop choosing for a large one, and seven pairs are identical at any size.

 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. How it is known

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. How it is known

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.

All themes · All essays