xraytools.

Direct Methods and the Sign Relation

The phases are not measured, but they are not free either: the amplitudes constrain them. For three reflections whose indices add up, the product of their signs is +1 more often the stronger they are — and that one relation is what solved small-molecule crystallography. This counts how often it holds on a structure whose answer is already known.

You supply
One of the named structures, and how far the series should run. Everything else — the amplitudes, the true signs, the predicted probability — is computed from the atoms.
Reading it
This counts the relation, it does not solve anything: the signs it checks against come from the published coordinates. A real program fixes an origin, propagates symbols through relations like these and ranks the results — a search whose answer nothing here could check.

Worked examples: caesium chloride — every relation holds · zirconia — the strong ones hold, the weak ones do not · quartz at 12 terms · quartz at 36 terms — more relations, same rule · albite — where no relation is certain

Input

up to h =

How far the series runs. Blank means 24. A triplet needs h + k inside the range, so the number of relations grows roughly as the square of this.

How often the relation holds

0.50.60.70.80.91.0Triple products from 0.56 to 0.65: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 75 per cent.Triple products from 0.56 to 0.65: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 75 per cent.Triple products from 0.65 to 0.70: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 77 per cent.Triple products from 0.65 to 0.70: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 77 per cent.Triple products from 0.70 to 0.76: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 79 per cent.Triple products from 0.70 to 0.76: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 79 per cent.Triple products from 0.76 to 1.63: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 85 per cent.Triple products from 0.76 to 1.63: the relation holds for 29 of 29 triplets, which is 100 per cent. Cochran predicts 85 per cent.Triple products from 1.63 to 1.82: the relation holds for 28 of 28 triplets, which is 100 per cent. Cochran predicts 95 per cent.Triple products from 1.63 to 1.82: the relation holds for 28 of 28 triplets, which is 100 per cent. Cochran predicts 95 per cent.0.00.51.01.5|E(h) E(k) E(h+k)|how often it holds

Point at a marker for that band’s count and what the theory predicted for it. Click it to keep it; click it again, click empty space, or press Escape to let go.

caesium chloride, Pm3m — 144 triplets from 24 non-zero coefficients, carrying detail to 0.17 Å.

Each marker is a band of triplets of similar strength, and the curve is Cochran’s prediction for this structure — not a fit to the points. The line at 0.5 is what guessing would give.

The markers sit above the curve, and that is the formula being conservative rather than wrong. Cochran’s derivation assumes the atoms are randomly placed; this projection is far more regular than that, which the statistic ⟨|E2 − 1|⟩ measures directly — 0.398 here against 0.968 for atoms at random. Symmetry makes the signs agree more often than chance, so the prediction is a floor and not an estimate. The same fact is why the Wilson plot refuses a verdict for several of these structures.

What that says

The relation holds for 144 of 144 triplets here, 100 per cent against the 50 you would get by guessing. Restrict it to the 36 whose predicted probability is above 0.9 and it holds for every one of them. That is the regime a program works in: it does not use the weak relations, it starts from the strong ones.

This does not solve anything. A direct-methods program fixes an origin, gives a few strong reflections symbolic signs, propagates them through relations like these, and ranks the resulting sign sets by a figure of merit — a search whose answer nothing here could check. What this page does instead is measure the relation the search rests on: the signs in the table are computed from the published coordinates, so every prediction can be marked right or wrong, which is the one thing a solver’s output would not allow.

The strongest relations

h k h+k |E1E2E3| predicted signs holds?
2 2 4 1.82 96% + + + yes
2 4 6 1.81 96% + + + yes
2 6 8 1.79 96% + + + yes
4 4 8 1.77 96% + + + yes
4 6 10 1.74 96% + + + yes
2 8 10 1.73 96% + + + yes
6 6 12 1.71 96% + + + yes
2 22 24 1.70 95% + + + yes
2 20 22 1.69 95% + + + yes
4 8 12 1.69 95% + + + yes
6 8 14 1.68 95% + + + yes
4 20 24 1.68 95% + + + yes
2 18 20 1.67 95% + + + yes
6 18 24 1.67 95% + + + yes
4 18 22 1.66 95% + + + yes
2 10 12 1.66 95% + + + yes
8 8 16 1.66 95% + + + yes
6 16 22 1.66 95% + + + yes
6 10 16 1.65 95% + + + yes
4 16 20 1.64 95% + + + yes

Strongest first, which is the order a program uses them in: it starts from the relations it can trust and never reaches the weak ones. Showing 20 of 144; the counts above are over all of them.

Why a product of three

A single sign is not a fact about the crystal. This projection has a centre of symmetry at x = 0 and another at x = ½, and either will do as an origin; move it to the second and F(h) is multiplied by (−1)h, so the sign of every odd reflection reverses — all 12 of them here. The triple product survives it, because (−1)h(−1)k(−1)h+k = +1 whatever h and k are. That is what makes it a structure invariant, and it is the reason the amplitudes can predict it: an amplitude does not move when the origin does, so no function of amplitudes could ever fix a quantity that does.

The reflections

h d / Å |F| E sign sign at the other origin
1 4.123 34.7 0.68 +
2 2.062 52.0 1.22 + +
3 1.374 26.2 0.73 +
4 1.031 37.4 1.21 + +
5 0.825 19.3 0.73 +
6 0.687 28.3 1.22 + +
7 0.589 15.2 0.74 +
8 0.515 22.5 1.20 + +
9 0.458 13.3 0.77 +
10 0.412 18.6 1.17 + +
11 0.375 11.8 0.81 +
12 0.344 15.5 1.16 + +
13 0.317 10.2 0.83 +
14 0.295 13.0 1.15 + +
15 0.275 8.5 0.83 +
16 0.258 10.8 1.15 + +
17 0.243 7.0 0.82 +
18 0.229 9.1 1.16 + +
19 0.217 5.7 0.80 +
20 0.206 7.6 1.17 + +
21 0.196 4.7 0.79 +
22 0.187 6.5 1.18 + +
23 0.179 4.0 0.78 +
24 0.172 5.6 1.18 + +

The E values are normalised in two steps, the same two the Wilson plot uses: divide |F|2 by Σg2 at each reflection’s own resolution, which removes the form-factor decay exactly, then by the mean over the reflections that are present. The second step makes ⟨E2⟩ equal to 1 by construction rather than by fit, and it is doing real work: without it copper comes out at exactly 2, because half its h00 are systematically absent and the survivors carry their share.

What the prediction was built from

Atoms in the cell2
Distinct projected scatterers 2 — atoms sharing an x project onto one another and act as one
Σ2 = Σg2 3,312.7
Σ3 = Σg3 171,178.7
Σ3Σ2−3/2 0.8978 — the whole prediction is tanh of this times the triple product
Equivalent equal scatterers 1.2, since that coefficient is 1/√N for N equal ones
E2 1.0000 — 1 by construction, which is what the second normalisation step buys
⟨|E2 − 1|⟩ 0.398 against 0.968 for a centrosymmetric structure of randomly placed atoms — which is the assumption behind the curve
Negative coefficients 0 of 24

The equivalent scatterer count is the number that decides everything: the relation reads 1/√N, so a structure with one heavy atom among light ones behaves like a much smaller one and gives up its signs easily. That is why the heavy-atom structures fell first, and why a page of equal light atoms is where these methods run out.

Try it

caesium chloride — every relation holds · zirconia — the strong ones hold, the weak ones do not · quartz at 12 terms · quartz at 36 terms — more relations, same rule · albite — where no relation is certain

Where this comes from