Multiscale fast summation of long‐range charge and dipolar interactions

Here we present a linear order multiscale method for the fast summation of long range forces in a system consisting of a large number of charge and dipolar particles. For a N‐body system, our algorithm requires an order of work that is proportional to O(N), in comparison to order O(N2) of the direct pairwise computation. Our method is demonstrated on two‐dimensional homogeneous point‐charge and dipolar systems, and a combined heterogeneous particle system, for the calculation of the induced electrostatic potential and energy. The electrostatic interaction is decomposed into a local part and a smooth part. The method thus, has several potential advantages over other O(N log N) or O(N) techniques, especially for calculation with moving particles or implicit charges locations. This approach is beneficial to large‐scale problems such as molecular statics, molecular dynamics, equilibrium statistics (Monte‐Carlo simulations), molecular docking, and in areas such as magnetism and astrophysics. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 717–731, 2001

[1]  Andrew W. Appel,et al.  An Efficient Program for Many-Body Simulation , 1983 .

[2]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[3]  Achi Brandt,et al.  Optimal multigrid algorithms for calculating thermodynamic limits , 1994 .

[4]  Ramzi Kutteh,et al.  Implementing the cell multipole method for dipolar and charged dipolar systems , 1995 .

[5]  A. Brandt Multilevel computations of integral transforms and particle interactions with oscillatory kernels , 1991 .

[6]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[7]  A. Brandt,et al.  Multilevel matrix multiplication and fast solution of integral equations , 1990 .

[8]  Achi Brandt,et al.  Multilevel Evaluation of Integral Transforms with Asymptotically Smooth Kernels , 1998, SIAM J. Sci. Comput..

[9]  D. Ron,et al.  Renormalization Multigrid (RMG): Statistically Optimal Renormalization Group Flow and Coarse-to-Fine Monte Carlo Acceleration , 2001 .

[10]  L Greengard,et al.  Fast Algorithms for Classical Physics , 1994, Science.

[11]  Vittorio Rosato,et al.  Simulation of a Two-Dimensional Dipolar System on a APE100/Quadrics Simd Architecture , 1997 .

[12]  Ruth Nussinov,et al.  A Method for Biomolecular Structural Recognition and Docking Allowing Conformational Flexibility , 1998, J. Comput. Biol..

[13]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

[14]  R Nussinov,et al.  Flexible docking allowing induced fit in proteins: Insights from an open to closed conformational isomers , 1998, Proteins.

[15]  A. B. MacIsaac,et al.  SIMULATION OF THE REORIENTATION TRANSITION IN ULTRATHIN MAGNETIC FILMS WITH STRIPED AND TETRAGONAL PHASES , 1998 .