Ornstein–Zernike equation

In statistical mechanics the Ornstein–Zernike (OZ) equation is an integral equation introduced by Leonard Ornstein and Frits Zernike that relates different correlation functions with each other. Together with a closure relation, it is used to compute the structure factor and thermodynamic state functions of amorphous matter like liquids or colloids.

Context

The OZ equation has practical importance as a foundation for approximations for computing the pair correlation function of molecules or ions in liquids, or of colloidal particles. The pair correlation function is related via Fourier transform to the static structure factor, which can be determined experimentally using X-ray diffraction or neutron diffraction.

The OZ equation relates the pair correlation function to the direct correlation function. The direct correlation function is only used in connection with the OZ equation, which can actually be seen as its definition.

Besides the OZ equation, other methods for the computation of the pair correlation function include the virial expansion at low densities, and the Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) hierarchy. Any of these methods must be combined with a physical approximation: truncation in the case of the virial expansion, a closure relation for OZ or BBGKY.

The equation

To keep notation simple, we only consider homogeneous fluids. Thus the pair correlation function only depends on distance, and therefore is also called the radial distribution function. It can be written

[math]\displaystyle{ g(\mathbf{r}_1,\mathbf{r}_2) = g(\mathbf{r}_1 - \mathbf{r}_2) \equiv g(\mathbf{r}_{12}) = g(|\mathbf{r}_{12}|) \equiv g(r_{12}) \equiv g(12), }[/math]

where the first equality comes from homogeneity, the second from isotropy, and the equivalences introduce new notation.

It is convenient to define the total correlation function as:

[math]\displaystyle{ h(12)\equiv g(12)-1 }[/math]

which expresses the influence of molecule 1 on molecule 2 at distance [math]\displaystyle{ \,r_{12}\, }[/math]. The OZ equation

[math]\displaystyle{ h(12) \; = \; c(12) \; + \; \rho \, \int \text{d}^3 \mathbf{r}_{3}\, c(13) \, h(32) }[/math]

splits this influence into two contributions, a direct and indirect one. The direct contribution defines the direct correlation function, [math]\displaystyle{ c(r). }[/math] The indirect part is due to the influence of molecule 1 on a third, labeled molecule 3, which in turn affects molecule 2, directly and indirectly. This indirect effect is weighted by the density and averaged over all the possible positions of molecule 3.

By eliminating the indirect influence, [math]\displaystyle{ \,c(r)\, }[/math] is shorter-ranged than [math]\displaystyle{ h(r) }[/math] and can be more easily modelled and approximated. The radius of [math]\displaystyle{ \,c(r)\, }[/math] is determined by the radius of intermolecular forces, whereas the radius of [math]\displaystyle{ \,h(r)\, }[/math] is of the order of the correlation length.

Fourier transform

The integral in the OZ equation is a convolution. Therefore, the OZ equation can be resolved by Fourier transform. If we denote the Fourier transforms of [math]\displaystyle{ h(\mathbf{r}) }[/math] and [math]\displaystyle{ c(\mathbf{r}) }[/math] by [math]\displaystyle{ \hat{h}(\mathbf{k}) }[/math] and [math]\displaystyle{ \hat{c}(\mathbf{k}) }[/math], respectively, and use the convolution theorem, we obtain

[math]\displaystyle{ \hat{h}(\mathbf{k}) \; = \; \hat{c}(\mathbf{k}) \; + \; \rho \; \hat{h}(\mathbf{k})\;\hat{c}(\mathbf{k})~ , }[/math]

which yields

[math]\displaystyle{ \hat{c}(\mathbf{k}) \; = \; \frac{\hat{h}(\mathbf{k})}{\;1 \;+\;\rho \;\hat{h}(\mathbf{k})\;} \qquad \text{ and } \qquad \hat{h}(\mathbf{k}) \; = \; \frac{\hat{c}(\mathbf{k})}{\; 1 \; - \; \rho \; \hat{c}(\mathbf{k}) \;} ~. }[/math]

Closure relations

As both functions, [math]\displaystyle{ \,h \, }[/math] and [math]\displaystyle{ \,c \, }[/math], are unknown, one needs an additional equation, known as a closure relation. While the OZ equation is purely formal, the closure must introduce some physically motivated approximation.

In the low-density limit, the pair correlation function is given by the Boltzmann factor,

[math]\displaystyle{ g(12)=\text{e}^{-\beta u(12)},\quad \rho\to 0 }[/math]

with [math]\displaystyle{ \beta=1/k_\text{B} T }[/math] and with the pair potential [math]\displaystyle{ u(r) }[/math].

Closure relations for higher densities modify this simple relation in different ways. The best known closure approximations are:

• The Percus–Yevick approximation for particles with impenetrable ("hard") core,
• the hypernetted-chain approximation, for particles with soft cores and attractive potential tails,
• the mean spherical approximation,
• the Rogers-Young approximation.

The latter two interpolate in different ways between the former two, and thereby achieve a satisfactory description of particles that have a hard core and attractive forces.

External links

• "The Ornstein–Zernike equation and integral equations". http://cbp.tnw.utwente.nl/PolymeerDictaat/node15.html. 
• "Multilevel wavelet solver for the Ornstein–Zernike equation". http://www4.ncsu.edu/~ctk/PAPERS/OZwavelet4.pdf. 
• "Analytical solution of the Ornstein–Zernike equation for a multicomponent fluid". http://www.iop.org/EJ/S/UNREG/lLP4nnFLwybFbk9aWg47cQ/article/0953-8984/12/38/101/c038l1.pdf. 
• "The Ornstein–Zernike equation in the canonical ensemble". http://www.iop.org/EJ/article/0295-5075/54/4/475/6545.html.