Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA.
J Chem Phys. 2012 Jun 28;136(24):244107. doi: 10.1063/1.4729535.
A method for calculating the dispersion energy between molecules modeled with the general effective fragment potential (EFP2) method and those modeled using a full quantum mechanics (QM) method, e.g., Hartree-Fock (HF) or second-order perturbation theory, is presented. C(6) dispersion coefficients are calculated for pairs of orbitals using dynamic polarizabilities from the EFP2 portion, and dipole integrals and orbital energies from the QM portion of the system. Dividing by the sixth power of the distance between localized molecular orbital centroids yields the first term in the commonly employed London series expansion. A C(8) term is estimated from the C(6) term to achieve closer agreement with symmetry adapted perturbation theory values. Two damping functions for the dispersion energy are evaluated. By using terms that are already computed during an ordinary HF or EFP2 calculation, the new method enables accurate and extremely rapid evaluation of the dispersion interaction between EFP2 and QM molecules.
本文提出了一种计算用广义有效片段势能(EFP2)方法建模的分子与用全量子力学(QM)方法建模的分子(例如 Hartree-Fock(HF)或二级微扰理论)之间的色散能的方法。使用 EFP2 部分的动态极化率以及系统 QM 部分的偶极积分和轨道能来计算轨道对的 C(6)色散系数。除以局域分子轨道质心之间距离的六次方得到通常采用的伦敦级数展开的第一项。从 C(6)项估计 C(8)项,以实现与自适应对称微扰理论值更接近的吻合。评估了两种用于色散能的阻尼函数。通过使用在普通 HF 或 EFP2 计算过程中已经计算的项,新方法能够精确且极其快速地评估 EFP2 和 QM 分子之间的色散相互作用。