English   español  
Please use this identifier to cite or link to this item: http://hdl.handle.net/10261/45988
logo share SHARE logo core CORE   Add this article to your Mendeley library MendeleyBASE

Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL
Exportar a otros formatos:


Particle-particle particle-mesh method for dipolar interactions:on error estimates and efficiency of schemes with analytical differentiation and mesh interlacing

AuthorsCerdà, Joan J. ; Ballenegger, V.; Holm, Christian
Issue Date2011
PublisherAmerican Institute of Physics
CitationJournal of Chemical Physics 135: 184110 (1-13) (2011)
AbstractThe interlaced and non-interlaced versions of the dipolar particle-particle particle-mesh (P3M) method implemented using the analytic differentiation scheme (AD-P3M) are presented together with their respective error estimates for the calculation of the forces, torques, and energies. Expressions for the optimized lattice Green functions, and for the Madelung self-forces, self-torques and self-energies are given. The applicability of the theoretical error estimates are thoroughly tested and confirmed in several numerical examples. Our results show that the accuracy of the calculations can be improved substantially when the approximate (mesh computed) Madelung self-interactions are subtracted. Furthermore, we show that the interlaced dipolar AD-P3M method delivers a significantly higher accuracy (which corresponds approximately to using a twice finer mesh) than the conventional method, allowing thereby to reduce the mesh size with respect to the non-interlaced version for a given accuracy. In addition, we present similar expressions for the dipolar ik-differentiation interlaced scheme, and we perform a comparison with the AD interlaced scheme. Rough tests for the relative speed of the dipolar P3M method using ik-differentiation and the interlaced/non-interlaced AD schemes show that when FFT computing time is the bottleneck, usually when working at high precisions, the interlaced AD-scheme can be several times faster than the other two schemes. For calculations with a low accuracy requirement, the interlaced version can perform worse than the ik and the non-interlaced AD schemes.
Publisher version (URL)http://dx.doi.org/10.1063/1.3657407
Appears in Collections:(IFISC) Artículos
Files in This Item:
File Description SizeFormat 
JChemPhys_135_184110.pdf1 MBAdobe PDFThumbnail
Show full item record
Review this work

Related articles:

WARNING: Items in Digital.CSIC are protected by copyright, with all rights reserved, unless otherwise indicated.