![]()
|
5. Phase determination by multiple isomorphous replacementA simple statement of the problem is that given |FP|, |FPH| (measured experimentally), and |FH|, fH (calculated from the site coordinates), the unknown fP is to be determined. Consider first the hypothetical case where all quantities are error-free. 5.1. Representation of the phase problem in the complex plane5.1.1. Centric caseThe structure will always appear centrosymmetric when projected down an evenfold axis. So for example P21 (unique b) has one centric zone: h0l; P21212 has three: hk0, h0l and 0kl; P1 and R3 have none. The centric zone reflections give a ready advantage in MIR studies, because fP can take only 1 of 2 possible values differing by 180° (e.g. 0/180° or 90/270°).
For each possibility of fP, |FPH|(calc) is compared with |FPH|(obs). 5.1.2. Acentric case
Traditionally, the determination of the phase in the acentric case has been taught by "Harker construction" (phase circles). However this is not a practical method for a computer program, and in any case proper account has to be taken of the experimental and other errors. Consequently all software uses probability density functions (see section 5.2).
By the cosine rule: |FPH|(calc) = Ö(|FP|2 + |FH|2 + 2|FP||FH|cos(fH-fP)) 5.1.3. SIRAS caseAtoms show anomalous scattering when the incident x-rays have a wavelength near to the absorption edge of the atoms. Part of the radiation is absorbed by the atom and re-emitted with a change in phase. The scattering factor for the atom can be written in terms of real and imaginary parts: f = f0 + Df' + i Df" where Df' is the dispersion component of the anomalous scattering and Df" the absorption component. The imaginary component scatters out of phase with the primary wave, leading to absorption. Both Df' and Df" vary with the wavelength especially near the absorption edge. The absorption coefficients for C, H, O and N are small enough to be ignored. However, for the heavy metal atoms of a derivative, these coefficients cannot be ignored, especially as the absorption edge for these elements is close to the wavelength of X-rays commonly used. The anomalous differences are largest when |FP| and |FH| are approximately perpendicular.
5.2. Probability density function5.2.1. SIR caseThe probability density function (PDF) is defined as: Pj(fP) = Kj exp(-|Dj|2 / (2 s2(|FP|) + 2 s2(|FPH|) + s2(|Dj|))) where Kj is a normalisation factor and s(|Dj|) is the r.m.s. lack of closure weighted by the phase probability. Note that this formulation ignores zone enhancement, centricity and non-isomorphism factors (see Read, R. in 1991 CCP4 Study Weekend Proceedings, pp 69-79). 5.2.2. SIRAS/MIR caseThe PDF is defined as P(fP) = Pj Pj(fP) Most programs use this, but the equation assumes independent information from each derivative, and therefore over-emphasises the contribution of FP (see R. Read, 1991 for rigorous derivation). 5.2.3. Visualisation of the PDF'sIf the PDF's P(fP) and Pj(fP) are plotted as functions of fP, either in Cartesian coordinates or in polar coordinates, P will have maxima at the most probable phase(s). Also in the polar plots the centroid of the probability-weighted unit circle is shown. The phase of the centroid is the best phase (minimum r.m.s. error in electron density). The radius of the centroid is the figure of merit m » cos(phase error). For electron density maps the Fourier coefficient m|FP|exp(i fP) is used; for difference maps m|FPH-FP|exp(i fP) is used. Example 1For better resolution, click on the pictures.
The figures show the individual PDF's (P1, P2, P3) and product PDF's (P123, P12) for one reflection from 3 simulated derivatives, for |FP| = 10.0.
Here, P1 and P2 are bimodal, whereas P3 is unimodal because |FPH3| = |FP| + |FH3|. The phasing power |FH|/s(|D|) determines the sharpness of the density function at the most probable phase, whereas the figure of merit determines the precision of the best phase. Multiplication of probabilities improves both the figure of merit and the phasing power. Addition of the low phasing power derivative 3 in this case makes little difference (compare P123 and P12). Example 2For better resolution, click on the pictures.
The figures show the effect of changing fH to 340° for derivative 3 (P3), while keeping all the other parameters the same. Even though derivative 3 has low phasing power, this considerably reduces the overall figure of merit, so even one poorly phased derivative can negate the accumulated effect of several good derivatives. Example 3For better resolution, click on the pictures.
The figures show the effect of varying the phasing power in the SIR case. fH has also been varied in order to separate the centroids in order to make the effect clear.
From this it can be seen that a phasing power below 0.5 does not contribute significantly to the phasing; often a cutoff of 1.0 is used, but at this level the contribution to the phasing is clearly significant.
|