Blake Thomas A, Xantheas Sotiris S
William R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, PO Box 999, MS K8-88, Richland, Washington 99352, USA.
J Phys Chem A. 2006 Sep 7;110(35):10487-94. doi: 10.1021/jp062472r.
We present the results of high level ab initio calculations for the structure, harmonic and anharmonic spectroscopic constants, and ring puckering barrier of cyclobutane (C4H8) in an effort to establish the minimum theoretical requirements needed for their accurate description. We have found that accurate estimates for the barrier between the minimum (D(2d)) and transition state (D(4h)) configurations require both higher levels of electron correlation [MP4, CCSD(T)] and orbital basis sets of quadruple-zeta quality or larger. By performing CCSD(T) calculations with basis sets as large as cc-pV5Z, we were able to obtain, for the first time, a value for the puckering barrier that lies within 10 cm(-1) (or 2%) from experiment, whereas the best previously calculated values were in errors exceeding 40% of experiment. Our best estimate of 498 cm(-1) for the puckering barrier is within 10 cm(-1) of the experimental value proposed originally, but it lies approximately 50 cm(-1) higher than the revisited value, which was obtained more recently using different assumptions regarding the coupling between the various modes. It is therefore suggested that revisiting the analysis of the experimental data might be warranted. Our best computed values (at the CCSD(T)/aug-cc-pVTZ level of theory) for the equilibrium structural parameters of C4H8 are r(C-C) = 1.554 A, r(C-H(alpha)) = 1.093 A, r(C-H(beta)) = 1.091 A, phi(C-C-C) = 88.1 degrees , alpha(H(alpha)-C-H(beta)) = 109.15 degrees , and theta = 29.68 degrees for the puckering angle. We have found that the puckering angle theta is more sensitive to the level of electron correlation than to the size of the basis set for a given method. We furthermore present anharmonic calculations that are based on a second-order perturbative evaluation of rovibrational parameters and their effects on the vibrational spectra and average structure. We have found that the anharmonic calculations predict the experimentally measured fundamental band origins within 1% (< or =30 cm(-1)) for most vibrations. The results of the current study can serve as a guide for future calculations on the substituted four-member ring hydrocarbon compounds. To this end we present a method for estimating the puckering barrier height at higher levels of electron correlation [MP4, CCSD(T)] from the MP2 results that can be used in chemically similar compounds.
我们展示了对环丁烷(C₄H₈)的结构、谐波和非谐波光谱常数以及环皱曲势垒进行的高水平从头算计算结果,旨在确定准确描述它们所需的最低理论要求。我们发现,要准确估计最低构型(D₂d)和过渡态构型(D₄h)之间的势垒,既需要更高水平的电子相关方法[MP4、CCSD(T)],也需要四重zeta质量或更大的轨道基组。通过使用高达cc-pV5Z的基组进行CCSD(T)计算,我们首次获得了一个皱曲势垒值,该值与实验值的偏差在10 cm⁻¹(或2%)以内,而此前计算的最佳值误差超过实验值的40%。我们对皱曲势垒的最佳估计值为498 cm⁻¹,与最初提出的实验值偏差在10 cm⁻¹以内,但比最近使用关于各种模式间耦合的不同假设重新获得的值大约高50 cm⁻¹。因此,建议重新审视实验数据的分析。我们对C₄H₈平衡结构参数的最佳计算值(在CCSD(T)/aug-cc-pVTZ理论水平下)为:r(C - C) = 1.554 Å,r(C - H(α)) = 1.093 Å,r(C - H(β)) = 1.091 Å,φ(C - C - C) = 88.1°,α(H(α) - C - H(β)) = 109.15°,皱曲角θ = 29.68°。我们发现,对于给定方法,皱曲角θ对电子相关水平的敏感度高于对基组大小的敏感度。此外,我们还展示了基于振动转动参数的二阶微扰评估及其对振动光谱和平均结构影响的非谐波计算。我们发现,对于大多数振动,非谐波计算预测的实验测量基频带起源误差在1%以内(≤30 cm⁻¹)。本研究结果可为未来对取代四元环烃化合物的计算提供指导。为此,我们提出了一种从MP2结果估计更高电子相关水平[MP4、CCSD(T)]下皱曲势垒高度的方法,该方法可用于化学性质相似的化合物。