Wolf method

The Wolf method . In order to calculate a volumetric property, the interaction between an ion and all the rest of the ions up to an infinite distance should be taken into account. However, in practice, summation has to be truncated somewhere. If the truncation is spherical, the sphere considered is not always electrically neutral, as the number of positive and negative ions at a distance no higher than cut radium is not the same. Of course, the charge excess is located in a shell whose thickness equals to the distance between neighbor ions. D. Wolf proposed the approximation that all the charge excess can be considered to be located exactly on the surface of truncation sphere: that is to say to consider the shell infinitely thin. Under this approximation, the potential energy of an ion can be calculated as the summation of the interaction of that ion with the rest of the ions contained in that sphere plus a correction term due to the interaction with the charged surface.

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{i} (R_{c})\approx E_{i}^{tot}(R_{c}) + E_{i}^{neutr}(R_{c}) }

where

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{i}^{tot}(R_{c}) = \sum_{j\neq i, r_{ij} < R_{c}}^{} \frac{q_{i}q_{j}}{r_{ij}} }

According to electrostatic rules, the energy of the interaction between an ion with charge Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle q_{i}} and a spherical charged surface of radium Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R{c}} centered in the ion Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle i} is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{i}^{neutr}(R_{c})\approx \frac{q_{i}\Delta q_{i}(R_{c})}{R_{c}}}

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta q_{i}(R_{c})} is the charge of the sphere of radium Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R{c}} centered in the ion $i$:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta q_{i}(R_{c})=\sum_{j, r_{ij}\leq R_{c}} q_{j}}

Thus, the energy in the Wolf method is calculated as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{i} (R_{c})\approx \sum_{j\neq i, r_{ij} < R_{c}}^{} \frac{q_{i}q_{j}}{r_{ij}} - \frac{q_{i}\Delta q_{i}(R_{c})}{R_{c}} }

The total potential can be split into two parts:

$U_{i}=\sum _{j\neq i}^{R_{c}}u_{ij}+u_{self}$

where

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle u_{ij}=q_{i}q_{j}\left(\frac{erfc(\alpha r_{ij})}{r_{ij}} - \frac{erfc(\alpha R_{c})}{R_{c}}\right)}

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle u_{self}= -q_{i}^{2}\left[\frac{1}{2} \frac{erfc(\alpha R_{c})}{R_{c}} + \frac{\alpha}{\sqrt{\pi}} \right] }

In the Wolf method it is suggested that the reciprocal space term of the Ewald Sum may be essentially ignored, provided one compensates the real space contribution by means of a truncation scheme.

Inhomogeneous systems

It appears to be the case (Ref. 3) that the Wolf method has problems for inhomogeneous systems.

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle U = U_{real} + U_{self}}

See also

• Ewald sum

References

Related reading

• D. Wolf, P. Keblinski, S. R. Phillpot, and J. Eggebrecht "Exact method for the simulation of Coulombic systems by spherically truncated, pairwise r^-1 summation", Journal of Chemical Physics 110 pp. 8254- (1999)
• Francisco Noé Mendoza, Jorge López-Lemus, Gustavo A. Chapela, and José Alejandre "The Wolf method applied to the liquid-vapor interface of water", Journal of Chemical Physics 129 024706 (2008)
• Björn Stenqvist, Martin Trulsson, Alexei I. Abrikosov and Mikael Lund "Direct summation of dipole-dipole interactions using the Wolf formalism", Journal of Chemical Physics 143 014109 (2015)