xraytools.

Fourier Synthesis and the Phase Problem

Single-crystal data

A diffraction pattern gives intensities, and correcting those gives |F|2, so amplitudes. The phases are never recorded at all, and the page opens by saying what that costs. This sums a real structure both ways, so you can see what the missing half was carrying.

Before this
A structure factor has an amplitude and a phase, and only the first survives a measurement. If that is not yet concrete, build one first.
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

Earlier on the path: Structure Factor Calculator Next on the path: The Patterson Function On From intensities to a structure, step 2 of 6

Notation here: F, |F|what each one means here

See also: Structure Factor Calculator · The Patterson Function · Difference Map

What each input changes
Terms
How many reflections the sum runs over — the resolution of the map. Too few and two atoms merge into one peak, which is a limit of the DATA and not of the method.
Teaching with this page
Objective
After this page a learner can state the phase problem precisely and show what a correct set of amplitudes with wrong phases produces.
Start from
this worked example
Ask first
A synthesis from six terms shows a maximum where the structure has no atom. Is the structure wrong?
Watch for
“Yes — density appears where the electrons are”
Then
The Patterson Function
Check yourself: A synthesis from six terms shows a maximum where the structure has no atom. Is the structure wrong?

No — a series cut short has ripples of its own Yes — density appears where the electrons are

The density is a sum over all reflections and any measurement supplies a finite number of them. Stopping at six is multiplying the true transform by a box, and a box in one space is a ripple in the other — so features appear beside real atoms and between them. Adding terms shrinks them and no number of terms removes them. Add some here and watch it happen.

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.252.03 Å0.504.07 Å0.756.10 Å1.008.14 Å04559111,366with the phasesevery phase set to zeroSi1 at x = 0.0037; nearest maximum 0.0037 awaySi1 at x = 0.0037; nearest maximum 0.0037 awayO1 at x = 0.0047; nearest maximum 0.0047 awayO1 at x = 0.0047; nearest maximum 0.0047 awayAl at x = 0.0088; nearest maximum 0.0088 awayAl at x = 0.0088; nearest maximum 0.0088 awayO5 at x = 0.0132; nearest maximum 0.0132 awayO5 at x = 0.0132; nearest maximum 0.0132 awayO6 at x = 0.0234; nearest maximum 0.0234 awayO6 at x = 0.0234; nearest maximum 0.0234 awayO2 at x = 0.0920; nearest maximum 0.0031 awayO2 at x = 0.0920; nearest maximum 0.0031 awayO4 at x = 0.1802; nearest maximum 0.0087 awayO4 at x = 0.1802; nearest maximum 0.0087 awaySi3 at x = 0.1810; nearest maximum 0.0079 awaySi3 at x = 0.1810; nearest maximum 0.0079 awayO8 at x = 0.1841; nearest maximum 0.0048 awayO8 at x = 0.1841; nearest maximum 0.0048 awayO3 at x = 0.1875; nearest maximum 0.0014 awayO3 at x = 0.1875; nearest maximum 0.0014 awaySi2 at x = 0.1919; nearest maximum 0.0030 awaySi2 at x = 0.1919; nearest maximum 0.0030 awayO7 at x = 0.2064; nearest maximum 0.0175 awayO7 at x = 0.2064; nearest maximum 0.0175 awayNa at x = 0.2317; nearest maximum 0.0428 awayNa at x = 0.2317; nearest maximum 0.0428 awayNa at x = 0.2683; nearest maximum 0.0428 awayNa at x = 0.2683; nearest maximum 0.0428 awayO7 at x = 0.2936; nearest maximum 0.0175 awayO7 at x = 0.2936; nearest maximum 0.0175 awaySi2 at x = 0.3081; nearest maximum 0.0030 awaySi2 at x = 0.3081; nearest maximum 0.0030 awayO3 at x = 0.3125; nearest maximum 0.0014 awayO3 at x = 0.3125; nearest maximum 0.0014 awayO8 at x = 0.3159; nearest maximum 0.0048 awayO8 at x = 0.3159; nearest maximum 0.0048 awaySi3 at x = 0.3190; nearest maximum 0.0079 awaySi3 at x = 0.3190; nearest maximum 0.0079 awayO4 at x = 0.3198; nearest maximum 0.0087 awayO4 at x = 0.3198; nearest maximum 0.0087 awayO2 at x = 0.4080; nearest maximum 0.0031 awayO2 at x = 0.4080; nearest maximum 0.0031 awayO6 at x = 0.4766; nearest maximum 0.0234 awayO6 at x = 0.4766; nearest maximum 0.0234 awayO5 at x = 0.4868; nearest maximum 0.0132 awayO5 at x = 0.4868; nearest maximum 0.0132 awayAl at x = 0.4912; nearest maximum 0.0088 awayAl at x = 0.4912; nearest maximum 0.0088 awayO1 at x = 0.4953; nearest maximum 0.0047 awayO1 at x = 0.4953; nearest maximum 0.0047 awaySi1 at x = 0.4963; nearest maximum 0.0037 awaySi1 at x = 0.4963; nearest maximum 0.0037 awaySi1 at x = 0.5037; nearest maximum 0.0037 awaySi1 at x = 0.5037; nearest maximum 0.0037 awayO1 at x = 0.5047; nearest maximum 0.0047 awayO1 at x = 0.5047; nearest maximum 0.0047 awayAl at x = 0.5088; nearest maximum 0.0088 awayAl at x = 0.5088; nearest maximum 0.0088 awayO5 at x = 0.5132; nearest maximum 0.0132 awayO5 at x = 0.5132; nearest maximum 0.0132 awayO6 at x = 0.5234; nearest maximum 0.0234 awayO6 at x = 0.5234; nearest maximum 0.0234 awayO2 at x = 0.5920; nearest maximum 0.0031 awayO2 at x = 0.5920; nearest maximum 0.0031 awayO4 at x = 0.6802; nearest maximum 0.0087 awayO4 at x = 0.6802; nearest maximum 0.0087 awaySi3 at x = 0.6810; nearest maximum 0.0079 awaySi3 at x = 0.6810; nearest maximum 0.0079 awayO8 at x = 0.6841; nearest maximum 0.0048 awayO8 at x = 0.6841; nearest maximum 0.0048 awayO3 at x = 0.6875; nearest maximum 0.0014 awayO3 at x = 0.6875; nearest maximum 0.0014 awaySi2 at x = 0.6919; nearest maximum 0.0030 awaySi2 at x = 0.6919; nearest maximum 0.0030 awayO7 at x = 0.7064; nearest maximum 0.0175 awayO7 at x = 0.7064; nearest maximum 0.0175 awayNa at x = 0.7317; nearest maximum 0.0428 awayNa at x = 0.7317; nearest maximum 0.0428 awayNa at x = 0.7683; nearest maximum 0.0428 awayNa at x = 0.7683; nearest maximum 0.0428 awayO7 at x = 0.7936; nearest maximum 0.0175 awayO7 at x = 0.7936; nearest maximum 0.0175 awaySi2 at x = 0.8081; nearest maximum 0.0030 awaySi2 at x = 0.8081; nearest maximum 0.0030 awayO3 at x = 0.8125; nearest maximum 0.0014 awayO3 at x = 0.8125; nearest maximum 0.0014 awayO8 at x = 0.8159; nearest maximum 0.0048 awayO8 at x = 0.8159; nearest maximum 0.0048 awaySi3 at x = 0.8190; nearest maximum 0.0079 awaySi3 at x = 0.8190; nearest maximum 0.0079 awayO4 at x = 0.8198; nearest maximum 0.0087 awayO4 at x = 0.8198; nearest maximum 0.0087 awayO2 at x = 0.9080; nearest maximum 0.0031 awayO2 at x = 0.9080; nearest maximum 0.0031 awayO6 at x = 0.9766; nearest maximum 0.0234 awayO6 at x = 0.9766; nearest maximum 0.0234 awayO5 at x = 0.9868; nearest maximum 0.0132 awayO5 at x = 0.9868; nearest maximum 0.0132 awayAl at x = 0.9912; nearest maximum 0.0088 awayAl at x = 0.9912; nearest maximum 0.0088 awayO1 at x = 0.9953; nearest maximum 0.0047 awayO1 at x = 0.9953; nearest maximum 0.0047 awaySi1 at x = 0.9963; nearest maximum 0.0037 awaySi1 at x = 0.9963; nearest maximum 0.0037 awayatoms, projectedx, in fractions of a and in Åelectrons per unit x

