Crystallography & Symmetry
patterns derived from their groups
Theme

The theme: The same arithmetic, renamed — page 8

A crystal form is an orbit. A twin law is a coset. The domain states left by a phase transition are the cosets of the low-symmetry group in the high-symmetry one. Four subjects that grew up in different centuries and different departments, doing one piece of arithmetic under four names.
The centrings six four-dimensional lattice groups allow, prime by prime. For each of the six four-dimensional lattices whose symmetry groups are maximal, with the group's order, and each prime from 2 to 7: the subspaces of the lattice modulo that prime the group keeps, sorted by what centring by each one makes. New Bravais types appear only at 2, three of them, all over the cube-and-line lattice, and at 5, where A₄ centred by fifths is A₄*. Both centrings of the hypercubic lattice by halves make D₄, which has three times the symmetry. Lattices

The body-centred hypercube is the checkerboard lattice

In three dimensions the centring letters come out as the subspaces a holohedry keeps modulo a prime, and only two and three make new lattices. In four dimensions the same census, run on the six lattices whose symmetry groups are maximal, finds a prime that the three-dimensional argument never could: five centres A₄ into its dual. It also finds that the most familiar centring does not survive the extra dimension. The body-centred hypercubic lattice is not a centred hypercube at all: its body centre is exactly as far away as an edge, and it is D₄, with three times the symmetry.

 Ninety-six four-colour designs. For each of the seventeen plane groups, the number of its four-colourings, as subgroups of index four, and the number of distinct designs they fall into under the affine normaliser, sorted by the group that permutes the four colours. The 281 colourings make 96 designs: 13 ℤ₄, 39 Klein, 39 D₄, 2 A₄, 3 S₄. The classification

Ninety-six four-colour designs

The seventeen plane groups have 281 four-colourings, counted as subgroups of index four, and many of them are the same design described from a different origin, on different axes or from a neighbouring cell. Identified under every affine map that carries the pattern's group to itself, they make ninety-six designs. The same computation run at two and three colours gives forty-six and twenty-three, and all three numbers match the classical counts of coloured plane groups.

 Which way a fibre of achiral crystals rotates light. For three achiral classes that permit a gyration tensor, every crystal direction that could be aligned along a fibre texture's axis, drawn on a disc as seen from above the z axis, and shaded by the sign of the rotation the texture then shows along its axis. In each class the rotation takes both signs, in regions separated by a cone on which it vanishes. For 4̅2m and mm2 that cone is exactly the class's mirror planes, where the texture inherits a mirror and must be inactive; for 4̅, which has no mirrors, the cone is where the rotation passes through zero between a right-handed texture and a left-handed one. What symmetry decides

A fibre of achiral crystals can rotate light

Four achiral crystal classes permit optical rotation, of both signs in different directions. Grind such a crystal into a random powder and the rotation averages away completely. Align the grains along one direction and it comes back: the texture rotates light along its axis exactly as one grain does, half as much the other way across it — and the direction chosen for the alignment decides the texture's hand.

 Which cells the flipping solver finds. For four cells of twenty atoms — all carbon, a quarter sulphur, one sulphur, one lead — the share of random starting phase sets that charge flipping takes to the right structure within 250 cycles, against the threshold below which density has its sign reversed. Each point pools several random arrangements of atoms with several starts each, and the atoms fall off with angle as real atoms do. The cells with sulphur are solved more often than carbon alone, the quarter-sulphur cell most of all, each at a threshold near one standard deviation; the lead cell is never solved at any threshold. How it is known

The heavy atom the flipping solver cannot see past

