MTZFIX

NAME

mtzfix - Fixes up Fourier map coefficients.

SYNOPSIS

mtzfix  [FLABEL <string>]  HKLIN in.mtz  HKLOUT out.mtz

DESCRIPTION

The program MTZFIX 'fixes up' the moduli of the map coefficients in an MTZ file produced by a structure-factor calculation or refinement program (typically Buster, Refmac or Sfall/Sigmaa), according to the formulae in the literature (Main, 1979; Read, 1986).  The program is able to detect which program originally created the file purely from the column contents (although this information can also obviously be deduced from the column labels!).  The table below compares the coefficients output by the programs Buster, Refmac & Sigmaa at the time of writing:

SourceCoefficientAcentricsCentrics
LiteratureFM2mFo - DFc mFo
FD2(mFo - DFc) mFo - DFc
FCDFcDFc
BUSTERFM2mFo - DFc 2mFo - DFc
FD2(mFo - DFc) 2(mFo - DFc)
FCDFcDFc
REFMACFM2mFo - DFc 2mFo - DFc
FD mFo - DFc mFo - DFc
FCDFcDFc
SIGMAAFM2mFo - DFc mFo
FD2(mFo - DFc) mFo - DFc
FCDFcDFc

where the 'coefficient' is:

FM for the '2Fo-Fc' map, FD for the 'Fo-Fc' map and FC for the 'Fc' map.

Differences from the literature values are indicated in red.  Note that for Refmac, if the purpose is only to compute the map, halving the acentric FD coefficient has the same effect as doubling the centric one (so Buster & Refmac are consistent in this respect).  However if the statistical analysis program EDSTATS is to be used it is essential to correct this inconsistency, otherwise incorrect map coefficients will be computed internally by EDSTATS.

The consistency checks performed by this program are (A=acentric, C=centric, else both):

A1FC = 2mFo - FM
A2FC >= 2mFo - FM
A3FD = FM - mFo
A4FD = 2(FM - mFo)
C5FC = 2mFo - FM
C6FM = mFo
C7FD = mFo-FC
C8FD = 2(mFo-FC)
C9FD <= mFo-FC
10FM = FC (unobserved only)
11FD = 0 (unobserved only)

The program makes a first pass through the input MTZ file to perform the above checks and work out what needs to be fixed.  Then if required, a second pass makes the necessary corrections, writes out a new MTZ file and repeats all the checks on the output MTZ file.  It is essential to check the log file to ensure that the output file passes all the tests; if it doesn't the faulty columns indicated by the program should not be used.

INPUT

Optional command line input: FLABEL <string>.  This specifies the column labels for Fo and sigma(Fo) in case they are not the assumed defaults.  The file label for sigma(Fo) is assumed to be 'SIG' followed by the Fo label specified; e.g. if the file labels for Fo and sigma(Fo) are 'FP' and 'SIGFP' then specify 'FLABEL FP'.

The assumed default file labels are as follows:

For Fo & sigma(Fo): 'FOSC' & 'SIGFOSC' (Buster), or 'F' & 'SIGF' (Refmac & Sigmaa).
The Fo and sigma(Fo) file labels may also be 'F_' and 'SIGF_' followed by the same arbitrary string (up to 24 characters).  In this case it's not necessary to specify an FLABEL string.

For Fc & phic: 'FC_ALL' & 'PHIC_ALL' (Refmac only), or 'FC' & 'PHIC' (all).  However, note that use of the standard FC and PHIC columns from Refmac will cause most of the tests to fail; it is essential that the program uses the default 'FC_ALL' & 'PHIC_ALL' columns (i.e. containing the total contribution to Fc including that from the bulk solvent).

For '2Fo-Fc' and associated phase: '2FOFCWT' & 'PH2FOFCWT' (Buster), or 'FWT' & 'PHWT' (Refmac & Sigmaa).

For 'Fo-Fc' and associated phase: 'FOFCWT' & 'PHFOFCWT' (Buster), or 'DELFWT' & 'PHDELWT' (Refmac & Sigmaa).

For FOM: 'FOM' (Buster & Refmac) or 'WCMB' (Sigmaa).

For free R flag: 'FREE' or 'FreeR_flag'.

INPUT FILES

MTZ file containing standard columns of Fo, sigma(Fo), Fc, phic, '2Fo-Fc' & 'Fo-Fc' map coefficients (moduli & phases), FOM and free R flag.  All columns with the exception of the free R flag are required.

OUTPUT FILES

'Fixed' MTZ file containing the same columns as input, with the map coefficients 'fixed'.  If the MTZ file has already been 'fixed' or 'fixing' is not necessary, then no MTZ output file is written.

REFERENCES

Main, P., Acta Cryst. (1979). A35, 779-785.
Read, R.J., Acta Cryst. (1986). A42, 140-149.

AUTHOR

EXAMPLES

This example assumes default labels (e.g. F/SIGF or FOSC/SIGFOSC for Fo):
mtzfix  HKLIN input.mtz  HKLOUT output.mtz
This example assumes labels FP/SIGFP for Fo & defaults for everything else:
mtzfix  FLABEL FP  HKLIN input.mtz  HKLOUT output.mtz