Abstract:We present an efficient polarizable electrostatic model, utilizing typed, atom-centered, polarizabilities and the fast direct approximation, designed for efficient use in molecular dynamics (MD) simulations. The model provides two convenient approaches to assigning partial charges in the context of the atomic polarizabilities. One is a generalization of RESP, called RESP-dPol, and the other, AM1-BCC-dPol, is an adaptation of the widely used AM1-BCC method. Both are designed to accurately replicate gas-phase QM electrostatic potentials. Benchmarks of this polarizable electrostatic model against gas-phase dipole moments, molecular polarizabilities, bulk liquid densities, and static dielectric constants of organic liquids, show good agreement with the reference values. Of note, the model yields markedly more accurate dielectric constants of organic liquids, relative to a matched non-polarizable force field. MD simulations with this method, which is currently parameterized for molecules containing elements C, N, O, and H, run about only 3.6-fold slower than fixed charge force fields, while simulations with the self-consistent mutual polarization average 4.5-fold slower. Our results suggest that RESP-dPol and AM1-BCC-dPol afford improved accuracy, relative to fixed charge force fields, and are good starting points for developing general, affordable, and transferable polarizable force fields. The software implementing these approaches has been designed to utilize the force field fitting frameworks developed and maintained by Open Force Field Initiative, setting the stage for further exploration of this approach to polarizable force field development.

Response Theory for Static and Dynamic Polarizabilities of Excited States

Excited State Polarizabilities in Solution Obtained by Cubic Response Theory: Calculations on Para-, Ortho-, and Meta-Nitroaniline

Ground- and Excited-State Hyperpolarizabilities of Cis-, Trans- and Diphenyl-Polyenes

Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study

Standard Response Calculation of Electric Dipole Polarizability and Specific Optical Rotation Power Densities

Excess Polarizabilities Upon The First Dipole-Allowed Excitation Of Some Conjugated Oligomers

Excess Polarizabilities Upon Excitation from the Ground State to the First Dipole‐allowed Excited State of Diphenylpolyenes

Linear Response Properties of Solvated Systems: A Computational Study

Accurate Molecular Polarizabilities Based on Continuum Electrostatics

Calculation of Dynamical Response Functions Using a Bound-state Method

A fast, convenient, polarizable electrostatic model for molecular dynamics

Pseudo linear-dependence and long-range interaction effects on the polarizability and hyperpolarizabilities of stereoregular polymers

Analytic response theory for the density matrix renormalization group

Robust and accurate computational estimation of the polarizability tensors of macromolecules

High-precision nonadiabatic calculations of dynamic polarizabilities and hyperpolarizabilities for the lowlying vibrational-rotational states of hydrogen molecular ions

Electrodynamic Response And Stability Of Molecular Crystals

Polaritonic response theory for exact and approximate wave functions

Dielectric Relaxation of Interacting/polarizable Polar Molecules with Linear Reaction Dynamics in a Weak Alternating Field

Nonlinear dielectric relaxation of polar molecules with linear reaction dynamics in self-consistent field approximation

Electric Polarization from Many-Body Neural Network Ansatz

Variational Perturbation Theory for Dynamic Polarizabilities and Dispersion Coefficients.