Introduction

Various pair distribution functions

There is a plethora of different pair distribution functions and various definitions. In this code, the definitions according to the book Underneath the Bragg Peaks by Takeshi Egami and Simon J.L. Billinge are used. 1 These are also the definitions used in the code PDFgetX3. 2 Let us start with the following function:

\[R_{ij}(r) = \frac{1}{N\mathrm{d}r}\ \frac{K_i K_j}{\langle K \rangle^2}\ n_{ij}(r)\]
\[R(r) = \displaystyle\sum_{i,j=1}^{n} R_{ij}(r)\]
_images/NaCl_R_r.png

Computed G(r) function of NaCl.

Herein, \(n_{ij}(r)\) is the number of atoms of type j between distances \(r\) and \(r+\mathrm{d}r\) from atoms of type i, \(N\) is the number of atoms in the unit cell, \(K_i\) and \(K_i\) are the atomic numbers if atoms i and j, and \(\langle K \rangle\) is the average atomic number of all atoms in the unit cell. While \(R(r)\) is a very intuitive function, it is not commonly used because it oscillates around \(4\pi\rho r^2\) (with the number density \(\rho\)), making it unpractical for plotting. An alternative function is:

\[g_{ij}(r) = \frac{R_{ij}(r)}{4\pi\rho r^2}\]
\[g(r) = \displaystyle\sum_{i,j=1}^{n} g_{ij}(r)\]
_images/NaCl_g_r.png

Computed g(r) function of NaCl.

\(g(r)\) oscillates around 1 and is therefore much nicer to plot than \(R(r)\). This is the default function used in xPDFsim. However, to obtain one of these functions from experimental data, one needs to make assumptions about the number density of the material. The function most commonly obtained from experimental data is:

\[G(r) = \frac{R(r)}{r}-4\pi\rho r\]
_images/NaCl_G_r-2.png

Computed G(r) function of NaCl.

In contrast to the functions above, \(G(r)\) cannot be expressed as a linear combination of the underlying element-element correlations.

Sources

  1. Egami, T. & Billinge, S. J. L. Underneath the Bragg Peaks: Structural Analysis of Complex Materials (Pergamon, 2003).

  2. Juhás, P., Davis, T., Farrow, C. L. & Billinge, S. J. L. PDFgetX3: a rapid and highly automatable program for processing powder diffraction data into total scattering pair distribution functions. J. Appl. Crystallogr. 46, 560–566 (2013).