xraytools.

Fourier Synthesis and the Phase Problem

A diffraction pattern gives |F|, and the electron density is the Fourier transform of the structure factors — so with amplitudes and phases the structure follows by summation. The phases are not measured. This sums a real structure both ways, so you can see what the missing half was carrying.

You supply
One of the named structures, and how far the series should run. Everything else — the coefficients, their signs, the electron count — is computed from the atoms.
Reading it
Summing over h alone gives the density projected down b and c, not a section — atoms sharing an x land on top of one another. And a truncated series rings: it merges neighbours the resolution cannot separate, and can dip below zero between the atoms.

Worked examples: zirconia — 7 of 12 signs negative · aragonite, alternating signs · quartz at 6 terms — the silicons merge · quartz at 30 terms — and separate · caesium chloride — where the demonstration is empty

Input

up to h =

How far the series runs. Blank means 12. Fewer terms is lower resolution — the peaks broaden and neighbouring atoms merge.

Gaps in the row of bars are systematic absences: a centred lattice or a glide plane makes whole classes of h00 vanish, so the series has fewer terms than its length suggests and the projection repeats more often than the cell does.

The projected electron density

0.000.00 Å0.251.23 Å0.502.46 Å0.753.69 Å1.004.91 Å075151226Si at x = 0.0000; nearest maximum 0.0000 awayO at x = 0.1466; nearest maximum 0.0590 awayO at x = 0.2669; nearest maximum 0.0613 awayO at x = 0.4135; nearest maximum 0.0865 awaySi at x = 0.4697; nearest maximum 0.0303 awaySi at x = 0.5303; nearest maximum 0.0303 awayO at x = 0.5865; nearest maximum 0.0865 awayO at x = 0.7331; nearest maximum 0.0613 awayO at x = 0.8534; nearest maximum 0.0590 awayatoms, projectedx, in fractions of a and in Åelectrons per unit x

&alpha;-quartz, SiO<sub>2</sub>, P3221 — 6 non-zero coefficients out of 6, carrying detail to 0.71 Å. The second curve is the same sum with the phases discarded.

Summing over h alone gives the density projected down b and c onto the a axis — not a section through the cell. So a peak sits at an atom’s x whatever its y and z are, and two atoms sharing an x project on top of one another: in rock salt the sodium and the chlorine both land at 0 and at ½, and the projection cannot tell them apart. A true section along x would need every hkl.

The phase problem

2 of the 6 coefficients are negative — a phase of 180°. A diffractometer measures the height of each bar and nothing else, so those signs are exactly what a measurement does not give you. Take every coefficient positive and the sum collapses: at x = 0 every cosine is 1, so every term reaches its maximum together and the map has one large peak at the origin. That peak is a property of having thrown the phases away, not of this crystal — it appears at the origin for any structure treated this way.

F(100) = -16.2, d = 4.255 A, phase 180 degF(200) = 18.1, d = 2.128 A, phase 0 degF(300) = -9.1, d = 1.418 A, phase 180 degF(400) = 13.4, d = 1.064 A, phase 0 degF(500) = 0.3, d = 0.851 A, phase 0 degF(600) = 4.5, d = 0.709 A, phase 0 deg123456h, in F(h00)F(h00)

Each bar is one coefficient, drawn with its sign. A measurement gives the height of every bar and not which side of the line it is on.

What was summed

Cell edge a 4.9134 Å
Atoms in the cell9
Distinct projected sites 9
F(000) 89.99 electrons
Mean of the curve 89.99 — which is F(000) again, and it never entered the series
Coefficients used 6 of 6
Negative coefficients 2
Highest value, phases discarded 213.1 at x = 0
Correlation of the two curves 0.2078

The series stops, and a stopped series rings. Between the atoms the sum overshoots and undershoots, leaving small maxima where the crystal has nothing. Here it stays above zero, which is not a general property: a few light atoms on a large F(000) baseline do not ripple that far, while rock salt at the same term count goes visibly negative. Atoms closer together in x than the series can resolve merge into one maximum: quartz’s two silicons, 0.056 apart, come out as a single peak at ½ until the series is long enough to separate them.

Where the atoms are

xAtoms there Nearest maximum
0.0000 Si 0.0000 away
0.1466 O 0.0590 away
0.2669 O 0.0613 away
0.4135 O 0.0865 away
0.4697 Si 0.0303 away
0.5303 Si 0.0303 away
0.5865 O 0.0865 away
0.7331 O 0.0613 away
0.8534 O 0.0590 away

The distance is to the nearest maximum in the curve above, and it is an annotation rather than a test: not every atom has a peak of its own once the projection has merged it with a neighbour, and not every maximum is an atom.

Try it

zirconia — 7 of 12 signs negative · aragonite, alternating signs · quartz at 6 terms — the silicons merge · quartz at 30 terms — and separate · caesium chloride — where the demonstration is empty