[CCP14 Home: (Frames | No Frames)]
CCP14 Mirrors: [UK] | [CA] | [US] | [AU]

(This Webpage Page in No Frames Mode)

CCP14

Tutorials and Examples

CRYSFIRE Powder Indexing System for DOS/Windows by Robin Shirley et al

For those who "work on the edge of the abyss of impossibility"

Using CRYSFIRE Powder Indexing System with less than 20 peaks (15 peaks)

The CCP14 Homepage is at http://www.ccp14.ac.uk

[Tutorials page] | [CRYSFIRE Tutorial Index] | [Crysfire What's New] | [LMGP Suite Chekcell Index] | [Powder Diffraction Indexing Resources]

Remember: "99 out of a 100 people who do indexing on less than 20 peaks are chasing rainbows"

This run through will use the peak list file obtained by using the corklz example file provided with Stefan Krumm's Winfit for Windows software. The peaks were fitted with Winfit and saved inside the winfit.dat peak result file. (Winfit is such an easy software package to use, no formal tutorial is required).


Obtaining and installing CRYSFIRE

NOTE: A graphical program for MS-Windows that links into CRYSFIRE to help check trial cells and determine spacegroups is the Chekcell software - part of the LMGP suite for Windows by Jean Laugier and Bernard Bochu.

Indexing strategy of indexing with less than 20 peaks

  1. First index on the native 15 peaks using the programs that work with less than 20 peaks: Dicvol, Taup, Kohl and Treor

  2. Then use Crysfire's Pattern Extension to extend the number of peaks out to 20, then run All the programs again.

    As noted by Robin:

    From: Robin Shirley [[email protected]]
    Date: Wed, 10 Jul 2002 16:31:23 GMT
    
    Jon and fellow Rietvelders/sdpd-ers
    
    > Does this just mean dividing the observed d-spacings by 2 (or 3,4..)
    > to get more "observed" values, or am I missing a more cunning trick?
    
    In essence that's exactly what lies behind it, though perhaps a bit more
    systematically than most humans would consider worth grinding through.
    
    For example, it works iteratively, making overlap-avoidance tests at each
    step, and produces the minimum extension set using the lowest available
    orders.
    
    But the bottom line is that no more information is injected into the
    pattern by this process - it would be impossible for things to be
    otherwise, since to win more degrees of freedom you have to make more
    actual observations (which is sometimes not experimentally possible).
    
    At least it would be impossible to conjure up more information (i.e.
    degrees of freedom) just through working from a fixed set of observed
    peak positions - in principle there might be more information to be
    found buried beneath the profile, though in practice it is usually very
    low-grade.
    
    In fact, not only does pattern extension *not* add more information, but
    if the original data contained systematic errors, then the extra
    higher-order lines will not reproduce these in a self-consistent way -
    all available correction and calibration should always be done *before*
    pattern extension.
    
    Having said all that, it is still a great convenience to be able to
    extend unavoidably sparse patterns from High-P/High-T experiments in a
    few seconds to give a format that none of the indexing programs will
    object to.  For good (though sparse) data, this can open the door to
    getting cell (or at least sub-cell) information that can otherwise be
    hard to lock onto.
    
    The authentic pattern that Lachlan supplied me for testing these issues
    (yes, it was Lachlan's turbo-charged lobbying on behalf of his
    High-P/High-T clientele that lay behind this) originally contained only
    7 lines.  As you can imagine, it was tough to interpret in that form,
    but becomes significantly more stable after pattern-extension, yielding
    some quite interesting possible solutions (or partial solutions).
    
    But it's certainly a case of caveat emptor, and is intended specifically
    for folk who know what they're doing.
    
    I'll try to put together a short paper on this when a bit more experience
    has been gained on how it usefully it performs.
    
    Robin Shirley
    

Note: The following is just a single example of using CRYSFIRE. There are other methods of interacting with the program and doing the following.


The following run-through shows how to import a WinFIT DAT peak listing file into Crysfire. However, extra example tutorials for importing data into Crysfire are:


A comment from Robin Shirley about filenames and descriptive text

From: Robin Shirley [[email protected]]
Organization: Psychology Dept, Surrey Univ. U.K.
To: Lachlan Cranswick [[email protected]]
Date: Mon, 22 Jul 2002 18:43:55 GMT

Incidentally, in this context, could you please emphasise in your
tutorials the importance of spending the few seconds that it takes to
give each new dataset a well-chosen name and some helpful brief
descriptive text in its description field, seeing that whatever is used
there will remain the default throughout the rest of the analysis?

Careful dataset naming really is important for keeping track of progress
in any serious study that can involve several dataset variants, and
such distinctive names have been assumed throughout Crysfire as the basis 
of its data organisation.

Users shouldn't rely just on having kept each problem separate within
its own data directory (folder).

While that's important, it's *not* sufficient, because (a) doing that
only identifies the directory, not the files within it, and (b) as the
study progresses, it's likely that several datasets and dataset variants
will be required, for example after recalibration, applying estimated
Z2theta or specimen-displacement corrections, rescaling, etc.

These will all need to be kept in the same directory so that, for
example, their trial cells can easily be loaded and examined (LC, M1,
etc), but if they aren't given different names (and preferably also
different description fields - see below) then complete confusion will
quickly result.