Sayre's identity degrades as a cell's atoms become unequal, most of all for a quarter of heavy atoms three times the weight of the rest, and the obvious prediction is that a solver built on the same sparseness should struggle there too. It does not. Charge flipping solves a cell of carbon with sulphur in it more often than carbon alone, and never solves a cell with one lead atom — because a threshold stated in standard deviations belongs to whichever atom owns the map's variance.

 Two hundred and twelve of the 219 are made of their own operations. The 219 affine space-group types divided by which of their own operations generate them, lattice translations included. 110 are generated both by their screws and glides alone and by their rotations, mirrors and inversions alone; 67 by their screws and glides only; 29 by their point operations only; 6 only by the two kinds together; and 7 by nothing they contain — P1, P2, Pm, Pmm2, P4, P3, P6 — which need a translation given outright. Into space

Seven space groups cannot make their own translations

Every space group contains a lattice of translations, and for all but seven of the 219 types those translations are products of the group's other operations. A screw done twice climbs its axis and a glide done twice slides along its plane, so the non-symmorphic operations make lattices by themselves. The seven that cannot are exactly the seven plane groups with no glide, repeated along a line they all fix, and the reason is a difference between a glide line and a glide plane.

 Diffuse scattering gathered under every Bragg peak. The one-phonon thermal diffuse scattering of a square lattice with springs to nearest and next-nearest neighbours, over the reciprocal plane from minus 2.5 to 2.5 in each index, with the Bragg reflections marked. The intensity rises as one over the squared distance to every reflection, because the acoustic phonons' frequencies vanish there, and it grows with the square of the scattering vector, so the lobes are strongest far from the origin. The lobes are not round: a small step from (1,0) along k gives 2.7 times the intensity of the same step along h, because it probes transverse phonons, which are softer, and that ratio is what an experiment reads the elastic constants from. How it is known

A warm plane crystal has no Bragg peaks

Thermal motion is disorder that moves, and in space it costs a crystal a little of each Bragg peak and gathers it into diffuse scattering underneath. In a plane it costs more than a little. Computed exactly from the phonons of a harmonic square lattice, an atom's mean-square wandering grows by the same amount every time the crystal doubles in size, so no Debye–Waller factor survives, and every Bragg peak fades as a power of the crystal's size — a power proportional to the temperature times the square of the reflection's distance from the origin.

 A looser fibre with more longitudinal response. The three piezoelectric coefficients of a polar fibre of grains with d33 = 90, d31 = −33 and d15 = 560 pC/N, as the tilt spread grows from nothing to a uniform one, with the degree-one and degree-three parts of d33 dashed. The degree-one part is 265 and the degree-three part -175; the negative octupole decays faster, so d33 rises from 90 to 168 at a spread of 30° before falling to nothing. What symmetry decides

A looser fibre, more piezoelectric

Spread the grains of a fibre about its axis and every property decays, but not at one rate: each spherical harmonic a property is built from is multiplied by its own average, and the higher the degree the faster it goes. A property made of two degrees therefore changes shape as the fibre loosens. For a polar grain whose shear coefficient dwarfs its longitudinal one, the fibre's longitudinal response nearly doubles before it falls; for another, a shear coefficient changes sign.

 Thirty-one band groups from three signs. Every operation of a ribbon keeps or reverses the direction along it, the direction across it and the sheet's two faces, so its point group is a subgroup of the eight sign triples. Each operation that keeps the along-direction can also carry half a period along the ribbon, a glide or a screw, and the choice must be consistent. The sixteen subgroups, weighted by their consistent choices, give 1 + 10 + 16 + 4 = 31, the number of band groups, and the enumeration finds all 31 as distinct groups. Into space

Every band group is a ribbon

A ribbon cut from a sheet has a symmetry group with translations along one direction only, one of 31 band groups. The 31 can be counted from three signs: whether an operation reverses the ribbon's direction, its width, or the sheet's faces, and whether it slides half a period. Cutting the eighty layers every way reaches all 31, seven of them as single edges and twenty-four only as ribbons cut symmetrically about a special line. Five of those twenty-four have no improper operation at all.

 The worst edge pair is the two most anomalous wavelengths. For each pair of the four wavelengths, and for the peak, inflection and high remote together, the median error of the phase difference between the whole structure and its anomalous atoms, recovered from both Friedel mates at each wavelength with a 1% error on every intensity. Beside each, the spread of f′ across the chosen wavelengths. The pairs that include a remote wavelength and one on the edge phase to seven or eight degrees; the peak and the inflection together, the two with the largest anomalous effect, to eighteen; the two remotes to fifty-five. How it is known

