DeYonker Nathan J, Cundari Thomas R, Wilson Angela K
Center for Advanced Scientific Computing and Modeling (CASCaM), Department of Chemistry, The University of North Texas, Denton, Texas 76203, USA.
J Chem Phys. 2006 Mar 21;124(11):114104. doi: 10.1063/1.2173988.
An alternative to the Gaussian-n (G1, G2, and G3) composite methods of computing molecular energies is proposed and is named the "correlation consistent composite approach" (ccCA, ccCA-CBS-1, ccCA-CBS-2). This approach uses the correlation consistent polarized valence (cc-pVXZ) basis sets. The G2-1 test set of 48 enthalpies of formation (DeltaHf), 38 adiabatic ionization potentials (IPs), 25 adiabatic electron affinities (EAs), and 8 adiabatic proton affinities (PAs) are computed using this approach, as well as the DeltaHf values of 30 more systems. Equilibrium molecular geometries and vibrational frequencies are obtained using B3LYP density functional theory. When applying the ccCA-CBS method with the cc-pVXZ series of basis sets augmented with diffuse functions, mean absolute deviations within the G2-1 test set compared to experiment are 1.33 kcal mol(-1) for DeltaHf,0.81 kcal mol(-1) for IPs, 1.02 kcal mol(-1) for EAs, and 1.51 kcal mol(-1) for PAs, without including the "high-level correction" (HLC) contained in the original Gn methods. Whereas the HLC originated in the Gaussian-1 method as an isogyric correction, it evolved into a fitted parameter that minimized the error of the composite methods, eliminating its physical meaning. Recomputing the G1 and G3 enthalpies of formation without the HLC reveals a systematic trend where most DeltaHf values are significantly higher than experimental values. By extrapolating electronic energies to the complete basis set (CBS) limit and adding G3-like corrections for the core-valence and infinite-order electron correlation effects, ccCA-CBS-2 often underestimates the experimental DeltaHf, especially for larger systems. This is desired as inclusion of relativistic and atomic spin-orbit effects subsequently improves theoretical DeltaHf values to give a 0.81 kcal mol(-1) mean absolute deviation with ccCA-CBS-2. The ccCA-CBS method is a viable "black box" method that can be used on systems with at least 10-15 heavy atoms.
本文提出了一种计算分子能量的方法,作为高斯-n(G1、G2和G3)复合方法的替代方案,并将其命名为“相关一致复合方法”(ccCA、ccCA-CBS-1、ccCA-CBS-2)。该方法使用相关一致极化价(cc-pVXZ)基组。使用此方法计算了包含48个生成焓(ΔHf)、38个绝热电离势(IP)、25个绝热电子亲和能(EA)和8个绝热质子亲和能(PA)的G2-1测试集,以及另外30个体系的ΔHf值。使用B3LYP密度泛函理论获得平衡分子几何结构和振动频率。当将ccCA-CBS方法与添加了弥散函数的cc-pVXZ系列基组一起应用时,与实验相比,G2-1测试集中ΔHf的平均绝对偏差为1.33 kcal mol⁻¹,IP为0.81 kcal mol⁻¹,EA为1.02 kcal mol⁻¹,PA为1.51 kcal mol⁻¹,不包括原始Gn方法中包含的“高级校正”(HLC)。虽然HLC最初在高斯-1方法中作为等旋校正出现,但它演变成一个拟合参数,使复合方法的误差最小化,从而消除了其物理意义。在没有HLC的情况下重新计算G1和G3生成焓,发现了一个系统趋势,即大多数ΔHf值明显高于实验值。通过将电子能量外推到完全基组(CBS)极限,并添加类似G3的核心价和无穷阶电子相关效应校正,ccCA-CBS-2通常会低估实验ΔHf,特别是对于较大的体系。这是理想的,因为随后包含相对论和原子自旋轨道效应会提高理论ΔHf值,使ccCA-CBS-2的平均绝对偏差为0.81 kcal mol⁻¹。ccCA-CBS方法是一种可行的“黑箱”方法,可用于至少含有10-15个重原子的体系。