|Simulation of neutrons and X-rays scattering|
|Partial structure factors|
The pair correlation function g(r) cannot be measured experimentally, therefore the g(r) obtained using simulations cannot be compared directly to experimental data. Nevertheless it is possible to compare the results of the numerical simulations to neutron and X-rays scattering experiments using the static neutron structure factor S(q):
where bj and represent respectively the neutron or the X-rays scattering length and the position of the atom j. N represents the total number of atoms in the system studied. The X-rays scattering is q point dependent with bj = fj (q) this is implemented in the rings-code from version 1.2.5, see "X-ray scattering factors computed from numerical Hartree-Fock wave functions", in [a] for details.
Equation [Eq. 1] defines for the atoms of the system the dispersion of a beam radiation at a given q vector. In neutron scattering experiment it is necessary to sum all possible orientations of q compared to the vector - .
This average on the orientations of the q vector lead to the Debye result:
Nevertheless the instantaneous individual atomic contributions introduced by this equation [Eq. 2] are not easy to interpret, It is more interesting to express these contributions introducing the radial distribution functions.
In order to achieve this goal it is first necessary to split the self-atomic contribution (j = k), from the contribution between distinct atoms:
with cj = .
4π cjbj2 represents the total scattering cross section of the material.
The function I(q) which describes the interaction between distinct atoms is related to the radial distribution functions thru the Fourier transform:
where the function G(r) is defined using the partial radial distribution functions
where cα = and bα represents the neutron scattering length of species α.
G(r) towards to the value - cαbα cβbβ for r = 0 and to 0 for r→∞.
Usually the self-contributions are substracted from equation [Eq. 3] and the structure factor is normalized using the relation:
It is therefore possible to write the structure factor [Eq. 4] in a more standard way:
where(r) is defined using the partial radial distribution functions :
In the case of a single atomic species system the normalization allows to obtain values of S(q) and(r) which are independent of the scattering factor/length and therefore independent of the measurement technique.
Nevertheless in the other cases the total functions are combinations of the partial functions weighted using the scattering factor and therefore depend on the measurement technique (Neutron, X-rays ...) used or simulated.
Figure [Fig. 1] presents a comparison bewteen the calculations of the total neutron structure factor done using on the one hand the atomic correlations [Eq. 2] and on the other hand the pair correlation functions [Eq. 7].
The material studied is a sample of glassy GeS2 at 300 K obtained using ab-initio molecular dynamics.
In several cases the neutron structure factor S(q) and the radial distribution function (r) [Eq. 8] can be compared to experimental data.
Nevertheless depending on the neutron scattering formalism used other function can be used to be compared to the experiment. Among these function the G(r) has already been defined [Eq. 5], the differential correlation function D(r), G(r) and the total correlation function T(r) are defined by:
(r) is equal to zero for r = 0 and towards to 1 for r→∞.
D(r) is equal to zero for r = 0 and towards to 0 for r→∞.
G(r) is equal to zero for r = 0 and towards to 0 for r→∞.
T(r) is equal to zero for r = 0 and towards to ∞ for r→∞.
The result of the evaluation of these different neutron weighted functions for a glassy GeS2 sample at 300 K obtained using ab-intio molecular dynamics is presented in figure [Fig. 2].
This allows to illustrate the differences between each of these functions.
|< Prev||Next >|