The two most anomalous wavelengths are the wrong pair

A MAD experiment measures both Friedel mates at two wavelengths near an absorption edge and solves for a phase. The obvious choice is the two wavelengths where the anomalous effects are largest, the peak of f″ and the dip of f′. Of all the pairs that use the edge, it is the worst. The sine of the phase is carried by f″, but the cosine is carried by how much f′ changes between the two wavelengths, and two points on the same edge barely change it.

 Every crystal class permits the Faraday effect. The thirty-two crystal classes, each with two counts: the number of independent components of the tensor that turns a magnetic field into an optical rotation, and the number of components of the natural gyration tensor. The first is never zero: 5 classes have 1, 11 classes have 2, 11 classes have 3, 3 classes have 5, 2 classes have 9. The second is zero in 17 classes, every one with an inversion centre or a mirror arrangement that cancels it. The Faraday count is the mean of the squared trace over the class, which is never below one because the isotropic part always survives. What symmetry decides

A field lets every class rotate light

Seventeen of the thirty-two classes forbid a crystal to rotate the plane of polarisation on its own. Put the same crystals in a magnetic field along the ray and not one of the thirty-two forbids it, the centrosymmetric ones included. The field supplies the one thing natural rotation needs a handed structure for, and time reversal is what tells the two rotations apart.

 Which plane groups can be coloured with five, seven, eleven or thirteen colours. For each of the seventeen plane groups and each of the primes 5, 7, 11 and 13, the number of its colourings with that many colours, as subgroups of that index, and the group that permutes the colours: cyclic (C), dihedral (D), or the affine group of the colours with a multiplier of order 3, 4 or 6 (written p:3, p:4, p:6). Totals 104, 172, 268 and 420. The square groups p4 appear at 5 and 13, the three-fold and six-fold groups p3 and p6 at 7 and 13, and p4m, p4g, p3m1, p31m and p6m never. The classification

Five colours, seven colours, and the primes that split

Four colours gave the seventeen plane groups 281 colourings, five different groups permuting the colours and ninety-six designs. Five colours give 104 colourings, three groups and fourteen designs, and from five upwards every prime behaves the same way. A prime number of colours that the point group cannot divide forces every colouring onto a sublattice the rotations keep, and whether the square and hexagonal patterns have one is decided by whether the prime splits in the Gaussian or the Eisenstein integers.

 Three difference maps of one substructure, and the map they approximate. Patterson syntheses over 1,681 reflections of a structure with 120 light atoms and two anomalous scatterers, from squared Bijvoet differences at one wavelength, squared dispersive differences between two wavelengths, and the two normalised and summed; the fourth panel is the Patterson of the anomalous atoms' own intensities. Circles mark the two vectors between the anomalous atoms. Each map's correlation with the exact one is printed beneath it: Bijvoet 0.775, dispersive 0.838, combined 0.999. How it is known

A second wavelength completes the map

The Bijvoet difference carries the anomalous substructure's intensity times the sine squared of a phase; the dispersive difference between two wavelengths carries it times the cosine squared. Normalised and added, the two make the substructure's own Patterson, correlation 0.999. It needs no fewer reflections than either half, because the threshold was never the average. What the dispersive half brings is a way of failing the Bijvoet half cannot have: two datasets, and a scale between them.

 Every fibre texture is a point in one plane. The two numbers that decide every coefficient of a polar fibre, the mean first and third Legendre polynomials of the grains' tilt, as a plane. A cone of grains at a single tilt is a point on the curve; any other texture is an average of such points and lies in the curve's convex hull. Spreads about the axis trace the inner path. For a grain with a large shear coefficient the largest d33 of any texture is at a point of the curve, the cone at 51°, 223 pC/N; the best spread about the axis reaches 168. What symmetry decides

