Atom-centered electric multipole moments dynamically generated from QM/MM MD simulations

Abstract

Atom-centered electric multipole moments can be extremely useful in chemistry as they enable the systematic mapping of a complex electrostatic problem to a simpler model. However, since they do not correspond to physical observables, there is no unique way to define them. In this work, we present an extension of the dynamically generated RESP charges (D-RESP) method, referred to as xDRESP, where atom-centered multipoles are computed from mixed quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. We compare the ability of xDRESP charges to reproduce the electrostatic potential, as well as molecular multipoles, against the performance of fixed point-charge models commonly used in force fields. Moreover, we highlight cases where DRESP atomic multipoles can provide valuable information about chemical systems, such as indicating when polarization plays a significant role, and chemical reactions, in which xDRESP atomic multipoles can be used as an on-the-fly analysis tool to track changes in electron density.

References in corpus (6)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Tinker-HP : Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields using GPUs and Multi-GPUs systems
• Analytic energy, gradient, and hessian of electrostatic embedding QM/MM based on electrostatic potential fitted atomic charges scaling linearly with the MM subsystem size
• MiMiC: A High-Performance Framework for Multiscale Molecular Dynamics Simulations
• Making Puzzle Pieces Fit or Reshaping MiMiC for Multiscale Simulations with CP2K and More
• OpenMM-MiMiC Interface for Efficient and Flexible Multiscale Simulations