MIR: Multiple Isomorphous Refinement
datset isomorphs to be included
limset heavy atoms to be included
datscl scale isomorph data to parent data
ref parameters to be refined
noref turn off parameter refinement
bins statistics grouping
stat statistics to be calculated and printed
nostat turn off statistics to be calculated and printed
mvfc1MVFC phase determination
fazinintroduce known phases (currently for
MVFC method only)
fazprb1modified Blow-Crick phase
determination
eiso isomorphous closure errors
eano anomalous closure errors
plthkl probability profile plotting
newtonadjustment of most probable phase by
Newton's method
cycle number of refinement cycles
1 If neither appears no phase determination will be
made.
MIR 1 phase refinement control
block block diagonal
full full matrix Bricogne method (not
implemented)
none no phase reference
2 list frequency on last cycle no listing
3 reflection listing control basic items
native basic + expanded list for parent, only
all basic + expanded list for all data sets
4 output fazprb closure errors on last cycle no
closure list closure errors
datset 1 data set name (12 characters)
2 data used for phasing
both use both isomorphs and anomalous
isom use isomorphs only
anom use anomalous only
none do not phase
3 data used for refinement (same controls as field 2)
atoms 1 data set name
2- atom labels of atoms to be used
limset 1 dataset name of isomorph
2 min sinθ/λ for hkl of this isomorph 0
3 max sinθ/λ for hkl of this isomorph bdf or
2
datscl 1 dataset name of isomorph
2 value to scale isomorph to parent bdf or 1
3 value of ΔB to scale isomorph to parent bdf or 0
4 observational errors (σ's) for each reflection
0.05
If an isomorph does not have these available, the program calculates them
as follows:
σ2= (1 +
C*F)2
where C is analogous to the instrumental instability factor used when
calculating observational errors from diffractometer data.
ref 1 atom label, dataset name or 'all' as
global signal
noref 2- referenced parameter as mnemonic
x y z atomic
coordinates
u u11 u22 u33 u12
u13 u23 thermal parameters
pp ppa pop. param
(real+anomalous)
k c overall scale and
ΔU
bins 1 number of statistics bins (max.
12)1 5
2 number of statistics bins (max. 12)1
3
3 statistics listing control for each cycle except last2
all list all categories of
statistics3
centric list centric statistics only
acentric list acentric statistics only
4 statistics listing control for last cycle3 (same
controls as in field 3)
centric
1 The total number of bins is the product of the number of
ranges and the bins. The total may be as few as 1
bin or as many as 50. To produce reasonable statistics a minimum of 100
reflections should appear in each bin and preferably several times that
number.
2 Statistics output can be quite copious. In addition to the
stat and nostat lines below, fields 3 and 4
offer an opportunity to limit the statistic listing to centric,
acentric or all.
3 As the isomorphous closure errors in
fazprb mode are updated from the centric values (if there were
any centric measurements), these are also printed under the all
option.
stat codes defining statistical analyses (see description
in Ref. Manual)
nostat fomhist figure-of-merit histogram
fom figure-of-merit by and
rcullis Cullis R-factor
rkraut Kraut R-factor
rfazpwr R phasing power; weighted closure error over
weighted heavy atom model
wteiso RMS weighted isomorphous closure error (on F)
wteano RMS weighted anomalous closure error (on F)
number number of reflections
fazdif phase differences between parent & heavy atom
model
all all statistics
lstcyc listed codes only apply to last cycle
If no fields follow this mnemonic, all statistics which have not
already been requested are printed on the last cycle. Only those codes listed
on the lstcyc line are printed on the last cycle.
for FAZPRB mode only:
parvar parent variance
eiso RMS isomorphous closure error
eisoadj RMS isomorph. closure errors adjusted for
variance
isodif RMS isomorphous differences
hamodel RMS heavy atom model
jntvar RMS joint variance of parent and isomorph
ermsano RMS anomalous closure errors
eanoadj RMS anom. closure errors adjusted for variances
anodif RMS anomalous differences
signs agreement of obs. and calc. anomalous diffrences
The statistical output can be copious. To save trees but not the movement
of electrons, the stat and nostat lines set the
print priority of the desired statistics to 3 and all other statistics to 5
(i.e. only output if reset pstat 5 is entered before MIR). A
very small group of statistics always has priority 2 for those who wish to
monitor the program while it is being run in the "interactive" mode. These
statistics include the average phase change from the previous cycle, a summary
providing the overall figure of merit, and for each isomorph the overall R
Cullis, R Kraut and R phasing power. In fazprb mode,
ermsiso and ermsano are also calculated.
The default is to print to all statistics on the last cycle and the
following statistics on all cycles: fom,
rcullis, rkraut, rfazpwr and in
fazprb mode ermsiso, eisoadj,
ermsano and eanoadj. Statistics relating to
anomalous information are printed only if there were observations.
mvfc 1 phase increment (max. 180o) 10
2 weak reflection cutoff ratio F/γ(δ) 2.5
3 phase refinement criterion selection 1
4 shift factor for phases (max. 10) 1
5 phase threshold to accept bimodal SIR soln. (max. 179o)
179
6 goodness-of-fit to initialize FOM bdf or
20
7 signal to not use bdf phases, enter newphases
use bdf
fazin Option Code
initial phase set init n 1
restraining phase set comb n
1
fazprb 1 phase incr. for integrating refln. probability profile
5
(min. 1o; max. 180o; recommended max.
12o)
2 signal to control refinement and statistics
most use the most probable phases
best use the 'best' phases
eiso and eano (for fazprb mode
only)
1 dataset ID
2 number of bins reserved for closure errors NO default
0 implies all bins used
3- closure errors
The number of closure errors must either be 1 or equal to the current
number of ranges (see bins line). If only 1
closure error is encountered, all ranges within the indicated
bins (field 2) are replaced with that value. Otherwise there is
corresponding replacement. The usual default value for the closure errors are
those on the bdf. If the bdf has not been used by MIR, default for isomorphous
closure errors is 1000 and anomalous closure errors is 592. If the number of
ranges or bins is changed (see bins
line) and there are closure errors on the binary data file, new closure errors
are expected to be supplied from the input stream.
plthkl (for fazprb mode only)
1-3 hkl indices of the first of a group of reflections to be plotted
4 number of reflections in group to be plotted (max. 100)
10
This line provides a facility for line printer plots of the phase
probability profile for a group of reflections. This operates only on the last
cycle of the program and prints at priority 4, and is intended primarily as a
teaching tool to study the effects of various parameters on the phase
distribution of selected reflections. Each plot takes a page of output so this
feature should be used sparingly. The program will not perform this operation
when the listing of every nth reflection has been requested on the
MIR line.
newton 1 centroid minimum level for centro reflections
0.5
2 centroid minimum level for noncentro reflections 0.8
cycle 1 number of refinement cycles (max. 10) 1