The best texture is a cone

Every fibre texture of a polar grain, however its orientations are distributed about the axis, reduces for piezoelectric purposes to two numbers, the averages of the first and third Legendre polynomials of the tilt. The textures fill a region of a plane bounded by one curve, and a coefficient is a straight line across that plane, so its best texture sits on the curve: a cone of grains at one angle. For a grain with a large shear coefficient the best cone lies at 51° and gives 223 pC/N, where the best spread about the axis gives 168.

 Whole stars of the hexagonal zone, and the boxes that hold them. For ten wavevectors of the primitive hexagonal zone, the number of cells in the smallest periodic box holding one arm of the star, in the smallest holding every arm, and, where it differs, in the smallest n by n grid of cells holding the star. The corner K is a single arm up to sign and is held whole by three cells, against nine for a grid; K raised to the top face, H, by six against eighteen; halfway to K, T, by twelve against thirty-six. M needs two cells for one arm and four for its three, and L two and eight. Every smallest box for a whole star keeps the full hexagonal symmetry. The classification

The corner of a hexagon fits in three cells

A periodic box for a hexagonal crystal holds a whole star of wavevectors only if it keeps the crystal's full symmetry, as a cubic box must. What the cube never offered is a box turned through thirty degrees. It holds the corner of the hexagonal zone in three cells where every grid needs nine, and the saving is always exactly a factor of three and always at the same stars. Turn a box through any other angle that the six-fold axis allows and it loses the mirrors, and comes with a mirror image.

 Every crystal class, as a reflection class and its relative invariants. The thirty-two crystal classes, each with the smallest reflection class containing it after a relabelling of the axes, the index, and the degrees of the secondary invariants that make its invariant ring a free module over the reflection class's polynomial ring. Twelve classes are reflection classes themselves; seventeen have index two and one secondary invariant besides 1; three — 1̄, 4̄ and 3̄ — have index four and three. What symmetry decides

What the twenty lack is a product of mirrors

Twelve crystal classes have invariant rings generated by three independent polynomials, and they are the twelve generated by their own mirrors. The other twenty were left as the classes whose counts need a correction. The correction has a simple shape. Each of the twenty sits inside one of the twelve, at index two or four, and its invariants are the larger class's invariants times one or three extra polynomials. Each extra polynomial is a product of the mirror planes the smaller class lacks. For the rotations of the cube it is a single polynomial of degree nine, the product of all nine mirrors of the full cubic class.

 Whole stars of the rhombohedral zone, and the boxes that hold them. For eight wavevectors of the rhombohedral zone, the number of cells in the smallest periodic box holding one arm of the star, in the smallest holding every arm, and, where it costs more, in the better of the two regular grids: the rhombohedral cell repeated n times along each edge, or the hexagonal triple cell repeated n by n by m. The axis end T is held whole in two cells against six, the face centre F in four against eight, and the in-plane third K₃ in nine against twenty-seven; L needs eight either way. Every smallest box for a whole star keeps all twelve operations of 3̅m. The classification

A rhombohedral box saves a factor of two the hexagon never offered

A rhombohedral crystal is drawn on two cells, the small rhombohedron and the hexagonal cell three times its size, and a periodic box is usually one of them repeated. The smallest box that holds a whole star of wavevectors is the star's own lattice instead, and across every star with denominators up to eight it beats the better of the two grids by one, two or three, never more. Three was the hexagon's number. Two is new, and it is the face centre's. And the boxes that keep the three-fold axis and lose the mirrors still come in mirror pairs, now at nine cells as well.

 Which plane groups can be coloured with five, twenty-five, seven and forty-nine colours. For each of the seventeen plane groups, its colourings with 5, 25, 7 and 49 colours, as subgroups of that index. At p² colours the count splits into those whose translations meet the lattice in the whole plane pL and those that meet it in a cyclic sublattice. The five groups with a mirror across a rotation, shaded, have none at 5 or 7 and p² of them at 25 and 49, all on pL. In all, 104, 1449, 172 and 4197 colourings. The classification

