Abstract:A set of atomic polarizability parameters for a new polarizable Gaussian model (pGM) has been developed with the goal to accurately reproduce the polarizability anisotropy, taking advantage of its ability to attenuate all short-range electrostatic interactions, by fitting the ab initio molecular polarizability tensors ( A <sub>pq</sub>) calculated at the B3LYP/aug-cc-pVTZ level. For comparison, we also rederived the parameters for three Thole models in which the 1-2 (bonded), 1-3 (separated by two bonds), and 1-4 (separated by three bonds) interactions are fully included. The average percent errors (APEs) of molecular polarizability tensors for 4842 molecules or dimers are 2.98, 3.76, 3.28, and 3.82% for the pGM, Thole linear, Thole exponential, and Thole Amoeba models, respectively, with atom-type independent, universal screening factors (USF). The APEs are reduced further to 2.30, 2.69, 2.25, and 2.48% for the four corresponding polarizable models with atom-type dependent, variable screening factors (VSF). It is encouraging that the pGM with variable screening factors achieved APEs of 1.83 for 1155 amino acid analogs, dipeptides, and tetrapeptides, 1.39 for 28 nucleic acid bases, 0.708 for 1464 water clusters, and 1.99 for 85 dimers of water and biological building blocks. Compared to the new set of models, the APEs of the old Thole models that were fitted to isotropic molecular polarizabilities are 8.7% for set A (without the 1-2 and 1-3 interactions) and 6.3% for set D (with the 1-2 and 1-3 interactions) models, respectively. MPAD, a metric of molecular polarization anisotropy difference based on the diagonal terms of molecular polarizability tensors was defined and applied to assess the polarizable models in reproducing the ab initio molecular polarization anisotropy. The MPADs are 3.71, 4.70, 4.11, and 4.77% for the pGM, Thole linear, Thole exponential, and Thole Amoeba USF models, respectively. The APEs are reduced further to 2.85, 3.58, 2.90, and 3.15% for the four corresponding VSF models. Thus, the new pGM and Thole models notably improve molecular polarization anisotropy. Since pGM attenuates all short-range electrostatic interactions, its application is expected to improve stability in charge fitting, energy, and force calculations and the accuracy of multibody polarization.

On combining Thole's induced point dipole model with fixed charge distributions in molecular mechanics force fields

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Atomic Polarizabilities for Interactive Dipole Induction Models

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A fast, convenient, polarizable electrostatic model for molecular dynamics

Gaussian induced dipole polarization model

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

An optimized charge penetration model for use with the AMOEBA force field.

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations

Accurate transferable polarization model derived from the monomer electron density

Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model

Development of Polarizable Gaussian Model for Molecular Mechanical Calculations I: Atomic Polarizability Parameterization To Reproduce ab Initio Anisotropy

From bonds to interactions: comprehensive molecular characterization via polarizable bond-dipole approach

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

Polarizable six-point water models from computational and empirical optimization

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

Accounting for electronic polarization in nonpolarizable force fields

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

Efficient polarizable QM/MM using the direct reaction field Hamiltonian with electrostatic potential fitted multipole operators