next up previous contents
Next: Parameters Up: Pair Distribution Function Previous: Pair Distribution Function   Contents


Theory and implementation

The PDF is an example of a pair correlation function, which describes how, on average, the atoms in a system are radially packed around each other. This proves to be a particularly effective way of describing the average structure of disordered molecular systems such as liquids. Also in systems like liquids, where there is continual movement of the atoms and a single snapshot of the system shows only the instantaneous disorder, it is extremely useful to be able to deal with the average structure.

The PDF is useful in other ways. For example, it is something that can be deduced experimentally from x-ray or neutron diffraction studies, thus providing a direct comparison between experiment and simulation. It can also be used in conjunction with the interatomic pair potential function to calculate the internal energy of the system, usually quite accurately.

Mathematically, the PDF can be computed using the following formula:

\begin{displaymath}
PDF(r)=\sum_{I = 1,J\ge I}^{N_{species}}n_In_J \omega_I \omega_J g_{IJ}(r)
\end{displaymath} (4.186)

where $N_{species}$ is the number of selected species, $n_I$ and $n_J$ are respectively the numbers of atoms of species I and J, $\omega_I$ and $\omega_J$ respectively the weights for species I and J (see Section 4.2.1 for more details) and $PDF_{\alpha\beta}(r)$ is the partial PDF for $I$ and $J$ species that can be defined as:
\begin{displaymath}
PDF_{IJ}(r) = \frac{\left\langle\sum_{\alpha = 1}^{n_I} n_{\alpha J}(r)\right\rangle}{n_I\rho_J 4\pi r^2dr}
\end{displaymath} (4.187)

where $\rho_J$ is the density of atom of specie J and $n_{\alpha J}(r)$ is the mean number of atoms of specie J in a shell of width dr at distance r of the atom $\alpha$ of specie I.

From the computation of PDF, two related quantities are computed in nMOLDYN, the Radial-Distribution Function (RDF) defined as:

\begin{displaymath}
RDF(r) = 4 \pi r^2 \rho_0 PDF(r)
\end{displaymath} (4.188)

and the Total-Correlation Function (TCF) defined as:
\begin{displaymath}
TCF(r) = 4\pi r \rho_0 (PDF(r) - 1.0)
\end{displaymath} (4.189)

where $\rho_0$ is the average atomic density defined as:
\begin{displaymath}
\rho_0 = \frac{N}{V}
\end{displaymath} (4.190)

where N is the total number of atoms of the system and $V$ the volume of the simulation box.

In nMOLDYN, the PDF, the RDF and the TCF are further splitted into an intra-and inter-molecular parts which added together give respectively the total PDF, RDF and TCF.


next up previous contents
Next: Parameters Up: Pair Distribution Function Previous: Pair Distribution Function   Contents
pellegrini eric 2009-10-06