Incremental update of electrostatic interactions in adaptively restrained particle simulations

Semeho Prince A. Edorh 1 Stephane Redon 1
1 NANO-D - Algorithms for Modeling and Simulation of Nanosystems
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, INPG - Institut National Polytechnique de Grenoble
Abstract : The computation of long-range potentials is one of the demanding tasks in Molecular Dynamics. During the last decades, an inventive panoply of methods was developed in order to reduce the CPU time of this task. In this work, we propose a fast method dedicated to the computation of the elec-trostatic potential in adaptively restrained systems. We exploit the fact that, in such systems, only some particles are allowed to move at each timestep. We developed an incremental algorithm derived from a multigrid-based alternative to traditional Fourier-based methods. Our algorithm was implemented inside LAMMPS, a popular molecular dynamics simulation package. We evaluated the method on different systems. We showed that the new algorithm's computational complexity scales with the number of active particles in the simulated system, and is able to outperform the well-established Particle Particle Particle Mesh (P3M) for adaptively restrained simulations.
Document type :
Journal articles
Complete list of metadatas

Cited literature [20 references]  Display  Hide  Download
Contributor : Nano-D Equipe <>
Submitted on : Monday, April 9, 2018 - 3:35:50 PM
Last modification on : Friday, November 9, 2018 - 2:16:04 PM


main_JCC (1).pdf
Files produced by the author(s)




Semeho Prince A. Edorh, Stephane Redon. Incremental update of electrostatic interactions in adaptively restrained particle simulations. Journal of Computational Chemistry, Wiley, 2018, 39 (20), pp.1455-1469. ⟨10.1002/jcc.25215⟩. ⟨hal-01761906⟩



Record views


Files downloads