Twenty-five colours, and the mirrors come back

With a prime number of colours from five up, five of the seventeen plane groups can never be coloured at all: the ones with a mirror across a rotation. With the square of a prime they all can. A colouring with twenty-five colours may meet the lattice in the whole plane modulo five, which every group keeps, and the colours then form a plane the point group acts on by its matrices. The designs come out as the prime's count again, with one more for each of the seventeen groups.

 Where the third wavelength should go. The median Cramér–Rao bound on the phase difference for two hundred reflections, measured at the inflection and the high remote and at a third wavelength whose energy runs across the edge, twice: with f′ known everywhere, and with f′ at the edge energies refined from the data. Known, the best third sits 1 eV above the edge at 6.26°, and a repeat of either wavelength is nearly as good. Refined, both repeats determine nothing, and the best third moves to +2 eV at 6.30°. How it is known

The third wavelength is there to measure the edge

A three-wavelength MAD experiment is usually justified as more anomalous signal: a third look at the edge, placed where the first two left the most out. Counted with the scattering factors known, a third wavelength buys nothing a longer look at the best pair would not buy better. Its real job appears only when f′ at the edge has to be found from the data, which it always does. Two wavelengths then carry one equation too few, and any third energy that is not a repeat supplies it, for two hundredths of a degree.

 The two forty-sixes over the ten plane point groups. The forty-six two-colour plane groups by the point group of the uncoloured group, and the forty-six rods that turn their axis over by the plane point group their operations show end on. The colourings pile onto the rectangular mirror groups — seventeen on 2mm, ten on m, none on 3 — and the rods onto the hexagonal and square groups, eleven on 6mm and nine on 4mm. Seventeen of the forty-six would have to change point group for the two lists to match. The classification

Two forty-sixes that do not match

The layer groups contain forty-six two-coloured plane patterns, and the rod groups contain forty-six threads that turn their axis over. Both lists are lists of a group with a half of index two, and both halves have a point group in the plane, so they can be laid side by side. They do not line up on any invariant tried: the plane spends its forty-six on rectangular mirrors, which halve a lattice easily, and the rods spend theirs on high rotations, which carry the most kinds of screw. The equal totals are two sums that happen to meet.

 The window in which A₄ and its dual are both minima. The exact second derivative of the Gaussian sum at A₄ and at A₄ along their five-dimensional family of changes of shape, for widths within six hundred-thousandths of one. A₄'s turns positive at t = 0.99998110 and A₄'s turns negative at 1.00001890, the reciprocal. Between them both are positive, so both lattices are local minima; at t = 1 itself both equal 8.59e-6. Lattices

The self-dual width is not flat

A search over four-dimensional lattices found that A₄ and its dual stop being traps for the Gaussian sum almost exactly at the self-dual width, and read that as a flat of the sum sitting on the width itself. Computed exactly rather than by finite differences, the curvature there is not zero. It is six parts in a million of the other curvature, positive, and the same for both lattices, because duality makes their Hessians one matrix at that width. Both are minima, in a window four parts in a hundred thousand wide, with basins a few millionths across.

 A scale error between two wavelengths has a sign. The contrast of the substructure's vectors in a dispersive difference Patterson, the weaker of its two peaks over the strongest other peak, as the second dataset is scaled wrongly, with no metal atom and with two metal atoms that do not scatter anomalously. An error that weakens the second dataset loses the map at 3.7% with no metal; one that strengthens it at 50%. The dashed line is where the substructure's peak is predicted to reach zero, a quarter of the fractional change in f′ between the wavelengths. The metals shorten the strong side and hardly move the weak one. How it is known

