Thomas Heike B, Hennemann Matthias, Kibies Patrick, Hoffgaard Franziska, Güssregen Stefan, Hessler Gerhard, Kast Stefan M, Clark Timothy
Computer-Chemie-Centrum, Department Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg , Nägelsbachstr. 25, 91052 Erlangen, Germany.
Physikalische Chemie III, Technische Universität Dortmund , Otto-Hahn-Str. 4a, 44227 Dortmund, Germany.
J Chem Inf Model. 2017 Aug 28;57(8):1907-1922. doi: 10.1021/acs.jcim.7b00080. Epub 2017 Aug 2.
A neglect of diatomic differential overlap (NDDO) Hamiltonian has been parametrized as an electronic component of a polarizable force field. Coulomb and exchange potentials derived directly from the NDDO Hamiltonian in principle can be used with classical potentials, thus forming the basis for a new generation of efficiently applicable multipolar polarizable force fields. The new hpCADD Hamiltonian uses force-field-like atom types and reproduces the electrostatic properties (dipole moment, molecular electrostatic potential) and Koopmans' theorem ionization potentials closely, as demonstrated for a large training set and an independent test set of small molecules. The Hamiltonian is not intended to reproduce geometries or total energies well, as these will be controlled by the classical force-field potentials. In order to establish the hpCADD Hamiltonian as an electronic component in force-field-based calculations, we tested its performance in combination with the 3D reference interaction site model (3D RISM) for aqueous solutions. Comparison of the resulting solvation free energies for the training and test sets to atomic charges derived from standard procedures, exact solute-solvent electrostatics based on high-level quantum-chemical reference data, and established semiempirical Hamiltonians demonstrates the advantages of the hpCADD parametrization.
忽略双原子微分重叠(NDDO)哈密顿量已被参数化为可极化力场的电子成分。原则上,直接从NDDO哈密顿量导出的库仑势和交换势可与经典势一起使用,从而构成新一代高效适用的多极可极化力场的基础。新的hpCADD哈密顿量使用类似力场的原子类型,并能紧密再现静电性质(偶极矩、分子静电势)和库普曼定理电离势,这在大量小分子训练集和独立测试集上得到了证明。该哈密顿量并非旨在很好地再现几何结构或总能量,因为这些将由经典力场势控制。为了将hpCADD哈密顿量确立为基于力场计算中的电子成分,我们测试了它与用于水溶液的三维参考相互作用位点模型(3D RISM)相结合时的性能。将训练集和测试集所得的溶剂化自由能与从标准程序导出的原子电荷、基于高水平量子化学参考数据的精确溶质 - 溶剂静电作用以及已确立的半经验哈密顿量进行比较,证明了hpCADD参数化的优势。