next up previous contents
Next: Parameters Up: Elastic Incoherent Structure Factor Previous: Elastic Incoherent Structure Factor   Contents


Theory and implementation

The Elastic Incoherent Structure Factor (EISF) is defined as the limit of the incoherent intermediate scattering function for infinite time,
\begin{displaymath}
EISF({\bf q}) \doteq \lim_{t \to \infty} {\cal F}_{\mathrm{inc}}({\bf q},t).
\end{displaymath} (4.174)

Using the above definition of the EISF one can decompose the incoherent intermediate scattering function as follows:
\begin{displaymath}
{\cal F}_{\mathrm{inc}}({\bf q},t) = EISF({\bf q}) + {\cal F}_{\mathrm{inc}}'({\bf q},t),
\end{displaymath} (4.175)

where ${\cal F}_{\mathrm{inc}}'({\bf q},t)$ decays to zero for infinite time. Taking now the Fourier transform it follows immediately that
\begin{displaymath}
{\cal S}_{\mathrm{inc}}({\bf q},\omega) = EISF({\bf q})\delta(\omega)
+ {\cal S}'_{\mathrm{inc}}({\bf q},\omega).
\end{displaymath} (4.176)

The EISF appears as the amplitude of the elastic line in the neutron scattering spectrum. Elastic scattering is only present for sytems in which the atomic motion is confined in space, as for solids. To understand which information is contained in the EISF we consider for simplicity a system where only one sort of atoms is visible to the neutrons. To a very good approximation this is the case for all systems containing a large amount of hydrogen atoms, as biological systems. Incoherent scattering from hydrogen dominates by far all other contributions. Using the definition of the van Hove self-correlation function $G_s({\bf r},t)$ [7],
\begin{displaymath}
b_{\mathrm{inc}}^2 G_s({\bf r},t) \doteq \frac{1}{2\pi^3}
\int d^3q \exp[-i{\bf q}\cdot{\bf r}] {\cal F}_{inc}({\bf q},t),
\end{displaymath} (4.177)

which can be interpreted as the conditional probability to find a tagged particle at the position ${\bf r}$ at time $t$, given it started at ${\bf r} = {\bf0}$, one can write:
\begin{displaymath}
EISF({\bf q}) = b_{\mathrm{inc}}^2\int d^3r \exp[i{\bf q}\cdot{\bf r}] 
G_s({\bf r},t=\infty).
\end{displaymath} (4.178)

The EISF gives the sampling distribution of the points in space in the limit of infinite time. In a real experiment this means times longer than the time which is observable with a given instrument. The EISF vanishes for all systems in which the particles can access an infinite volume since $G_s({\bf r},t)$ approaches $1/V$ for large times. This is the case for molecules in liquids and gases.

For computational purposes it is convenient to use the following representation of the EISF [14]:

\begin{displaymath}
EISF({\bf q}) = \sum^{N_{species}}_{I = 1} n_I \omega_{I,\mathrm{inc}} EISF_I(q)
\end{displaymath} (4.179)

where $N_{species}$ is the number of selected species, $n_I$ the number of atoms of species I, $\omega_{I,\mathrm{inc}}$ the weight for specie I (see Section 4.2.1 for more details) and for each specie the following expression for the elastic incoherent scattering function is
\begin{displaymath}
EISF_I({\bf q}) = \frac{1}{n_I}\sum^{n_I}_{\alpha} \langle\vert\exp[i{\bf q}\cdot{\bf R}_\alpha]\vert^2\rangle.
\end{displaymath} (4.180)

This expression is derived from definition (4.174) of the EISF and expression (4.142) for the intermediate scattering function, using that for infinite time the relation
\begin{displaymath}
\langle\exp[-i{\bf q}\cdot{\bf R}_\alpha(0)]
\exp[i{\bf q}\c...
...
\langle\vert\exp[i{\bf q}\cdot{\bf R}_\alpha]\vert^2 \rangle
\end{displaymath} (4.181)

holds. In this way the computation of the EISF is reduced to the computation of a static thermal average. We remark at this point that the length of the MD trajectory from which the EISF is computed should be long enough to allow for a representative sampling of the conformational space.

nMOLDYN allows one to compute the elastic incoherent structure factor on a grid of equidistantly spaced points along the q-axis:

\begin{displaymath}
EISF(q_m) \doteq \sum^{N_{species}}_{I = 1} n_I \omega_I EISF_I(q_m), m = 0\ldots N_q - 1.
\end{displaymath} (4.182)

where $N_q$ is a user-defined number of q-shells, the values for $q_m$ are defined as $q_m = q_{min} + m\cdot\Delta q$, and for each specie the following expression for the elastic incoherent scattering function is:
\begin{displaymath}
EISF_I(q_m) = \frac{1}{n_I}\sum^{n_I}_{\alpha} \overline{ \langle\vert\exp[i{\bf q}\cdot{\bf R}_\alpha]\vert^2\rangle}^q.
\end{displaymath} (4.183)

Here the symbol $\overline{\rule{0pt}{5pt}\ldots}^{q}$ denotes an average over the q-vectors having the same modulus $q_m$. The program corrects the atomic input trajectories for jumps due to periodic boundary conditions.


next up previous contents
Next: Parameters Up: Elastic Incoherent Structure Factor Previous: Elastic Incoherent Structure Factor   Contents
pellegrini eric 2009-10-06