Dataset Information

Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential.

ABSTRACT: Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF in the current work explicitly depends on the local geometry (bond and angle) of the molecule. To our knowledge, this is the first study that derives energy and force expressions due to GDCF in a multipole-based polarizable FF framework. Due to the inclusion of GDCF, the AMOEBA+ water model is noticeably improved in terms of describing the monomer properties, cluster binding/interaction energy, and a variety of liquid properties, including the infrared spectra that previous flexible water models were not able to capture.

SUBMITTER: Liu C

PROVIDER: S-EPMC7384396 | biostudies-literature | 2020 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential.

Liu Chengwen C Piquemal Jean-Philip JP Ren Pengyu P

The journal of physical chemistry letters 20191230 2

Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF ...[more]

PMID: 31865706

Similar Datasets

Project description:The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.

Dataset Information

Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential.

Publications

Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure