Program N_TREOR This is a Stand-alone version of N-TREOR (Version November 2001). REFERENCE: New techniques for indexing: N-TREOR in EXPO J.Appl.Cryst. (2000) 33, 1180-1186 by Angela Altomare, Carmelo Giacovazzo, Antonietta Guagliardi, Anna Grazia Giuseppina Moliterni, Rosanna Rizzi and Per-Erik Werner In EXPO N-TREOR is combined with a peak search procedure and a a program set for structure determination. If you are only interested in the description of how to run the stand-alone version of the program go to the headline: N_TREOR the stand-alone version of N-TREOR N-TREOR is an updated version of TREOR90 described below. The main changes characterizing N-TREOR with respect to TREOR90 are: 1) Several decisions normally made by the user of TREOR90 when a default run fails or when a plausible solution is found are automatically made by N-TREOR a) If a plausible solution is found, the program will look for a better solution in a new trial. This may be achieved by requiring a higher M20 value or all the first 20 peaks to be indexed. b) If N-TREOR does not give satisfactory results in a default run it repeats the unit cell search with wider error limits. c) If still no solution is found, with increased error limits, the maximum (h,k,l) Miller indices of the orthorhombic and monoclinic base lines are increased and the error limits are narrowed to 50 per cent of the default values. ( The reason for narrowing the error limits is that a large number of big erroneous unit cells may otherwise be generated in the trial phase, and correct solutions may be eliminated before refinement. There is always a need for more accurate data if the unit cell is big. The figure-of-merit tests will otherwise fail.) 2) The maximum cell volume and the maximum axis values investigated by N-TREOR are preset to 4000 A**3 and 35 A, respectively, if the maximum observed d value is greater than 10 A. If the maximum observed d value is less than 10 A the default values are (as shown in the key-word list below) 2000 A**3 and 25 A, respectively. 3) The dominant-zone tests in TREOR90 are improved, and if more than the first 6 peaks can be indexed by a two-dimensional rectangular cell,but no three-dimensional solution is found, the first non-indexed peak is reported as a possible impurity line that may be deleted in a second test. 4) Wavelength dependence of D1 and D2 ( error tolerances, see the key-word list below) has been introduced. DX = DXCU*(WAVE/1.560598)**2 where DXCU = D1 and D2, respectively for CuK-alpha1 radiation 5) Monoclinic solutions are checked for possible rhombohedral symmetry. If the symmetry is likely to be rhombohedral, N-TREOR calculates the possible hexagonal axes. In general, however, the user has to check for geometrical ambiguities. In N_TREOR version November 2001, the rhombohedral cell is also calculated and refined by PIRUM. 6) Possible zero shifts in the 2theta Bragg angles are tested. This is probably the most serious systematic error in powder diffractometry. Although well known internal as well as external standard techniques exist, that can be used to completely (or almost completely) eliminate such errors, the experience from a large number of indexing problems shows that zero point errors are frequently present in data sets used for powder indexing. ( Suprisingly enough even modern synchrotron data, with exceptionally high RELATIVE accuracy, have sometimes been found to have large zero point errors.) 7) N-TREOR can suggest more candidate cells. In order to select the best one, a modified de Wolff figure of merit M'20 has been introduced, with M'20 = (7-N-par)*M20 where N-par is the number of cell parameters to be determined ( 6 for triclinic system, 4 for monoclinic,....), and M20 is the original de Wolff figure of merit. Even if M20 is effective as a tool for assessing the reliability of the cell, the heuristic relation (M'20=.....) allows the program to prefer a higher symmetry solution. E.g., of two possible unit cells with similar M20 values, one monoclinic and the other triclinic, the monoclinic is preferred by M'20. 8) The program PIRUM, originally an interactive program, has been suitably modified in order to automatically refine the unit cell parameters, without any user intervention. N-TREOR performs a statistical study of index parity of the assigned reflections in order to detect the presence of doubled axes or of additional lattice points (A, B, C, I, R or F centered cell). In these cases the reflection conditions are suitably modified and the unit cell parameters are again refined. It may be necessary, however, to check the analysis of systematic extinctions in the condensed output file. The reason is that mutiple indexed reflections are not analysed. Additional refinement cycles are performed by PIRUM after applying two supplementary origin shifts of +- one quarter of the shift value selected by N-TREOR according to point 6 above. See also the keyword LIMIT. 9) It is not possible to interrupt the calculations by pressing the key...I....This option is not available any more. 10) It is still possible to use the key-word system as described below. Normally, however, only WAVE, LIMIT, CHOICE and END* are needed. 11) Last corrections: September 2004 N_TREOR the stand-alone version of N-TREOR Data input file: LINE ONE. TITLE Any text in col.2-80 LINE SET TWO. One d-value and optionally INTENSITY on each line. (Free format) It is also possible to use other types of input data. See keyword CHOICE in the key- word list below. Note that the default value of CHOICE is 4 in N_TREOR (i.e. d-values). The data must be given in order, starting with the low order lines.Generally the first 20-25 lines should be used. Remaining lines (if any) may be used in later final refine- ments. (The dialogue version of PIRUM). INTENSITY data is optional. It is never used by the program. If present, it will only be transferred to the output lists as (max.) 4 digits integers. STOP LINE FOR LINE SET TWO IS A BLANK LINE The program will ask for file names and Theta-shift. The normal Theta-shift is the step used in the data collection; (Usually =0.01 for Wave=1.5406 A) i.e. 0.02 deg. in 2theta. LINE SET THREE.(See general description of this line set below) Below 4 examples of line set three for N_TREOR are given: Example 1. end* Comment:This is the only key-word needed for a normal N_TREOR run if CuK-alpha1 radiation has been used and d-values are given. Example 2. wave=1.529, end* Comment: Error limits are wave-length dependent. It is therefore important to give the wave-lenght, although the d-values are not wave-lenght dependent. Do not forget the , before next key-word. Example 3. limit=1, end* Comment: LIMIT is a new key-word, not present earlier in TREOR. LIMIT can be any integer between 1 and 10 (default value is 10). The parameter LIMIT can be given a value less than 10 in order to reduce the maximum number of zero shifts tested. ( The number of zero shifts applied will be dependent on the result of the the calculations and is therefore usually less than 10). If you have good reason to belive that the zero point error is negligible, you may save a lot of computer time by limit=1 Example 4. choice=0, limit=1, end* Comment: Input data are sine square theta values derived from an instrument with strictly monochromatic CuK-alpha1 radiation. About limit=1, see example 3 above. WARNINGS: A) In the present version of the program the automatic tests of index parity are NOT complete. See comments about VISUAL INSPECTION in the condensed output list. ( May be improved in a later program version.) B) On a PC the calculations may be very time consuming. It can be recommended to use LIMIT=1 ( see Example 3 above) in a first test run, at least if you have good reason to believe that the zero point error is negligible. C) OBS ! The best solution may be found ANYWHERE in the condesed output file. The reason is that the automatic Pirum run may fail and give more unindexed lines. It is recommended to always use a separate Pirum run to check the possible solution/s. END OF INSTRUCTIONS FOR THE STAND-ALONE VERSION OF N-TREOR TTTTTTT RRRR EEEEEE OOOOO RRRR 9999 0000 T R R E O O R R 9 9 0 0 T R R E O O R R 9 9 0 0 T RRRR EEEEEE O O RRRR 99999 0 0 T R R E O O R R 9 0 0 T R R E O O R R 9 0 0 T R R EEEEEE OOOOO R R 9 0000 This is an autoindexing program and should be used to find a physically plausible solution to an indexing problem. It is not a program for final refinement of the unit cell dimensions. (c.f. PIRUM or NBS*AIDS83) Information for the programmer: The TREOR99 version is written in Microsoft Fortran for PC. The interrupt function (..press the letter i) is generally not needed and if you want to use other compilers and/or computers whithout any possibility to sense keyboard oper- ations during the calculations, you may skip this in the following way Skip the include 'flib.fi' and 'flib.fd' Replace the original subroutine chkkey by the following dummy routine: Subroutine chkkey(code) integer*2 code code=0 return end End of programmers information FEBRUARY 1999 Two new keywords, ZEROSET and KREDUCT, are introduced. See keyword list below. KREDUCT may be used in order to save only the best trial solutions.It is presupposed then, that input data is very accurate and a cell axis is expected to be very long (>20 A) In order to avoid run time errors, a large number of small changes have been made. The earlier program versions were somewhat dependent on the wave-length used. The parameters: D1, D2 and SSQTL are now automatically ( i.e. if SELECT=-1 ,the default value) adjusted for WAVE. As in earlier program versions one can forget about WAVE if CHOICE=4, i.e. if input data are d-values. If correct WAVE is given, also the correct diffraction angles are reported in the output lists. Monoclinic solutions are checked for possible rhombohedral symmetry. If several high symmetry solutions are found, it is recommended to decrease the parameter VOL in order to find the smallest monoclinic cell (1/6 of the 'true' hexa- gonal cell). Then treor will check if the true symmetry is likely to be rhombohedral and calculate the possible hexagonal axes. SEPTEMBER 1997 (Below 1997 is changed to 1999 (changed february 1999) Keywords (see below) may be written in upper or lower case. The program will run in 1999 mode if the keyword SELECT has the value -1 (default value) The normal TREOR99 run is: Title line: Any text in col. 2-80 Data lines: Less than 99 lines (Normally 20-25 lines). Free format. Only one value on each line is necessary. Optionally an intensity may be given as a second data on the line. Intensity data is not used by the program and only transferred to the output as max. 4 digits integers A blank line: choice=4, (if D-values were used on the data lines) Do not forget the , before next key-word Key-words may be given in arbitrary order end* and on arbitrary lines but the end* must be the last one. Usually no other key-words than CHOICE need to be changed from the default values. If a solution is found the program will automatically decrease NIX or increase MERIT and rerun the problem in order to check if a better solution can be found. However, triclinic tests will only be made if the first solution was triclinic. If no solution is found the program will automatically change the parameters IDIV and the error limits D1 and D2 i.e. add 0.0001*(Wave/1.5405981)**2)and rerun the problem. Of course the TREOR99 mode may take a minute more on a Pentium but it is a strongly recommended, simple, safe and powerful procedure. Powder indexing depends to 95 per cent on the data quality If you are using a Pentium processor the CPU time will probably never exceed 3 minutes. ( A PC 486 may be almost ten times slower.) It is possible to run the program in TREOR4, TREOR90 or TREOR99 mode TREOR4 mode if VOL=positive value and SELECT is zero (or more.) TREOR90 mode if VOL=negative value and SELECT is zero (or more) TREOR99 mode if SELECT=-1 (Default) (Normally preferred mode) If the TREOR99 mode is used, you normally only need to choose the correct CHOICE parameter. In all modes most of the key-word parameters may be changed from the default values ( c.f. the key-word list below). Do not change default values without good reasons. For example do not change VOL from the default value, unless you have reason to believe that the unit cell volume is more than 2000 A**3 ( You will not reduce CPU time much if you are using a VOL that is less than 2000 A**3) If TREOR99 finds a plausible solution ( according to the NIX and MERIT parameters) this will be announced on the display and the user may ( although not recommended ) interrupt the calculations and look at the condensed output list. Interrupt...Press the letter i on the keyboard. If the program is running on a Pentium there is no reason for doing so. (The interrupt function is at present only available in the PC version.) Problem notes: A solution that is very close to a higher symmetry may be lost in the condensed output list, although an acceptable lower symmetry solution may be reported in the long output list. This may happen if two cell dimensions are almost equal. JUNE 1995 Within this section JUNE 1995 to JUNE 1992 a very short description for the lazy user of TREOR90 (who does not want to read the complete documentation file) is given. The normal TREOR90 run is: Title line Data lines (20-25 lines, only one d-value, col.1-16, on each line) A blank line CHOICE=4, (if D-values were used on the data lines) SELECT=0, (to avoid the TREOR99 mode) VOL=-2000, (note the minus sign. --- a treor90 run) END* (stop card) It is now possible to interrupt the calculations. This can be done by the letter i on the keyboard. Then the program will stop soon. The reason for this option is that one should never hesitate to give VOL eq. a negative value i.e. the most efficient test of all symmetries. One can always stop for example triclinic tests from the keyboard. (This option is only available on the PC version of the program) If no solution is found rerun the problem, but include the following keyword line: D1=0.0003, D2=0.0005, (note the , after each keyword value) Another test is to set IDIV=0, which means that the first seven lines will not be changed by ( what treor may erronously judge as ) higher order lines The parameters MERIT and NIX may also be changed. ( See below.) The first accepted solution may not be the best one. Therefore, follow the instructions on the condensed output file. Do not print the long output file. Use your editor and look for the most promising M-TEST lines. This program should be used to find a physical plausible solution of the indexing problem. The refinement is only preliminar ( not Hess weighted). Especially lines at high diffraction angles may be unindexed by treor. In final refinements all lines, all extinction conditions, Hess weights ( i.e. all diffraction angles should be given equal weights)and all knowledge about intensity distribution between overlaps ( if available) should be used. From this department a flexible dialouge program,(version 930101 of PIRUM) may be distributed. Program NBS*AIDS83 is (probably ?) recommended by ICDD. JUNE 1992 It is strongly recommended to run TREOR90 on a PC/AT using a 486 CPU. Otherwise the VAX version TREOR90V may be used. AUGUST 1990 OBS. In the PC/AT and VAX versions the subroutines ORTAL, MAEG and COUNT are not vectorized. Vectorized versions of these subroutines are available for CONVEX computers. The original TREOR90 has been written for a CONVEX vector processor. This should be kept in mind when the comments below are read. The program may be very time-consuming on a PC (unless a 486 processor is used). 1) Dominant zone test is added for the orthorhombic symmetry. 2) Dominant zone test is added for the triclinic symmetry. 3) Higher order lines among the first seven lines ( used in the base line sets) are automatically excluded from the trial phase of the calculations. 4) If a monoclinic or triclinic solution is found, the program will end with a unit cell reduction followed by a conversion of the reduced cell to a conventional cell according to the metric symmetry. The reduction should be valid unless syste- matic extinctions are found in the trial cell. 5) If a satisfactory solution is found, only the condensed out- put file is needed. It contains all relevant information and only one indexed list. 6) The general output list (that is normally not needed, cf. 5) will only list trials where M20 ( or Mxx if less lines are available) is 6 or more and not more than 3 lines among the first 20 (or xx) lines are unindexed. 7) If the parameter VOL is given with a negative sign all symme- tries are tested until a final solution is found- if possible. OBS. This is the normal procedure for TREOR90 8) An algoritm for successive reduction of trial-cell volumes is used in monoclinic and triclinic tests if a negative VOL parameter is given. It is based on the input cell volume limit and the number of trial cells found with IQ ( See keyword IQ) or more than IQ indexable lines. 9) It is strongly recommended to give only the first ( well checked and accurately measured ) 25 lines in the diffraction data list (See LINE SET TWO). 10) It is expected that more than 95 per cent of monoclinic and higher symmetry patterns and probably more than 50 per cent of triclinic patterns will be indexed PRESUPPOSED the DATA QUALITY is high ( i.e. average differences between calculated and observed diffraction angles less than 0.02 deg. and also the weak lines included in the data). The experience of tri- clinic patterns is limited, however. 11) Obs. It is important to check cubic, tetragonal and hexagonal solutions by a second run with KS=0 and THS=0 ( See key-word list.) Do not trust cubic, tetragonal or hexagonal solutions without an orthorhombic test. 12) The reason for testing the symmetries in correct order ( from cubic to triclinic) and to START the orthrhombic, monoclinic and triclinic tests with dominant zone tests is that by this procedure false solutions are avoided. 13) For a normal TREOR90 run only the keywords CHOICE=X, (see key-word list) VOL=-2000, (OBS. The negative sign.) END* should be given after the diffraction data list. Computing times of more than 1 minute is rare for monoclinic or higher symmetries on a CONVEX computer. Computing times of more than 5 minutes (on a CONVEX) for a triclinic pattern has not yet been found. For a VAX (Micro VAX II) computing times may be more than 50 times longer. The PC/AT 486 is faster than the Micro Vax II (but slower than CONVEX. The source code for VAX is not exactly the same as for CONVEX. There are very small differences between the PC/AT and the VAX versions, however. The input of file names and OPEN statements must be changed if you want to run this program on a VAX. Furthermore, the VAX version uses a subroutine to measure the CPU-time. 14) The input format for LINE SET TWO ( See below) is changed in agreement with the output format of the diffraction data file from the Guinier-H{gg film scanner system ( at Stockholm University). The change is mainly of interest for output of intensities. 15) The original key-word instructions given below are relevant as long as a positive VOL parameter is given. 16) If VOL is given a negative value (see 13 above) the following key-words are fixed: MONO=135 and MONOSET=7. Other key-words may be used as in the description below. 17) On the output lists M-TEST= xx UNINDEXED IN THE TEST= y usually means that xx is identical with M(20) and y is the number of unindexed lines whithin the first 20 lines ( i.e. used for the MERIT test). If less than 20 lines are available xx and y refer to the number of lines used. November 1988 29 11 88 Trial-and-error program for indexing of unknown powder patterns. Cubic, Tetragonal, Hexagonal, Orthorhombic, Monoclinic and Triclinic symmetries. Version 2 1/9-75 = Version 26/4 plus DENS,EDENS and MOLW. See Keyword list below. Version 3 8/5-80 New output form Version 4 2/10-84 = Version 3 plus The following new options.... 1. IDIV. See keyword IDIV below. 2. Monoclinic (020)-test Ref: Smith,G.S. and Kahara,E J.Appl.Cryst. 8 (1975) 681 3. SHORT. See keyword SHORT below. Short axis test. (Indexing of dominant zones.) 4. TRIC. See keyword TRIC below. Indexing of triclinic patterns. The source code was modified in order to decrease the CPU-times in September 1988. The changes have no influence on input or out- put from the program, but CPU-time reductions of 20-50 per cent have been observed. Version 5. (=Version November 1988) 29/11 1988 Dominant zone test introduced also for orthorhombic symmetry. In version 4 high symmetry short axis solutions were only found indirectly from the monoclinic tests. Condensed output file. A complete list of observed and calculated lines is only given for the solution (if it is found) i.e. for an indexing where the stop limits ( See keywords MERIT and NIX ) are fullfilled. Normally only the condensed output file is needed. If the stop limits are fullfilled the unit cell is refined three cycles more. OBS. Final least-squares refinement should be made by a separate program (for example by PIRUM). The TREOR program is written in order to FIND a plausible cell, not to produce the ultimate refinement. Only the first part of the difference analysis table is printed if no solution is found. (Usually it is not needed as you should rerun the problem after modifications of the input data.) If you have any questions, write to.... P.-E.Werner Dept. of Structural Chemistry Arrhenius Laboratory Stockholm University S-106 91 Stockholm, SWEDEN TEL: 08 / 16 23 93 FAX: 46-8-15 21 87 EMAIL: PEW@STRUC.SU.SE It is believed, however, that the following documentation should be sufficient for all careful readers. GOOD LUCK! R E F E R E N C E S Basic principles. Werner,P.-E., Z.Kristallogr. 120 (1964) 375-387 TREOR, a semi-exhaustive trial-and-error powder indexing program for all symmetries. Werner,P.-E., Eriksson,L. and Westdahl,M., J. Appl. Crystallogr. 18 (1985) 367-370 Refinement of unit cell. Werner,P.-E.,Arkiv Kemi 31(1969) 513-516 Figure of merit. De Wolff,P.M.,J.Appl.Crystallogr. 1(1968)108-113 Geometrical ambiguities. Mighell, A.D. and Santoro, A., J. Appl. Crystallogr. 8 (1975) 372 G E N E R A L C O M M E N T S This is a general trial-and-error indexing program for X-ray diffraction powder patterns (i.e. all symmetries included). Historical information ---- In order to reduce computing times on computers without hardware floating point processers, parts of the program have been written for integer calculations. The parameters given as normal values in the keyword list below should be regarded as an important part of the program. They are based on experience from many successful runs on structures con- firmed by single crystal data. The parameters VOL and CEM, however, may be selected for the actual data set and the symmetry tried. ...For a monoclinic trial the parameter MONO must be non-zero. ...For a triclinic trial the parameter TRIC must be 1. A TREOR99 run (i.e. SELECT=-1) or a TREOR90 run (i.e. VOL= a nega- tive number) will automatically check all symmetries. Most of the powder patterns used to check the program have been obtained by focusing Guinier-Hagg cameras. The photographs have been measured by.... 1. The method described by Hagg,G., Rev.Sci.Instr.18 (1947) 371 and Westman,S. and Magneli,A., Acta Chem. Scand. 11 (1957) 1587 2. The method described by Malmros, G. and Werner, P.-E., Acta Chem. Scand. 27 (1973) 493 3. The film scanner system SCANPI ( written for the Guinier film scanners LS18 and LS20) The program has also been tested on a large number of NBS-data sets. (JCPDS data sets.) The accurate data obtained by NBS,National Bureau of Standards, is clearly sufficient for successful indexing (in spite of the fact that they are now usually obtained by powder diffracto- meters. Unfortunately, however, many diffractometer data sets found in the litterature show parabolic deviations between observed and calculated diffraction angles.) The following citations, however, should be emphasized.... 'The paramount importance of resolution for indexing work explains the high success rate for focussing camera data, especially from Guinier-Hagg instruments, whose resolution can omly be described as superb. It is rather less common ( and considerably more expensive) to obtain as good resolution with diffractometer data.' 'Powder indexing is not like structure analysis, which works well on good data, and will usually get by on poor data given a little more time and attention. Powder indexing works beautifully on good data, but with poor data it will usually not work at all' Ref: Data accuracy for powder indexing.Shirley,R.NBS Spec. Publ. 567 (1980) P.370 and P.362 respectively. WARNING! A zero point error is much more serious than statistical errors of the same magnitude. Sigma(Two theta) should be less than 0.02 deg. ******************************************* * DO NOT WASTE COMPUTER TIME ON BAD DATA. * ******************************************* An indexing algorithm cannot be statet rigorously because of the unpredictable distribution of unobserved lines and the errors of measurements. One (or two) impurity lines may cause severe problems.More impurity lines may rule out the possibility to index the pattern. Therefore, it is expected that various methods may be useful for various powder patterns. For example, a multitude of non-systematic extinctions may not appreciably affect the power of trial-and-error methods. The least-squares refinement of the unit cell dimensions should normally not be considered as the ultimate refinement. The main purpose of this program is to FIND the unit cell. The program PIRUM (version 930101) may be used for ultimate refinements. PIRUM (version 930101) is a dialouge version of the old PIRUM ( cf. ref. /Refinement of unit cell/ given above. ) Extinction conditions, Hess weights and max. accepted deviations between observed and calculated 2theta in degrees are normally used in version 930101 of PIRUM. In old PIRUM versions, parameters like D1, D2, and SSQTL (See keyword list below) were used. (cf. also the NBS*AIDS83 program. PIRUM is designed to be more user-friendly, however.) A limited number of nonsense cells may be printed on the output file. You should look for max. De Wolff figure of merit ( not F- index) and min. number of unindexed lines. WARNING. You should not accept unindexed lines unless you are able to explain them. On the other hand, you should not put in uncertain (doubtful) lines in this program. They may be tested later by a refinement program (ex. PIRUM). I N P U T D A T A LINE ONE. TITLE Any text in col.2-80 LINE SET TWO. One SQ and optionally INTENSITY on each line. (Free format) SQ (=Sine square theta). It is also possible to use other types of input data. See keyword CHOICE The SQ data must be given in order, starting with the low order lines.Generally the first 20-25 lines should be used. Remainging lines (if any) may be used in later final refine- ments. (Program PIRUM). INTENSITY data is optional. It is never used by the program. If present, it will only be transferred to the output lists as (max.) 4 digits integers. STOP LINE FOR LINE SET TWO IS A BLANK LINE (OR A NEGATIVE SQ) LINE SET THREE. GENERAL INSTRUCTIONS. All parameters in line set three have preset values. A preset value is denoted 'NORMAL VALUE 'below. Any 'NORMAL VALUE' may be changed in the following way: KEYWORD1=VALUE1, KEYWORD2 = VALUE2, KEYWORD3=VALUE3, ......., END* 1. The keywords are listed below 2. You must not forget = 3. The value may be given in free format ( integer or real ). 4. You must not forget , You may use arbitrary positions on the lines. All blanks are irrelevant. The number of lines is arbitrary. You may give one or more key-word(s) on each line. Line set three must end with the keyword END* (OBS. asterisk) S T R A T E G Y Unless SELECT=-1 (i.e. TREOR99 mode) or VOL=(a negative value) (i.e. TREOR90 mode), the program will not differ much from earlier program versions. See the comments on the top of this list. Then (if VOL=negative value) only parameters such as NIX, MERIT, IDIV and in exeptional cases D1, SSQTL and/or D2 may be changed if indexing is not successful. Usually the main problem, however, is the quality of your diffraction data. Therefore, if the first run does not give a satisfactory solution, it may be recommended to increse D1 and D2 to 0.0003 and 0.0005, respectively. If the TREOR99 mode is used (i.e. SELECT=-1) the program will, regardless a plausible solution is found or not, rerun the problem and automatically adjust some parameters. Therefore it is generally not necessary for the user to change the default values in MERIT, NIX, D1 and/or D2. If you are not using the normal TREOR99 procedure (i.e. SELECT= -1) or the TREOR90 procedure(i.e. negative VOL) the standard procedure is to start with the higher symmetries: cubic, tetragonal, hexagonal and orthorhombic ( in one run). Next the monoclinic symmetry may be tried. More than one job may be needed..successively increasing the number of base line sets, and cell volume (See keywords: VOL, CEM and MONOSET) If formula weight and density are known, they may be used. (See keywords: DENS, EDENS and MOLW). The CPU-time needed will then usually be strongly reduced. (Unfortunately they are usually not known and therefore they have not been used very much.) LINE SET THREE EXAMPLES: (TREOR4, not normal TREOR99 or TREOR90 examples.) EXAMPLE 1.Next line (except C in col.1) represents a line set 3 END* Cubic, tetragonal, hexagonal and orthorhombic symmetries are tried. It may be recommended to try a smaller VOL limit even if a solution with acceptable figure of merit has been obtained. Sometimes it is difficult to find the necessary transformations between a high symmetry unit cell of too large dimensions and the primitive one. EXAMPLE 2. Next two line is a line set 3. KS=0,THS=0,OS1=0, CEM=20, V O L = 1000 , MONO=130,END* This is an example of a first monoclinic trial. ( See keyword MONO). Note that it is irrelevant if you give 'CEM=20.0' or 'CEM=20' etc. EXAMPLE 3. Next.....etc. KS=0,THS=0,OS1=0, CEM=20, VOL=1500, MONO=130, END* If example 2 is unsuccessful you may increase the VOL parameter to 500 EXAMPLE 4. Next.....etc. KS=0,THS=0,OS1=0,CEM=20, MONOSET=7,LIST=1, DENS=3.123,EDENS=0.2,MOLW=234, END* If you have any possibility to put in density and formula weight, the CPU-time will be much reduced. This may also be tried if you expect that the lattice contains a dominant zone i.e. if in a test run you get a large number of trial cells when using the keyword SHORT=1. EXAMPLE 5. Next....etc. CEM=20,VOL=700,TRIC=1,MERIT=20,END* This is a triclinic test ( OBS. time-consuming) ( See. keyword TRIC). A de Wolff figure merit of 20 may (sometimes) be needed for a triclinic cell The examples given above illustrate a step-wise strategy for indexing. However, the VOL parameter may be estimated from the D-value of the 20th line. (cf. keyword TRIC) WARNING. If the unit cell has a small volume, for example 250 A**3 and VOL=2000 is used, the correct solution may be lost in the trial process. The reason is that a large number of large trial cells may erroneously index more lines than the correct cell. The problem is less severe in TREOR90 as a negative VOL parameter will cause the program to test ( for all but the triclinic symmetry) half the maximum volume in a first step. WARNING. Estimation of the unit cell volume from the relations VOL(monoclinic cell)= 20*D(20)**3 where D(20)= the D- value for line number 20, and VOL(orthorhombic) = 31*D(20)**3 are much less reliable than the correspon- ding relation for the triclinic symmetry. VOL(triclinic)=13.39*D(20)**3 Triclinic structures have no systematic extinctions! For structures containing atoms with large differences in scattering factors ( eg. metal-organic structures) the general rule may fail also in a triclinic case. Ref: Smith,G.S. J Appl. Crystallogr. 10 (1977) 252 It is usually easy to put in a known ( or expected ) cell edge into the program. Example: A monoclinic trial with the restric- tion that one cell axis is X.XX A. Add this D-value in line set two. Suppose it will be line number 2. Then set MH2=1, MK2=1, ML2=0 and MS2=1. Then the line will be used as A-axis or ( the unique) B-axis in the monoclinic test. Conclusion: It is usually easy to put in prior knowledge and constraints -for example density- into the program. ( This statement is made because of some misunderstandings in the literature.) H O W T O I N T E R P R E T T H E O U T P U T. As in all good detective stories, the solution of the problem will usually be given on the last page..... i.e. the output list will be interrupted as soon as a unit cell that will satisfy the criteria set by the keywords NIX and MERIT are fullfilled. The main rule is that if all the first 20 lines are indexed and the De Wolff figure of merit M(20) is greater than 9, then the indexing problem is in principle solved. This does not mean that the cell is reduced, that a cell axis may not be double etc., UNIT CELLS OBTAINED BY THE PROGRAM SHOULD BE CAREFULLY CHECKED A. If M(20) is less than 10 or more than one line is unindexed within the 20 first observed lines the solution is probably meaningless. Is any low-order line wrong ? B. Check for common factors in the quadratic forms. Example: A teragonal pattern may have H*H + K*K = 5*N i.e. the A-axis is 2.3607 ( square root of 5) times shorter than given on the output list. Example: If all H, K or L are even, the corresponding cell axis should be divided by 2. C. If the unit cell obtained is centered, derive a primitive cell. ( Run program MODCELL or a corresponding NBS program) D. Reduce the primitive cell and derive the conventional cell. (Run program REDUCT or a corresponding NBS program) E. Hexagonal and tetragonal cells are sometimes indexed as orthorhombic. Example: A=B*1.7321 i.e. a possible hexagonal cell. F. Check for geometrical ambiguities. See reference above. It is also strongly recommended to chech cubic, tetragonal and hexagonal solutions by an orthorhombic test. Put KS=0 and THS=0 and re-run the problem. There are two reasons for this procedure.... 1. It may help you to identify geometrical ambiguities. 2. It has been found that sometimes very small orthorhombic unit cells can be indexed in an acceptable way ( i.e. fullfill the De Wolff criteria) by a larger unit cell of higher symmetry. Although the unit cells are sometimes related to each other, the relations are often difficult to detect, and therefore it is often convenient to let the program derive both solutions. G. The De Wolff figure of merits are derived from the assump- tion that no systematic extinctions are present and that all lines are indexed. A high figure of merit has no meaning unless all lines are indexed. The De Wolff figure of merit will increase in the final refinement made with program PIRUM, where the systematic extinctions can be taken into account. H. If possible, use the density and formula weight to check that the unit cell contains an integral number of formula units. I. If a cell axis is more than 20 A....be suspicious! It has been found that the De Wolff figure of merit may fail in such cases. (Require M(20) > 20) J. If one cell edge is much shorter than the others..........be suspicious! It may be a dominant zone problem and the De Wolff test may fail. (This problem is usually not severe in TREOR90, where short axis tests are made prior to the general tests.) K. If a table starts with...NOT REFINED UNIT CELL... two parameters are probably almost identical ( the symmetry may be higher) and the trial cell parameters are used to print the list. L. If no satisfactory solution is found ( See the keywords NIX and MERIT), the program may end with a small table con- taining a difference analysis. The program is described in Z. Kristallogr.120 (1964) p.381-382 (Werner,P.-E.) where it is named I1. The most interesting differences are those that have high multiplicities (on the top of the list) and are not too small (to the right of the list). In the present state of the program, the difference table is usually not needed. M. Why not solve the crystal structure from your powder data ? This is the ultimate way to prove the unit cell! K E Y W O R D L I S T KEYWORD. NORMAL COMMENT. VALUE. KH =4 Max H for cubic base line. KK =4 Max K for cubic base line. KL =4 Max L for cubic base line. OBS. The program will only generate H greater than or equal to K and K greater than or equal to L for this line. KS =6 Max H+K+L for this line. OBS. If KS=0 cubic test omitted. OBS. The cubic base lines are (1) and (2). * * * * * * * * * * * * * * * * * * * * * * * * * * * THH =4 Max H for tetragonal and hexagonal base lines. THK =4 Max K for tetragonal and hexagonal base lines. THL =4 Max L for tetragonal and hexagonal base lines. OBS. The program will only generate H greater than or equal to K for these lines. THS =4 Max H+K+L for these lines. OBS. If THS=0 tetragonal and hexagonal tests omitted. OBS. The tetragonal and hexagonal base lines are (1,2),(1,3) and (2,3) * * * * * * * * * * * * * * * * * * * * * * * * * * * OH1 =2 Max H for the first orthorhombic base line. OK1 =2 Max K for the first orthorhombic base line. OL1 =2 Max L for the first orthorhombic base line. OBS. The program will only generate H greater than or equal to K, and K greater than or equal to L for this line. This is also valid if the SELECT parameter is used. (See below). OS1 =3 Max H+K+L for this line. OBS. If OS1=0 orthorhombic test omitted. OH2 =2 Max H for the second orthorhombic base line. OK2 =2 Max K for the second orthorhombic base line. OL2 =2 Max L for the second orthorhombic base line. OS2 =4 Max H+K+L for this line. OH3 =2 Max H for the third orthorhombic base line. OK3 =2 Max K for the third orthorhombic base line. OL3 =2 Max L for the third orthorhombic base line. OS3 =4 Max H+K+L for this line. OBS. The orthorhombic base lines are (1,2,3) (1,2,4) (1,2,5) (1,3,4) (2,3,4) and (1,2,6) if SELECT0=0 (See SELECT below) * * * * * * * * * * * * * * * * * * * * * * * * * * * MH1 =2 Max Abs(H) for the first monoclinic base line. MK1 =2 Max K for the first monoclinic base line. ML1 =2 Max L for the first monoclinic base line. OBS. The program will only generate H greater than or equal to L for this line. EQ. TO L FOR THIS LINE. This is also valid if SELECT is used. (See SELECT below) MS1 =2 Max Abs(H)+K+L for this line The normal (and fast) way to test an expected cell axis is to put it in as SQ number one (in card set two) and set MH1=1, MK1=1, ML1=0 and MS1=1 MH2 =2 Max Abs(H) for the second monoclinic base line. MK2 =2 Max K for the second monoclinic base line. ML2 =2 Max L for the second monoclinic base line. MS2 =3 Max Abs(H)+K+L for this line. MH3 =2 Max Abs(H) for the third monoclinic base line. MK3 =2 Max K for the third monoclinic base line. ML3 =2 Max L for the third monoclinic base line. MS3 =3 Max Abs(H)+K+L for this line. MH4 =2 Max Abs(H) for the fourth monoclinic base line. MK4 =2 Max K for the fourth monoclinic base line. ML4 =2 Max L for the fourth monoclinic base line. MS4 =4 Max Abs(H)+K+L for this line. OBS. The monoclinic base lines are (1,2,3,4) (1,2,3,5) and (1,2,4,5) If SELECT is less than 6. (See SELECT below) MONOSET =0 This parameter makes it possible to use more than 3 base line sets in the monoclinic trials. If MONOSET is: Greater than 3, base line set (1,3,4,5) will be used Greater than 4, base line set (1,2,3,6) will be used Greater than 5, base line set (2,3,4,5) will be used Greater than 6, base line set (1,2,3,7) will be used Thus max 7 base line sets can be used. In N-TREOR (=STREOR), TREOR99 and TREOR90 modes MONOSET is fixed = 7 MONOGAM=1 The best 5 (in N-TREOR 10) trial parameter sets stored (See IQ) for each base line set will be refined before next base line set is tested. If MONOGAM=0 all base line sets are tried before any refinement is made. MONOGAM is only used in monoclinic tests. It is recommended to use MONOGAM=1 because a refined cell parameter set is always tested for the stop limits NIX and MERIT. Thus CPU-time may be saved. MONO =0 Max beta angle allowed in a cell. OBS. No monoclinic test if MONO=0 (See also SHORT) In N-TREOR (=STREOR), TREOR99 and TREOR90 modes MONO is fixed = 135 SHORT =1 Short axis test. The parameter is only used for monoclinic tests. The first six lines are tested for the occurrence of a common zero index in the six first lines. If SHORT=0 no short axis test. If you want to make this test without repeating other monoclinic tests, you may give MONO a negative sign. * * * * * * * * * * * * * * * * * * * * * * * * * * * USE =19 -or equal to the number of input lines if there are less than 19 lines, -or equal to the number of lines with sine square thetas less than 0.327 -USE is the number of lines used in the trial-indexing part of the calculations. OBS. Max USE=20 OBS. If you want to change USE, you should also change IQ. (See IQ). IQ =USE-3 The number of indexable lines required in the trial- indexing procedure if the cell should be stored for ev. least-squares refinement. These reciprocal cell parameters are printed if LIST=1 LIST =0 See IQ above. SELECT =-1 In N-TREOR (=STREOR) and TREOR99 SELECT must be -1 If you want to avoid the TREOR99 mode (but why ?..) set SELECT=0 or... if SELECT is positive the orthorhombic base lines are (SELECT,1,2) (SELECT,1,3) and (SELECT,2,3) or.. if SELECT is greater than 5 the monoclinic base lines are (SELECT,1,2,3) (SELECT,1,2,4) and (SELECT,1,3,4) (Positive SELECT can not be used in TREOR99 mode.) MERIT =10 The De Wolff figure of merit required as stop limit. Ref: De Wolff,P.M. J. Appl. Crystallogr. 1 (1968) 108-113 ( For cubic, tetragonal and hexagonal symmetries are the different quadratic forms as given in Int. Tabl. of X-Ray Crystallogr. (1968) Vol.2 p.109-145 used in the calculation of the number of theoretical lines.) OBS. The figure of merit calculations are not strictly valid unless all 20 first lines are indexed. NIX =1 If a cell after least squares refinemnet has a figure of merit equal to or greater than MERIT and the number of not indexable lines among the USE first lines is less than or equal to NIX, the calculations are stopped. OBS. Otherwise he calculations will end with a difference analysis (Program I1. Werner,P.-E. Z.Kristallogr. 120 (1964) 375-378) (Now not very useful list. --'historical') IDIV =1 The 7 first lines are adjusted by (eventually occurring) higher order lines. If IDIV=0 no corrections. Usually the default value 1 is o.k. There are exeptions, however. If indexing is not successful, you may try IDIV=0 WAVE =1.5405981 Wave length. (in Angstroem) As a rule one should not change WAVE If D-values are used in the input data file (See CHOICE=4) one can always pretend that WAVE was 1.5405981 A. WAVE is then a formal parameter only related to D1, SSQTL and D2 (See below). VOL =2000 Max cell volume (in Angstroem**3) A new option available in TREOR90 is to give a negative value of VOL, ex. VOL=-2000. See comments number 16 on the top of this list. CEM =25 Max cell edge (in Angstroem) The CPU-time depends on VOL and CEM D1 =0.0002 (if WAVE=1.5405981) See D2 below. SSQTL =0.05 (if WAVE=1.5405981) See D2 below. D2 =0.0004 (if WAVE=1.5405981) A line is regarded as indexed if.. sine square theta is less than SSQTL and Abs(sine square theta observed minus sine square theta calculated) is less than D1 or... if sine square theta is greater than SSQTL and the corresponding difference is less than D2. D1, SSQTL and D2 are used in the trial indexing part as well as in the least squares refinements If SELECT=-1 are D1, SSQTL and D2 automatically multiplied by the factor (WAVE/1.5405981)**2 (i.e. D1, SSQTL and D2 depend on WAVE) CHOICE =4 Indicator defining SQ on card set two.. CHOICE=0 SQ=Sine square theta =1 SQ=1/(D*D) (D-spaceing in Angstroem) =2 SQ=Theta (Theta=Bragg angle in deg.) =3 SQ=2*Theta =4 SQ=D OBS next line... Default value in earlier TREOR-versions was 0 DENS =0 Density. (DENS=0 density not used.) If only an integral number of molecules in the unit cell is accepted DENS, EDENS and MOLW may be used. (On your own responsibility) DENS = density in gram per cm**3 EDENS =0 Not used unless DENS equals non zero. EDENS= Max deviation in DENS. OBS. DENS and EDENS are used in trial calculations i.e. they are used on non refined unit cells. Therefore, do not use too small EDENS MOLW =0 Not used unless DENS ( and EDENS ) are non zero. Mol. weight in A.U. (OBS. Crystal water included.) It is not recommended to use DENS, EDENS and MOLW TRIC =0 No triclinic test. If TRIC=1 all higher symmetry tests are omitted and a triclinic test is made. It is presupposed that all higher symmetries have been tried in earlier runs. Although it is in principle possible to index any pattern as triclinic, the indexing algorithm used here is not effective for higher symmetries. OBS. See comment 7 on the top of this list. TRIC is not used in N-TREOR (=STREOR), TREOR99 and TREOR90 modes. ZEROSET =0 This keyword may be given as an integer. ZEROSET*0.005 degrees are added to all theta- values (i.e. ZEROSET*0.01 deg. to all 2theta). Example: If 2theta-zero at 0.05 deg., then set ZEROSET=-5 In N-TREOR different zero point positions are automatically tested as explained under the key-word LIMIT KREDUCT =0 If KREDUCT is given, only trial solutions with IQ=USE-(KREDUCT-1) number of lines indexed may be saved for trial-refinements. OBS. If KREDUCT=0 (or is not given) IQ=USE-3 (See keyword IQ above) These reciprocal cell parameters are printed if LIST=1 Thus, if KREDUCT=1 only trial-solutions where all lines are indexed will be saved for further tests. LIMIT =10 New key-word for program versions STREOR and NTREOR. Max. number of origin shifts. The program will ask for the theta-shift. Normal theta-shift is 0.01 deg. if wave=1.5406 A (i.e. 0.02 deg. in 2 theta) Tests will be made for origin theta-shifts until (LIMIT-1)*(theta-shift) Thus, if LIMIT=1, no origin shift will be tested. If you can expect that zero shift is negligible you may save a lot of computing time if LIMIT=1 is used. END* This keyword denotes the end of the parameter list. (i.e. end of card set three) C O M M E N T S F O R T H E P R O G R A M M E R THE FILES ARE OPENED IN THE MAIN PROGRAM (THE FIRST PROG). THE LOGICAL UNITS ARE.. NUIT=9 THE CONDENSED OUTPUT FILE. IIN=8 THE DATA INPUT FILE. IOUT=7 THE OUTPUT FILE. NDISP=6 OUTPUT (ON DISPLAY) OF TRIAL PARAMETERS IF KEYWORD LIST=1 (SEE KEYWORDS IQ AND LIST) LKEY=5 KEY-BOARD. THE LOGICAL UNIT NUMBERS 5,6,7,8 AND 9 ARE GIVEN IN THE MAIN PROGRAM AND MAY BE CHANGED FOR YOUR COMPUTER. THEY NEED NOT BE CHANGED IN ANY OTHER PLACE OF THE PROGRAM, HOWEVER. THE PROGRAM IS MAINLY WRITTEN IN FORTRAN (II) AND (IV), BUT FORTRAN 77 HAS BEEN USED TO SOME EXTENT. (SEE FOR EXAMPLE SUBROUTINE TWODIM.)-IT IS THE INTENTION, HOWEVER, THAT IT SHOULD NOT BE DIFFICULT TO REWRITE THE FORTRAN 77 STATEMENTS IF ONLY FORTRAN(IV) IS AVAILABLE. It may be recommended to use the WATCOM FORTRAN 77 compiler. VERSION 4 OF THE PROGRAM HAS BEEN DEVELOPED AT STOCKHOLM UNIVERSITY USING A VAX 11/750 COMPUTER. VERSION 5 WAS DEVELOPED FOR CONVEX 210, VAX 11/750 AND IBM PC/AT. VERSION TREOR90 IS WRITTEN FOR CONVEX 210. A NON-VECTORIZED VERSION IS ALSO AVAILABLE. TRICLINIC TESTS MAY BE VERY TIMECONSUMING ON A VAX, HOWEVER. CALLS FROM THE MAIN PROGRAM ARE TO... PWINL.....THE DATA INPUT ROUTINE. TREOB.....THE TRIAL MODULE (THE MOST TIME-CONSUMING PART). TREOC.....PROG. FOR DIFFERENCE ANALYSIS AND ORGANISATION FOR TREOD. TREOD.....LEAST SQUARES REFINEMENTS OF THE BEST TRIAL CELLS.