Physics:Shifted force method

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. 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. 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, although an adaption of the SF method for interfaces has been recently suggested. 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. 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.