previous (interaction potentials)

theory index

next (thermodynamic properties)

Fourier transformation is problematic for long range potentials (for example Coulomb interaction) since they usually do not become zero on the discretization interval. Different approaches have to be used in such cases. Please note, that this text concerns only the Hypernetted Chain closure relation at the moment. Percus-Yevick related discussion will be added later (PY shouldn’t be used with long ranged (Coulomb) potential, since it doesn’t perform well!).

Ng renormalization (*J. Chem. Phys.* **1974**, *61*, 2680)

As suggested by Ng, for potentials in the form U(r) = A/r we use modified functions

where is some suitable well-behaved function. The expression for the direct correlation function can be easily proven. c(r) is given within the Hypernetted Chain approximation by

The second expression is used as a new (short-ranged) version of the HNC closure relation for the Ornstein-Zernike equation.

Regarding the suitable functions to be used in this procedure, there are some possibilities. We are looking for a function that is capable of counteracting the long range tail of the interaction potential. One of such functions is

with being adjustable parameter (pyOZ uses 1.08). This function agrees with the A/r potential overall with the exception of small distances. Therefore, subtracting this function from the interaction potential leads to a short ranged function (see figure for comparison, click for bigger version).

For the purpose of pyOZ, we are interested in the Fourier-Bessel transform of this function, which is performed using the sine transform (see the links if you are not familiar with these transforms!).

Integration per parts leads to

Since the first term vanishes at integration limits, we are left with

Solution of the integral is known, since it is equivalent to the Fourier transform of a gaussian

Since gaussian is an even function, the Fourier transform involves only cosine terms and the integral from to can be replaced by twice the integral from zero to infinity. The used FT was unitary (having the same normalization factor for forward and inverse transform). In order to have non-unitary transform (used i pyOZ), we need to divide by the normalization factor , yielding

Taking this expression and substituting the integral discussed above yields

Last expression is the sought Fourier-Bessel transform of the function .

Modified solution procedure can be described as follows (see Ornstein-Zernike equation theory and the pyOZ algorithm description to get more information). The O-Z equation is written by means of short-ranged functions

Fourier transformation on short-range functions (the indirect correlation function uses modified closure defined previously) is performed discretely without problems. Recall, that Fourier transform of is known analytically.

Then, uncorrected (original) functions in the Fourier space are restored (FT is linear, i.e., FT(f+g) = FT(f)+FT(g))

This equation is solved for and indirect **short-ranged** correlation function is calculated.

After transforming back to real space, short-ranged Gamma is used in the closure relation, new **short-ranged** direct correlation function is calculated and used for next iteration (until convergence). In reality, pyOZ uses this algorithm for convenience even in case of short-range only potentials. In such a case, the correction (and its Fourier transform as well) are set to zero, restoring thereby the non-modified Ornstein-Zernike equation.

previous (interaction potentials)

theory index

next (thermodynamic properties)