International
Tables for Crystallography Volume B Reciprocal space Edited by U. Shmueli © International Union of Crystallography 2010 
International Tables for Crystallography (2010). Vol. B, ch. 2.2, pp. 215243
https://doi.org/10.1107/97809553602060000764 Chapter 2.2. Direct methods ^{a}Dipartimento Geomineralogico, Campus Universitario, 70125 Bari, Italy, and Institute of Crystallography, Via G. Amendola, 122/O, 70125 Bari, Italy Direct methods are essentially reciprocalspace techniques, developed historically to solve smallmolecule crystal structures. Their success in this area (in practice, they solve the phase problem) is based on numerous theoretical achievements which concern origin specification (Section 2.2.3), the concepts of structure invariants and seminvariants, normalization of the structure factors (Section 2.2.4), inequalities among structure factors, and probabilistic phase relationships (Section 2.2.5). Probabilistic phase relationships are at the core of direct methods: triplet (via Cochran and via the P10 formula), quartet (according to Hauptman and to Giacovazzo) and quintet invariant phase estimates are described, along with estimates of onephase and of twophase structure seminvariants (via representation theory). Determinantal formulas are also quoted. Techniques for estimating invariants by using some prior information and their realspace counterparts are discussed. The success of direct methods is intimately connected with the phasing procedures. The most important tools of the procedures (e.g. the tangent formula, magic integers, randomstart approaches and figures of merit for recognizing the correct solution) are analysed. The integration of direct methods with macromolecular crystallography is discussed in Section 2.2.10: in particular we refer to the ab initio methods for solving protein structures (in combination with directspace techniques like electrondensity modification, envelope determination, histogram matching etc.) as well as to the combination of direct methods with SIR–MIR and SAD–MAD techniques (both in `onestep' procedures and in the twostep method, the latter requiring the prior determination of the heavyatom or anomalousscatterer substructure). 
Direct methods are today the most widely used tool for solving small crystal structures. They work well both for equalatom molecules and when a few heavy atoms exist in the structure. In recent years the theoretical background of direct methods has been improved to take into account a large variety of prior information (the form of the molecule, its orientation, a partial structure, the presence of pseudosymmetry or of a superstructure, the availability of isomorphous data or of data affected by anomalousdispersion effects, …). Owing to this progress and to the increasing availability of powerful computers, the phase problem for small molecules has been solved in practice: a number of effective, highly automated packages are today available to the scientific community.
The combination of direct methods with socalled directspace methods have recently allowed the ab initio crystal structure solution of proteins. The present limit of complexity is about 2500 nonhydrogen atoms in the asymmetric unit, but diffraction data at atomic resolution (~1 Å) are required. Trials are under way to bring this limit to 1.5 Å and have shown some success.
The theoretical background and tables useful for origin specification are given in Section 2.2.3; in Section 2.2.4 the procedures for normalizing structure factors are summarized. Phasedetermining formulae (inequalities, probabilistic formulae for triplet, quartet and quintet invariants, and for one and twophase s.s.'s, determinantal formulae) are given in Section 2.2.5. In Section 2.2.6 the connection between direct methods and related techniques in real space is discussed. Practical procedures for solving smallmolecule crystal structures are described in Sections 2.2.7 and 2.2.8, and references to the most extensively used packages are given in Section 2.2.9. The integration of direct methods, isomorphous replacement and anomalousdispersion techniques is briefly discussed in Section 2.2.10.
The reader interested in a more detailed description of the topic is referred to a recent textbook (Giacovazzo, 1998).

The normalized structure factors E (see also Chapter 2.1 ) are calculated according to (Hauptman & Karle, 1953) where is the squared observed structurefactor magnitude on the absolute scale and is the expected value of .
depends on the available a priori information. Often, but not always, this may be considered as a combination of several typical situations. We mention:

When probability theory is not used, the quasinormalized structure factors and the unitary structure factors are often used. and are defined according to Since is the largest possible value for represents the fraction of with respect to its largest possible value. Therefore If atoms are equal, then .
N.s.f.'s cannot be calculated by applying (2.2.4.1) to observed s.f.'s because: (a) the observed magnitudes (already corrected for Lp factor, absorption, …) are on a relative scale; (b) cannot be calculated without having estimated the vibrational motion of the atoms.
This is usually obtained by the well known Wilson plot (Wilson, 1942), according to which observed data are divided into ranges of and averages of intensity are taken in each shell. Reflection multiplicities and other effects of spacegroup symmetry on intensities must be taken into account when such averages are calculated. The shells are symmetrically overlapped in order to reduce statistical fluctuations and are restricted so that the number of reflections in each shell is reasonably large. For each shell should be obtained, where K is the scale factor needed to place Xray intensities on the absolute scale, B is the overall thermal parameter and is the expected value of in which it is assumed that all the atoms are at rest. depends upon the structural information that is available (see Section 2.2.4.1 for some examples).
Equation (2.2.4.3) may be rewritten as which plotted at various should be a straight line of which the slope (2B) and intercept (ln K) on the logarithmic axis can be obtained by applying a linear leastsquares procedure.
Very often molecular geometries produce perceptible departures from linearity in the logarithmic Wilson plot. However, the more extensive the available a priori information on the structure is, the closer, on the average, are the Wilsonplot curves to their leastsquares straight lines.
Accurate estimates of B and K require good strategies (Rogers & Wilson, 1953) for:

Once K and B have been estimated, values can be obtained from experimental data by where is the expected value of for the reflection h on the basis of the available a priori information.
Under some fairly general assumptions (see Chapter 2.1 ) probability distribution functions for the variable for cs. and ncs. structures are (see Fig. 2.2.4.1) andrespectively. Corresponding cumulative functions are (see Fig. 2.2.4.2)
Some moments of the distributions (2.2.4.4) and (2.2.4.5) are listed in Table 2.2.4.1. In the absence of other indications for a given crystal structure, a cs. or an ncs. space group will be preferred according to whether the statistical tests yield values closer to column 2 or to column 3 of Table 2.2.4.1.

For further details about the distribution of intensities see Chapter 2.1 .
From the earliest periods of Xray structure analysis several authors (Ott, 1927; Banerjee, 1933; Avrami, 1938) have tried to determine atomic positions directly from diffraction intensities. Significant developments are the derivation of inequalities and the introduction of probabilistic techniques via the use of joint probability distribution methods (Hauptman & Karle, 1953).
An extensive system of inequalities exists for the coefficients of a Fourier series which represents a positive function. This can restrict the allowed values for the phases of the s.f.'s in terms of measured structurefactor magnitudes. Harker & Kasper (1948) derived two types of inequalities:
Type 1. A modulus is bound by a combination of structure factors: where m is the order of the point group and .
Applied to loworder space groups, (2.2.5.1) gives The meaning of each inequality is easily understandable: in , for example, must be positive if is large enough.
Type 2. The modulus of the sum or of the difference of two structure factors is bound by a combination of structure factors: where stands for `real part of'. Equation (2.2.5.2) applied to P1 gives
A variant of (2.2.5.2) valid for cs. space groups is After Harker & Kasper's contributions, several other inequalities were discovered (Gillis, 1948; Goedkoop, 1950; Okaya & Nitta, 1952; de Wolff & Bouman, 1954; Bouman, 1956; Oda et al., 1961). The most general are the Karle–Hauptman inequalities (Karle & Hauptman, 1950): The determinant can be of any order but the leading column (or row) must consist of U's with different indices, although, within the column, symmetryrelated U's may occur. For and , equation (2.2.5.3) reduces to which, for cs. structures, gives the Harker & Kasper inequality For , equation (2.2.5.3) becomes from which where If the moduli , , are large enough, (2.2.5.4) is not satisfied for all values of . In cs. structures the eventual check that one of the two values of does not satisfy (2.2.5.4) brings about the unambiguous identification of the sign of the product .
It was observed (Gillis, 1948) that `there was a number of cases in which both signs satisfied the inequality, one of them by a comfortable margin and the other by only a relatively small margin. In almost all such cases it was the former sign which was the correct one. That suggests that the method may have some power in reserve in the sense that there are still fundamentally stronger inequalities to be discovered'. Today we identify this power in reserve in the use of probability theory.
For any space group (see Section 2.2.3) there are linear combinations of phases with cosines that are, in principle, fixed by the magnitudes alone (s.i.'s) or by the values and the trigonometric form of the structure factor (s.s.'s). This result greatly stimulated the calculation of conditional distribution functions where , is an s.i. or an s.s. and is a suitable set of diffraction magnitudes. The method was first proposed by Hauptman & Karle (1953) and was developed further by several authors (Bertaut, 1955a,b, 1960; Klug, 1958; Naya et al., 1964, 1965; Giacovazzo, 1980a). From a probabilistic point of view the crystallographic problem is clear: the joint distribution , from which the conditional distributions (2.2.5.5) can be derived, involves a number of normalized structure factors each of which is a linear sum of random variables (the atomic contributions to the structure factors). So, for the probabilistic interpretation of the phase problem, the atomic positions and the reciprocal vectors may be considered as random variables. A further problem is that of identifying, for a given Φ, a suitable set of magnitudes on which Φ primarily depends. The formulation of the nested neighbourhood principle first (Hauptman, 1975) fixed the idea of defining a sequence of sets of reflections each contained in the succeeding one and having the property that any s.i. or s.s. may be estimated via the magnitudes constituting the various neighbourhoods. A subsequent more general theory, the representation method (Giacovazzo, 1977a, 1980b), arranges for any Φ the set of intensities in a sequence of subsets in order of their expected effectiveness (in the statistical sense) for the estimation of Φ.
In the following sections the main formulae estimating loworder invariants and seminvariants or relating phases to other phases and diffraction magnitudes are given.
The basic formula for the estimation of the triplet phase given the parameter is Cochran's (1955) formula where , is the atomic number of the jth atom and is the modified Bessel function of order n. In Fig. 2.2.5.1 the distribution is shown for different values of G.
The conditional probability distribution for , given a set of and , is given (Karle & Hauptman, 1956; Karle & Karle, 1966) bywhere is the most probable value for . The variance of may be obtained from (2.2.5.7) and is given by which is plotted in Fig. 2.2.5.2.
Equation (2.2.5.9) is the socalled tangent formula. According to (2.2.5.10), the larger is α the more reliable is the relation .
For an equalatom structure .
The basic conditional formula for sign determination of in cs. crystals is Cochran & Woolfson's (1955) formula where is the probability that is positive and k ranges over the set of known values . The larger the absolute value of the argument of tanh, the more reliable is the phase indication.
An auxiliary formula exploiting all the 's in reciprocal space in order to estimate a single Φ is the formula (Hauptman & Karle, 1958; Karle & Hauptman, 1958) given by where C is a constant which differs for cs. and ncs. crystals, is the average value of and p is normally chosen to be some small number. Several modifications of (2.2.5.12) have been proposed (Hauptman, 1964, 1970; Karle, 1970a; Giacovazzo, 1977b).
A recent formula (Cascarano, Giacovazzo, Camalli et al., 1984) exploits information contained within the second representation of Φ, that is to say, within the collection of special quintets (see Section 2.2.5.6): where k is a free vector. The formula retains the same algebraic form as (2.2.5.6), but where , is assumed to be zero if it is experimentally negative. The prime to the summation warns the reader that precautions have to be taken in order to avoid duplications in the contributions.
G may be positive or negative. In particular, if the triplet is estimated negative.
The accuracy with which the value of Φ is estimated strongly depends on . Thus, in practice, only a subset of reciprocal space (the reflections k with large values of ) may be used for estimating Φ.
(2.2.5.13) proved to be quite useful in practice. Positive triplet cosines are ranked in order of reliability by (2.2.5.13) markedly better than by Cochran's parameters. Negative estimated triplet cosines may be excluded from the phasing process and may be used as a figure of merit for finding the correct solution in a multisolution procedure.
A strength of direct methods is that no knowledge of structure is required for their application. However, when some a priori information is available, it should certainly be a weakness of the methods not to make use of this knowledge. The conditional distribution of Φ given and the first three of the five kinds of a priori information described in Section 2.2.4.1 is (Main, 1976; Heinermann, 1977a) where stand for h, , , and for . The quantities have been calculated in Section 2.2.4.1 according to different categories: is a suitable average of the product of three scattering factors for the ith atomic group, p is the number of atomic groups in the cell including those related by symmetry elements. We have the following categories.

In early papers (Hauptman & Karle, 1953; Simerska, 1956) the phase was always expected to be zero. Schenk (1973a,b) [see also Hauptman (1974)] suggested that Φ primarily depends on the seven magnitudes: , called basis magnitudes, and , called cross magnitudes.
The conditional probability of Φ in P1 given seven magnitudes according to Hauptman (1975) is where L is a suitable normalizing constant which can be derived numerically, For equal atoms . Denoting gives Fig. 2.2.5.3 shows the distribution (2.2.5.18) for three typical cases. It is clear from the figure that the cosine estimated near π or in the middle range will be in poorer agreement with the true values than the cosine near 0 because of the relatively larger values of the variance. In principle, however, the formula is able to estimate negative or enantiomorphsensitive quartet cosines from the seven magnitudes.

Distributions (2.2.5.18) (solid curve) and (2.2.5.20) (dashed curve) for the indicated values in three typical cases. 
In the cs. case (2.2.5.18) is replaced (Hauptman & Green, 1976) by where is the probability that the sign of is positive or negative, and The normalized probability may be derived by . More simple probabilistic formulae were derived independently by Giacovazzo (1975, 1976): where and . Q is never allowed to be negative.
According to (2.2.5.20) is expected to be positive or negative according to whether is positive or negative: the larger is C, the more reliable is the phase indication. For , (2.2.5.18) and (2.2.5.20) are practically equivalent in all cases. If N is small, (2.2.5.20) is in good agreement with (2.2.5.18) for quartets strongly defined as positive or negative, but in poor agreement for enantiomorphsensitive quartets (see Fig. 2.2.5.3).
In cs. cases the sign probability for is where G is defined by (2.2.5.21).
All three cross magnitudes are not always in the set of measured reflections. From marginal distributions the following formulae arise (Giacovazzo, 1977c; Heinermann, 1977b):
Equations (2.2.5.20) and (2.2.5.23) are easily modifiable when some cross magnitudes are not in the measurements. If is not measured then (2.2.5.20) or (2.2.5.23) are still valid provided that in G it is assumed that . For example, if and are not in the data then (2.2.5.21) and (2.2.5.22) become In space groups with symmetry higher than more symmetryequivalent quartets can exist of the type where are rotation matrices of the space group. The set is called the first representation of Φ. In this case Φ primarily depends on more than seven magnitudes. For example, let us consider in Pmmm the quartet Quartets symmetry equivalent to Φ and respective cross terms are given in Table 2.2.5.1.

Experimental tests on the application of the representation concept to quartets have been made (Busetta et al., 1980). It was shown that quartets with more than three cross magnitudes are more accurately estimated than other quartets. Also, quartets with a cross reflection which is systematically absent were shown to be of significant importance in direct methods. In this context it is noted that systematically absent reflections are not usually included in the set of diffraction data. This custom, not exceptionable when only triplet relations are used, can give rise to a loss of information when quartets are used. In fact the usual programs of direct methods discard quartets as soon as one of the cross reflections is not measured, so that systematic absences are dealt with in the same manner as those reflections which are outside the sphere of measurements.
A quintet phase may be considered as the sum of three suitable triplets or the sum of a triplet and a quartet, i.e. or It depends primarily on 15 magnitudes: the five basis magnitudesand the ten cross magnitudesIn the following we will denoteConditional distributions of Φ in P1 and given the 15 magnitudes have been derived by several authors and allow in favourable circumstances in ncs. space groups the quintets having Φ near 0 or near π or near to be identified. Among others, we remember:

For cs. cases (2.2.5.24) reduces to Positive or negative quintets may be identified according to whether G is larger or smaller than zero.
If is not measured then (2.2.5.24) and (2.2.5.25) are still valid provided that in (2.2.5.25) .
If the symmetry is higher than in then more symmetryequivalent quintets can exist of the type where are rotation matrices of the space groups. The set is called the first representation of Φ. In this case Φ primarily depends on more than 15 magnitudes which all have to be taken into account for a careful estimation of Φ (Giacovazzo, 1980a).
A wide use of quintet invariants in directmethods procedures is prevented for two reasons: (a) the large correlation of positive quintet cosines with positive triplets; (b) the large computing time necessary for their estimation [quintets are phase relationships of order , so a large number of quintets have to be estimated in order to pick up a sufficient percentage of reliable ones].
In a crystal structure with N identical atoms the joint probability distribution of n normalized s.f.'s under the following conditions:
is given (Tsoucaris, 1970) [see also Castellano et al. (1973) and Heinermann et al. (1979)] by for cs. structures and for ncs. structures. In (2.2.5.27) and (2.2.5.28) we have denoted is an element of , and is the covariance matrix with elements is a K–H determinant: therefore . Let us call the K–H determinant obtained by adding to the last column and line formed by , and , respectively. Then (2.2.5.27) and (2.2.5.28) may be written and respectively. Because is a constant, the maximum values of the conditional joint probabilities (2.2.5.29) and (2.2.5.30) are obtained when is a maximum. Thus the maximum determinant rule may be stated (Tsoucaris, 1970; Lajzérowicz & Lajzérowicz, 1966): among all sets of phases which are compatible with the inequality the most probable one is that which leads to a maximum value of .
If only one phase, i.e. , is unknown whereas all other phases and moduli are known then (de Rango et al., 1974; Podjarny et al., 1976) for cs. crystals and for ncs. crystals where Equations (2.2.5.31) and (2.2.5.32) generalize (2.2.5.11) and (2.2.5.7), respectively, and reduce to them for . Fourthorder determinantal formulae estimating triplet invariants in cs. and ncs. crystals, and making use of the entire data set, have recently been secured (Karle, 1979, 1980a).
Advantages, limitations and applications of determinantal formulae can be found in the literature (Heinermann et al., 1979; de Rango et al., 1975, 1985). Taylor et al. (1978) combined K–H determinants with a magicinteger approach. The computing time, however, was larger than that required by standard computing techniques. The use of K–H matrices has been made faster and more effective by de Gelder et al. (1990) (see also de Gelder, 1992). They developed a phasing procedure (CRUNCH) which uses random phases as starting points for the maximization of the K–H determinants.
According to the representations method (Giacovazzo, 1977a, 1980a,b):

The more general expressions for the s.s.'s of first rank are
In other words:

The set of special quartets (2.2.5.35a) and (2.2.5.35b) constitutes the first representations of Φ.
Structure seminvariants of the second rank can be characterized as follows: suppose that, for a given seminvariant Φ, it is not possible to find a vectorial index h and a rotation matrix such that is a structure invariant. Then Φ is a structure seminvariant of the second rank and a set of structure invariants ψ can certainly be formed, of type by means of suitable indices h and l and rotation matrices and . As an example, for symmetry class 222, or or are s.s.'s of the first rank while is an s.s. of the second rank.
The procedure may easily be generalized to s.s.'s of any order of the first and of the second rank. So far only the role of onephase and twophase s.s.'s of the first rank in direct procedures is well documented (see references quoted in Sections 2.2.5.9 and 2.2.5.10).
Let be our onephase s.s. of the first rank, where In general, more than one rotation matrix and more than one vector h are compatible with (2.2.5.36). The set of special triplets is the first representation of . In cs. space groups the probability that , given and the set , may be estimated (Hauptman & Karle, 1953; Naya et al., 1964; Cochran & Woolfson, 1955) by where In (2.2.5.37), the summation over n goes within the set of matrices for which (2.2.5.35a,b) is compatible, and h varies within the set of vectors which satisfy (2.2.5.36) for each . Equation (2.2.5.36) is actually a generalized way of writing the socalled relationships (Hauptman & Karle, 1953).
If is a phase restricted by symmetry to and in an ncs. space group then (Giacovazzo, 1978) If is a general phase then is distributed according to where with a reliability measured by The second representation of is the set of special quintets provided that h and vary over the vectors and matrices for which (2.2.5.36) is compatible, k over the asymmetric region of the reciprocal space, and over the rotation matrices in the space group. Formulae estimating via the second representation in all the space groups [all the base and cross magnitudes of the quintets (2.2.5.40) now constitute the a priori information] have been secured (Giacovazzo, 1978; Cascarano & Giacovazzo, 1983; Cascarano, Giacovazzo, Calabrese et al., 1984). Such formulae contain, besides the contribution of order provided by the first representation, a supplementary (not negligible) contribution of order arising from quintets.
Denoting formulae (2.2.5.37), (2.2.5.38), (2.2.5.39) still hold provided that is replaced by where m is the number of symmetry operators and is the Hermite polynomial of order four.
is assumed to be zero if it is computed negative. The prime to the summation warns the reader that precautions have to be taken in order to avoid duplication in the contributions.
Twophase s.s.'s of the first rank were first evaluated in some cs. space groups by the method of coincidence by Grant et al. (1957); the idea was extended to ncs. space groups by Debaerdemaeker & Woolfson (1972), and in a more general way by Giacovazzo (1977e,f).
The technique was based on the combination of the two triplets which, subtracted from one another, give If all four 's are sufficiently large, an estimate of the twophase seminvariant is available.
Probability distributions valid in according to the neighbourhood principle have been given by Hauptman & Green (1978). Finally, the theory of representations was combined by Giacovazzo (1979a) with the joint probability distribution method in order to estimate twophase s.s.'s in all the space groups.
According to representation theory, the problem is that of evaluating via the special quartets (2.2.5.35a) and (2.2.5.35b). Thus, contributions of order will appear in the probabilistic formulae, which will be functions of the basis and of the cross magnitudes of the quartets (2.2.5.35) . Since more pairs of matrices and can be compatible with (2.2.5.34), and for each pair more pairs of vectors and may satisfy (2.2.5.34), several quartets can in general be exploited for estimating Φ. The simplest case occurs in where the two quartets (2.2.5.35) suggest the calculation of the sixvariate distribution function which leads to the probability formula where is the probability that the product is positive, and It may be seen that in favourable cases .
For the sake of brevity, the probabilistic formulae for the general case are not given and the reader is referred to the original papers.
The statistical treatment suggested by Wilson for scaling observed intensities corresponds, in direct space, to the origin peak of the Patterson function, so it is not surprising that a general correspondence exists between probabilistic formulation in reciprocal space and algebraic properties in direct space.
For a structure containing atoms which are fully resolved from one another, the operation of raising to the nth power retains the condition of resolved atoms but changes the shape of each atom. Let where is an atomic function and is the coordinate of the `centre' of the atom. Then the Fourier transform of the electron density can be written as If the atoms do not overlap and its Fourier transform gives is the scattering factor for the jth peak of :
We now introduce the condition that all atoms are equal, so that and for any j. From (2.2.6.1) and (2.2.6.2) we may write where is a function which corrects for the difference of shape of the atoms with electron distributions and . Since the Fourier transform of both sides gives from which the following relation arises: For , equation (2.2.6.4) reduces to Sayre's (1952) equation [but see also Hughes (1953)] If the structure contains resolved isotropic atoms of two types, P and Q, it is impossible to find a factor such that the relation holds, since this would imply values of such that and simultaneously. However, the following relationship can be stated (Woolfson, 1958): where and are adjustable parameters of . Equation (2.2.6.6) can easily be generalized to the case of structures containing resolved atoms of more than two types (von Eller, 1973).
Besides the algebraic properties of the electron density, Patterson methods also can be developed so that they provide phase indications. For example, it is possible to find the reciprocal counterpart of the function For the function (2.2.6.7) coincides with the usual Patterson function ; for , (2.2.6.7) reduces to the double Patterson function introduced by Sayre (1953). Expansion of as a Fourier series yields Vice versa, the value of a triplet invariant may be considered as the Fourier transform of the double Patterson.
Among the main results relating direct and reciprocalspace properties it may be remembered:

A traditional procedure for phase assignment may be schematically presented as follows:

In very complex structures a large initial set of known phases seems to be a basic requirement for a structure to be determined. This aim can be achieved, for example, by introducing a large number of permutable phases into the initial set. However, the introduction of every new symbol implies a fourfold increase in computing time, which, even in fast computers, quickly leads to computingtime limitations. On the other hand, a relatively large starting set is not in itself enough to ensure a successful structure determination. This is the case, for example, when the triplet invariants used in the initial steps differ significantly from zero. New strategies have therefore been devised to solve more complex structures.

Some references for directmethods packages are given below. Other useful packages using symbolic addition or multisolution procedures do exist but are not well documented.
CRUNCH: Gelder, R. de, de Graaff, R. A. G. & Schenk, H. (1993). Automatic determination of crystal structures using Karle–Hauptman matrices. Acta Cryst. A49, 287–293.
DIRDIF: Beurskens, P. T., Beurskens G., de Gelder, R., GarciaGranda, S., Gould, R. O., Israel, R. & Smits, J. M. M. (1999). The DIRDIF99 program system. Crystallography Laboratory, University of Nijmegen, The Netherlands.
MITHRIL: Gilmore, C. J. (1984). MITHRIL. An integrated directmethods computer program. J. Appl. Cryst. 17, 42–46.
MULTAN88: Main, P., Fiske, S. J., Germain, G., Hull, S. E., Declercq, J.P., Lessinger, L. & Woolfson, M. M. (1999). Crystallographic software: teXsan for Windows. http://www.rigaku.com/downloads/journal/Vol15.1.1998/texsan.pdf .
PATSEE: Egert, E. & Sheldrick, G. M. (1985). Search for a fragment of known geometry by integrated Patterson and direct methods. Acta Cryst. A41, 262–268.
SAPI: Fan, H.F. (1999). Crystallographic software: teXsan for Windows. http://www.rigaku.com/downloads/journal/Vol15.1.1998/texsan.pdf .
SnB: Weeks, C. M. & Miller, R. (1999). The design and implementation of SnB version 2.0. J. Appl. Cryst. 32, 120–124.
SHELX97 and SHELXS: Sheldrick, G. M. (2000). The SHELX home page. http://shelx.uniac.gwdg.de/SHELX/ .
SHELXD: Sheldrick, G. M. (1998). SHELX: applications to macromolecules. In Direct methods for solving macromolecular structures, edited by S. Fortier, pp. 401–411. Dordrecht: Kluwer Academic Publishers.
SIR97: Altomare, A., Burla, M. C., Camalli, M., Cascarano, G. L., Giacovazzo, C., Guagliardi, A., Moliterni, A. G. G., Polidori, G. & Spagna, R. (1999). SIR97: a new tool for crystal structure determination and refinement. J. Appl. Cryst. 32, 115–119.
SIR2004: Burla, M. C., Caliandro, R., Camalli, M., Carrozzini, B., Cascarano, G. L., De Caro, L., Giacovazzo, C., Polidori, G. & Spagna, R. (2005). SIR2004: an improved tool for crystal structure determination and refinement. J. Appl. Cryst. 38, 381–388.
XTAL3.6.1: Hall, S. R., du Boulay, D. J. & OlthofHazekamp, R. (1999). Xtal3.6 crystallographic software. http://xtal.sourceforge.net/ .
The smallest protein molecules contain about 400 nonhydrogen atoms, so they cannot be solved ab initio by the algorithms specified in Sections 2.2.7 and 2.2.8. However, traditional direct methods are applied for:
The application of standard tangent techniques to (a) and (b) has not been found to be very satisfactory (Coulter & Dewar, 1971; Hendrickson et al., 1973; Weinzierl et al., 1969). Tangent methods, in fact, require atomicity and nonnegativity of the electron density. Both these properties are not satisfied if data do not extend to atomic resolution (d > 1.2 Å). Because of series termination and other errors the electrondensity map at d > 1.2 Å presents large negative regions which will appear as false peaks in the squared structure. However, tangent methods use only a part of the information given by the Sayre equation (2.2.6.5). In fact, (2.2.6.5) express two equations relating the radial and angular parts of the two sides, so obtaining a large degree of overdetermination of the phases. To achieve this Sayre (1972) [see also Sayre & Toupin (1975)] suggested minimizing (2.2.10.1) by least squares as a function of the phases: Even if tests on rubredoxin (extensions of phases from 2.5 to 1.5 Å resolution) and insulin (Cutfield et al., 1975) (from 1.9 to 1.5 Å resolution) were successful, the limitations of the method are its high cost and, especially, the higher efficiency of the leastsquares method. Equivalent considerations hold for the application of determinantal methods to proteins [see Podjarny et al. (1981); de Rango et al. (1985) and literature cited therein].
A question now arises: why is the tangent formula unable to solve protein structures? Fan et al. (1991) considered the question from a firstprinciple approach and concluded that:
Sheldrick (1990) suggested that direct methods are not expected to succeed if fewer than half of the reflections in the range 1.1–1.2 Å are observed with (a condition seldom satisfied by protein data).
The most complete analysis of the problem has been made by Giacovazzo, Guagliardi et al. (1994). They observed that the expected value of α (see Section 2.2.7) suggested by the tangent formula for proteins is comparable with the variance of the α parameter. In other words, for proteins the signal determining the phase is comparable with the noise, and therefore the phase indication is expected to be unreliable.
Quite relevant results have recently been obtained by integrating direct methods with some additional experimental information. In particular, we will describe the combination of direct methods with:
Point (d) will not be treated here, as it is described extensively in IT F, Part 13 .
Ab initio techniques do not require prior information of any atomic positions. The recent tremendous increase in computing speed led to direct methods evolving towards the rapid development of multisolution techniques. The new algorithms of the program ShakeandBake (Weeks et al., 1994; Weeks & Miller, 1999; Hauptman et al., 1999) allowed an impressive extension of the structural complexity amenable to direct phasing. In particular we mention: (a) the minimal principle (De Titta et al., 1994), according to which the phase problem is considered as a constrained global optimization problem; (b) the refinement procedure, which alternately uses direct and reciprocalspace techniques; and (c) the parametershift optimization technique (Bhuiya & Stanley, 1963), which aims at reducing the value of the minimal function (Hauptman, 1991; De Titta et al., 1994). An effective variant of Shakeand Bake is SHELXD (Sheldrick, 1998) which cyclically alternates tangent refinement in reciprocal space with peaklist optimisation procedures in real space (Sheldrick & Gould, 1995). Detailed information on these programs is available in IT F (2001), Part 16 .
A different approach is used by ACORN (Foadi et al., 2000), which first locates a small fragment of the molecule (eventually by molecularreplacement techniques) to obtain a useful nonrandom starting set of phases, and then refines them by means of solventflattening techniques.
The program SIR2004 (Burla et al., 2005) uses the tangent formula as well as automatic Patterson techniques to obtain a first imperfect structural model; then directspace techniques are used to refine the model. The Patterson approach is based on the use of the superposition minimum function (Buerger, 1959; Richardson & Jacobson, 1987; Sheldrick, 1992; Pavelcík, 1988; Pavelcík et al., 1992; Burla et al., 2004). It may be worth noting that even this approach is of multisolution type: up to 20 trial solutions are provided by using as pivots the highest maxima in the superposition minimum function.
It is today possible to solve structures up to 2500 nonhydrogen atoms in the asymmetric unit provided data at atomic (about 1 Å) resolution are available. Proteins with data at quasiatomic resolution (say up to 1.5–1.6 Å) can also be solved, but with greater difficulties (Burla et al., 2005). A simple evaluation of the potential of the ab initio techniques suggests that the structural complexity range and the resolution limits amenable to the ab initio approach could be larger in the near future. The approach will profit by general technical advances like the increasing speed of computers and by the greater efficiency of informatic tools (e.g. faster Fouriertransform techniques). It could also profit from new specific crystallographic algorithms (for example, Oszlányi & Süto, 2004). It is of particular interest that extrapolating moduli and phases of nonmeasured reflections beyond the experimental resolution limit makes the ab initio phasing process more efficient, and leads to crystal structure solution even in cases in which the standard programs do not succeed (Caliandro et al., 2005a,b). Moreover, the use of the extrapolated values improves the quality of the final electrondensity maps and makes it easier to recognize the correct one among several trial structures.
SIR–MIR cases are characterized by a situation in which there is one native protein and one or more heavyatom substructures. In this situation the phasing procedure may be a twostep process: in the first stage the heavyatom positions are identified by Patterson techniques (Rossmann, 1961; Okaya et al., 1955) or by direct methods (Mukherjee et al., 1989). In the second step the protein phases are estimated by exploiting the substructure information. Direct methods are able to contribute to both steps (see Sections 2.2.10.5 and 2.2.10.6). In Section 2.2.10.4 we show that direct methods are also able to suggest alternative onestep procedures by estimating structure invariants from isomorphous data.
The theoretical basis was established by Hauptman (1982a): his primary interest was to establish the twophase and threephase structure invariants by exploiting the experimental information provided by isomorphous data. The protein phases could be directly assigned via a tangent procedure.
Let us denote the modulus of the isomorphous difference aswhere the subscripts d and p denote the derivative and the protein, respectively.
Denote also by and atomic scattering factors for the atom labelled j in a pair of isomorphous structures, and let and denote corresponding normalized structure factors. Then where The conditional probability of the twophase structure invariant given and is (Hauptman, 1982a)where Threephase structure invariants were evaluated by considering that eight invariants exist for a given triple of indices h, k, l : So, for the estimation of any , the joint probability distribution has to be studied, from which eight conditional probability densities can be obtained: for .
The analytical expressions of are too intricate and are not given here (the reader is referred to the original paper). We only say that may be positive or negative, so that reliable triplet phase estimates near 0 or near π are possible: the larger , the more reliable the phase estimate.
A useful interpretation of the formulae in terms of experimental parameters was suggested by Fortier et al. (1984): according to them, distributions do not depend, as in the case of the traditional threephase invariants, on the total number of atoms per unit cell but rather on the scattering difference between the native protein and the derivative (that is, on the scattering of the heavy atoms in the derivative).
Hauptman's formulae were generalized by Giacovazzo et al. (1988): the new expressions were able to take into account the resolution effects on distribution parameters. The formulae are completely general and include as special cases native protein and heavyatom isomorphous derivatives as well as Xray and neutron diffraction data. Their complicated algebraic forms are easily reduced to a simple expression in the case of a native protein heavyatom derivative: in particular, the reliability parameter for is where indices p and H warn that parameters have to be calculated over protein atoms and over heavy atoms, respectively, and Δ is a pseudonormalized difference (with respect to the heavyatom structure) between moduli of structure factors.
Equation (2.2.10.2) may be compared with Karle's (1983) algebraic rule: if the sign of is plus then the value of is estimated to be zero; if its sign is minus then the expected value of is close to π. In practice Karle's rule agrees with (2.2.10.2) only if the Cochrantype term in (2.2.10.2) may be neglected. Furthermore, (2.2.10.2) shows that large reliability values do not depend on the triple product of structurefactor differences, but on the triple product of pseudonormalized differences.
A similar mathematical approach has been applied to estimate quartet invariants via isomorphous data. The result may be summarized as follows: a quartet is a phase relationship of order (Giacovazzo & Siliqi, 1996a,b; see also Kyriakidis et al., 1996), with reliability factor equal to where Q_{4} is a suitable normalizing factor.
As previously stressed, equations (2.2.10.2) and (2.2.10.3) are valid if the lack of isomorphism and the errors in the measurements are assumed to be negligible. At first sight this approach seems more appealing than the traditional twostep procedures, however it did not prove to be competitive with them. The main reason is the absence in the Hauptman and Giacovazzo approaches of a probabilistic treatment of the errors: such a treatment, on the contrary, is basic for the traditional SIR–MIR techniques [see Blow & Crick (1959) and Terwilliger & Eisenberg (1987) for two related approaches].
The problem of the errors in the probabilistic scenario defined by the joint probability distribution functions approach has recently been overcome by Giacovazzo et al. (2001). In their probabilistic calculations the following assumptions were made:where j refers to the jth derivative. is the error, which can include model as well as measurement errors.
A more realistic expression for the reliability factor G of triplet invariants is obtained by including the expression (2.2.10.4) in the probabilistic approach. Then the reliability parameter of the triplet invariants is transformed into (Giacovazzo et al., 2001)where .
Equation (2.2.10.5) suggests how the error influences the reliability of the triplet estimate: even quite a small value of may be critical if the scattering power of the heavyatom substructure is a very small percentage of the derivative scattering power.
A onestep procedure has been implemented in a computer program (Giacovazzo et al., 2002): it has been shown that the method is able to derive automatically, from the experimental data and without any user intervention, good quality (i.e. perfectly interpretable) electrondensity maps.
2.2.10.5. SIR–MIR case: the twostep procedure. Finding the heavyatom substructure by direct methods
The first trials for finding the heavyatom substructure were based on the following assumption: the modulus of the isomorphous difference,is assumed at a first approximation as an estimate of the heavyatom structure factor F_{H}. Perutz (1956) approximated F_{H}^{2} with the difference . Blow (1958) and Rossmann (1960) suggested a better approximation: . A deeper analysis was performed by Phillips (1966), Dodson & Vijayan (1971), Blessing & Smith (1999) and GrosseKunstleve & Brunger (1999). The use of direct methods requires the normalization of and application of the tangent formula (Wilson, 1978).
A sounder procedure has been suggested by Giacovazzo et al. (2004): they studied, for the SIR case, the joint probability distribution functionunder the following assumptions:

Let us suppose that the various heavyatom substructures have been determined. They may be used as additional prior information for a more accurate estimate of the values. To this purpose the distributionsmay be used under the assumption (2.2.10.6). and , for j = 1, …, n, are the structure factors of the jth derivative and of the jth heavyatom substructure, respectively, both normalized with respect to the protein. Any joint probability density (2.2.10.8) may be reliably approximated by a multidimensional Gaussian distribution (Giacovazzo & Siliqi, 2002), from which the following conditional distribution is obtained:where , the expected value of , is given byand .
is the reliability factor of the phase estimate. A robust phasing procedure has been established which, starting from the observed moduli , is able to automatically provide, without any user intervention, a highquality electrondensity map of the protein (Giacovazzo et al., 2002).
If the frequency of the radiation is close to an absorption edge of an atom, then that atom will scatter the Xrays anomalously (see Chapter 2.4 ) according to . This results in the breakdown of Friedel's law. It was soon realized that the Bijvoet difference could also be used in the determination of phases (Peerdeman & Bijvoet, 1956; Ramachandran & Raman, 1956; Okaya & Pepinsky, 1956). Since then, a great deal of work has been done both from algebraic (see Chapter 2.4 ) and from probabilistic points of view. In this section we are only interested in the second.
SAD (single anomalous dispersion) and MAD (multiple anomalous dispersion) techniques can be used. Both are characterized by one protein structure and one anomalousscatterer substructure. The experimental diffraction data differ only because of the different anomalous scattering (not because of different anomalousscatterer substructures). In the MAD case the anomalousscatterer substructure is in some way `overdetermined' by the data and, therefore, it is more convenient to use a twostep procedure: first define the positions of the anomalous scatterers, and then estimate the protein phase values. For completeness, we describe the onestep procedures in Section 2.2.10.8. These are based on the estimation of the structure invariants and on the application of the tangent formula. The twostep procedures are described in the Sections 2.2.10.9 and 2.2.10.10.
Probability distributions of diffraction intensities and of selected functions of diffraction intensities for dispersive structures have been given by various authors [Parthasarathy & Srinivasan (1964), see also Srinivasan & Parthasarathy (1976) and relevant literature cited therein]. We describe here some probabilistic formulae for estimating invariants of low order.

GivenHauptman and Giacovazzo found the following conditional distribution:
The definitions of Ω and ω are rather extensive and so the reader is referred to the published papers. We only add that Ω is always positive and that ω, the expected value of Φ, may lie anywhere between 0 and 2π. Understanding the role of the various parameters in equation (2.2.10.11) is not easy. Giacovazzo et al. (2003) found an equivalent simpler expression from which interpretable estimates of the parameters were obtained. In the same paper the limitations of the approach (versus the twostep procedures) were clarified.
2.2.10.9. SAD–MAD case: the twostep procedures. Finding the anomalousscatterer substructure by direct methods
The anomalousscatterer substructure is traditionally determined by the techniques suggested by Karle and Hendrickson (Karle, 1980b; Hendrickson, 1985; Pähler et al., 1990; Terwilliger, 1994). The introduction of selenium into proteins as selenomethionine encouraged the secondgeneration direct methods programs [Shake and Bake by Miller et al. (1994); Half bake by Sheldrick (1998); SIR2000N by Burla et al. (2001); ACORN by Foadi et al. (2000)] to locate Se atoms. Since the number of Se atoms may be quite large (up to 200), direct methods rather than Patterson techniques seem to be preferable. Shake and Bake, Half Bake and ACORN obtain the coordinates of the anomalous scatterers from a singlewavelength set of data. When more sets of diffraction data are available the solutions obtained by the other sets are used to confirm the correct solution.
A different approach has been suggested in two recent papers (Burla et al., 2002; Burla, Carrozzini et al., 2003): the estimates of the amplitudes of the structure factors of the anomalously scattering substructure are derived, via the rigorous method of the joint probability distribution functions, from the experimental diffraction moduli relative to n wavelengths. To do that, first the joint distributionis calculated, where A_{oa}, B_{oa}, E_{oa}, , , , are the real and imaginary components of E_{oa}, , , respectively, K is a symmetric square matrix of order (4n + 2), K^{−1} = {λ_{ij}} is its inverse, and T is a suitable vector with components defined in terms of the variables . E_{oa} is the normalized structure factor of the anomalous scatterer substructure calculated by neglecting anomalous scattering components. Then the conditional distributionis derived, from whichis obtained, whereThe standard deviation of the estimate is also calculated:from which
The advantage of the above approach is that the estimates can simultaneously exploit both the anomalous and the dispersive differences. The computing procedure proposed by Burla, Carrozzini et al. (2003) is the following:

The application of the above procedure to several MAD cases showed that the various wavelength combinations are not equally informative. A criterion based on the correlation among the various Δ_{ano} values was also provided (see also Schneider & Sheldrick, 2002) for predicting the most informative combinations.
Once the anomalousscatterer substructure has been found, the corresponding structure factors are known in modulus and phase. Then the conditional joint probability distributionmay be calculated (Giacovazzo & Siliqi, 2004), from which the conditional distributionmay be derived.
It has been shown that the most probable phase of , say , is the phase of the vectorand the reliability parameter of the phase estimate is nothing other than the modulus of (2.2.10.14). The first term in (2.2.10.14) is a Simlike contribution; the other terms, through the weights w, take into account the errors and the experimental differences , and .
References
Allegra, G. (1979). Derivation of threephase invariants from the Patterson function. Acta Cryst. A35, 213–220.Altomare, A., Burla, M. C., Camalli, M., Cascarano, G. L., Giacovazzo, C., Guagliardi, A., Moliterni, A. G. G., Polidori, G. & Spagna, R. (1999). SIR97: a new tool for crystal structure determination and refinement. J. Appl. Cryst. 32, 115–119.
Anzenhofer, K. & Hoppe, W. (1962). Phys. Verh. Mosbach. 13, 119.
Ardito, G., Cascarano, G., Giacovazzo, C. & Luić, M. (1985). 1Phase seminvariants and Harker sections. Z. Kristallogr. 172, 25–34.
Argos, P. & Rossmann, M. G. (1980). Molecular replacement method. In Theory and Practice of Direct Methods in Crystallography, edited by M. F. C. Ladd & R. A. Palmer, pp. 381–389. New York: Plenum.
Avrami, M. (1938). Direct determination of crystal structure from Xray data. Phys. Rev. 54, 300–303.
Baggio, R., Woolfson, M. M., Declercq, J.P. & Germain, G. (1978). On the application of phase relationships to complex structures. XVI. A random approach to structure determination. Acta Cryst. A34, 883–892.
Banerjee, K. (1933). Determination of the signs of the Fourier terms in complete crystal structure analysis. Proc. R. Soc. London Ser. A, 141, 188–193.
Bertaut, E. F. (1955a). La méthode statistique en cristallographie. I. Acta Cryst. 8, 537–543.
Bertaut, E. F. (1955b). La méthode statistique en cristallographie. II. Quelques applications. Acta Cryst. 8, 544–548.
Bertaut, E. F. (1960). Ordre logarithmique des densités de répartition. I. Acta Cryst. 13, 546–552.
Beurskens, P. T., Beurskens, G., de Gelder, R., GarciaGranda, S., Gould, R. O., Israel, R. & Smits, J. M. M. (1999). The DIRDIF99 program system. Crystallography Laboratory, University of Nijmegen, The Netherlands.
Beurskens, P. T., Gould, R. O., Bruins Slot, H. J. & Bosman, W. P. (1987). Translation functions for the positioning of a well oriented molecular fragment. Z. Kristallogr. 179, 127–159.
Beurskens, P. T., Prick, A. J., Doesburg, H. M. & Gould, R. O. (1979). Statistical properties of normalized differencestructure factors for noncentrosymmetric structures. Acta Cryst. A35, 765–772.
Bhuiya, A. K. & Stanley, E. (1963). The refinement of atomic parameters by direct calculation of the minimum residual. Acta Cryst. 16, 981–984.
Blessing, R. H. & Smith, G. D. (1999). Difference structurefactor normalization for heavyatom or anomalousscattering substructure determinations. J. Appl. Cryst. 32, 664–670.
Blow, D. M. (1958). An Xray examination of some crystal forms of pig and rabbit haemoglobin. Acta Cryst. 11, 125–126.
Blow, D. M. & Crick, F. H. C. (1959). The treatment of errors in the isomorphous replacement method. Acta Cryst. 12, 794–802.
Böhme, R. (1982). Direkte Methoden für Strukturen mit Uberstruktureffekten. Acta Cryst. A38, 318–326.
Bouman, J. (1956). A general theory of inequalities. Acta Cryst. 9, 777–780.
Bricogne, G. (1984). Maximum entropy and the foundation of direct methods. Acta Cryst. A40, 410–415.
Britten, P. L. & Collins, D. M. (1982). Information theory as a basis for the maximum determinant. Acta Cryst. A38, 129–132.
Buerger, M. J. (1959). Vector Space and its Applications in Crystal Structure Investigation. New York: John Wiley.
Burla, M. C., Caliandro, R., Camalli, M., Carrozzini, B., Cascarano, G. L., De Caro, L., Giacovazzo, C., Polidori, G. & Spagna, R. (2005). SIR2004: an improved tool for crystal structure determination and refinement. J. Appl. Cryst. 38, 381–388.
Burla, M. C., Caliandro, R., Carrozzini, B., Cascarano, G. L., De Caro, L., Giacovazzo, C. & Polidori, G. (2004). Ab initio protein phasing: the Patterson deconvolution method in SIR2002. J. Appl. Cryst. 37, 258–264.
Burla, M. C., Camalli, M., Carrozzini, B., Cascarano, G. L., Giacovazzo, C., Polidori, G. & Spagna, R. (2001). SIR2000N, a program for large and small crystal structures. J. Appl. Cryst. 34, 523–526.
Burla, M. C., Camalli, M., Carrozzini, B., Cascarano, G. L., Giacovazzo, C., Polidori, G. & Spagna, R. (2003). SIR2002: the program. J. Appl. Cryst. 36, 1103.
Burla, M. C., Carrozzini, B., Cascarano, G. L., Giacovazzo, C. & Polidori, G. (2003). SAD or MAD phasing: location of the anomalous scatterers. Acta Cryst. D59, 662–669.
Burla, M. C., Carrozzini, B., Cascarano, G. L., Giacovazzo, C., Polidori, G. & Siliqi, D. (2002). MAD phasing: probabilistic estimate of F_{oa}. Acta Cryst. D58, 928–935.
Burla, M. C., Cascarano, G., Giacovazzo, C., Nunzi, A. & Polidori, G. (1987). A weighting scheme for tangent formula development. III. The weighting scheme of the SIR program. Acta Cryst. A43, 370–374.
Busetta, B. (1976). An example of the use of quartet and triplet structure invariants when enantiomorph discrimination is difficult. Acta Cryst. A32, 139–143.
Busetta, B., Giacovazzo, C., Burla, M. C., Nunzi, A., Polidori, G. & Viterbo, D. (1980). The SIR program. I. Use of negative quartets. Acta Cryst. A36, 68–74.
Caliandro, R., Carrozzini, B., Cascarano, G. L., De Caro, L., Giacovazzo, C. & Siliqi, D. (2005a). Phasing at resolution higher than the experimental resolution. Acta Cryst. D61, 556–565.
Caliandro, R., Carrozzini, B., Cascarano, G. L., De Caro, L., Giacovazzo, C. & Siliqi, D. (2005b). Ab initio phasing at resolution higher than experimental resolution. Acta Cryst. D61, 1080–1087.
Camalli, M., Giacovazzo, C. & Spagna, R. (1985). From a partial to the complete crystal structure. II. The procedure and its applications. Acta Cryst. A41, 605–613.
Cascarano, G. & Giacovazzo, C. (1983). Onephase seminvariants of first rank. I. Algebraic considerations. Z. Kristallogr. 165, 169–174.
Cascarano, G., Giacovazzo, C., Burla, M. C., Nunzi, A. & Polidori, G. (1984). The distribution of . Acta Cryst. A40, 389–394.
Cascarano, G., Giacovazzo, C., Calabrese, G., Burla, M. C., Nunzi, A., Polidori, G. & Viterbo, D. (1984). Onephase seminvariants of first rank. II. Probabilistic considerations. Z. Kristallogr. 167, 37–47.
Cascarano, G., Giacovazzo, C., Camalli, M., Spagna, R., Burla, M. C., Nunzi, A. & Polidori, G. (1984). The method of representations of structure seminvariants. The strengthening of triplet relationships. Acta Cryst. A40, 278–283.
Cascarano, G., Giacovazzo, C. & Luić, M. (1985a). Noncrystallographic translational symmetry: effects on diffractionintensity statistics. In Structure and Statistics in Crystallography, edited by A. J. C. Wilson, pp. 67–77. Guilderland, USA: Adenine Press.
Cascarano, G., Giacovazzo, C. & Luić, M. (1985b). Direct methods and superstructures. I. Effects of the pseudotranslation on the reciprocal space. Acta Cryst. A41, 544–551.
Cascarano, G., Giacovazzo, C. & Luić, M. (1987). Direct methods and structures showing superstructure effects. II. A probabilistic theory of triplet invariants. Acta Cryst. A43, 14–22.
Cascarano, G., Giacovazzo, C. & Luić, M. (1988a). Direct methods and structures showing superstructure effects. III. A general mathematical model. Acta Cryst. A44, 176–183.
Cascarano, G., Giacovazzo, C. & Luić, M. (1988b). Direct methods and structures showing superstructure effects. IV. A new approach for phase solution. Acta Cryst. A44, 183–188.
Cascarano, G., Giacovazzo, C., Luić, M., Pifferi, A. & Spagna, R. (1987). 1Phase seminvariants and Harker sections. II. A new procedure. Z. Kristallogr. 179, 113–125.
Cascarano, G., Giacovazzo, C. & Viterbo, D. (1987). Figures of merit in direct methods: a new point of view. Acta Cryst. A43, 22–29.
Castellano, E. E., Podjarny, A. D. & Navaza, J. (1973). A multivariate joint probability distribution of phase determination. Acta Cryst. A29, 609–615.
Cochran, W. (1955). Relations between the phases of structure factors. Acta Cryst. 8, 473–478.
Cochran, W. & Douglas, A. S. (1957). The use of a highspeed digital computer for the direct determination of crystal structure. II. Proc. R. Soc. London Ser. A, 243, 281–288.
Cochran, W. & Woolfson, M. M. (1955). The theory of sign relations between structure factors. Acta Cryst. 8, 1–12.
Coulter, C. L. & Dewar, R. B. K. (1971). Tangent formula applications in protein crystallography: an evaluation. Acta Cryst. B27, 1730–1740.
Crowther, R. A. & Blow, D. M. (1967). A method of positioning a known molecule in an unknown crystal structure. Acta Cryst. 23, 544–548.
Cutfield, J. F., Dodson, E. J., Dodson, G. G., Hodgkin, D. C., Isaacs, N. W., Sakabe, K. & Sakabe, N. (1975). The high resolution structure of insulin: a comparison of results obtained from leastsquares phase refinement and difference Fourier refinement. Acta Cryst. A31, S21.
Debaerdemaeker, T., Tate, C. & Woolfson, M. M. (1985). On the application of phase relationships to complex structures. XXIV. The Sayre tangent formula. Acta Cryst. A41, 286–290.
Debaerdemaeker, T. & Woolfson, M. M. (1972). On the application of phase relationships to complex structures. IV. The coincidence method applied to general phases. Acta Cryst. A28, 477–481.
Debaerdemaeker, T. & Woolfson, M. M. (1983). On the application of phase relationships to complex structures. XXII. Techniques for random refinement. Acta Cryst. A39, 193–196.
Declercq, J.P., Germain, G., Main, P. & Woolfson, M. M. (1973). On the application of phase relationships to complex structures. V. Finding the solution. Acta Cryst. A29, 231–234.
Declercq, J.P., Germain, G. & Woolfson, M. M. (1975). On the application of phase relationships to complex structures. VIII. Extension of the magicinteger approach. Acta Cryst. A31, 367–372.
De Titta, G. T., Edmonds, J. W., Langs, D. A. & Hauptman, H. (1975). Use of the negative quartet cosine invariants as a phasing figure of merit: NQEST. Acta Cryst. A31, 472–479.
DeTitta, G. T., Weeks, C. M., Thuman, P., Miller, R. & Hauptman, H. A. (1994). Structure solution by minimalfunction phase refinement and Fourier filtering. I. Theoretical basis. Acta Cryst. A50, 203–210.
Dodson, E. & Vijayan, M. (1971). The determination and refinement of heavyatom parameters in protein heavyatom derivatives. Some model calculations using acentric reflexions. Acta Cryst. B27, 2402–2411.
Egert, E. (1983). Patterson search – an alternative to direct methods. Acta Cryst. A39, 936–940.
Egert, E. & Sheldrick, G. M. (1985). Search for a fragment of known geometry by integrated Patterson and direct methods. Acta Cryst. A41, 262–268.
Eller, G. von (1973). Génération de formules statistiques entre facteurs de structure: la méthode du polynome. Acta Cryst. A29, 63–67.
Fan, H.F. (1999). Crystallographic software: teXsan for Windows. http://www.rigaku.com/downloads/journal/Vol15.1.1998/texsan.pdf .
Fan, H. F., Hao, Q. & Woolfson, M. M. (1991). Proteins and direct methods. Z. Kristallogr. 197, 197–208.
Fan, H.F., Yao, J.X., Main, P. & Woolfson, M. M. (1983). On the application of phase relationships to complex structures. XXIII. Automatic determination of crystal structures having pseudotranslational symmetry by a modified MULTAN procedure. Acta Cryst. A39, 566–569.
Foadi, J., Woolfson, M. M., Dodson, E. J., Wilson, K. S., Jiaxing, Y. & Chaode, Z. (2000). A flexible and efficient procedure for the solution and phase refinement of protein structures. Acta Cryst. D56, 1137–1147.
Fortier, S. & Hauptman, H. (1977). Quintets in : probabilistic theory of the fivephase structure invariant in the space group . Acta Cryst. A33, 829–833.
Fortier, S., Weeks, C. M. & Hauptman, H. (1984). On integrating the techniques of direct methods and isomorphous replacement. III. The threephase invariant for the native and twoderivative case. Acta Cryst. A40, 646–651.
Freer, A. A. & Gilmore, C. J. (1980). The use of higher invariants in MULTAN. Acta Cryst. A36, 470–475.
French, S. & Wilson, K. (1978). On the treatment of negative intensity observations. Acta Cryst. A34, 517–525.
Gelder, R. de (1992). Thesis. University of Leiden, The Netherlands.
Gelder, R. de, de Graaff, R. A. G. & Schenk, H. (1990). On the construction of Karle–Hauptman matrices. Acta Cryst. A46, 688–692.
Gelder, R. de, de Graaff, R. A. G. & Schenk, H. (1993). Automatic determination of crystal structures using Karle–Hauptman matrices. Acta Cryst. A49, 287–293.
Germain, G., Main, P. & Woolfson, M. M. (1970). On the application of phase relationships to complex structures. II. Getting a good start. Acta Cryst. B26, 274–285.
Germain, G., Main, P. & Woolfson, M. M. (1971). The application of phase relationships to complex structures. III. The optimum use of phase relationships. Acta Cryst. A27, 368–376.
Giacovazzo, C. (1974). A new scheme for seminvariant tables in all space groups. Acta Cryst. A30, 390–395.
Giacovazzo, C. (1975). A probabilistic theory in of the invariant . Acta Cryst. A31, 252–259.
Giacovazzo, C. (1976). A probabilistic theory of the cosine invariant . Acta Cryst. A32, 91–99.
Giacovazzo, C. (1977a). A general approach to phase relationships: the method of representations. Acta Cryst. A33, 933–944.
Giacovazzo, C. (1977b). Strengthening of the triplet relationships. II. A new probabilistic approach in . Acta Cryst. A33, 527–531.
Giacovazzo, C. (1977c). On different probabilistic approaches to quartet theory. Acta Cryst. A33, 50–54.
Giacovazzo, C. (1977d). Quintets in and related phase relationships: a probabilistic approach. Acta Cryst. A33, 944–948.
Giacovazzo, C. (1977e). A probabilistic theory of the coincidence method. I. Centrosymmetric space groups. Acta Cryst. A33, 531–538.
Giacovazzo, C. (1977f). A probabilistic theory of the coincidence method. II. Noncentrosymmetric space groups. Acta Cryst. A33, 539–547.
Giacovazzo, C. (1978). The estimation of the onephase structure seminvariants of first rank by means of their first and second representation. Acta Cryst. A34, 562–574.
Giacovazzo, C. (1979a). A probabilistic theory of twophase seminvariants of first rank via the method of representations. III. Acta Cryst. A35, 296–305.
Giacovazzo, C. (1979b). A theoretical weighting scheme for tangentformula development and refinement and Fourier synthesis. Acta Cryst. A35, 757–764.
Giacovazzo, C. (1980a). Direct Methods in Crystallography. London: Academic Press.
Giacovazzo, C. (1980b). The method of representations of structure seminvariants. II. New theoretical and practical aspects. Acta Cryst. A36, 362–372.
Giacovazzo, C. (1980c). Triplet and quartet relations: their use in direct procedures. Acta Cryst. A36, 74–82.
Giacovazzo, C. (1983a). From a partial to the complete crystal structure. Acta Cryst. A39, 685–692.
Giacovazzo, C. (1983b). The estimation of twophase invariants in when anomalous scatterers are present. Acta Cryst. A39, 585–592.
Giacovazzo, C. (1987). One wavelength technique: estimation of centrosymmetrical twophase invariants in dispersive structures. Acta Cryst. A43, 73–75.
Giacovazzo, C. (1988). New probabilistic formulas for finding the positions of correctly oriented atomic groups. Acta Cryst. A44, 294–300.
Giacovazzo, C. (1998). Direct Phasing in Crystallography. New York: IUCr, Oxford University Press.
Giacovazzo, C., Cascarano, G. & Zheng, C.D. (1988). On integrating the techniques of direct methods and isomorphous replacement. A new probabilistic formula for triplet invariants. Acta Cryst. A44, 45–51.
Giacovazzo, C., Guagliardi, A., Ravelli, R. & Siliqi, D. (1994). Ab initio direct phasing of proteins: the limits. Z. Kristallogr. 209, 136–142.
Giacovazzo, C., Ladisa, M. & Siliqi, D. (2002). Crystal structure solution of proteins by direct methods: an automatic procedure for SIR–MIR and SIRAS–MIRAS cases. Acta Cryst. A58, 598–604.
Giacovazzo, C., Ladisa, M. & Siliqi, D. (2003). The estimation of threephase invariants when anomalous scatterers are present: the limits. Acta Cryst. A59, 569–576.
Giacovazzo, C., Moustiakimov, M., Siliqi, D. & Pifferi, A. (2004). Locating heavy atoms by integrating direct methods and SIR techniques. Acta Cryst. A60, 233–238.
Giacovazzo, C. & Siliqi, D. (1996a). On integrating direct methods and isomorphous replacement techniques. I. A distribution function for quartet invariants. Acta Cryst. A52, 133–142.
Giacovazzo, C. & Siliqi, D. (1996b). On integrating direct methods and isomorphous replacement techniques. II. The quartet invariant estimate. Acta Cryst. A52, 143–151.
Giacovazzo, C. & Siliqi, D. (2002). The method of joint probability distribution functions applied to SIR–MIR and to SIRAS–MIRAS cases. Acta Cryst. A58, 590–597.
Giacovazzo, C. & Siliqi, D. (2004). Phasing via SAD/MAD data: the method of the joint probability distribution functions. Acta Cryst. D60, 73–82.
Giacovazzo, C., Siliqi, D. & GarcíaRodríguez, L. (2001). On integrating direct methods and isomorphousreplacement techniques: triplet estimation and treatment of errors. Acta Cryst. A57, 571–575.
Gillis, J. (1948). Structure factor relations and phase determination. Acta Cryst. 1, 76–80.
Gilmore, C. J. (1984). MITHRIL. An integrated directmethods computer program. J. Appl. Cryst. 17, 42–46.
Goedkoop, J. A. (1950). Remarks on the theory of phase limiting inequalities and equalities. Acta Cryst. 3, 374–378.
Gramlich, V. (1984). The influence of rational dependence on the probability distribution of structure factors. Acta Cryst. A40, 610–616.
Grant, D. F., Howells, R. G. & Rogers, D. (1957). A method for the systematic application of sign relations. Acta Cryst. 10, 489–497.
GrosseKunstleve, R. W. & Brunger, A. T. (1999). A highly automated heavyatom search procedure for macromolecular structures. Acta Cryst. D55, 1568–1577.
Hall, S. R., du Boulay, D. J. & OlthofHazekamp, R. (1999). Xtal3.6 crystallographic software. http://xtal.sourceforge.net/ .
Harker, D. & Kasper, J. S. (1948). Phases of Fourier coefficients directly from crystal diffraction data. Acta Cryst. 1, 70–75.
Hauptman, H. (1964). The role of molecular structure in the direct determination of phase. Acta Cryst. 17, 1421–1433.
Hauptman, H. (1965). The average value of . Z. Kristallogr. 121, 1–8.
Hauptman, H. (1970). Communication at New Orleans Meeting of Am. Crystallogr. Assoc.
Hauptman, H. (1974). On the identity and estimation of those cosine invariants, , which are probably negative. Acta Cryst. A30, 472–476.
Hauptman, H. (1975). A new method in the probabilistic theory of the structure invariants. Acta Cryst. A31, 680–687.
Hauptman, H. (1982a). On integrating the techniques of direct methods and isomorphous replacement. I. The theoretical basis. Acta Cryst. A38, 289–294.
Hauptman, H. (1982b). On integrating the techniques of direct methods with anomalous dispersion. I. The theoretical basis. Acta Cryst. A38, 632–641.
Hauptman, H., Fisher, J., Hancock, H. & Norton, D. A. (1969). Phase determination for the estriol structure. Acta Cryst. B25, 811–814.
Hauptman, H. & Green, E. A. (1976). Conditional probability distributions of the fourphase structure invariant in . Acta Cryst. A32, 45–49.
Hauptman, H. & Green, E. A. (1978). Pairs in : probability distributions which lead to estimates of the twophase structure seminvariants in the vicinity of π/2. Acta Cryst. A34, 224–229.
Hauptman, H. & Karle, J. (1953). Solution of the Phase Problem. I. The Centrosymmetric Crystal. Am. Crystallogr. Assoc. Monograph No. 3. Dayton, Ohio: Polycrystal Book Service.
Hauptman, H. & Karle, J. (1956). Structure invariants and seminvariants for noncentrosymmetric space groups. Acta Cryst. 9, 45–55.
Hauptman, H. & Karle, J. (1958). Phase determination from new joint probability distributions: space group . Acta Cryst. 11, 149–157.
Hauptman, H. & Karle, J. (1959). Table 2. Equivalence classes, seminvariant vectors and seminvariant moduli for the centered centrosymmetric space groups, referred to a primitive unit cell. Acta Cryst. 12, 93–97.
Hauptman, H. A. (1991). In Crystallographic Computing 5: from Chemistry to Biology, edited by D. Moras, A. D. Podjarny & J. C. Thierry. IUCr/Oxford University Press.
Hauptman, H. A., Xu, H., Weeks, C. M. & Miller, R. (1999). Exponential ShakeandBake: theoretical basis and applications. Acta Cryst. A55, 891–900.
Heinermann, J. J. L. (1977a). The use of structural information in the phase probability of a triple product. Acta Cryst. A33, 100–106.
Heinermann, J. J. L. (1977b). Thesis. University of Utrecht.
Heinermann, J. J. L., Krabbendam, H. & Kroon, J. (1979). The joint probability distribution of the structure factors in a Karle–Hauptman matrix. Acta Cryst. A35, 101–105.
Heinermann, J. J. L., Krabbendam, H., Kroon, J. & Spek, A. L. (1978). Direct phase determination of triple products from Bijvoet inequalities. II. A probabilistic approach. Acta Cryst. A34, 447–450.
Hendrickson, W. A. (1985). Analysis of protein structure from diffraction measurement at multiple wavelengths. Trans. Am. Crystallogr. Assoc. 21, 11–21.
Hendrickson, W. A., Love, W. E. & Karle, J. (1973). Crystal structure analysis of sea lamprey hemoglobin at 2 Å resolution. J. Mol. Biol. 74, 331–361.
Hoppe, W. (1963). Phase determination and zero points in the Patterson function. Acta Cryst. 16, 1056–1057.
Hughes, E. W. (1953). The signs of products of structure factors. Acta Cryst. 6, 871.
Hull, S. E. & Irwin, M. J. (1978). On the application of phase relationships to complex structures. XIV. The additional use of statistical information in tangentformula refinement. Acta Cryst. A34, 863–870.
Hull, S. E., Viterbo, D., Woolfson, M. M. & ShaoHui, Z. (1981). On the application of phase relationships to complex structures. XIX. Magicinteger representation of a large set of phases: the MAGEX procedure. Acta Cryst. A37, 566–572.
International Tables for Crystallography (2001). Vol. F, Macromolecular Crystallography, edited by M. G. Rossmann & E. Arnold. Dordrecht: Kluwer Academic Publishers.
Jaynes, E. T. (1957). Information theory and statistical mechanics. Phys. Rev. 106, 620–630.
Karle, J. (1970a). An alternative form for B_{3.0}, a phase determining formula. Acta Cryst. B26, 1614–1617.
Karle, J. (1970b). Partial structures and use of the tangent formula and translation functions. In Crystallographic Computing, pp. 155–164. Copenhagen: Munksgaard.
Karle, J. (1972). Translation functions and direct methods. Acta Cryst. B28, 820–824.
Karle, J. (1979). Triple phase invariants: formula for centric case from fourthorder determinantal joint probability distributions. Proc. Natl Acad. Sci. USA, 76, 2089–2093.
Karle, J. (1980a). Triplet phase invariants: formula for acentric case from fourthorder determinantal joint probability distributions. Proc. Natl Acad. Sci. USA, 77, 5–9.
Karle, J. (1980b). Some developments in anomalous dispersion for the structural investigation of macromolecular systems in biology. Int. J. Quantum Chem. Quantum Biol. Symp. 7, 357–367.
Karle, J. (1983). A simple rule for finding and distinguishing triplet phase invariants with values near 0 or π with isomorphous replacement data. Acta Cryst. A39, 800–805.
Karle, J. & Hauptman, H. (1950). The phases and magnitudes of the structure factors. Acta Cryst. 3, 181–187.
Karle, J. & Hauptman, H. (1956). A theory of phase determination for the four types of noncentrosymmetric space groups , , , . Acta Cryst. 9, 635–651.
Karle, J. & Hauptman, H. (1958). Phase determination from new joint probability distributions: space group . Acta Cryst. 11, 264–269.
Karle, J. & Hauptman, H. (1961). Seminvariants for noncentrosymmetric space groups with conventional centered cells. Acta Cryst. 14, 217–223.
Karle, J. & Karle, I. L. (1966). The symbolic addition procedure for phase determination for centrosymmetric and noncentrosymmetric crystals. Acta Cryst. 21, 849–859.
Klug, A. (1958). Joint probability distributions of structure factors and the phase problem. Acta Cryst. 11, 515–543.
Koch, M. H. J. (1974). On the application of phase relationships to complex structures. IV. Automatic interpretation of electrondensity maps for organic structures. Acta Cryst. A30, 67–70.
Krabbendam, H. & Kroon, J. (1971). A relation between structure factor, triple products and a single Patterson vector, and its application to sign determination. Acta Cryst. A27, 362–367.
Kroon, J., Spek, A. L. & Krabbendam, H. (1977). Direct phase determination of triple products from Bijvoet inequalities. Acta Cryst. A33, 382–385.
Kyriakidis, C. E., Peschar, R. & Schenk, H. (1996). The estimation of fourphase structure invariants using the single difference of isomorphous structure factors. Acta Cryst. A52, 77–87.
Lajzérowicz, J. & Lajzérowicz, J. (1966). Loi de distribution des facteurs de structure pour un répartition non uniforme des atomes. Acta Cryst. 21, 8–12.
Langs, D. A. (1985). Translation functions: the elimination of structuredependent spurious maxima. Acta Cryst. A41, 305–308.
Lessinger, L. & Wondratschek, H. (1975). Seminvariants for space groups and . Acta Cryst. A31, 521.
Mackay, A. L. (1953). A statistical treatment of superlattice reflexions. Acta Cryst. 6, 214–215.
Main, P. (1976). Recent developments in the MULTAN system. The use of molecular structure. In Crystallographic Computing Techniques, edited by F. R. Ahmed, pp. 97–105. Copenhagen: Munksgaard.
Main, P. (1977). On the application of phase relationships to complex structures. XI. A theory of magic integers. Acta Cryst. A33, 750–757.
Main, P., Fiske, S. J., Germain, G., Hull, S. E., Declercq, J.P., Lessinger, L. & Woolfson, M. M. (1999). Crystallographic software: teXsan for Windows. http://www.rigaku.com/downloads/journal/Vol15.1.1998/texsan.pdf .
Main, P. & Hull, S. E. (1978). The recognition of molecular fragments in E maps and electron density maps. Acta Cryst. A34, 353–361.
Miller, R., Gallo, S. M., Khalak, H. G. & Weeks, C. M. (1994). SnB: crystal structure determination via shakeandbake. J. Appl. Cryst. 27, 613–621.
Mukherjee, A. K., Helliwell, J. R. & Main, P. (1989). The use of MULTAN to locate the positions of anomalous scatterers. Acta Cryst. A45, 715–718.
Narayan, R. & Nityananda, R. (1982). The maximum determinant method and the maximum entropy method. Acta Cryst. A38, 122–128.
Navaza, J. (1985). On the maximumentropy estimate of the electron density function. Acta Cryst. A41, 232–244.
Navaza, J., Castellano, E. E. & Tsoucaris, G. (1983). Constrained density modifications by variational techniques. Acta Cryst. A39, 622–631.
Naya, S., Nitta, I. & Oda, T. (1964). A study on the statistical method for determination of signs of structure factors. Acta Cryst. 17, 421–433.
Naya, S., Nitta, I. & Oda, T. (1965). Affinement tridimensional du sulfanilamide β. Acta Cryst. 19, 734–747.
Nordman, C. E. (1985). Introduction to Patterson search methods. In Crystallographic Computing 3. Data Collection, Structure Determination, Proteins and Databases, edited by G. M. Sheldrick, G. Kruger & R. Goddard, pp. 232–244. Oxford: Clarendon Press.
Oda, T., Naya, S. & Taguchi, I. (1961). Matrix theoretical derivation of inequalities. II. Acta Cryst. 14, 456–458.
Okaya, J. & Nitta, I. (1952). Linear structure factor inequalities and the application to the structure determination of tetragonal ethylenediamine sulphate. Acta Cryst. 5, 564–570.
Okaya, Y. & Pepinsky, R. (1956). New formulation and solution of the phase problem in Xray analysis of noncentric crystals containing anomalous scatterers. Phys. Rev. 103, 1645–1647.
Okaya, Y., Saito, Y. & Pepinsky, R. (1955). New method in Xray crystal structure determination involving the use of anomalous dispersion. Phys. Rev. 98, 1857–1858.
Oszlányi, G. & Süto, A. (2004). Ab initio structure solution by charge flipping. Acta Cryst. A60, 134–141.
Ott, H. (1927). Zur Methodik der Struckturanalyse. Z. Kristallogr. 66, 136–153.
Pähler, A., Smith, J. L. & Hendrickson, W. A. (1990). A probability representation for phase information from multiwavelength anomalous dispersion. Acta Cryst. A46, 537–540.
Parthasarathy, S. & Srinivasan, R. (1964). The probability distribution of Bijvoet differences. Acta Cryst. 17, 1400–1407.
Pavelcík, F. (1988). Pattersonoriented automatic structure determination: getting a good start. Acta Cryst. A44, 724–729.
Pavelcík, F., Kuchta, L. & Sivý, J. (1992). Pattersonoriented automatic structure determination. Utilizing Patterson peaks. Acta Cryst. A48, 791–796.
Peerdeman, A. F. & Bijvoet, J. M. (1956). The indexing of reflexions in investigations involving the use of the anomalous scattering effect. Acta Cryst. 9, 1012–1015.
Perutz, M. F. (1956). Isomorphous replacement and phase determination in noncentrosymmetric space groups. Acta Cryst. 9, 867–873.
Phillips, D. C. (1966). Advances in protein crystallography. In Advances in Structure Research by Diffraction Methods, Vol. 2, edited by R. Brill & R. Mason, pp. 75–140. New York: John Wiley.
Piro, O. E. (1983). Information theory and the phase problem in crystallography. Acta Cryst. A39, 61–68.
Podjarny, A. D., Schevitz, R. W. & Sigler, P. B. (1981). Phasing lowresolution macromolecular structure factors by matricial direct methods. Acta Cryst. A37, 662–668.
Podjarny, A. D., Yonath, A. & Traub, W. (1976). Application of multivariate distribution theory to phase extension for a crystalline protein. Acta Cryst. A32, 281–292.
Rae, A. D. (1977). The use of structure factors to find the origin of an oriented molecular fragment. Acta Cryst. A33, 423–425.
Ramachandran, G. N. & Raman, S. (1956). A new method for the structure analysis of noncentrosymmetric crystals. Curr. Sci. (India), 25, 348.
Rango, C. de (1969). Thesis. Paris.
Rango, C. de, Mauguen, Y. & Tsoucaris, G. (1975). Use of highorder probability laws in phase refinement and extension of protein structures. Acta Cryst. A31, 227–233.
Rango, C. de, Mauguen, Y., Tsoucaris, G., Dodson, E. J., Dodson, G. G. & Taylor, D. J. (1985). The extension and refinement of the 1.9 Å spacing isomorphous phases to 1.5 Å spacing in 2Zn insulin by determinantal methods. Acta Cryst. A41, 3–17.
Rango, C. de, Tsoucaris, G. & Zelwer, C. (1974). Phase determination from the Karle–Hauptman determinant. II. Connexion between inequalities and probabilities. Acta Cryst. A30, 342–353.
Richardson, J. W. & Jacobson, R. A. (1987). Patterson and Pattersons, edited by J. P. Glusker, B. K. Patterson & M. Rossi, pp. 310–317. Oxford University Press.
Rogers, D., Stanley, E. & Wilson, A. J. C. (1955). The probability distribution of intensities. VI. The influence of intensity errors on the statistical tests. Acta Cryst. 8, 383–393.
Rogers, D. & Wilson, A. J. C. (1953). The probability distribution of Xray intensities. V. A note on some hypersymmetric distributions. Acta Cryst. 6, 439–449.
Rossmann, M. G. (1960). The accurate determination of the position and shape of heavyatom replacement groups in proteins. Acta Cryst. 13, 221–226.
Rossmann, M. G. (1961). The position of anomalous scatterers in protein crystals. Acta Cryst. 14, 383–388.
Rossmann, M. G., Blow, D. M., Harding, M. M. & Coller, E. (1964). The relative positions of independent molecules within the same asymmetric unit. Acta Cryst. 17, 338–342.
Sayre, D. (1952). The squaring method: a new method for phase determination. Acta Cryst. 5, 60–65.
Sayre, D. (1953). Double Patterson function. Acta Cryst. 6, 430–431.
Sayre, D. (1972). On leastsquares refinement of the phases of crystallographic structure factors. Acta Cryst. A28, 210–212.
Sayre, D. & Toupin, R. (1975). Major increase in speed of leastsquares phase refinement. Acta Cryst. A31, S20.
Schenk, H. (1973a). Direct structure determination in and other noncentrosymmetric symmorphic space groups. Acta Cryst. A29, 480–481.
Schenk, H. (1973b). The use of phase relationships between quartets of reflexions. Acta Cryst. A29, 77–82.
Schneider, T. R. & Sheldrick, G. M. (2002). Substructure solution with SHELXD. Acta Cryst. D58, 1772–1779.
Sheldrick, G. M. (1990). Phase annealing in SHELX90: direct methods for larger structures. Acta Cryst. A46, 467–473.
Sheldrick, G. M. (1992). Crystallographic Computing 5, edited by D. Moras, A. D. Podjarny & J. C. Thierry, pp. 145–157. Oxford University Press.
Sheldrick, G. M. (1998). SHELX: applications to macromolecules. In Direct Methods for Solving Macromolecular Structures, edited by S. Fortier, pp. 401–411. Dordrecht: Kluwer Academic Publishers.
Sheldrick, G. M. (2000). The SHELX home page. http://shelx.uniac.gwdg.de/SHELX/ .
Sheldrick, G. M. & Gould, R. O. (1995). Structure solution by iterative peaklist optimization and tangent expansion in space group . Acta Cryst. B51, 423–431.
Sim, G. A. (1959). The distribution of phase angles for structures containing heavy atoms. II. A modification of the normal heavyatoms method for noncentrosymmetrical structures. Acta Cryst. 12, 813–815.
Simerska, M. (1956). Czech. J. Phys. 6, 1.
Simonov, V. I. & Weissberg, A. M. (1970). Calculation of the signs of structure amplitudes by a binary function section of interatomic vectors. Sov. Phys. Dokl. 15, 321–323.
Sint, L. & Schenk, H. (1975). Phase extension and refinement in noncentrosymmetric structures containing large molecules. Acta Cryst. A31, S22.
Srinivasan, R. & Parthasarathy, S. (1976). Some Statistical Applications in Xray Crystallography. Oxford: Pergamon Press.
Taylor, D. J., Woolfson, M. M. & Main, P. (1978). On the application of phase relationships to complex structures. XV. Magic determinants. Acta Cryst. A34, 870–883.
Terwilliger, T. C. (1994). MAD phasing: Bayesian estimates of F_{A}. Acta Cryst. D50, 11–16.
Terwilliger, T. C. & Eisenberg, D. (1987). Isomorphous replecement: effect of errors on the phase probability distribution. Acta Cryst. A43, 6–13.
Tsoucaris, G. (1970). A new method for phase determination. The maximum determinant rule. Acta Cryst. A26, 492–499.
Van der Putten, N. & Schenk, H. (1977). On the conditional probability of quintets. Acta Cryst. A33, 856–858.
Vaughan, P. A. (1958). A phasedetermining procedure related to the vectorcoincidence method. Acta Cryst. 11, 111–115.
Vermin, W. J. & de Graaff, R. A. G. (1978). The use of Karle–Hauptman determinants in smallstructure determinations. Acta Cryst. A34, 892–894.
Vicković, I. & Viterbo, D. (1979). A simple statistical treatment of unobserved reflexions. Application to two organic substances. Acta Cryst. A35, 500–501.
Weeks, C. M., DeTitta, G. T., Hauptman, H. A., Thuman, P. & Miller, R. (1994). Structure solution by minimalfunction phase refinement and Fourier filtering. II. Implementation and applications. Acta Cryst. A50, 210–220.
Weeks, C. M. & Miller, R. (1999). The design and implementation of SnB version 2.0. J. Appl. Cryst. 32, 120–124.
Weinzierl, J. E., Eisenberg, D. & Dickerson, R. E. (1969). Refinement of protein phases with the Karle–Hauptman tangent fomula. Acta Cryst. B25, 380–387.
White, P. & Woolfson, M. M. (1975). The application of phase relationships to complex structures. VII. Magic integers. Acta Cryst. A31, 53–56.
Wilkins, S. W., Varghese, J. N. & Lehmann, M. S. (1983). Statistical geometry. I. A selfconsistent approach to the crystallographic inversion problem based on information theory. Acta Cryst. A39, 47–60.
Wilson, A. J. C. (1942). Determination of absolute from relative Xray intensity data. Nature (London), 150, 151–152.
Wilson, K. S. (1978). The application of MULTAN to the analysis of isomorphous derivatives in protein crystallography. Acta Cryst. B34, 1599–1608.
Wolff, P. M. de & Bouman, J. (1954). A fundamental set of structure factor inequalities. Acta Cryst. 7, 328–333.
Woolfson, M. M. (1958). Crystal and molecular structure of p,p′dimethoxybenzophenone by the direct probability method. Acta Cryst. 11, 277–283.
Woolfson, M. M. (1977). On the application of phase relationships to complex structures. X. MAGLIN – a successor to MULTAN. Acta Cryst. A33, 219–225.
Yao, J.X. (1981). On the application of phase relationships to complex structures. XVIII. RANTAN – random MULTAN. Acta Cryst. A37, 642–664.