P3M algorithm for dipolar interactions.

An extension to the P(3)M algorithm for electrostatic interactions is presented that allows to efficiently compute dipolar interactions in periodic boundary conditions. Theoretical estimates for the root-mean-square error of the forces, torques, and the energy are derived. The applicability of the estimates is tested and confirmed in several numerical examples. A comparison of the computational performance of the new algorithm to a standard dipolar-Ewald summation methods shows a performance crossover from the Ewald method to the dipolar P(3)M method for as few as 300 dipolar particles. In larger systems, the new algorithm represents a substantial improvement in performance with respect to the dipolar standard Ewald method. Finally, a test comparing point-dipole-based and charged-pair based models shows that point-dipole-based models exhibit a better performance than charged-pair based models.

[1]  Christian Holm,et al.  How to Mesh up Ewald Sums , 2000 .

[2]  K. Kremer,et al.  Advanced Computer Simulation Approaches for Soft Matter Sciences III , 2005 .

[3]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[4]  J. W. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. II. Equivalence of boundary conditions , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[5]  S. Kantorovich,et al.  Aggregate formation in ferrofluid monolayers: simulations and theory , 2008, Journal of physics. Condensed matter : an Institute of Physics journal.

[6]  Tim N. Heinz,et al.  Combining the lattice-sum and reaction-field approaches for evaluating long-range electrostatic interactions in molecular simulations. , 2005, The Journal of chemical physics.

[7]  Gerhard Hummer,et al.  Simulation and Theory of Electrostatic Interactions in Solution: Computational Chemistry, Biophysics and Aqueous Solutions, Santa Fe, New Mexico, U. S. A., 23-25 June 1999 , 1999 .

[8]  H. G. Petersen,et al.  An algorithm for the simulation of condensed matter which grows as the 3/2 power of the number of particles , 1988 .

[9]  O. Lenz,et al.  The optimal P3M algorithm for computing electrostatic energies in periodic systems. , 2007, The Journal of chemical physics.

[10]  Molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[12]  M. Deserno,et al.  HOW TO MESH UP EWALD SUMS. II. AN ACCURATE ERROR ESTIMATE FOR THE PARTICLE-PARTICLE-PARTICLE-MESH ALGORITHM , 1998, cond-mat/9807100.

[13]  S. J. Singer,et al.  Computer simulations of a two-dimensional system with competing interactions. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  B. M. Berkovskiĭ,et al.  Magnetic Fluids and Applications Handbook , 1996 .

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

[16]  Celeste Sagui,et al.  P3M and PME: A comparison of the two methods , 1999 .

[17]  Axel Arnold,et al.  Efficient methods to compute long-range interactions for soft matter systems , 2005 .

[18]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[19]  Richard J. Needs,et al.  An Optimized Ewald Method for Long-Ranged Potentials , 1994 .

[20]  T. Darden,et al.  Molecular dynamics simulations of biomolecules: long-range electrostatic effects. , 1999, Annual review of biophysics and biomolecular structure.

[21]  Philippe H. Hünenberger,et al.  Calculation of the group-based pressure in molecular simulations. I. A general formulation including Ewald and particle-particle-particle-mesh electrostatics , 2002 .

[22]  H. G. Petersen.  Accuracy and efficiency of the particle mesh Ewald method , 1995 .

[23]  T. Schlick,et al.  New splitting formulations for lattice summations , 2001 .

[24]  Christian Holm,et al.  Estimate of the Cutoff Errors in the Ewald Summation for Dipolar Systems , 2001 .

[25]  Ramzi Kutteh,et al.  Efficient dipole iteration in polarizable charged systems using the cell multipole method and application to polarizable water , 1995 .

[26]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[27]  Hans-Jörg Limbach,et al.  ESPResSo - an extensible simulation package for research on soft matter systems , 2006, Comput. Phys. Commun..

[28]  J. Weis,et al.  The structure of ferrofluids: A status report , 2005 .

[29]  T. Darden,et al.  Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions , 2000 .

[30]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[31]  Henrik Gordon Petersen,et al.  On the fast multipole method for computing the energy of periodic assemblies of charged and dipolar particles , 1993 .

[32]  Christian Holm,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998 .

[33]  Philippe H. Hünenberger,et al.  Optimal charge-shaping functions for the particle–particle—particle–mesh (P3M) method for computing electrostatic interactions in molecular simulations , 2000 .

[34]  J. Perram,et al.  Cutoff Errors in the Ewald Summation Formulae for Point Charge Systems , 1992 .

[35]  Stefan Odenbach,et al.  Ferrofluids: Magnetically Controllable Fluids And Their Applications , 2010 .

[36]  C. Holm,et al.  Microstructure analysis of monodisperse ferrofluid monolayers: theory and simulation. , 2008, Physical chemistry chemical physics : PCCP.

[37]  Stefan Odenbach,et al.  Magnetoviscous effects in ferrofluids , 2002 .

[38]  V. F. Medvedev,et al.  Magnetic Fluids: Engineering Applications , 1993 .