ANALYZE_MIR: Analyzing an MIR dataset
ANALYZE_MIR is a routine to calculate difference Pattersons and anomalous difference Pattersons for an MIR dataset, to calculate correlation coefficients among them all, and to set up a "solve_mir.script" file that can be used to run SOLVE on the dataset. You usually do not have to worry about this routine at all because it is ordinarily called right after running SCALE_MIR for automatic structure determination. It is ordinarily followed by running SOLVE with the solve_mir.script file written out by this routine or by using the keyword "SOLVE" after running this routine.
ANALYZE_MIR assumes that you have a datafile ("mir_fbar.scl") that contains Fnat, sigma, and (Fbar,sigma,DelAno,sigma) for each derivative of MIR data. That is, there are exactly 2 columns of data for the native and 4 data columns for each derivative. This is a dorgbn file. It is ordinarily created by SCALE_MIR, but you can create your own if you like.
To run ANALYZE_MIR, you need to input your standard setup file ("solve.setup"), and the name of the mir datafile. ANALYZE_MIR will calculate origin-removed difference Pattersons for all isomorphous and anomalous differences. These maps are all compared to each other and the correlation coefficients are displayed in a table. In a typical MIR experiment the anomalous and isomorphous difference Pattersons have correlations with each other on the order of 0.0 to 0.2 or so (pretty low, so don't be worried).
ANALYZE_MIR writes out a script file "solve_mir.script" that you can use to go on with the SOLVE routine to solve this mir structure. You can edit this script file to modify it if you like or you can run it as is. You can also go right on with SOLVE by adding the command "SOLVE" after "ANALYZE_MIR".
Please note that there is a slight difference between running ANALYZE_MIR right after SCALE_MIR and in a separate session. If you run them right after each other, you should input all the information about each derivative at once, as in the sample script for automatic analysis of MIR data. (If you wanted to input information about derivative #1 in two places in your script file, you would need to use the keyword "gotoderiv 1" the second time. You can avoid this by putting all the information for derivative #1 together.)
A typical input script for ANALYZE_MIR follows:
!--------------------------Run ANALYZE_MIR---------------------------------
@solve.setup ! standard information for this crystal
mirfbarfile mir_fbar.scl ! input dorgbn file with Fnat,sig, and
! (Fbar,sig,Delano,Sig) for each wavelength
logfile analyze_mir.logfile ! write out most information to this file
derivative 1 ! derivative #1 information is to follow
label deriv #1 Hg ! label for deriv 1
atomname HG ! this is an HG derivative
derivative 2
label deriv #2 Iodine
atomname I-
ANALYZE_MIR !analyze all the Pattersons and
!write solve_mir.script
!-----------------------------------------------------------------------------
Keywords for ANALYZE_MIR:
LOGFILE xx.logfile log file for output will be xx.logfile
mirfbarfile xx.scl input dorgbn file with Fnat, sig, and
(Fbar,sig,Delano,Sig)
for each wavelength will be xx.scl
[Default="mir_fbar.scl"]
SCRIPTFILE xxx Output script file containing instructions
for running SOLVE written to xxx
[default="solve_mir.script"]. Starting SOLVE with
this script is equivalent to following
ANALYZE_MIR with SOLVE.
derivative n derivative number (n=1,2,3...) for values
to be entered next of atomname and label
label xxxx xxx is a label for this derivative
ATOMNAME XXXX XXXX is the atom name of the atoms in
this derivative. This name can be in
SOLVE's database or you can
enter it using the NEWATOMTYPE keyword.
PLEASE NOTE: the f' and f" values assumed
by SOLVE are for Cu Kalpha radiation (1.54 A).
If you collected your MIR data at a synchrotron
then you will want to define a new atom type
with the appropriate values of f' (fprimv) and
f" (fprprv). See the HEAVY or SOLVE keyword
lists for entering a NEWATOMTYPE.
The atom types recognized by SOLVE are:
H, H-1, He, Li, Li+1, Be, Be+2, B, C, Cv, N, O, O-1,
F, F-1, Ne, Na, Na+1, Mg, Mg+2, Al, Al+3, Si, Siv, Si+4,
P, S, Cl, Cl-1, Ar, K, K+1, Ca, Ca+2, Sc, Sc+3, Ti, Ti+2,
Ti+3, Ti+4, V, V+2, V+3, V+5, Cr, Cr+2, Cr+3, Mn, Mn+2, Mn+3,
Mn+4, Fe, Fe+2, Fe+3, Co, Co+2, Co+3, Ni, Ni+2, Ni+3, Cu,
Cu+1, Cu+2, Zn, Zn+2, Ga, Ga+3, Ge, Ge+4, As, Se, Br,
Br-1, Kr, Rb, Rb+1, Sr, Sr+2, Y, Y+3, Zr, Zr+4, Nb, Nb+3,
Nb+5, Mo, Mo+3, Mo+5, Mo+6, Tc, Ru, Ru+3, Ru+4, Rh, Rh+3,
Rh+4, Pd, Pd+2, Pd+4, Ag, Ag+1, Ag+2, Cd, Cd+2, In, In+3,
Sn, Sn+2, Sn+4, Sb, Sb+3, Sb+5, Te, I, I-1, Xe, Cs, Cs+1,
Ba, Ba+2, La, La+3, Ce, Ce+3, Ce+4, Pr, Pr+3, Pr+4, Nd,
Nd+3, Pm, Pm+3, Sm, Sm+3, Eu, Eu+2, Eu+3, Gd, Gd+3, Tb,
Tb+3, Dy, Dy+3, Ho, Ho+3, Er, Er+3, Tm, Tm+3, Yb, Yb+2,
Yb+3, Lu, Lu+3, Hf, Hf+4, Ta, Ta+5, W, W+6, Re, Os, Os+4,
Ir, Ir+3, Ir+4, Pt, Pt+2, Pt+4, Au, Au+1, Au+3, Hg, Hg+1,
Hg+2, Tl, Tl+1, Tl+3, Pb, Pb+2, Pb+4, Bi, Bi+3, Bi+5, Po,
At, Rn, Fr, Ra, Ra+2, Ac, Ac+3, Th, Th+4, Pa, U, U+3, U+4,
U+6, Np, Np+3, Np+4, Np+6, Pu, Pu+3, Pu+4, Pu+6, Am, Cm, Bk,
Cf
NSHELLS n Number of shells for analysis is n
[default=10]