A scale error between two wavelengths has a sign

A dispersive difference compares two datasets, so a wrong scale factor between them adds the whole structure to it. On a structure whose own Patterson is a haze, the substructure's vectors survived a scale error of fifty per cent, and heavy atoms that do not scatter anomalously were expected to take over much sooner. They do, but only in one direction. Scale the second dataset the other way and the substructure is erased at four per cent whatever else is in the crystal, because the error's linear term cancels the signal. Scaling the datasets onto each other is itself an error of that sign.

 Every four-dimensional lattice with a five-fold rotation, on one interval. The lattices with a five-fold rotation in four dimensions are one ring, ℤ[ζ₅], with the form Tr(a·x·x̄), and the number t = log(a/σ(a)) / 2 log φ, folded into the interval from 0 to 1, says which lattice a is. The dots are the members with a = u + v√5 for u up to 9, 8 distinct lattices. The two ends are A₄* and A₄, with 240 symmetries each; every point between has 20. The arcs join each member to the one a centring by fifths makes of it, its reflection t ↦ 1 − t. Lattices

A five-fold lattice has one free number

In four dimensions a lattice can have a five-fold rotation, and A₄ and its dual both do. Every other lattice with one turns out to be the same ring, the integers of the fifth roots of unity, with a single number deciding its shape. Folded onto an interval, the number runs from A₄* at one end to A₄ at the other, and every lattice strictly between has twenty symmetries. Centring by fifths does not make a new lattice. It reflects the interval end to end, exchanging A₄ with A₄* and pairing every other lattice with a partner, and no lattice sits at the middle to be its own.

 Every six-dimensional lattice with a seven-fold rotation, on one torus. The lattices with a rotation of order seven in six dimensions are the ring ℤ[ζ₇] with three positive weights, one on each complex line it spans. Up to scale and isometry the weights live on a torus, drawn as a sixty-degree rhombus with opposite sides glued, shaded by how well each lattice packs spheres. The Galois automorphism turns the torus a third of a turn about three points, on the long diagonal: A₆* with 10080 symmetries, A₆ with 10080, and the ideal (1 − ζ)² with 672, the densest lattice in the family. The small dots are integral members, each with 28 symmetries. Lattices

A seven-fold lattice lives on a folded torus

In six dimensions a lattice can have a seven-fold rotation, and every one that does is the integers of the seventh roots of unity with three weights on it. Up to scale the weights live on a hexagonal torus that a Galois symmetry folds by a third-turn, and the fold has three fixed points: A₆*, A₆, and a third lattice that packs better than either. Centring by sevenths is no longer a reflection, as fifths was in four dimensions. It moves every lattice a third of the way along the torus and carries the three special ones round in a cycle.

 No triple of wavelengths is best for both the phase and the edge. Every triple of energies on a grid across a model selenium edge, placed by two Cramér–Rao bounds with the edge's position and white line refined from the data: the median phase error across, and the error in the edge's position up, on a logarithmic scale. The front of triples that no other beats on both runs from +1, +2, +150 eV (6.15°, ±0.066 eV) to -120, -1, +3 eV (7.54°, ±0.022 eV). The classic peak, inflection and high remote sits on it, at 6.42° and ±0.026 eV. How it is known

The wavelengths that phase best do not find the edge

A MAD experiment refines the anomalous atom's scattering factors so that it can solve for the phases, and those same factors report where the absorption edge sits and how tall its white line stands, which is chemistry. Scored for the edge instead of the phase, the best triple of wavelengths moves: two straddling the step and one far below it, where the phase-best triple put two on the step and one far above. No triple is best at both, and the textbook peak, inflection and remote turns out to be the compromise between them.

All themes · All essays