next up previous contents
Next: Parameters Up: Spatial Density Previous: Spatial Density   Contents


Theory and implementation

The SD can be seen as an generalization of the pair distribution function. Indeed, pair distribution functions are defined as orientionally averaged distribution functions. are in the sense that, Altough these correlation functions reflects many key features of the short-range order in molecular systems, it should be realized that an average spatial assembly of non-spherical particles can not be uniquely characterized from these one-dimensionals functions. So, structural models postulated for the molecular ordering in nonsimple systems based only on one-dimensional PDF will always be somewhat ambiguous. The goal of SD analysis is to provide greater clarity in the structual analysis of molecular systems by utilizing distribution function which span both the radial and angular coordinates of the separation vector. This can provide useful information about the average local structure in a complex system.

nMOLDYN allows one to compute the SD in spherical coordinates on a set of concentrics shells surrounding the centers of mass of selected triplets of atoms using the formula:

\begin{displaymath}
SD(r_l,\theta_m ,\phi_n) \doteq \frac{1}{N_{triplets N_{grou...
...N_r - 1, m = 0\ldots N_{\theta} - 1, n = 0\ldots N_{\phi} - 1.
\end{displaymath} (4.194)

where $N_{triplets}$ and $N_{groups}$ are respectively the number of triplets and groups, $r_l$, $\theta_m$ and $\phi_n$ are the spherical coordinates at which the SD is evaluated, $N_r$, $N_{\theta}$ and $N_{\phi}$ are respectively the number of discrete r, $\theta$ and $\phi$ values and $n_{tg}(r_l,\theta_m,\phi_n)$ is the number of group of atoms of type g whose centers of mass is found to be in the volume element defined by $[r, r+dr]$, $[\theta , \theta + d\theta]$ and $[\phi , \phi + d\phi]$ in the spherical coordinates basis centered on the center of mass of triplet t.

So technically, nMOLDYN proceeds more or less on the following way:


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