Logo image
Scalable evaluation of polarization energy and associated forces in polarizable molecular dynamics: II. Toward massively parallel computations using smooth particle mesh Ewald
Journal article   Open access   Peer reviewed

Scalable evaluation of polarization energy and associated forces in polarizable molecular dynamics: II. Toward massively parallel computations using smooth particle mesh Ewald

Louis Lagardère, Filippo Lipparini, Étienne Polack, Benjamin Stamm, Éric Cancès, Michael Schnieders, Pengyu Ren, Yvon Maday and Jean-Philip Piquemal
Journal of chemical theory and computation, Vol.11(6), pp.2589-2599
06/09/2015
DOI: 10.1021/acs.jctc.5b00171
PMID: 26575557
url
http://hdl.handle.net/11568/872932View
Open Access

Abstract

In this article, we present a parallel implementation of point dipole-based polarizable force fields for molecular dynamics (MD) simulations with periodic boundary conditions (PBC). The smooth particle mesh Ewald technique is combined with two optimal iterative strategies, namely, a preconditioned conjugate gradient solver and a Jacobi solver in conjunction with the direct inversion in the iterative subspace for convergence acceleration, to solve the polarization equations. We show that both solvers exhibit very good parallel performances and overall very competitive timings in an energy and force computation needed to perform a MD step. Various tests on large systems are provided in the context of the polarizable AMOEBA force field as implemented in the newly developed Tinker-HP package, which is the first implementation of a polarizable model that makes large-scale experiments for massively parallel PBC point dipole models possible. We show that using a large number of cores offers a significant acceleration of the overall process involving the iterative methods within the context of SPME and a noticeable improvement of the memory management, giving access to very large systems (hundreds of thousands of atoms) as the algorithm naturally distributes the data on different cores. Coupled with advanced MD techniques, gains ranging from 2 to 3 orders of magnitude in time are now possible compared to nonoptimized, sequential implementations, giving new directions for polarizable molecular dynamics with periodic boundary conditions using massively parallel implementations.
Particle Size Surface Properties Molecular Dynamics Simulation

Details

Metrics

55 readers on Mendeley
1 readers on CiteULike
Logo image