IDXREF.LP: Difference between revisions
| (2 intermediate revisions by the same user not shown) | |||
| Line 87: | Line 87: | ||
From the difference vectors, the "reduced cell" (essentially a P1 cell, with a<b<c) has been established, together with its axes and angles. Furthermore the relation of the reciprocal cell axes (found in the beginning) with respect to the reduced cell is given. If the user supplies UNIT_CELL_CONSTANTS (and SPACE_GROUP_NUMBER >0 ) then these are converted to a reduced cell and given here. | From the difference vectors, the "reduced cell" (essentially a P1 cell, with a<b<c) has been established, together with its axes and angles. Furthermore the relation of the reciprocal cell axes (found in the beginning) with respect to the reduced cell is given. If the user supplies UNIT_CELL_CONSTANTS (and SPACE_GROUP_NUMBER >0 ) then these are converted to a reduced cell and given here. | ||
'''If the difference vectors are not (close to) integers, something is wrong''' - see [[Problems]]. | |||
== Results from local indexing == | == Results from local indexing == | ||
| Line 282: | Line 283: | ||
0.54 1.00 | 0.54 1.00 | ||
This is just a little help to tell the user how big the OSCILLATION_RANGE can be without producing overlap. The maximum oscillation range is less than that given by the table, since the crystal mosaicity has to be subtracted from the table value. Please see [[ | This is just a little help to tell the user how big the OSCILLATION_RANGE can be without producing overlap. The maximum oscillation range is less than that given by the table, since the crystal mosaicity has to be subtracted from the table value. Please see [[Choice of OSCILLATION RANGE]]. | ||
cpu time used 2.8 sec | cpu time used 2.8 sec | ||
| Line 304: | Line 305: | ||
JOB= DEFPIX INTEGRATE CORRECT | JOB= DEFPIX INTEGRATE CORRECT | ||
This is a feature (not a bug) to make the user aware of a possible problem. | This is a feature (not a bug) to make the user aware of a possible problem. | ||
== See also == | |||
[[Indexing]] | |||
Latest revision as of 19:20, 9 September 2016
IDXREF.LP is the logfile of the IDXREF "job". It reports on aspects of the analysis of difference vectors, the lattice(s) encountered, their interpretation in terms of 44 Bravais lattices, and the refinement of all geometric parameters of the experiment based on the strong reflections identified in the SPOT_RANGE by COLSPOT.
In the following, extracts from a typical IDXREF.LP are shown, together with comments, in the order as they occur.
Listing of parameter values used for indexing
INPUT PARAMETER VALUES ---------------------- NAME_TEMPLATE_OF_DATA_FRAMES=../kb/G3-16mer-A2-1_1_???.img TIFF STARTING_ANGLE= 0.0000 STARTING_FRAME= 1 NX= 3072 NY= 3072 QX= 0.073242 QY= 0.073242 ROTATION_AXIS= 1.000000 0.000000 0.000000 OSCILLATION_RANGE= 1.0000 X-RAY_WAVELENGTH= 1.000000 INCIDENT_BEAM_DIRECTION= 0.000000 0.000000 1.000000 DIRECTION_OF_DETECTOR_X-AXIS= 1.00000 0.00000 0.00000 DIRECTION_OF_DETECTOR_Y-AXIS= 0.00000 1.00000 0.00000 DETECTOR_DISTANCE= 349.994 ORGX= 1532.00 ORGY= 1566.00 INDEX_MAGNITUDE= 8 INDEX_ERROR= 0.050 INDEX_QUALITY= 0.80 SEPMIN= 6.00 CLUSTER_RADIUS= 3 MAXIMUM_ERROR_OF_SPOT_POSITION= 3.0 MAXIMUM_ERROR_OF_SPINDLE_POSITION= 2.0 SPACE GROUP AND CELL PARAMETERS ARE UNKNOWN
The values shown are taken from XDS.INP, or if not given there, the defaults are output. "SPACE GROUP AND CELL PARAMETERS ARE UNKNOWN" means that XDS.INP had SPACE_GROUP_NUMBER=0
Which data were used for indexing?
AUTOINDEXING IS BASED ON 4043 SPOTS LOCATED IN THE TRUSTED_REGION= 0.00 1.20 OF THE DETECTOR OF THE FOLLOWING DATA IMAGES: SPOT_RANGE= 1 90
The "SPOTS" mentioned above are read from SPOT.XDS, and are collected by COLSPOT. Nota bene: SPOT.XDS is a text file and could be written by a custom program, or a SPOT.XDS written by COLSPOT could be modified afterwards - see Indexing!
Determination of the reciprocal lattice basis
NUMBER OF DIFFERENCE VECTOR CLUSTERS USED 198 MAXIMUM RADIUS OF DIFFERENCE VECTOR CLUSTERS (pixels) 3 MINIMUM DISTANCE BETWEEN DIFFRACTION SPOTS (pixel) 6.0 MINIMUM ALLOWED DISTANCE BETWEEN REC. LATTICE POINTS 0.1256E-02 OBSERVED BASIS CELL VOLUME 0.1080E+07 DIMENSION OF SPACE SPANNED BY DIFFERENCE VECTOR CLUSTERS 3 # COORDINATES OF REC. BASIS VECTOR LENGTH 1/LENGTH 1 0.0040197-0.0034658 0.0044763 0.0069432 144.03 2 0.0060960 0.0063989-0.0005531 0.0088551 112.93 3 -0.0064850 0.0072590 0.0114902 0.0150590 66.41
Above, the values of CLUSTER_RADIUS (3) and SEPMIN (6.0) are repeated. The reciprocal cell axis lengths (given as 1/LENGTH) are derived from the difference vectors between reciprocal lattice points. If the user supplies UNIT_CELL_CONSTANTS (and SPACE_GROUP_NUMBER is >0), then the supplied unit cell axes are matched here against the observed difference vectors.
The output continues with showing the difference vectors in h,k,l units - ideally these should be close to integral numbers (1, 2 3, ...), which they are in this example:
CLUSTER COORDINATES AND INDICES WITH RESPECT TO REC. LATTICE BASIS VECTORS
# COORDINATES OF VECTOR CLUSTER FREQUENCY CLUSTER INDICES
1 -0.0040208 0.0034733-0.0044142 991. -0.99 0.00 0.00
2 -0.0101816-0.0029195-0.0037781 945. -0.99 -1.01 0.01
3 -0.0122581-0.0127817 0.0012276 937. 0.00 -2.00 0.01
4 -0.0061055-0.0064076 0.0006063 921. 0.00 -1.00 0.00
5 -0.0020733-0.0098616 0.0050198 904. 1.00 -1.00 0.00
6 -0.0080465 0.0069232-0.0089221 848. -2.00 0.00 0.00
7 -0.0141979 0.0005603-0.0082059 844. -1.99 -1.00 0.01
8 0.0064218-0.0072694-0.0114905 801. 0.00 -0.01 -1.00
9 -0.0019301 0.0133409-0.0094462 794. -1.99 1.00 0.00
10 -0.0206410 0.0078000 0.0032924 757. -1.98 -1.00 1.01
11 -0.0203128-0.0058325-0.0075707 744. -1.98 -2.01 0.02
12 -0.0023531 0.0037679 0.0159009 740. 1.01 0.01 0.99
...
58 0.0000713 0.0231797-0.0145601 486. -3.01 1.99 0.00
59 -0.0119425-0.0263910-0.0096567 482. 0.00 -3.01 -0.99
60 -0.0117272-0.0032032-0.0242257 480. -3.00 -1.00 -0.99
PARAMETERS OF THE REDUCED CELL (ANGSTROEM & DEGREES)
66.41 112.93 144.03 89.86 89.89 89.76
# COORDINATES OF REC. BASIS VECTOR REDUCED CELL INDICES
1 0.0040197-0.0034658 0.0044763 0.00 0.00 1.00
2 0.0060960 0.0063989-0.0005531 0.00 1.00 0.00
3 -0.0064850 0.0072590 0.0114902 -1.00 0.00 0.00
From the difference vectors, the "reduced cell" (essentially a P1 cell, with a<b<c) has been established, together with its axes and angles. Furthermore the relation of the reciprocal cell axes (found in the beginning) with respect to the reduced cell is given. If the user supplies UNIT_CELL_CONSTANTS (and SPACE_GROUP_NUMBER >0 ) then these are converted to a reduced cell and given here.
If the difference vectors are not (close to) integers, something is wrong - see Problems.
Results from local indexing
RESULTS FROM LOCAL INDEXING OF 3000 OBSERVED SPOTS *****
MAXIMUM MAGNITUDE OF INDEX DIFFERENCES ALLOWED 8
MAXIMUM ALLOWED DEVIATION FROM INTEGERAL INDICES 0.050
MIMINUM QUALITY OF INDICES FOR EACH SPOT IN A SUBTREE 0.80
QUALITY OF INDICES REQUIRED TO INCLUDE SECOND SUBTREE 0.00
NUMBER OF SUBTREES 118
SUBTREE POPULATION
1 2873
2 6
3 3
4 2
5 2
6 2
...
The "subtrees" each refer to their own lattice. The list above indicates that 2873 out of the strongest 3000 reflections can be indexed with a single lattice. 6 reflections correspond to the second-best lattice. If the diffraction pattern arises from split crystals, or there are two (or more) non-equivalent lattices because e.g. ORGX ORGY (in XDS.INP) denote a position right in the middle between two reflections, then several lattices are listed here that have a substantial number of reflections. In such a case IDXREF will choose the lattice with most reflections, but the user should be aware that other lattices exist!
Finding the origin of the reciprocal lattice (=direct beam position, if normal geometry)
***** SELECTION OF THE INDEX ORIGIN OF THE REFLECTIONS *****
The origin of the reflection indices determined so far is
0,0,0 by default which is usually correct. In certain critical
cases it may happen that this automatic choice is wrong which
leads to misindexing of the reflections by a constant offset.
You may replace the default by specifying INDEX_ORIGIN= h k l
in the input file "XDS.INP" and rerun the IDXREF step.
Below you find a list of possible alternatives together with a
measure of their likelihood.
QUALITY small values mean a high likelihood for this offset
DELTA distance between given 1532.00 1566.00 and computed
direct beam position (pixels) on the detector
XD, YD computed direct beam position (pixels) on detector
X,Y,Z computed coordinates of the direct beam wave vector
DH,DK,DL mean absolute difference between observed and
fitted indices
INDEX_ QUALITY DELTA XD YD X Y Z DH DK DL
ORIGIN
0 0 0 1.0 0.5 1531.7 1566.3 -0.0001 0.0001 1.0000 0.06 0.06 0.10
0 0 -1 4.4 30.4 1550.9 1542.1 0.0039 -0.0050 1.0000 0.13 0.11 0.16
0 0 1 4.4 31.4 1512.4 1590.6 -0.0041 0.0051 1.0000 0.13 0.11 0.15
SELECTED: INDEX_ORIGIN= 0 0 0
Errors in the values of ORGX, ORGY (as supplied in XDS.INP) are the most common single source of indexing failure. XDS tries several possible origins (3 in the example above) around the supplied values and gives an estimate of indexing quality ("QUALITY"; 1.0 is best) for each of them. DH,HK,DL should ideally be 0; they correspond to the deviation of H,K,L from being integer. It is advisable to always use INDEX_ORIGIN= 0 0 0 (the default); if the QUALITY and DH,DK,DL indicators say that the supplied ORGX, ORGY are wrong then the latter should be fixed.
It is important to realize that ORGX and ORGY are the coordinates of the point of the detector which is closest to the crystal; this is not the same as the direct beam coordinates! However, in practice, at synchrotron beamlines the detector is perpendicular to the beam, in which case taking the direct beam position as ORGY ORGY is accurate enough; IDXREF (and INTEGRATE, CORRECT) refines the beam direction (and other geometric parameters) anyway. This is also why there is no need to specify highly accurate values for BEAM_DIRECTION, ROTATION_AXIS, ORGX, ORGY.
The word "SELECTED:" may be a bit misleading - XDS does not select the INDEX_ORIGIN, it's the user who does this.
First refinement of geometry parameters
***** REFINED SOLUTION BASED ON INDEXED REFLECTIONS IN SUBTREE # 1 ***** REFINED VALUES OF DIFFRACTION PARAMETERS DERIVED FROM 2873 INDEXED SPOTS REFINED PARAMETERS: DISTANCE BEAM AXIS CELL ORIENTATION STANDARD DEVIATION OF SPOT POSITION (PIXELS) 5.16 STANDARD DEVIATION OF SPINDLE POSITION (DEGREES) 4.60 CRYSTAL MOSAICITY (DEGREES) 0.200 DIRECT BEAM COORDINATES (REC. ANGSTROEM) -0.002368 0.008316 0.999963 DETECTOR COORDINATES (PIXELS) OF DIRECT BEAM 1531.82 1565.97 DETECTOR ORIGIN (PIXELS) AT 1541.74 1531.14 CRYSTAL TO DETECTOR DISTANCE (mm) 306.77 LAB COORDINATES OF DETECTOR X-AXIS 1.000000 0.000000 0.000000 LAB COORDINATES OF DETECTOR Y-AXIS 0.000000 1.000000 0.000000 LAB COORDINATES OF ROTATION AXIS 0.999916 -0.012791 -0.002291 COORDINATES OF UNIT CELL A-AXIS 24.214 -28.552 -44.997 COORDINATES OF UNIT CELL B-AXIS 68.543 70.993 -7.718 COORDINATES OF UNIT CELL C-AXIS 74.857 -62.751 80.778 REC. CELL PARAMETERS 0.017084 0.010103 0.007889 90.239 89.766 89.803 UNIT CELL PARAMETERS 58.534 98.983 126.753 89.760 90.235 90.198 SPACE GROUP NUMBER 1
Based on the parameters known at this point, and all reflections assigned to the strongest lattice, a first refinement is done.
Second refinement, after rejecting reflections that do not fit well
***** INDEXING OF OBSERVED SPOTS IN SPACE GROUP # 1 *****
1103 OUT OF 4043 SPOTS INDEXED.
0 REJECTED REFLECTIONS (REASON: OVERLAP)
2940 REJECTED REFLECTIONS (REASON: TOO FAR FROM IDEAL POSITION)
EXPECTED ERROR IN SPINDLE POSITION 0.809 DEGREES
EXPECTED ERROR IN DETECTOR POSITION 1.90 PIXELS
***** DIFFRACTION PARAMETERS USED AT START OF INTEGRATION *****
REFINED VALUES OF DIFFRACTION PARAMETERS DERIVED FROM 1103 INDEXED SPOTS
REFINED PARAMETERS: DISTANCE BEAM AXIS CELL ORIENTATION
STANDARD DEVIATION OF SPOT POSITION (PIXELS) 1.86
STANDARD DEVIATION OF SPINDLE POSITION (DEGREES) 0.79
CRYSTAL MOSAICITY (DEGREES) 0.200
DIRECT BEAM COORDINATES (REC. ANGSTROEM) 0.002112 0.014343 0.999895
DETECTOR COORDINATES (PIXELS) OF DIRECT BEAM 1531.92 1565.77
DETECTOR ORIGIN (PIXELS) AT 1522.67 1502.92
CRYSTAL TO DETECTOR DISTANCE (mm) 320.89
LAB COORDINATES OF DETECTOR X-AXIS 1.000000 0.000000 0.000000
LAB COORDINATES OF DETECTOR Y-AXIS 0.000000 1.000000 0.000000
LAB COORDINATES OF ROTATION AXIS 0.999918 -0.009393 -0.008744
COORDINATES OF UNIT CELL A-AXIS 25.289 -29.852 -46.602
COORDINATES OF UNIT CELL B-AXIS 71.309 74.418 -8.739
COORDINATES OF UNIT CELL C-AXIS 78.373 -65.138 84.444
REC. CELL PARAMETERS 0.016435 0.009668 0.007556 90.014 89.938 89.901
UNIT CELL PARAMETERS 60.847 103.437 132.348 89.986 90.062 90.099
SPACE GROUP NUMBER 1
Based on the results from the first refinement, all reflections found by COLSPOT are indexed. In this case, a bit more than 1/4 of these are indexed with low error. This leads to the message "!!! ERROR !!! INSUFFICIENT PERCENTAGE (< 50%) OF INDEXED REFLECTIONS" at the bottom of IDXREF.LP , since that fraction is less than MINIMUM_FRACTION_OF_INDEXED_SPOTS (default 0.50).
Determination of Bravais lattices consistent with the observed spot positions
*********** DETERMINATION OF LATTICE CHARACTER AND BRAVAIS LATTICE *********** The CHARACTER OF A LATTICE is defined by the metrical parameters of its reduced cell as described in the INTERNATIONAL TABLES FOR CRYSTALLOGRAPHY Volume A, p. 746 (KLUWER ACADEMIC PUBLISHERS, DORDRECHT/BOSTON/LONDON, 1989). Note that more than one lattice character may have the same BRAVAIS LATTICE. A lattice character is marked "*" to indicate a lattice consistent with the observed locations of the diffraction spots. These marked lattices must have low values for the QUALITY OF FIT and their implicated UNIT CELL CONSTANTS should not violate the ideal values by more than MAXIMUM_ALLOWED_CELL_AXIS_RELATIVE_ERROR= 0.03 MAXIMUM_ALLOWED_CELL_ANGLE_ERROR= 3.0 (Degrees) LATTICE- BRAVAIS- QUALITY UNIT CELL CONSTANTS (ANGSTROEM & DEGREES) CHARACTER LATTICE OF FIT a b c alpha beta gamma * 31 aP 0.0 60.8 103.4 132.3 90.0 89.9 89.9 * 44 aP 0.4 60.8 103.4 132.3 90.0 90.1 90.1 * 34 mP 2.1 60.8 132.3 103.4 90.0 90.1 90.1 * 33 mP 2.4 60.8 103.4 132.3 90.0 90.1 90.1 * 35 mP 3.1 103.4 60.8 132.3 90.1 90.0 90.1 * 32 oP 3.5 60.8 103.4 132.3 90.0 90.1 90.1 29 mC 248.7 60.8 215.5 132.3 90.0 90.1 73.7 28 mC 249.1 60.8 271.5 103.4 90.0 90.1 77.1 39 mC 250.2 215.5 60.8 132.3 90.1 90.0 73.7 ... ...
The above list is sorted by the "Quality of fit" - good values are below 10. Triclinic (Bravais lattice "aP") is always the best since it has no restrictions and can thus most easily fit the reduced cell. The unit cell constants are not cleaned to obey the restrictions, e.g. orthorhombic does not necessarily have alpha=beta=gamma=90°. (Please note that, when specifying unit cell constants in XDS.INP, all restrictions have to be met.)
For protein crystals the possible space group numbers corresponding to
each Bravais-type are given below for your convenience. Note, that
reflection integration is based only on orientation and metric of the
lattice. It does not require knowledge of the correct space group!
Thus, if no such information is provided by the user in XDS.INP,
reflections are integrated assuming a triclinic reduced cell lattice;
the space group is assigned automatically or by the user in the last
step (CORRECT) when integrated intensities are available.
****** LATTICE SYMMETRY IMPLICATED BY SPACE GROUP SYMMETRY ******
BRAVAIS- POSSIBLE SPACE-GROUPS FOR PROTEIN CRYSTALS
TYPE [SPACE GROUP NUMBER,SYMBOL]
aP [1,P1]
mP [3,P2] [4,P2(1)]
mC,mI [5,C2]
oP [16,P222] [17,P222(1)] [18,P2(1)2(1)2] [19,P2(1)2(1)2(1)]
oC [21,C222] [20,C222(1)]
oF [22,F222]
oI [23,I222] [24,I2(1)2(1)2(1)]
tP [75,P4] [76,P4(1)] [77,P4(2)] [78,P4(3)] [89,P422] [90,P42(1)2]
[91,P4(1)22] [92,P4(1)2(1)2] [93,P4(2)22] [94,P4(2)2(1)2]
[95,P4(3)22] [96,P4(3)2(1)2]
tI [79,I4] [80,I4(1)] [97,I422] [98,I4(1)22]
hP [143,P3] [144,P3(1)] [145,P3(2)] [149,P312] [150,P321] [151,P3(1)12]
[152,P3(1)21] [153,P3(2)12] [154,P3(2)21] [168,P6] [169,P6(1)]
[170,P6(5)] [171,P6(2)] [172,P6(4)] [173,P6(3)] [177,P622]
[178,P6(1)22] [179,P6(5)22] [180,P6(2)22] [181,P6(4)22] [182,P6(3)22]
hR [146,R3] [155,R32]
cP [195,P23] [198,P2(1)3] [207,P432] [208,P4(2)32] [212,P4(3)32]
[213,P4(1)32]
cF [196,F23] [209,F432] [210,F4(1)32]
cI [197,I23] [199,I2(1)3] [211,I432] [214,I4(1)32]
This is just the mapping from Bravais lattice to possible spacegroups.
Maximum oscillation range to prevent angular overlap at high resolution limit
assuming zero (!) mosaicity.
Maximum oscillation range High resolution limit
(degrees) (Angstrom)
2.15 4.00
1.61 3.00
1.07 2.00
0.54 1.00
This is just a little help to tell the user how big the OSCILLATION_RANGE can be without producing overlap. The maximum oscillation range is less than that given by the table, since the crystal mosaicity has to be subtracted from the table value. Please see Choice of OSCILLATION RANGE.
cpu time used 2.8 sec elapsed wall-clock time 1.7 sec !!! ERROR !!! INSUFFICIENT PERCENTAGE (< 50%) OF INDEXED REFLECTIONS AUTOMATIC DATA PROCESSING STOPPED. AS THE CRITERIA FOR A GOOD SOLUTION ARE RATHER STRICT, YOU MAY CHOOSE TO CONTINUE DATA PROCESSING AFTER CHANGING THE "JOB="-CARD IN "XDS.INP" TO "JOB= DEFPIX INTEGRATE CORRECT". IF THE BEST SOLUTION IS REALLY NONSENSE YOU SHOULD FIRST HAVE A LOOK AT THE ASCII-FILE "SPOT.XDS". THIS FILE CONTAINS THE INITIAL SPOT LIST SORTED IN DECREASING SPOT INTENSITY. SPOTS NEAR THE END OF THE FILE MAY BE ARTEFACTS AND SHOULD BE ERASED. ALTERNATIVELY YOU MAY TRY DIFFERENT VALUES FOR "INDEX_ORIGIN" AS SUGGESTED IN THE ABOVE LISTING. IF THE CRYSTAL HAS SLIPPED AT THE BEGINNING OF DATA COLLECTION YOU MAY CHOOSE TO SKIP SOME OF THE FIRST FRAMES BY CHANGING THE "DATA_RANGE=" IN FILE "XDS.INP" AND START ALL OVER AGAIN.
End of IDXREF.LP. In this case XDS would not automatically continue with the DEFPIX step. Rather, the user has to explicitly state that s/he wants to do this, by changing the JOB line in XDS.INP to
JOB= DEFPIX INTEGRATE CORRECT
This is a feature (not a bug) to make the user aware of a possible problem.