| Latest revision |
Your text |
| Line 1: |
Line 1: |
| The '''Wolf method''' <ref>[http://dx.doi.org/10.1103/PhysRevLett.68.3315 Dieter Wolf "Reconstruction of NaCl surfaces from a dipolar solution to the Madelung problem", Physical Review Letters '''68''' pp. 3315-3318 (1992)]</ref>. 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 the cut radius 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.
| | {{stub-general}} |
| | |
| :<math>E_{i} (R_{c})\approx E_{i}^{tot}(R_{c}) + E_{i}^{neutr}(R_{c})
| |
| </math>
| |
| | |
| where
| |
| | |
| :<math>E_{i}^{tot}(R_{c}) = \sum_{j\neq i, r_{ij} < R_{c}}^{} \frac{q_{i}q_{j}}{r_{ij}}
| |
| </math>
| |
| | |
| According to electrostatic rules, the energy of the interaction between an ion with charge <math>q_{i}</math> and a spherical charged surface of radium <math>R{c}</math> centered in the ion <math>i</math> is
| |
| | |
| :<math>E_{i}^{neutr}(R_{c})\approx \frac{q_{i}\Delta q_{i}(R_{c})}{R_{c}}</math>
| |
| | |
| <math>\Delta q_{i}(R_{c})</math> is the charge of the sphere of radium <math>R{c}</math> centered in the ion <math>i</math>:
| |
| | |
| :<math>\Delta q_{i}(R_{c})=\sum_{j, r_{ij}\leq R_{c}} q_{j}</math>
| |
| | |
| Thus, the energy in the Wolf method is calculated as
| |
| | |
| :<math>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}}
| |
| </math>
| |
| | |
| The total potential can be split into two parts:
| |
| | |
| :<math>U_{i}=\sum_{j\neq i}^{R_{c}}u_{ij} + u_{self}</math>
| |
| | |
| where
| |
| | |
| :<math>u_{ij}=q_{i}q_{j}\left(\frac{erfc(\alpha r_{ij})}{r_{ij}} - \frac{erfc(\alpha R_{c})}{R_{c}}\right)</math>
| |
| | |
| | |
| :<math>u_{self}= -q_{i}^{2}\left[\frac{1}{2} \frac{erfc(\alpha R_{c})}{R_{c}} + \frac{\alpha}{\sqrt{\pi}} \right]
| |
| </math>
| |
| | |
| | |
| | |
| In the Wolf method it is suggested that the reciprocal space term of the Ewald Sum may be essentially ignored, | | 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. | | provided one compensates the real space contribution by means of a truncation scheme. |
| Line 41: |
Line 5: |
| ==Inhomogeneous systems== | | ==Inhomogeneous systems== |
| It appears to be the case (Ref. 3) that the Wolf method has problems for inhomogeneous systems. | | It appears to be the case (Ref. 3) that the Wolf method has problems for inhomogeneous systems. |
|
| |
| :<math>U = U_{real} + U_{self}</math>
| |
|
| |
|
| ==See also== | | ==See also== |
| *[[Ewald sum]] | | *[[Ewald sum]] |
| ==References== | | ==References== |
| <references/>
| | #[http://dx.doi.org/10.1103/PhysRevLett.68.3315 Dieter Wolf "Reconstruction of NaCl surfaces from a dipolar solution to the Madelung problem", Physical Review Letters '''68''' pp. 3315-3318 (1992)] |
| ;Related reading
| | #[http://dx.doi.org/10.1063/1.478738 D. Wolf, P. Keblinski, S. R. Phillpot, and J. Eggebrecht "Exact method for the simulation of Coulombic systems by spherically truncated, pairwise r<sup>-1</sup> summation", Journal of Chemical Physics '''110''' pp. 8254- (1999)] |
| *[http://dx.doi.org/10.1063/1.478738 D. Wolf, P. Keblinski, S. R. Phillpot, and J. Eggebrecht "Exact method for the simulation of Coulombic systems by spherically truncated, pairwise r<sup>-1</sup> summation", Journal of Chemical Physics '''110''' pp. 8254- (1999)]
| | #[http://dx.doi.org/10.1063/1.2948951 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)] |
| *[http://dx.doi.org/10.1063/1.2948951 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)]
| |
| *[http://dx.doi.org/10.1063/1.4923001 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)]
| |
| | |
| | |
| [[Category: Computer simulation techniques]] | | [[Category: Computer simulation techniques]] |
| [[Category: Electrostatics]] | | [[Category: Electrostatics]] |