DRAFT DICTIONARY CBF/imgCIF Extensions Dictionary |
Draft version 1.5 for comment |
Description:
Data items in the AXIS category record the information required to describe the various goniometer, detector, source and other axes needed to specify a data collection or the axes defining the coordinate system of an image. The location of each axis is specified by two vectors: the axis itself, given by a unit vector in the direction of the axis, and an offset to the base of the unit vector. The vectors defining an axis are referenced to an appropriate coordinate system. The axis vector, itself, is a dimensionless unit vector. Where meaningful, the offset vector is given in millimetres. In coordinate systems not measured in metres, the offset is not specified and is taken as zero. The available coordinate systems are: The imgCIF standard laboratory coordinate system The direct lattice (fractional atomic coordinates) The orthogonal Cartesian coordinate system (real space) The reciprocal lattice An abstract orthogonal Cartesian coordinate frame For consistency in this discussion, we call the three coordinate system axes X, Y and Z. This is appropriate for the imgCIF standard laboratory coordinate system, and last two Cartesian coordinate systems, but for the direct lattice, X corresponds to a, Y to b and Z to c, while for the reciprocal lattice, X corresponds to a*, Y to b* and Z to c*. For purposes of visualization, all the coordinate systems are taken as right-handed, i.e., using the convention that the extended thumb of a right hand could point along the first (X) axis, the straightened pointer finger could point along the second (Y) axis and the middle finger folded inward could point along the third (Z) axis. THE IMGCIF STANDARD LABORATORY COORDINATE SYSTEM The imgCIF standard laboratory coordinate system is a right-handed orthogonal coordinate similar to the MOSFLM coordinate system, but imgCIF puts Z along the X-ray beam, rather than putting X along the X-ray beam as in MOSFLM. The vectors for the imgCIF standard laboratory coordinate system form a right-handed Cartesian coordinate system with its origin in the sample or specimen. The origin of the axis system should, if possible, be defined in terms of mechanically stable axes to be be both in the sample and in the beam. If the sample goniometer or other sample positioner has two axes the intersection of which defines a unique point at which the sample should be mounted to be bathed by the beam, that will be the origin of the axis system. If no such point is defined, then the midpoint of the line of intersection between the sample and the center of the beam will define the origin. For this definition the sample positioning system will be set at its initial reference position for the experiment. | Y (to complete right-handed system) | | | | | |________________X / principal goniometer axis / / / / /Z (to source) Axis 1 (X): The X-axis is aligned to the mechanical axis pointing from the sample or specimen along the principal axis of the goniometer or sample positioning system if the sample positioning system has an axis that intersects the origin and which form an angle of more than 22.5 degrees with the beam axis. Axis 2 (Y): The Y-axis completes an orthogonal right-handed system defined by the X-axis and the Z-axis (see below). Axis 3 (Z): The Z-axis is derived from the source axis which goes from the sample to the source. The Z-axis is the component of the source axis in the direction of the source orthogonal to the X-axis in the plane defined by the X-axis and the source axis. If the conditions for the X-axis can be met, the coordinate system will be based on the goniometer or other sample positioning system and the beam and not on the orientation of the detector, gravity etc. The vectors necessary to specify all other axes are given by sets of three components in the order (X, Y, Z). If the axis involved is a rotation axis, it is right-handed, i.e. as one views the object to be rotated from the origin (the tail) of the unit vector, the rotation is clockwise. If a translation axis is specified, the direction of the unit vector specifies the sense of positive translation. Note: This choice of coordinate system is similar to but significantly different from the choice in MOSFLM (Leslie & Powell, 2004). In MOSFLM, X is along the X-ray beam (the CBF/imgCIF Z axis) and Z is along the rotation axis. In some experimental techniques, there is no goniometer or the principal axis of the goniometer is at a small acute angle with respect to the source axis. In such cases, other reference axes are needed to define a useful coordinate system. The order of priority in defining directions in such cases is to use the detector, then gravity, then north. If the X-axis cannot be defined as above, then the direction (not the origin) of the X-axis should be parallel to the axis of the primary detector element corresponding to the most rapidly varying dimension of that detector element's data array, with its positive sense corresponding to increasing values of the index for that dimension. If the detector is such that such a direction cannot be defined (as with a point detector) or that direction forms an angle of less than 22.5 degrees with respect to the source axis, then the X-axis should be chosen so that if the Y-axis is chosen in the direction of gravity, and the Z-axis is chosen to be along the source axis, a right-handed orthogonal coordinate system is chosen. In the case of a vertical source axis, as a last resort, the X-axis should be chosen to point North. All rotations are given in degrees and all translations are given in mm. Axes may be dependent on one another. The X-axis is the only goniometer axis the direction of which is strictly connected to the hardware. All other axes are specified by the positions they would assume when the axes upon which they depend are at their zero points. When specifying detector axes, the axis is given to the beam centre. The location of the beam centre on the detector should be given in the DIFFRN_DETECTOR category in distortion-corrected millimetres from the (0,0) corner of the detector. It should be noted that many different origins arise in the definition of an experiment. In particular, as noted above, it is necessary to specify the location of the beam centre on the detector in terms of the origin of the detector, which is, of course, not coincident with the centre of the sample. The unit cell, reciprocal cell and crystallographic orthogonal Cartesian coordinate system are defined by the CELL and the matrices in the ATOM_SITES category. THE DIRECT LATTICE (FRACTIONAL COORDINATES) The direct lattice coordinate system is a system of fractional coordinates aligned to the crystal, rather than to the laboratory. This is a natural coordinate system for maps and atomic coordinates. It is the simplest coordinate system in which to apply symmetry. The axes are determined by the cell edges, and are not necessarily othogonal. This coordinate system is not uniquely defined and depends on the cell parameters in the CELL category and the settings chosen to index the crystal. Molecules in a crystal studied by X-ray diffracraction are organized into a repeating regular array of unit cells. Each unit cell is defined by three vectors, a, b and c. To quote from Drenth, "The choice of the unit cell is not unique and therefore, guidelines have been established for selecting the standard basis vectors and the origin. They are based on symmetry and metric considerations: "(1) The axial system should be right handed. (2) The basis vectors should coincide as much as possible with directions of highest symmetry." (3) The cell taken should be the smallest one that satisfies condition (2) (4) Of all the lattice vectors, none is shorter than a. (5) Of those not directed along a, none is shorter than b. (6) Of those not lying in the ab plane, none is shorter than c. (7) The three angles between the basis vectors a, b and c are either all acute (<90\%) or all obtuse (>=90\%)." These rules do not produce a unique result that is stable under the assumption of experimental errors, and the the resulting cell may not be primitive. In this coordinate system, the vector (.5, .5, .5) is in the middle of the given unit cell. Grid coordinates are an important variation on fractional coordinates used when working with maps. In imgCIF, the conversion from fractional to grid coordinates is implicit in the array indexing specified by _array_structure_list.dimension. Note that this implicit grid-coordinate scheme is 1-based, not zero-based, i.e. the origin of the cell for axes along the cell edges with no specified _array_structure_list_axis.displacement will have grid coordinates of (1,1,1), i.e. array indices of (1,1,1). THE ORTHOGONAL CARTESIAN COORDINATE SYSTEM (REAL SPACE) The orthogonal Cartesian coordinate system is a transformation of the direct lattice to the actual physical coordinates of atoms in space. It is similar to the laboratory coordinate system, but is anchored to and moves with the crystal, rather than being schored to the laboratory. The transformation from fractional to orthogonal cartesian coordinates is given by the _atom_sites.Cartn_transf_matrix[i][j] and _atom_sites.Cartn_transf_vector[i] tags. A common choice for the matrix of the transformation is given in the 1992 PDB format document | a b cos(\g) c cos(\b) | | 0 b sin(\g) c (cos(\a) - cos(\b)cos(\g))/sin(\g) | | 0 0 V/(a b sin(\g)) | This is a convenient coordinate system in which to do fitting of models to maps and in which to understand the chemistry of a molecule. THE RECIPROCAL LATTICE The reciprocal lattice coordinate system is used for diffraction intensitities. It is based on the reciprocal cell, the dual of the cell, in which reciprocal cell edges are derived from direct cell faces: a* = bc sin(\a)/V b* = ac sin(\b)/V c* = ab sin(\g)/V cos(\a*) = (cos(\b) cos(\g) - cos(\a))/(sin(\b) sin(\g)) cos(\b*) = (cos(\a) cos(\g) - cos(\b))/(sin(\a) sin(\g)) cos(\g*) = (cos(\a) cos(\b) - cos(\g))/(sin(\a) sin(\b)) V = abc SQRT(1 - cos(\a)^2^ - cos(\b)^2^ - cos(\g)^2^ + 2 cos(\a) cos(\b) cos(\g) ) In this form the dimensions of the reciprocal lattice are in reciprocal \%Angstroms (\%A^-1^). A dimensionless form can be obtained by multiplying by the wavelength. Reflections are commonly indexed against this coordinate system as (h, k, l) triples. References: Drenth, J., "Introduction to basic crystallography." chapter 2.1 in Rossmann, M. G. and Arnold, E. "Crystallography of biological macromolecules", Volume F of the IUCr's "International tables for crystallography", Kluwer, Dordrecht 2001, pp 44 -- 63 Leslie, A. G. W. and Powell, H. (2004). MOSFLM v6.11. MRC Laboratory of Molecular Biology, Hills Road, Cambridge, England. http://www.CCP4.ac.uk/dist/X-windows/Mosflm/. Stout, G. H. and Jensen, L. H., "X-ray structure determination", 2nd ed., Wiley, New York, 1989, 453 pp. __, "PROTEIN DATA BANK ATOMIC COORDINATE AND BIBLIOGRAPHIC ENTRY FORMAT DESCRIPTION," Brookhaven National Laboratory, February 1992.Examples:
Example 1 -
This example shows the axis specification of the axes of a kappa-
geometry goniometer [see Stout, G. H. & Jensen, L. H. (1989). X-ray
structure determination. A practical
guide, 2nd ed. p. 134. New York: Wiley Interscience].
There are three axes specified, and no offsets. The outermost axis,
omega, is pointed along the X axis. The next innermost axis, kappa,
is at a 50 degree angle to the X axis, pointed away from the source.
The innermost axis, phi, aligns with the X axis when omega and
phi are at their zero points. If T-omega, T-kappa and T-phi
are the transformation matrices derived from the axis settings,
the complete transformation would be:
X' = (T-omega) (T-kappa) (T-phi) X
loop_ _axis.id _axis.type _axis.equipment _axis.depends_on _axis.vector[1] _axis.vector[2] _axis.vector[3] omega rotation goniometer . 1 0 0 kappa rotation goniometer omega -.64279 0 -.76604 phi rotation goniometer kappa 1 0 0 |
Example 2 -
This example shows the axis specification of the axes of a
detector, source and gravity. The order has been changed as a
reminder that the ordering of presentation of tokens is not
significant. The centre of rotation of the detector has been taken
to be 68 millimetres in the direction away from the source.
loop_ _axis.id _axis.type _axis.equipment _axis.depends_on _axis.vector[1] _axis.vector[2] _axis.vector[3] _axis.offset[1] _axis.offset[2] _axis.offset[3] source . source . 0 0 1 . . . gravity . gravity . 0 -1 0 . . . tranz translation detector rotz 0 0 1 0 0 -68 twotheta rotation detector . 1 0 0 . . . roty rotation detector twotheta 0 1 0 0 0 -68 rotz rotation detector roty 0 0 1 0 0 -68 |
Example 3 -
This example show the axis specification of the axes for a map,
using fractional coordinates. Each cell edge has been divided
into a grid of 50 divisions in the ARRAY_STRUCTURE_LIST_AXIS
category. The map is using only the first octant of the grid
in the ARRAY_STRUCTURE_LIST category.
The fastest changing axis is the gris along A, then along B,
and the slowest is along C.
The map sampling is being done in the middle of each grid
division
loop_ _axis.id _axis.system _axis.vector[1] _axis.vector[2] _axis.vector[3] CELL_A_AXIS fractional 1 0 0 CELL_B_AXIS fractional 0 1 0 CELL_C_AXIS fractional 0 0 1 loop_ _array_structure_list.array_id _array_structure_list.index _array_structure_list.dimension _array_structure_list.precedence _array_structure_list.direction _array_structure_list.axis_id MAP 1 25 1 increasing CELL_A_AXIS MAP 1 25 2 increasing CELL_B_AXIS MAP 1 25 3 increasing CELL_C_AXIS loop_ _array_structure_list_axis.axis_id _array_structure_list_axis.fract_displacement _array_structure_list_axis.fract_displacement_increment CELL_A_AXIS 0.01 0.02 CELL_B_AXIS 0.01 0.02 CELL_C_AXIS 0.01 0.02 |
Example 4 -
This example show the axis specification of the axes for a map,
this time as orthogonal \%Angstroms, using the same coordinate system
as for the atomic coordinates. The map is sampling every 1.5
\%Angstroms (1.5e-7 millimeters) in a map segment 37.5 \%Angstroms on
a side.
loop_ _axis.id _axis.system _axis.vector[1] _axis.vector[2] _axis.vector[3] X orthogonal 1 0 0 Y orthogonal 0 1 0 Z orthogonal 0 0 1 loop_ _array_structure_list.array_id _array_structure_list.index _array_structure_list.dimension _array_structure_list.precedence _array_structure_list.direction _array_structure_list.axis_id MAP 1 25 1 increasing X MAP 2 25 2 increasing Y MAP 3 25 3 increasing Z loop_ _array_structure_list_axis.axis_id _array_structure_list_axis.displacement _array_structure_list_axis.displacement_increment X 7.5e-8 1.5e-7 Y 7.5e-8 1.5e-7 Z 7.5e-8 1.5e-7 |
Mandatory category: no
Copyright © 2005 International Union of Crystallography