Since the description field gets appended to every summary-file solution 
line, an important opportunity is lost if its contents can't act as an 
aide memoire of the characteristics of the particular dataset or dataset
variant that was used for that trial solution.

This becomes particularly relevant when datasets are merged and/or
rescaling/unscaling is used, as was discussed above and as is
increasingly likely to happen in response to the new volume estimates
and rescaling prompts in CF2002.

I'm already regretting that I short-sightedly left it far too easy in
WF2crys and XF2crys for users to default to an undistinctive and
meaningless dataset name like Crys.cdt (I've noticed that you are 
yourself a frequent offender in this regard!), and I may well look at 
ways to prevent this in the next revisions of those programs.

As things stand, it's too easy for a moment's impatience on the part of
the user at this basic choice point to handicap all subsequent work on
that sample with a basically null name (and description) - I really
should have had the foresight to protect users from this pitfall.

(A bit of a rant, perhaps, but I do think the point needs making.)

With best wishes

Robin


From (in this case) the C:\crysfire\examples\15peak directory, run Crysfire.

Then type IM to go into the IMport area; followed by WF to import a Winfit file, type in the filename (default is winfit.dat), then answer the questions such as the value of the wavelength. When prompted, load the newly created Crysfire CDT file using the LO command. Crysfire will warn that with less than 20 peaks, not all the indexing programs will run - though it does advise the EP Pattern Extension command is available.

Be wary that if you are low on DOS RAM - the last few peaks may not get imported into the CDT file. Check this out after running the importing routine.


Indexing on the 15 peaks

Run all the relevant Crysfire based indexing programs as per the Crysfire Basic tutorial (via the IN (Indexing) command)

Due to there only being 15 peaks, only the following programs can be run.

  • Taup
  • Dicvol - high symmetry
  • Treor
  • Dicvol - low symmetry

Crysfire will explain this in the IN Indexing menu option, as showing in the following screen dump.

Crysfire summary file on the 15 peaks. - Taup, Dicvol - high symmetry, Treor, Dicvol - low symmetry


Performing Pattern Extension in Crysfire

To extend the pattern to 20 peaks, run the EP (Extend Pattern) command.

Crysfire recommends and prompts the saving of the modified data to a CDT file with EP to the filename (e.g., in this case, crys1ep.CDT).

Saving the pattern extended data


CDT file BEFORE pattern extension

""CRYS v9.30c dataset, saved as crys1 by WF2CRYS v1.7b on 27Jul02 at 21:21:40"
""
"To use with CRYS versions prior to v9.30, delete these first 3 lines"
"crys1",15,0,0,0,0,1.54056,0,1,"P"
0,0,0,0,23.077446,25.4587
0,0,0,0,25.578506,54.818542
0,0,0,0,29.387346,425.113647
0,0,0,0,35.154282,135.198227
0,0,0,0,35.955048,52.457371
0,0,0,0,37.777325,58.142704
0,0,0,0,39.366253,82.505386
0,0,0,0,43.118942,34.948811
0,0,0,0,43.352829,259.191711
0,0,0,0,47.085106,46.223732
0,0,0,0,47.502411,149.731659
0,0,0,0,48.560741,137.30661
0,0,0,0,52.561459,101.193619
0,0,0,0,56.56464,24.265675
0,0,0,0,57.508236,313.949036