Point at a tick along the foot to see which atoms project onto it. Click it to keep it; click it again, click empty space, or press Escape to let go.

low albite, NaAlSi3O8, C1 — 6 non-zero coefficients out of 12, carrying detail to 0.61 Å. 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 amplitudes this page sums are not what a detector records. On a single crystal the chain runs: an integrated count over the time a reflection passes through the diffracting condition, divided by the Lorentz factor for how long that took, divided by the polarisation factor for the beam, corrected for absorption through the crystal and for extinction where a reflection is strong, and scaled — which gives I ∝ |F|², so |F| is a square root away and the sign is already gone. The corrections page works that chain through in detail for a powder diffractometer; the factors have the same names in both geometries and different forms, because the Lorentz factor is a statement about how a reflection sweeps through the Ewald sphere and a spinning crystal does not sweep the way a powder ring does. Take the shape of the chain from there and not the formulae.

The phase problem

A structure factor is a wave, and a wave has a size and a starting point. |F| is the size — how strongly that set of planes scatters. The phase is the starting point: where the crests of that wave sit relative to the origin of the cell. Adding the waves up puts density where crests from many reflections coincide, so the phases are what decide where the atoms are, and the amplitudes only how much scattering there is to place. That is why moving the origin changes every phase and no amplitude: shift the cell by t and F(h) is multiplied by exp(2πi h·t), which turns each wave without resizing it. An individual phase is therefore partly a statement about a choice of origin, which is exactly why direct methods work with combinations — like the triplet on the sign-relation page — that survive that choice.

