CHARMM c42b2 rdfsol.doc



File: RDFSol -=- Node: Top
Up: (commands.doc) -=- Next: Syntax



                         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 consisting of
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 the special treatment of solvent molecules
(TIP3 water is supported by special routines, others need the use of the
PROTotype facility (see see *note Selection:(proto.doc))). 


* 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



File: RDFSol -=- Node: Syntax
Previous: Top -=- Up: Top -=- Next: General


                         Syntax for the RDFSOL command

[SYNTAX RadialDistributionFunctions_with_SOLvent]

Syntax:

RDFSOL  [ RDF int ]  [ DDIP int ]  [ QDIP int ]  [ HD int ] -

                  [ setA-spec ]  [ setB-spec ]  [ around-spec ] -

                  [ SAME ]  [ RMAX real ]  [ NBIN int ] -

                  [ VOLUme real ]  [ PRECise ]  [ BRUTe ]  [ MINI ]  -

                  [ traj-spec ]  [ SPFAc int]



      setA-spec:: [ SITE ]  [ atom-selection ]
                  [ XREF real ]  [ YREF real ]  [ ZREF real ]
                  [ WATEr ]
                  [ PROTo int ]


      setB-spec:: [ SITE ]  [ atom-selection ]
                  [ WATEr ]               
                  [ PROTo int ]


      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)



File: RDFSol -=- Node: General
Previous: Syntax -=- Up: Top -=- Next: Sets


                         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 keywords:

  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.

The integer after each function to be calculated gives the unit number
the respective function is to be written to.



File: RDFSol -=- Node: Sets
Previous: General -=- Up: Top -=- Next: Limit


                         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).
Finally for each of the two sets a previously defined prototype set
(see *note Protoypes:(proto.doc)) can be used. In this case
the center of geometry (or mass with keyword MASS) and dipole of each
individual set member will be used in the requested functions.
For both sets WATEr is the default.



File: RDFSol -=- Node: Limit
Previous: Sets -=- Up: Top -=- Next: Options


                         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.



File: RDFSol -=- Node: Options
Previous: Limit -=- Up: Top -=- Next: Traj


                         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)
----------------------------------------------------------------------
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. Currently only one
                function can be calculated if MINI is used. Its major use
                is the computation of the distance dependent Kirkwood
                G-factor (with DDIP, second column). Here, one needs to go
                'into the corners' (i.e. sqrt(3)/2 * L for a cubic box)
                without counting pairs twice.
		(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)
----------------------------------------------------------------------



File: RDFSol -=- Node: Traj
Previous: Options -=- Up: Top -=- Next: Caveats


                         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.



File: RDFSol -=- Node: Caveats
Previous: Traj -=- Up: Top -=- Next: Examples


                         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



File: RDFSol -=- Node: Examples
Previous: Caveats -=- Up: Top


                         Examples

(See also test/c30test/rdfsol.inp test/c30test/rdfsol2.inp testcases)

----------------------------------------------------------------------

RDFSOL RDF 10 SETA WATER SAME RMAX 7.5 NBIN 150 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 10 SETA WATER SAME RMAX 7.5 NBIN 150 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.

----------------------------------------------------------------------

RDFSOL RDF 10 QDIP 11 DDIP 12 SETA WATER SAME RMAX 7.5 NBIN 150 PRECISE -
       FIRSTUNIT 13 NUNIT 1

Same sets as in the first example but here all three functions are
calculated at once (i.e. the trajectory is only read once).

----------------------------------------------------------------------

NIH Helix/Biowulf Systems
charmm.org Homepage

CHARMM Documentation / Rick_Venable@nih.gov