Radial Correlation Functions
The RDFSOL command computes radially resolved correlation
functions, such as radial distribution functions or orientational
correlation functions. The function of interest is computed either between
pairs of atoms from two atom selections, or between a pair constituted by
one atom selection and a reference point, which can be either a fixed point
in space or the center of mass of a set of atoms. As the name rdfSOL
suggests, the routine allows for special treatment of solvent molecules
(currently only support for TIP3 is implemented).
* Menu:
* Syntax:: The syntax of the RDFSOL command
* General:: General overview
* Sets:: Set selections
* Limit:: Locally limiting sets in each frame
* Options:: other general options
* Traj:: Trajectory specifications
* Caveats:: Some limitations/todos to keep in mind
* Examples:: Just what it says
Syntax for the RDFSOL command
[SYNTAX RadialDistributionFunctions_with_SOLvent]
Syntax:
RDFSOL [ RDF ] [ setA-spec ] [ setB-spec ] [ around-spec ] -
[ DDIP ]
[ QDIP ]
[ HD ]
[ SAME ] [ RMAX real ] [ NBIN int ] [ IOUT int ] -
[ VOLUme real ] [ PRECise ] [ BRUTe ] [ MINI ] -
[ traj-spec ] [ SPFAc int]
setA-spec:: [ SITE ] [ atom-selection ]
[ XREF real ] [ YREF real ] [ ZREF real ]
[ WATEr ]
setB-spec:: [ SITE ] [ atom-selection ]
[ WATEr ]
around-spec:: AROUnd [ RAROund real ] [ LOCAl ] -
[ atom-selection ]
[ XREF real ] [ YREF real ] [ ZREF real ]
traj-spec:: [ FIRStu int ] [ NUNIt int ] [ BEGIn int ] -
[ STOP int ] [ SKIP int ]
atom-selection::= see *note Selection:(select.doc)
General overview and options
RDFSOL calculates radially resolved pair-distribution or angular
correlation functions between two sets of atoms (setA and setB).
The type of function computed is selected by keyword:
RDF - Radial Distribution Function
If one of the two sets is WATEr then distribution functions
for the oxygen and hydrogens are computed.
If both sets are WATEr then the hydrogen-hydrogen distribution
function will also be evaluated.
If the first set is not water and the keyword SITE is present,
the center of mass of the set is taken as single center. Else
the average over all points in setA is taken.
DDIP - Dipole-Dipole correlation function. If one or both sets are
not WATEr the center of mass and dipole moment of this set
is used (no matter whether the keyword SITE is present or
not). In this case setA must not be a fixed point in space
since the dipole moment is not defined in this case.
QDIP - Charge-Dipole correlation function. As with RDF setA can
either be a SITE or the average of all points.
SiteA/B specifications
SetA and SetB are two sets for which the chosen function is evaluated
for all pairs A-B. Both can be WATEr, in which all TIP3 residues present
will be included. In this case the oxygen positions will be used as the
centers of the molecules. In both cases the center of mass (and set
dipole moment if needed) can be used if the keyword SITE is present.
SetA can be a fixed point in space: (XREF/YREF/ZREF) (if SITE is present
but no atom selection (0/0/0) will be used as default).
For both sets WATEr is the default.
Limiting Sets
If only a subset which is localized around a certain point should be
used in each frame this can be achieved by the AROUnd keyword. If it is
present setA will be re-selected in each frame. If the keyword LOCAl is
present setB will also be re-selected. RAROund <real> is the radius
around the selected center within which an atom must lie to be available
for evaluation in this frame. The center itself can either be a fixed
point in space (XREF/YREF/ZREF) or an atom selection of which the center
of mass will be used.
Other options
----------------------------------------------------------------------
SAME - if this keyword is present, only setA is used for both
sets, thus calculating auto-functions (this algorithm
should be faster than the general one if setA and setB
use the same selection)
----------------------------------------------------------------------
RMAX <real> - the maximum distance up to which a pair A-B is evaluated
(default: 7.5A)
----------------------------------------------------------------------
NBIN <int> - the number of bins used to sample (each bin is RMAX/NBIN
wide)
(default: 150)
----------------------------------------------------------------------
IOUT <int> - unit number the function(s) should be written to
(default: stddout)
----------------------------------------------------------------------
VOLUme <real> - the volume of the total system, necessary for the
normalization. If not specified by the user and crystal
is in use, the resulting cell volume will be used.
Finally, if crystal is not used and no volume
is specified, and if both sets are localized (see AROUnd),
the volume of the limiting sphere will be used.
----------------------------------------------------------------------
PRECise - if RDFs are calculated and one or both sets contain
WATEr, some pairs including water hydrogens will be
missed since only oxygen distances are evaluated.
If PRECise is present, these pairs are also included
which results in a slightly diminished efficiency of the
cubing algorithm
----------------------------------------------------------------------
BRUTe - use a simple double loop algorithm rather than a cubing
algorithm
----------------------------------------------------------------------
MINI - use 'real' minimum image conventions. Only used if both
sets are water and useful only if distances farther than
the minimum image distance (when pairs could be counted
twice) are evaluated. Currently this is needed only if
the distance dependent Kirkwood G-factor is computed
(with DDIP, second column) as one needs to go 'into the
corners' (i.e. sqrt(3)/2 * L for a cubic box).
(caution: needs lots of memory)
----------------------------------------------------------------------
SPFAc - if images are present, the number of total
atoms/pairs/cubes may change from frame to frame.
So an estimate of the needed space needs to be made
before reading the trajectory so SPFAc times the actual
values is allocated
(default: 3)
----------------------------------------------------------------------
Trajectory specifications
These are the usual specs. The trajectory is read NUNIt units starting
with FIRSTu reading from frame BEGIn to STOP where SKIP frames are
skipped between reading.
Caveats and Comments
- When computing dipole-dipole correlations for a set which is not WATEr,
only its center of mass and dipole moment for primary atoms will be
evaluated. So if a part of a large molecule which is re-centered
bysegment (e.g. a protein) and "sticks out" of the primary box, some
pairs may not be sampled.
- Normalization of RDFs differs slightly from that used in COOR ANAL
- no excluded volume correction
- point-point (e.g. two SITEs or DDIP for two non-WATEr sets...) not yet
implemented
Examples
(See also test/c30test/rdfsol.inp test/c30test/rdfsol2.inp testcases)
----------------------------------------------------------------------
RDFSOL RDF SETA WATER SAME RMAX 7.5 NBIN 150 IOUT 10 PRECISE -
FIRSTUNIT 11 NUNIT 1
This will calculate g_OO, g_OH and g_HH for all waters in the simulated
system up to 7.5 A into 150 bins. One trajectory file will be read from
unit 11 and the result output to unit 10.
----------------------------------------------------------------------
RDFSOL RDF SETA WATER SAME RMAX 7.5 NBIN 150 IOUT 10 PRECISE -
AROUND RAROUND 7.5 LOCAL SELECT ATOM PROT 1 NH END -
FIRSTUNIT 11 NUNIT 1
The same as above but only waters around the NH of residue 1 of segment
PROT will be considered.
----------------------------------------------------------------------
NIH/DCRT/Laboratory for Structural Biology
FDA/CBER/OVRR Biophysics Laboratory
Modified, updated and generalized by C.L. Brooks, III
The Scripps Research Institute