radial | rdf

Calculate radial distribution function.

radial <outfilename> <spacing> <maximum> <solvent mask1>[<solute mask2>] [noimage] [density <density> | volume][center1 | center2 | nointramol] [<name>][intrdf <file>] [rawrdf <file>]

<outfilename> File to write RDF to, required.
<spacing> Bin spacing, required.
<maximum> Max bin value, required.
<solvent mask1> Atoms to calculate RDF for, required.
[<solute mask2>] (Optional) If specified calculate RDF of all atoms in <solvent mask1> to each atom in <solute mask2>.
[noimage] Do not image distances.
[density <density>] Use density value of <density> for normalization (default 0.033456 molecules Å^3).
[volume] Determine density for normalization from average volume of input frames.
[center1] Calculate RDF from geometric center of atoms in <solvent mask1> to all atoms in <solute mask2>.
[center2] Calculate RDF from geometric center of atoms in <solute mask2> to all atoms in <solvent mask1>.
[nointramol] Ignore intra-molecular distances.
[<name>] Name radial dataset.
[intrdf <file>] Calculate integral of RDF bin values (averaged over # of frames but otherwise not normalized) and write to <file> (can be same as <output_filename>).
[rawrdf <file>] Write raw (non-normalized) RDF values to <file>.

DataSet Aspects:
[int] (intrdf only) Integral of RDF bin values.
[raw] (rawrdf only) Raw (non-normalized) RDF values.

Calculate the radial distribution function (RDF, aka pair correlation function) of atoms in <solvent mask1> (note that this mask does not need to be solvent, but this nomenclature is used for clarity). If an optional second mask (<solute mask2>) is given, calculate the RDF of ALL atoms in <solvent mask1> to EACH atom in <solute mask2>. If desired, the geometric center of atoms in <solvent mask1> or <solute mask2> can be used by specifying the center1 or center2 keywords respectively, or alternatively intra-molecular distances can be ignored by specifying the nointramol keyword.

The RDF is calculated from the histogram of the number of particles found as a function of distance R, normalized by the expected number of particles at that distance. The normalization is calculated from:

\text{Density} \cdot \left( \left[  \frac{4 \pi}{3}  (R + dR^3 ) \right] - \left[ \frac{4 \pi}{3} dR^3 \right] \right)

where dR is equal to the bin spacing. Some care is required by the user in order to normalize the RDF correctly.
The default density value is 0.033456 molecules \AA^3 , which corresponds to a density of water approximately equal to 1.0 g mL^{-1}. To convert a standard density in g/mL, multiply the density by \frac{0.6022}{Mr}, where Mr is the mass of the molecule in atomic mass units. Alternatively, if the volume keyword is specified the density is determined from the average volume of the system over all Frames.

Note that correct normalization of the RDF depends on the number of atoms in each mask; if multiple topology files are being processed that result in changes in the number of atoms in each mask, the normalization will be off.