Logo image
Polarizable Molecular Dynamics in a Polarizable Continuum Solvent
Journal article   Open access   Peer reviewed

Polarizable Molecular Dynamics in a Polarizable Continuum Solvent

Filippo Lipparini, Louis Lagardère, Christophe Raynaud, Benjamin Stamm, Eric Cancès, Benedetta Mennucci, Michael Schnieders, Pengju Ren, Yvon Maday and Jean-Philip Piquemal
Journal of chemical theory and computation, Vol.11(2), pp.623–634-634
01/12/2015
DOI: 10.1021/ct500998q
PMCID: PMC4621208
PMID: 26516318
url
http://hdl.handle.net/11568/724064View
Open Access

Abstract

We present, for the first time, scalable polarizable molecular dynamics (MD) simulations within a polarizable continuum solvent with molecular shape cavities and exact solution of the mutual polarization. The key ingredients are a very efficient algorithm for solving the equations associated with the polarizable continuum, in particular, the domain decomposition Conductor-like Screening Model (ddCOSMO), which involves a rigorous coupling of the continuum with the polarizable force field achieved through a robust variational formulation and an effective strategy to solve the coupled equations. The coupling of ddCOSMO with nonvariational force fields, including AMOEBA, is also addressed. The MD simulations are feasible, for real-life systems, on standard cluster nodes; a scalable parallel implementation allows for further acceleration in the context of a newly developed module in Tinker, named Tinker-HP. NVE simulations are stable, and long-term energy conservation can be achieved. This paper is focused on the methodological developments, the analysis of the algorithm, and the stability of the simulations; a proof-of-concept application is also presented to attest to the possibilities of this newly developed technique.
Chemical Sciences Other

Details

Metrics

Logo image