1 of the 6 coefficients is 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.

Every coefficient in this sum is +|F| or −|F|, so “the phase problem” here is a choice between two signs rather than an angle anywhere on a circle. That is a property of these structures, not of the transform. F(h00) comes out real, and the projection centrosymmetric, exactly when the space group has an operation sending x to −x; every structure offered here has one, including the five with no centre of symmetry of their own — quartz, berlinite, cristobalite, zinc blende and urea. A structure in P1 has none, and its projection is not centrosymmetric: one carbon at 0.1, 0.2, 0.3 gives F(100) a phase of 36° on the structure-factor page. Restore the other two indices and the phases go back to being continuous even here: a general hkl of a non-centrosymmetric structure has a phase that is not 0 or 180°, and recovering it is a harder problem than choosing a sign.

F(200) = 0.3 electrons, d = 3.638 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 0.3. Nothing in it gives the side.F(200) = 0.3 electrons, d = 3.638 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 0.3. Nothing in it gives the side.F(400) = 130.6 electrons, d = 1.819 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 130.6. Nothing in it gives the side.F(400) = 130.6 electrons, d = 1.819 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 130.6. Nothing in it gives the side.F(600) = 123.9 electrons, d = 1.213 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 123.9. Nothing in it gives the side.F(600) = 123.9 electrons, d = 1.213 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 123.9. Nothing in it gives the side.F(800) = -5.9 electrons, d = 0.909 Å. Phase 180 degrees, which is what the side of the line means: this wave is subtracted. A measurement gives an intensity, which correcting and square-rooting turns into the height, 5.9. Nothing in it gives the side.F(800) = -5.9 electrons, d = 0.909 Å. Phase 180 degrees, which is what the side of the line means: this wave is subtracted. A measurement gives an intensity, which correcting and square-rooting turns into the height, 5.9. Nothing in it gives the side.F(1000) = 84.9 electrons, d = 0.728 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 84.9. Nothing in it gives the side.F(1000) = 84.9 electrons, d = 0.728 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 84.9. Nothing in it gives the side.F(1200) = 38.8 electrons, d = 0.606 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 38.8. Nothing in it gives the side.F(1200) = 38.8 electrons, d = 0.606 Å. Phase 0 degrees, which is what the side of the line means: this wave is added. A measurement gives an intensity, which correcting and square-rooting turns into the height, 38.8. Nothing in it gives the side.123456789101112h, in F(h00)F(h00)

Point at any bar for its value, its phase and what the side of the line means. Click it to keep it; click it again, click empty space, or press Escape to let go.

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 8.1372 Å
Atoms in the cell52
Distinct projected sites 52
F(000) 519.90 electrons
Mean of the curve 519.90 — which is F(000) again, and it never entered the series
Coefficients used 6 of 12
Negative coefficients 1
Highest value, phases discarded 1,288.6 at x = 0
Correlation of the two curves 0.9983

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.0606 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.0037 Si1 0.0037 away
0.0047 O1 0.0047 away
0.0088 Al 0.0088 away
0.0132 O5 0.0132 away
0.0234 O6 0.0234 away
0.0920 O2 0.0031 away
0.1802 O4 0.0087 away
0.1810 Si3 0.0079 away
0.1841 O8 0.0048 away
0.1875 O3 0.0014 away
0.1919 Si2 0.0030 away
0.2064 O7 0.0175 away
0.2317 Na 0.0428 away
0.2683 Na 0.0428 away
0.2936 O7 0.0175 away
0.3081 Si2 0.0030 away
0.3125 O3 0.0014 away
0.3159 O8 0.0048 away
0.3190 Si3 0.0079 away
0.3198 O4 0.0087 away
0.4080 O2 0.0031 away
0.4766 O6 0.0234 away
0.4868 O5 0.0132 away
0.4912 Al 0.0088 away
0.4953 O1 0.0047 away
0.4963 Si1 0.0037 away
0.5037 Si1 0.0037 away
0.5047 O1 0.0047 away
0.5088 Al 0.0088 away
0.5132 O5 0.0132 away
0.5234 O6 0.0234 away
0.5920 O2 0.0031 away
0.6802 O4 0.0087 away
0.6810 Si3 0.0079 away
0.6841 O8 0.0048 away
0.6875 O3 0.0014 away
0.6919 Si2 0.0030 away
0.7064 O7 0.0175 away
0.7317 Na 0.0428 away
0.7683 Na 0.0428 away
0.7936 O7 0.0175 away
0.8081 Si2 0.0030 away
0.8125 O3 0.0014 away
0.8159 O8 0.0048 away
0.8190 Si3 0.0079 away
0.8198 O4 0.0087 away
0.9080 O2 0.0031 away
0.9766 O6 0.0234 away
0.9868 O5 0.0132 away
0.9912 Al 0.0088 away
0.9953 O1 0.0047 away
0.9963 Si1 0.0037 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

Where this comes from