next up previous contents
Next: Parameters Up: Velocity Autocorrelation Function Previous: Velocity Autocorrelation Function   Contents


Theory and implementation

The Velocity AutoCorrelation Function (VACF) is another interesting property describing the dynamics of a molecular system. Indeed, it reveals the underlying nature of the forces acting on the system.

In a molecular system that would be made of non interacting particles, the velocities would be constant at any time triggering the VACF to be a constant value. Now, if we think about a system with small interactions such as in a gas-phase, the magnitude and direction of the velocity of a particle will change gradually over time due to its collision with the other particles of the molecular system. In such a system, the VACF will be represented by a decaying exponential.

In the case of solid phase, the interaction are much stronger and, as a results, the atoms are bound to a given position from which they will move backwards and forwards oscillating between positive and negative values of their velocity. The oscillations will not be of equal magnitude however, but will decay in time, because there are still perturbative forces acting on the atoms to disrupt the perfection of their oscillatory motion. So, in that case the VACF will look like a damped harmonic motion.

Finally, in the case of liquid phase, the atoms have more freedom than in solid phase and because of the diffusion process, the oscillatory motion seen in solid phase will be cancelled quite rapidly depending on the density of the system. So, the VACF will just have one very damped oscillation before decaying to zero. This decaying time can be considered as the average time for a collision between two atoms to occur before they diffuse away.

Mathematically, the VACF of atom $\alpha$ in an atomic or molecular system is usually defined as

\begin{displaymath}
C_{vv ; \alpha\alpha}(t) \doteq
\frac{1}{3}\langle {\bf v}_\alpha(t_0)\cdot{\bf v}_\alpha(t_0+t)\rangle_{t_0}.
\end{displaymath} (4.36)

In some cases, e.g. for non-isotropic systems, it is useful to define VACF along a given axis,
\begin{displaymath}
C_{vv ; \alpha\alpha}(t;{\bf n}) \doteq
\langle v_\alpha(t_0;{\bf n})v_\alpha(t_0+t;{\bf n})\rangle_{t_0},
\end{displaymath} (4.37)

where $v_\alpha(t;{\bf n})$ is given by
\begin{displaymath}
v_\alpha(t;{\bf n}) \doteq
{\bf n}\cdot{\bf v}_{\alpha}(t).
\end{displaymath} (4.38)

The vector n is a unit vector defining a space-fixed axis.

The VACF of the particles in a many body system can be related to the incoherent dynamic structure factor by the relation:

\begin{displaymath}
lim_{q\to 0} \frac{\omega^2}{q^2}{\cal S}({\bf q},\omega) =
G(\omega),
\end{displaymath} (4.39)

where $G(\omega)$ is the Density Of States (DOS). For an isotropic system it reads
$\displaystyle G(\omega)$ $\textstyle =$ $\displaystyle \sum_{\alpha}b^2_{\alpha,inc}
\tilde C_{vv ; \alpha\alpha}(\omega),$ (4.40)
$\displaystyle \tilde C_{vv ; \alpha\alpha}(\omega)$ $\textstyle =$ $\displaystyle \frac{1}{2\pi}\int_{-\infty}^{+\infty}dt  \exp[-i\omega t]
C_{vv ; \alpha\alpha}(t).$ (4.41)

For non-isotropic systems relation (4.39) holds if the DOS is computed from the atomic velocity autocorrelation functions $C_{vv ; \alpha\alpha}(t;{\bf n}_q)$, where ${\bf n}_q$ is the unit vector in the direction of q.


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