论文标题
从头算有效的单电子电位运营商。 I.充电转移能量在有效碎片电位中的应用
Ab Initio Effective One-Electron Potential Operators. I. Applications for Charge-Transfer Energy in Effective Fragment Potentials
论文作者
论文摘要
有效的单电子电位(OEP)的概念在有效地描述化学系统的电子结构(尤其是扩展的分子聚集体,例如凝结相中的相互作用分子)方面已经有用。在此,提出了一种基于碎片的计算方法(例如第二代有效片段势(EFP2)方法)调整为基于电子排斥积分(ERI)的有效消除电子排斥积分(ERI)的一般方法。 OEP运算符矩阵元件的两种一般类型通过分布式多极扩展或在本工作中开发的扩展密度拟合方案进行区分和处理。然后,将OEP技术应用于在许多应用程序中使用不完整的EFP2设置在相互作用能量和分子动力学模拟中使用不完整的问题,这是由于评估电荷转移(CT)效应的计算成本相对较高而与其他效果相比。基于OEP的替代CT能量模型是在与Hartree-Fock非相互作用的无相相的气相参考波函数的背景下提出的,与EFP2公式兼容。发现与EFP2方案中的CT能量方法相比,计算成本可以降低20倍,而不会损害各种弱相互作用的中性分子复合物的准确性。因此,人们认为所提出的模型可以在EFP2框架中使用,从而使CT能量项不再是基于EFP2的复杂系统模拟中的瓶颈。
The concept of the effective one-electron potentials (OEP) has been useful for many decades in efficient description of electronic structure of chemical systems, especially extended molecular aggregates such as interacting molecules in condensed phases. Here, a general method for effective OEP-based elimination of electron repulsion integrals (ERI), that is tuned towards the fragment-based calculation methodologies such as the second generation of the effective fragment potentials (EFP2) method, is presented. Two general types of the OEP operator matrix elements are distinguished and treated either via the distributed multipole expansion or the extended density fitting schemes developed in this work. The OEP technique is then applied to address the problem of using incomplete EFP2 settings in many applications in interaction energy and molecular dynamics simulations due to relatively high computational cost of evaluating the charge transfer (CT) effects as compared to other effects. The alternative OEP-based CT energy model is proposed in the context of the intermolecular perturbation theory with Hartree-Fock non-interacting gas-phase reference wavefunctions, compatible with the EFP2 formulation. It is found that the computational cost can be reduced up to 20 times as compared to the CT energy method within the EFP2 scheme without compromising the accuracy for a wide range of weakly interacting neutral molecular complexes. Therefore, it is believed that the proposed model can be used within the EFP2 framework, making the CT energy term no longer the bottleneck in EFP2-based simulations of complex systems.