P3M algorithm for dipolar interactions
Abstract
An extension to the P3M 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 P3M 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.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- December 2008
- DOI:
- 10.1063/1.3000389
- arXiv:
- arXiv:0805.4783
- Bibcode:
- 2008JChPh.129w4104C
- Keywords:
-
- 71.15.Pd;
- 41.20.Cv;
- 02.70.Ns;
- Molecular dynamics calculations and other numerical simulations;
- Electrostatics;
- Poisson and Laplace equations boundary-value problems;
- Molecular dynamics and particle methods;
- Physics - Computational Physics
- E-Print:
- 56 pages, 7 figures