NEWLAUE - LAUE SIMULATIONS AND INDEXING ======================================= INTRODUCTION The program NEWLAUE may be used to generate and display a Laue pattern given basic information such as cell parameters, crystal to film distance, wavelength range, minimum spacing, crystal type and setting. It is based on M. Elder's LGEN program written originally for a PERQ computer at the Daresbury Laboratory . A large amount of development was carried out by E.M. Pantos and T. Jhahjadi. The current version has been further developed by I.J. Clifton and J.W. Campbell. When the program starts the option menu is displayed on the left of the screen. Options may be selected by typing the first four characters of the displayed name. A prompt is output at the bottom of the screen requesting the selection of an option with the current default option being displayed in square brackets. Options are of two types, those which set parameter values and those which are program commands. When a parameter value is to be input, this may be done, in general, in a number of ways. If the default option is the required option then a carriage return, , may be input in which case the program will prompt for the required value. Alternatively the value itself may be given directly (or the parameter name followed by the value). If the default option is not the required option then the required parameter name is given. If this is not followed by a value then the program will prompt for the required value. If it is followed by a value then this value will become the new value for the parameter. CTRL/Z reprints the option menu. This may be necessary on some terminals if the alpha screen is erased after a plot has been displayed. CTRL/Z may also be used to cancel any other prompts throughout the program (but not cancel cursor options). Calculation of the pattern commences when the START option is selected. START has to be reused every time a pattern parameter is altered. List of sections: Parameters Commands Files PARAMETERS Introduction This section describes the parameters used by the NEWLAUE program. List of subsections in this section: CELL WMIN WMAX DMIN DISTANCE RADIUS ORIENT PX PY PZ SYSTEM SPINDLE ABSENCES CELL The CELL parameters, A, B, C, ALPHA, BETA and GAMMA are used to determine the shape of the unit cell. Not all these parameters are needed in more symmetric crystal systems ( see SYSTEM ). Each parameter is individually selected by name and the program prompts for input value (cell length in Angstroms, angles in degrees). For some crystal systems only some of the parameters need be defined, e.g., for Cubic systems only A is sufficient, B=A, C=A , ALPHA=90, BETA=ALPHA, GAMMA=ALPHA are implied. If the user insists in supplying values to surplus parameters, correct or incorrect ones, the program will ignore them, in fact it resets them when the calculation starts. In such cases the list in the Menu and the hard copy will leave the value field blank. WMIN The minimum wavelength value in Angstroms. WMAX The maximum wavelength value in Angstroms. The program produces statistics for the number of spots in 15 wavelength bins, W to W+dW, covering the range WMIN to WMAX. DMIN This sets the lower limit on the interplanar separations included. To have all spots plotted choose DMIN < WMIN/2. High order index planes (and hence in general weak reflections) will have small interplanar spacings. The value of DMIN cuts off the spots due to interplanar spacings smaller than that length (in Angstroms). DISTANCE Sets crystal to film distance (in mm). For transmission (forward reflection) Laue enter a positive number. For back reflection Laue enter a negative number. RADIUS The radius of the exposed area of the film (in mm). Topographers Note: The cassette used in topography stations is square with half-diagonal (equivalent to radius here) of 50mm. Notice that if the SCALE option is set to FULL all spots within this radius will be plotted. If SCALE is set for another device and if the value of radius is large compared with the physical dimensions of the plotting area for that device (defined by the square border) will not be plotted. ORIENT It prompts for orientation method option. A = Mis-setting angles: This the orientation method used by protein crystallographers. Prompts are given for: * Reciprocal axis along rotation axis * Reciprocal axis nearest x-ray beam axis in plane defined by beam axis and rotation axis. * Camera type or fiducial positions. The prompt is: >>> Camera type (C1=Arndt-Wonacott, C2=Huber) or fiducial values [current values]: For the camera type question the reply may be either C1, C2 or six numerical values. The fiducial values for the camera types C1 and C2 are as follows: C1 (Arndt-Wonacott): -5000 -4000 -5000 4000 5000 4000 C2 (Huber): -4000 5000 -4000 -5000 4000 -5000 These fiducial positions are used by NEWLAUE for plotting on predicted patterns as well as being made available for output to a .gen file via the DATSAV command. V = Vector (Poles or Direction): This option, for topographers, invokes code taken from program FSLP. There is little information at present on the options and method of calculation used as well as conventions assumed. One useful reference is the report "Stereographic Projections by Computer", R.A.Ploc and P.C.Barnett, Atomic Energy of Canada Ltd, AECL-3741 and papers by the same authors in Journal of Materials Science (1983), pp.1083-1088 and in Journal of Applied Crystallography 5 (1972) 135 and 11 (1978) 713. It prompts for Poles or Directions defined in the form PQ=xx,xx,xx or Dq=xx,xx,xx where Q is A,B,C,1 or 2 and xx is a value For example: DA=1,1,1, DB=0,-1,0 is valid input DA corresponds to a vector parallel to the X-ray beam DB is a vector which is horizontal DC is a vector which is vertical and perpendicular to both DA and DB Note that only two directions are needed, the third is implied. The program gives the option to define the third direction as well but it is not clear why. Thus for <111> (body diagonal) enter DA=1,1,1 and DB=any vector perpendicular to [111], for example 1,-1,0. The effect of this method of orienting the crystal is to set the mis-setting angles Phi(x), Phi(y), Phi(z) by the amount needed to correspond to Orientation +c*,+a*. It is not clear at present whether this means that a* is along DA and a* along the projection of DB. The values of the vectors are listed before returning to the command menu. Notice that subsequent change of the mis-setting angles using the PX,PY,PZ commands directly does not update the vector values. PX Prompts for rotation angle phi_x (Deg) about the axis X along the beam axis. LAUE assumes that the beam is normal to the film. Positive is clockwise, looking in the same direction as the axis. PY Prompts for rotation angle phi_y (Deg.) about axis Y perpendicular both to the beam and the spindle axis, forming a right-handed set with these axes. The Y axis corresponds to the x (horizontal) axis of the coordinate system of the film. Positive is clockwise, looking in the same direction as the axis. PZ Prompts for rotation angle phi_z (Deg.) about axis Z along the spindle axis. This axis corresponds to the y (vertical) axis of the coordinate system of the film. Positive is clockwise, looking in the same direction as the axis. SYSTEM The SYSTEM command is used to set the axes and centring type used and hence the lattice type. If any conflict arises between the system and the cell parameters the crystal system takes priority. The value of `crystal system' determines the cell parameters required as shown below. System Required params Implied relationship Triclinic a,b,c No relationships alpha,beta,gamma Monoclinic a,b,c beta a NE b NE c (NB Unique axis b only) alpha=gamma=90 Orthorhombic a,b,c a NE b NE c alpha=beta=gamma=90 Tetragonal a,c a = b NE c Hexagonal a,c a=b NE c alpha=beta=90 gamma=120 Rhombohedral a alpha a = b = c alpha=beta=gamma not 90 Cubic a a=b=c alpha=beta=gamma=90 Lattice centring types P Primitive R Rhombohedral centring A Face centering in which there is a translation vector b/2 c/2 in the A faces of the basic unit of the space lattice B Face centering in which there is a translation vector a/2 c/2 in the B faces of the basic unit of the space lattice C Face centering in which there is a translation vector a/2 b/2 in the C faces of the basic unit of the space lattice F Face centering equivalent to simultaneous ABC face centering I Centering - translation vector a/2 b/2 c/2 giving point of intersection of body diagonal of the basic unit of the space lattice NB There are two possible descriptions of a rhombohedral lattice: on hexagonal axes (HEXAG,R) or on rhombohedral axes (RHOMBO,P) SPINDLE Defines the value of Spindle angle setting (in Deg). It is used on the oscillation/Laue camera to rotate the crystal. The spindle axis is always assumed to be perpendicular to the beam. It defines the Z axis and hence the Y axis (see PY, PZ). Clockwise is positive. ABSENCES Some Laue spots are absent in particular crystal classes because of the symmetry. For example 110 is forbidden in the face centred cubic system. For each class there is a set of relations between Miller indices hkl and their absence e.g. hkl must be all odd or all even for the reflection to be present, if they are mixed then they are absent. The absences routine allows the user to add more conditions, eg for the silicon structure where reflections such as 222 are forbidden because of the two atom basis. You can either ADD to the rules for the crystal which you have, or START again and enter the rules completely. Consult a crystallography text for the rules of systematic absences. COMMANDS Introduction This section describes the commands which may be used in the NEWLAUE program. List of subsections in this section: VECTOR DATGET DATSAV START PLOT HARD SCALE EXPAND INDEX AUTO PICSAV OVERLAY HELP STOP VECTOR Tells you which orientation vectors have been set in ORIENT when the VECTOR method has been used. The values of vectors are reset to 0.0 when the ORIENT command is invoked and recalculated only if the VECTOR suboption is used. DATGET The DATGET command allows the user to retrieve data created by a previous run of the program or in user edited file. The format is that of a .gen file as described in the documentation of the program GENLAUE. Directories may be given as part of the file name. The filename indicated on the screen will have the directory and file extension removed. If no extension is specified then .gen will be added. DATSAV The DATSAV command allows the user to save the data currently in use concerning his crystal. The file format is that of a .gen file as described in the documentation of the program GENLAUE. The file may either be used as input to GENLAUE for the next stage of the processing or for re-input to the program NEWLAUE. Directories may be given as part of the file name. The filename indicated on the screen will have the directory and file extension removed. If no extension is specified then .gen will be added. When a .gen file is written, then the user has the option of writing the current values of these items (default values or values read in as indicated above) or changing some or all of these values. The initial prompts are: >>> Enter name of .gen file in which to save control data [default]: >>> Use default .gen file specific parameter values [y]: or >>> Use current .gen file specific parameter values [y]: The first form with the word 'default' is output if no .gen file has been read in and the second form with the word 'current' is used if a .gen file has been read. If a reply is no to the latter questions, the following questions are then asked. A carriage return on its own confirms the current values in each case. >>> Film (f) or image plate (i) [f]: (The default will be 'i' if an image-plate .gen file was read in) >>> xc, yc, wc [current values]: >>> twist, tilt, bulge [current values]: >>> yscale [current value]: >>> xrej, yrej [current values]: >>> Camera type (C1=Arndt-Wonacott, C2=Huber) or fiducial values [current values]: (Omitted for image-plate) >>> wfilm, fidbox [current values]: >>> n1od [current value]: >>> Film type (CEA or DEF) or baseod, g1od, curve [current values]: (Omitted for image-plate) >>> nxrast, nyrast, dscal [current values]: >>> xcen_in, ycen_in, rastsiz [current values]: >>> xlow, xhigh, ylow, yhigh [current values]: For the camera type question the reply may be either C1, C2 or six numerical values. The fiducial values for the camera types C1 and C2 are as follows: C1 (Arndt-Wonacott): -5000 -4000 -5000 4000 5000 4000 C2 (Huber): -4000 5000 -4000 -5000 4000 -5000 For the film type question the built in consants provided for CEA abd DEF film are as follows: CEA: baseod = 0.04, g1od = 2.3, curve = 0.07 DEF: baseod = 0.05, g1od = 2.4, curve = 0.06 START The START command prompts for the projection to be used, and then commences the business of checking the input data, calculating the orientation matrix and the coordinates of Laue spots. The wavelength and spatial multiplicities are displayed before the spot diagram can be generated using the PLOT command. NOTE:At present only up to 50000 reflections can be calculated and displayed. More development work is required before the program can operate correctly for more reflections. Projections Currently LAUE supports three projections: a "Laue projection" in which the positions of reflections are generated as in an actual Laue diffraction experiment; and inverted Stereographic and Gnomonic projections of the poles of lattice planes. The projections are inverted through the centre of the plot so that a pole in the projection is plotted on the same side of the centre as the corresponding Laue reflection. This is done for ease of comparison. But note that spots on the OUTSIDE of the Laue plot correspond to spots on the INSIDE of the projected plots. Stereographic projection In this option, the crystal-to-film distance is interpreted as the radius of the sphere of projection, and the radius of the plot is the "Primitive circle" (90 deg circle). Gnomonic projection In this option, the crystal-to-film distance is interpreted as the radius of the sphere of projection, and the radius of the plot corresponds to a "Minimum circle" whose radius in Laue plot space is supplied by the user. (ie. spots just outside the minimum circle are mapped to spots just inside the boundary circle of the plot.) The mimimum radius specified must be greater than or equal to 5.0. PLOT This command displays the Laue spot diagram. The terminal is reset to alpha mode at the end of the plotting. If you are using a PERICOM terminal use the blue graphics/alpha keys to alternate between graph and menu. Black, Blue, Green and Red colours are used for spots and annotation. (Follow the PLOT command with the HARD command to get a Postscript file for hard copy output). When plots are drawn, the user has the option of selecting points to be drawn for each reflection as an alternative to a variety of symbols indicating the reflection types. The following new prompt is output: >>> Plot reflection symbols (S) or points (P) [S]: When symbols are requested, crosses (+) are used for singlet spots, and stars (*) for multiplets. To produce plots scaled for a given plotting device see SCALE. HARD The current spot diagram and parameters used to produce it are saved in a Postscript file. The file's name is either the default name 'laue.ps' or a user specified name. The file name is restricted to 16 characters. SCALE This command allows you to set a scaling factor for the pattern so that when plotted on a hard copy device it would match a Laue photograph obtained experimentally. It is also possible to achieve approximate scaling for the VDU type you are using. The options are: Plotters : HP7550 HP7475 HP7470 BENSON VERSATEC Terminals: SELANAR PERICOM VT100 VC3100 Other : FULL USERSET If no scaling is chosen (device FULL) the pattern is scaled so that all the spots within the RADIUS value chosen will be displayed. If USERSET is selected you may define your own scaling factor. Note: Not all the devices specified are still available at Daresbury. EXPAND The EXPAND command allows the user to specify an area of the screen to be looked at more closely, by specifying opposite corners of a rectangle to be enlarged. Note that expansion operates only once, i.e., you cannot expand an already expanded version. Also, a square section is expanded from the largest side of the rectangle specified. INDEX The INDEX command allows the user to identify, by Miller indices, screen co-ordinates and wavelengths, the reflections on the screen. You may optionally select to have indexed spots annotated on the graph. Four options are offered: NONE : No spots are annotated on graph. Indices and wavelength multiplicities are displayed on VDU screen only. SINGLE : Every spot identified using the cursor is annotated. LOW INDEX: All spots with indices above given value are annotated. This includes high order spots which after dividing their Miller indices by the higest common denominator have absolute values => index value supplied by user. ALL : All spots are annotated (up to highest index 99). In all last three cases only the low order index hkl sets are displayed if the spot ha wavelength multiplicity. The spots are identified by positioning the cursor and pressing any alphanumeric key EXCEPT~the~RETURN~or~ENTER~keys. YOU SHOULD NOT ATTEMPT TO ABORT THE OPTION BY TYPING CTRL/Z WHEN THE CURSOR IS DISPLAYED, THIS WILL CRASH THE PROGRAM SINCE THE CURSOR ROUTINE DOES NOT TRAP CTRL/Z (end-of-file mark). The presence of a graph is assumed. INDEX may also be used on an expanded graph but not on an expansion of an expanded one. Note:When the indices are found, the terminal is switched to its ANSI (text) mode and the found indices are output. The terminal isthen switched immediately back to the graphics mode. If the terminal is incorrectly set then the text screen may be cleared when the graphics mode is re-entered. This has the effect that the displayed indices on the text screen are lost. On a SELANAR the terminal can be correctly set by toggling the appropriate bit in SETUP C (this is the 4'th flag from the left. This can be achieved by the following keystrokes: SET-UP up_arrow up_arrow right_arrow right_arrow right_arrow right_arrow right_arrow right_arrow SHIFT/T SET-UP. AUTO The AUTO command calculates the orientation phi_x, phi_y, phi_z of the setting based on the positions of nodal spots on the film input via the tablet (or optionally from a file or the screen). Cell parameters and crystal to film information must be provided. The program calculates reciprocal lattice vectors for low index spots (a maximum of 280 can be stored) and the angles between them. These angles are then compared with those based on the entered spot positions in order to index these spots and hence derive the crystal orientation. The user is asked to supply a maximum value for h**2 + k**2 + l**2 to determine how many reciprocal lattice vectors are calculated. A low value speeds up the calculation but may miss some spots. A relatively high value should usually be used for non-primitive cells to make up for the systematically absent spots. The user is also prompted to supply an estimate of the error in spot positions, 1 to 1.5 mms is usually sensible for this. If the setting is known very approximately, for example the crystal may be approximately set along axes, then the user can greatly enhance the chance of AUTO succeeding (and speed up the operation) by providing the indices of the spots input, as well as their (x,y) positions. The indices can be obtained by generating the approximate Laue pattern, then using the INDEX option to identify the Miller indices of the specific spots which are chosen for AUTO. When refining angles during the auto-indexing option the program will issue a prompt after every 50 iterations even if its current refinement cycle is not complete. This allows the user to decide whether or not to continue the refinement which might possibly take a long time to complete. The program saves up to 12 solutions which have the lowest rms values. The rms values in question are either the unrefined values, if no refinement has been requested by the user, or the latest values when the user stops making further attempts at refinement. After all solutions have been processed, these values are displayed in a table on the terminal sorted in order of increasing rms values. The follwing prompt is then output: >>> Enter w=write a data file, r=reset spindle and phis or c=continue: The replies are as follows: w (write) Save the solution in a .gen file. This is useful if more than one solution is to be examined further. The position, in the table, is used to identify the required solution and the name of an output .gen file is requested as in the DATSAV option. The above prompt is then repeated. r (reset) Replace the current orientation parameters with the parameters obtained for one of the solutions given in the table. The position, in the table, is used to identify the required solution. The above prompt is then repeated. c (continue) Return to the main program menu without any further alteration to the orientation parameters. PICSAV The PICSAV command allows the user to save the co-ordinate points of all the spots in this picture in a format suitable for reading on the VAX thus giving far greater possibilities for graphic output. The extension ".pic" is always added to the filename. Only files in current directory are accepted. OVERLAY This command sets or clears the OVERLAY MODE. If set, then subsequent graphs are superimposed on one another. Spots from graphs produced in overlay mode use different point symbols from the first graph drawn with overlay mode off. This mode may for instance be used to produce patterns from twins or mixed crystal systems. If the HARD option is used after two or more patterns have been superimposed the parameter values of the last pattern are printed only. The overlay mode is set to OFF after HARD is used. HELP The HELP command has now been withdrawn as its impementation was Vax specific. The information that was presented is now available in this documentation file.. STOP Ends the program and resets the screen. If you abort the program by using the CNTL/C or CNTL/Y keys you will need to reset the terminal manually with the SETUP, RESET keys. FILES NEWLAUE can read and writes .gen files. Thus NEWLAUE may be used as a means of preparing a .gen file for GENLAUE if required. There are a number of items in the .gen file which are not required by NEWLAUE. These items have default values written into the program and are as follows: Compound name: TEST SAMPLE Crystal identifier: CRYS xc, yc, wc: 0.0, 0.0, 0.0 twist, tilt, bulge: 0.0, 0.0, 0.0 yscale: 0.998 xrej, yrej: 0.0, 0.0 nfid fidx1 fidy1 fidx2 fidy2 fidx3 fidy3: 3 -5000 -4000 -5000 4000 5000 4000 wfilm fidbox n1od baseod g1od curve: 120 5 128 0.04 2.3 0.07 nrec ylen dscal: 2400 2400 2.0 If a .gen file is read into NEWLAUE then these default values are replaced by the corresponding values read fron the .gen file. When a .gen file is written, the current values for these parameters may be written to the output file or, with the exception of the compound name and crystal identifier, alternative values may be specified.