next up previous contents
Next: Parameters Up: Coordination number Previous: Coordination number   Contents


Theory and implementation

In chemistry, the CN is the total number of neighbours of a central atom in a molecule or ion. The definition used in nMOLDYN is somewhat different and can be seen as an extension of as the former definition. Indeed, in nMOLDYN, the CN is not defined over one defined central atom but around the centers of gravity of a set of group of atoms. So, if only one group made of only atom is selected for the analysis, then, the definition is the same as the original definition. In that context, the CN is defined as:
\begin{displaymath}
n(r, r+dr) = \frac{1}{N_{G}}\sum_{g=1}^{N_{G}}\sum_{I=1}^{N_{species}} n_{gI}(r, r+dr)
\end{displaymath} (4.191)

where $N_{G}$ is the number of groups of atoms, $N_{species}$ is the number of species found in the system and $n_{gI}(r)$ is the CN defined for specie I defined as the number of atoms of species I found in a shell of width dr at a distance r of the center of gravity of the group of atom g.

nMOLDYN allows one to compute the CN on a set of equidistantly spaced distances at different times:

\begin{displaymath}
CN(r_m) \doteq \frac{1}{N_{frames}}\frac{1}{N_{G}}\sum_{f=1}...
..._f),
\qquad m = 0\ldots N_r - 1,\; n = 0\ldots N_{frames} - 1.
\end{displaymath} (4.192)

where $N_r$ and $N_{frames}$ are respectively the number of distances and times at which the CN is evaluated and
\begin{displaymath}
CN_{gI}(r_m, t_f) = n_{gI}(r_m, t_f),
\end{displaymath} (4.193)

is the number of atoms of specie I found within $[r_m,r_m + dr]$ at frame f from the center of gravity of group g.

From these expression, several remarks can be done. Firstly, the Eqs 4.192 and 4.193 can be restricted to intramolecular and intermolecular distances only. Secondly, these equations can be averaged over the selected frames providing a time averaged intra and intermolecular CN. Finally, the same equations (time-dependent and time-averaged) can be integrated over $r$ to provide a cumulative CN. nMOLDYN computes all these variations.

The concept of CN is useful for structure-related analysis. It can reveal for instance some packing effects that may have occured during the simulation.


next up previous contents
Next: Parameters Up: Coordination number Previous: Coordination number   Contents
pellegrini eric 2009-10-06