CDT file AFTER pattern extension

"Crysfire 9.45.15 dataset, saved as crys1EP on 27Jul02 at 21:30:15"
",EP20"
"To use with CRYS versions prior to v9.30, delete these first 3 lines"
"crys1EP",20,0,0,0.0000000000000,0,1.54056000,0.0000000000000,0,"P"
  674.3630135640981,  0,  0,  0, 23.0774460000000, 25.45800000
  825.8914333229660,  0,  0,  0, 25.5785060000000, 54.81854000
 1084.3671725247637,  0,  0,  0, 29.3873460000000,425.11364000
 1537.0452911027076,  0,  0,  0, 35.1542820000000,135.19822000
 1605.5288895572519,  0,  0,  0, 35.9550480000000, 52.45737000
 1766.3192315616914,  0,  0,  0, 37.7773250000000, 58.14270000
 1912.0238138513860,  0,  0,  0, 39.3662530000000, 82.50538000
 2275.8251900599657,  0,  0,  0, 43.1189420000000, 34.94881000
 2299.3891744257908,  0,  0,  0, 43.3528290000000,259.19171000
 2688.9592859349787,  0,  0,  0, 47.0851060000000, 46.22373000
 2734.0613033657769,  0,  0,  0, 47.5024110000000,149.73165000
 2849.7934223408283,  0,  0,  0, 48.5607410000000,137.30660000
 3304.1399922362721,  0,  0,  0, 52.5614590000000,101.19361000
 3783.7558567225660,  0,  0,  0, 56.5646400000000, 24.26567000
 3900.1957912623138,  0,  0,  0, 57.5082360000000,313.94903000
 4337.4686900990546,  0,  0,  0, 60.9686023975848,  0
 6069.2671220768825,  0,  0,  0, 73.7527582619853,  0
 6148.1811644108302,  0,  0,  0, 74.3108366474561,  0
 6422.1155582290076,  0,  0,  0, 76.2366863689592,  0
 7065.2769262467655,  0,  0,  0, 80.7007999910034,  0


Indexing on the 20 Pattern Extended peaks

Now run all the Crysfire based indexing programs as per the Crysfire Basic tutorial (via the IN (Indexing) command)

  • Taup
  • Dicvol - high symmetry
  • ITO
  • Fjzn
  • Treor
  • Kohl
  • Dicvol - low symmetry
  • Lzon

Crysfire summary file on the extended peaks. - Taup, Dicvol - high symmetry, ITO, Fjzn, Treor, Kohl, Dicvol - low symmetry, Lzon


As per the example shown in the Crysfire Basic tutorial, it would be good to now run MMap and evaluate the trial solutions for both the raw data, and extended data using Chekcell to obtain the best cell and spacegroup. However:


While using the various Indexing programs incorporated inside CRYSFIRE may give reasonable results for the above data. It is not difficult to get reasonable results that are also wrong. For instance, if you were a bit too enthusiastic in assuming you had a new single phase, you may have forgotten to do a search match to check if the following could have been known phases. For example the above data gives a very nice fit with being a mixture of corundum (Aluminium Oxide) and calcite (Calcium Carbonate)

The peak labelled as being a 3.354 Angstrom peak is actually a Copper K beta (1.39217 Angstrom) satellite of the main Calcite peak.


[Tutorials page] | [CRYSFIRE Tutorial Index] | [Crysfire What's New] | [LMGP Suite Chekcell Index] | [Powder Diffraction Indexing Resources]

[CCP14 Home: (Frames | No Frames)]
CCP14 Mirrors: [UK] | [CA] | [US] | [AU]

(This Webpage Page in No Frames Mode)

If you have any queries or comments, please feel free to contact the CCP14