# IDXREF.LP: Difference between revisions

No edit summary |
|||

Line 200: | Line 200: | ||

SPACE GROUP NUMBER 1 | 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 . | 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 [http://xds.mpimf-heidelberg.mpg.de/html_doc/xds_parameters.html#MINIMUM_FRACTION_OF_INDEXED_SPOTS= MINIMUM_FRACTION_OF_INDEXED_SPOTS] (default 0.50). |

## Revision as of 22:58, 19 June 2012

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.

## 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).