Physics:Shifted force method

From HandWiki

The net electrostatic force acting on a charged particle with index [math]\displaystyle{ i }[/math] contained within a collection of particles is given as: [math]\displaystyle{ \mathbf{F}(\mathbf{r}) = \sum_{j \ne i}F(r)\mathbf{\hat{r}} \, \, \,; \, \, F(r) = \frac{q_{i}q_{j}}{4\pi\varepsilon_0 r^{2}} }[/math]

where [math]\displaystyle{ \mathbf{r} }[/math] is the spatial coordinate, [math]\displaystyle{ j }[/math] is a particle index, [math]\displaystyle{ r }[/math] is the separation distance between particles [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math], [math]\displaystyle{ \mathbf{\hat{r}} }[/math] is the unit vector from particle [math]\displaystyle{ j }[/math] to particle [math]\displaystyle{ i }[/math], [math]\displaystyle{ F(r) }[/math] is the force magnitude, and [math]\displaystyle{ q_{i} }[/math] and [math]\displaystyle{ q_{j} }[/math] are the charges of particles [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math], respectively. With the electrostatic force being proportional to [math]\displaystyle{ r^{-2} }[/math], individual particle-particle interactions are long-range in nature, presenting a challenging computational problem in the simulation of particulate systems. To determine the net forces acting on particles, the Ewald or Lekner summation methods are generally employed. One alternative and usually computationally faster technique based on the notion that interactions over large distances (e.g. > 1 nm) are insignificant to the net forces acting in certain systems is the method of spherical truncation.[1] The equations for basic truncation are:

[math]\displaystyle{ \displaystyle F_{CUT}(r) = \begin{cases} \frac{q_{i}q_{j}}{4\pi\varepsilon_0 r^{2}} & \text{for } r \le r_c \\ 0 & \text{for } r \gt r_c. \end{cases} }[/math]

where [math]\displaystyle{ r_c }[/math] is the cutoff distance. Simply applying this cutoff method introduces a discontinuity in the force at [math]\displaystyle{ r_c }[/math] that results in particles experiencing sudden impulses when other particles cross the boundary of their respective interaction spheres. In the particular case of electrostatic forces, as the force magnitude is large at the boundary, this unphysical feature can compromise simulation accuracy. A way to correct this problem is to shift the force to zero at [math]\displaystyle{ r_c }[/math], thus removing the discontinuity.[2] This can be accomplished with a variety of functions, but the most simple/computationally efficient approach is to simply subtract the value of the electrostatic force magnitude at the cutoff distance as such:

[math]\displaystyle{ \displaystyle F_{SF}(r) = \begin{cases} \frac{q_{i}q_{j}}{4\pi\varepsilon_0 r^{2}} - \frac{q_{i}q_{j}}{4\pi\varepsilon_0 r_c^{2}} & \text{for } r \le r_c \\ 0 & \text{for } r \gt r_c. \end{cases} }[/math]

As mentioned before, the shifted force (SF) method is generally suited for systems that do not have net electrostatic interactions that are long-range in nature. This is the case for condensed systems that show electric-field screening effects. Note that anisotropic systems (e.g. interfaces) may not be accurately simulated with the SF method,[3] although an adaption of the SF method for interfaces has been recently suggested.[4] Additionally, note that certain system properties (e.g. energy-dependent observables) will be more greatly influenced by the use of the SF method than others. It is not safe to assume, without reasonable argument, that the SF method can be used to accurately determine a certain property for a given system. If the accuracy of the SF method need be tested, this may be done by testing for convergence (i.e. showing that simulation results do not significantly change with increasing cutoff) or by comparing with results obtained through other electrostatics techniques (such as Ewald) that are known to perform well.[5] As a rough rule of thumb, results obtained with the SF method tend to be sufficiently accurate when the cutoff is at least five times larger than the distance of the near neighbor interactions.

With the SF method, a discontinuity is still present in the derivative of the force, and it may be preferable for ionic liquids to further alter the force equation as to remove that discontinuity.[6]

References

  1. Fennell, C. J.; Gezelter, J. D. (2006). "Is the Ewald summation still necessary? Pairwise alternatives to the accepted standard for long-range electrostatics". The Journal of Chemical Physics 124 (23): 234104. doi:10.1063/1.2206581. PMID 16821904. Bibcode2006JChPh.124w4104F. 
  2. Toxvaerd, S. R.; Dyre, J. C. (2011). "Communication: Shifted forces in molecular dynamics". The Journal of Chemical Physics 134 (8): 081102. doi:10.1063/1.3558787. PMID 21361519. Bibcode2011JChPh.134h1102T. 
  3. Spohr, E. (1997). "Effect of electrostatic boundary conditions and system size on the interfacial properties of water and aqueous solutions". The Journal of Chemical Physics 107 (16): 6342–6348. doi:10.1063/1.474295. Bibcode1997JChPh.107.6342S. 
  4. Welch, D. A.; Mehdi, B. L.; Hatchell, H. J.; Faller, R.; Evans, J. E.; Browning, N. D. (2015). "Using molecular dynamics to quantify the electrical double layer and examine the potential for its direct observation in the in-situ TEM". Advanced Structural and Chemical Imaging 1. doi:10.1186/s40679-014-0002-2. 
  5. Hansen, J. S.; Schrøder, T. B.; Dyre, J. C. (2012). "Simplistic Coulomb Forces in Molecular Dynamics: Comparing the Wolf and Shifted-Force Approximations". The Journal of Physical Chemistry B 116 (19): 5738–5743. doi:10.1021/jp300750g. PMID 22497264. 
  6. Kale, S.; Herzfeld, J. (2011). "Pairwise Long-Range Compensation for Strongly Ionic Systems". Journal of Chemical Theory and Computation 7 (11): 3620–3624. doi:10.1021/ct200392u. PMID 22247701.