The net electrostatic force acting on a charged particle with index formula_1 contained within a collection of particles is given as:
formula_2
where formula_3 is the spatial coordinate, formula_4 is a particle index, formula_5 is the separation distance between particles formula_1 and formula_4, formula_8 is the unit vector from particle formula_4 to particle formula_1, formula_11 is the force magnitude, and formula_12 and formula_13 are the charges of particles formula_1 and formula_4, respectively. With the electrostatic force being proportional to formula_16, 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:
formula_17
where formula_18 is the cutoff distance. Simply applying this cutoff method introduces a discontinuity in the force at formula_18 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 formula_18, 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:
formula_21
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.
References.
 ;templatestyles src="Reflist/styles.css" / ;