Computational implementation of integral equations

From SklogWiki
Jump to navigation Jump to search

Integral equations are solved numerically. One has the Ornstein-Zernike relation, and a closure relation, (which incorporates the bridge function ). The numerical solution is iterative;

  1. trial solution for
  2. calculate
  3. use the Ornstein-Zernike relation to generate a new etc.

Note that the value of is local, i.e. the value of at a given point is given by the value of at this point. However, the Ornstein-Zernike relation is non-local. The way to convert the Ornstein-Zernike relation into a local equation is to perform a (fast) Fourier transform (FFT). Note: convergence is poor for liquid densities. (See Ref.s 1 to 6).

Picard iteration[edit]

Picard iteration generates a solution of an initial value problem for an ordinary differential equation (ODE) using fixed-point iteration. Here are the four steps used to solve integral equations:

Closure relation [edit]

(Note: for linear fluids )

Perform the summation[edit]

where is the separation between molecular centers and the sets of Euler angles needed to specify the orientations of the two molecules, with

with .

Define the variables[edit]

Thus

.

Evaluate[edit]

Evaluations of are performed at the discrete points where the are the roots of the Legendre polynomial where are the roots of the Chebyshev polynomial and where are the roots of the Chebyshev polynomial thus

where



where is the angular, , part of the rotation matrix , and

For the limits in the summations

The above equation constitutes a separable five-dimensional transform. To rapidly evaluate this expression it is broken down into five one-dimensional transforms:

Operations involving the and basis functions are performed in complex arithmetic. The sum of these operations is asymptotically smaller than the previous expression and thus constitutes a ``fast separable transform". and are parameters; is the number of nodes in the Gauss integration, and the the max index in the truncated rotational invariants expansion.

Integrate over angles [edit]

Use Gauss-Legendre quadrature for and Use Gauss-Chebyshev quadrature for , and . Thus

where the Gauss-Legendre quadrature weights are given by

while the Gauss-Chebyshev quadrature has the constant weight

Perform FFT from Real to Fourier space [edit]

This is non-trivial and is undertaken in three steps:

Conversion from axial reference frame to spatial reference frame[edit]

this is done using the Blum transformation (Refs 7, 8 and 9):

Fourier-Bessel Transforms[edit]

(see Blum and Torruella Eq. 5.6 in Ref. 7 or Lado Eq. 39 in Ref. 3), where is a Bessel function of order . `step-down' operations can be performed by way of sin and cos operations of Fourier transforms, see Eqs. 49a, 49b, 50 of Lado Ref. 3. The Fourier-Bessel transform is also known as a Hankel transform. It is equivalent to a two-dimensional Fourier transform with a radially symmetric integral kernel.


Conversion from the spatial reference frame back to the axial reference frame[edit]

this is done using the Blum transformation

Ornstein-Zernike relation [edit]

For simple fluids:

For molecular fluids (see Eq. 19 of Lado Ref. 3)

where and are matrices with elements .

For mixtures of simple fluids (see Ref. 10 Juan Antonio Anta PhD thesis pp. 107--109):

Conversion back from Fourier space to Real space[edit]

(basically the inverse of step 2).

Axial reference frame to spatial reference frame[edit]

Inverse Fourier-Bessel transform[edit]

'Step-up' operations are given by Eq. 53 of Ref. 3. The inverse Hankel transform is

Change from spatial reference frame back to axial reference frame[edit]

.

Ng acceleration[edit]

Angular momentum coupling coefficients[edit]

References[edit]

  1. M. J. Gillan "A new method of solving the liquid structure integral equations" Molecular Physics 38 pp. 1781-1794 (1979)
  2. Stanislav Labík, Anatol Malijevský and Petr Voncaronka "A rapidly convergent method of solving the OZ equation", Molecular Physics 56 pp. 709-715 (1985)
  3. F. Lado "Integral equations for fluids of linear molecules I. General formulation", Molecular Physics 47 pp. 283-298 (1982)
  4. F. Lado "Integral equations for fluids of linear molecules II. Hard dumbell solutions", Molecular Physics 47 pp. 299-311 (1982)
  5. F. Lado "Integral equations for fluids of linear molecules III. Orientational ordering", Molecular Physics 47 pp. 313-317 (1982)
  6. Enrique Lomba "An efficient procedure for solving the reference hypernetted chain equation (RHNC) for simple fluids" Molecular Physics 68 pp. 87-95 (1989)
  7. L. Blum and A. J. Torruella "Invariant Expansion for Two-Body Correlations: Thermodynamic Functions, Scattering, and the Ornstein—Zernike Equation", Journal of Chemical Physics 56 pp. pp. 303-310 (1972)
  8. L. Blum "Invariant Expansion. II. The Ornstein-Zernike Equation for Nonspherical Molecules and an Extended Solution to the Mean Spherical Model", Journal of Chemical Physics 57 pp. 1862-1869 (1972)
  9. L. Blum "Invariant expansion III: The general solution of the mean spherical model for neutral spheres with electostatic interactions", Journal of Chemical Physics 58 pp. 3295-3303 (1973)
  10. P. G. Kusalik and G. N. Patey " On the molecular theory of aqueous electrolyte solutions. I. The solution of the RHNC approximation for models at finite concentration", Journal of Chemical Physics 88 pp. 7715-7738 (1988)