Huang Hsing-Hsiang, Wang Yi-Siang, Chao Sheng D
Institute of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan R.O.C.
School of Chemistry & Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States.
ACS Omega. 2022 May 31;7(23):20059-20080. doi: 10.1021/acsomega.2c01888. eCollection 2022 Jun 14.
We extend our previous quantum chemistry calculations of interaction energies for 31 homodimers of small organic functional groups (the SOFG-31 data set) by including 239 heterodimers with monomers selected within the SOFG-31 data set, thus resulting in the SOFG-31+239 data set. The minimum-level theoretical scheme contains (1) the basis set superposition error corrected supermolecule (BSSE-SM) approach for intermolecular interactions; (2) the second-order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets for the geometry optimization and correlation energy calculations; and (3) the single-point energy calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods for the MP2 energy calibrations. In addition, we have performed a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT) in order to gain chemical insights. That the above procedure cannot be further reduced has been proven to be very crucial for constructing reliable data sets of interaction energies. The calculated CCSD(T)/CBS interaction energy data can serve as a benchmark for testing or training less accurate but more efficient calculation methods, such as the electronic density functional theory. As an application, we employ a segmental SAPT model previously developed for the SOFG-31 data set to predict binding energies of large heterodimer complexes. These model energy "quanta" can be used in coarse-grained molecular dynamics simulations by avoiding large-scale calculations.
我们扩展了之前对31个小有机官能团同二聚体(SOFG - 31数据集)相互作用能的量子化学计算,纳入了239个异二聚体,其单体选自SOFG - 31数据集,从而得到了SOFG - 31 + 239数据集。最低水平的理论方案包括:(1) 用于分子间相互作用的基组叠加误差校正超分子(BSSE - SM)方法;(2) 采用邓宁的aug - cc - pVXZ(X = D、T、Q)基组的二阶莫勒 - 普莱斯特定理微扰理论(MP2),用于几何优化和相关能计算;以及(3) 使用经过充分测试的MP2能量校准外推方法,在完整基组极限下采用单、双和微扰三激发耦合簇方法[CCSD(T)/CBS]进行单点能量计算。此外,我们基于对称适配微扰理论(SAPT)进行了一系列并行的能量分解计算,以获得化学见解。已证明上述步骤无法进一步简化对于构建可靠的相互作用能数据集至关重要。计算得到的CCSD(T)/CBS相互作用能数据可作为测试或训练不太精确但更高效的计算方法(如电子密度泛函理论)的基准。作为应用,我们采用先前为SOFG - 31数据集开发的分段SAPT模型来预测大型异二聚体复合物的结合能。这些模型能量“量子”可通过避免大规模计算用于粗粒度分子动力学模拟。