Research Group of Theoretical Chemistry and Molecular Modeling, Hasselt University, Agoralaan Gebouw D, Diepenbeek B-3590, Belgium.
J Chem Phys. 2009 Dec 14;131(22):224321. doi: 10.1063/1.3270190.
A benchmark theoretical study of the electronic ground state and of the vertical and adiabatic singlet-triplet (ST) excitation energies of benzene (n=1) and n-acenes (C(4n+2)H(2n+4)) ranging from naphthalene (n=2) to heptacene (n=7) is presented, on the ground of single- and multireference calculations based on restricted or unrestricted zero-order wave functions. High-level and large scale treatments of electronic correlation in the ground state are found to be necessary for compensating giant but unphysical symmetry-breaking effects in unrestricted single-reference treatments. The composition of multiconfigurational wave functions, the topologies of natural orbitals in symmetry-unrestricted CASSCF calculations, the T1 diagnostics of coupled cluster theory, and further energy-based criteria demonstrate that all investigated systems exhibit a (1)A(g) singlet closed-shell electronic ground state. Singlet-triplet (S(0)-T(1)) energy gaps can therefore be very accurately determined by applying the principles of a focal point analysis onto the results of a series of single-point and symmetry-restricted calculations employing correlation consistent cc-pVXZ basis sets (X=D, T, Q, 5) and single-reference methods [HF, MP2, MP3, MP4SDQ, CCSD, CCSD(T)] of improving quality. According to our best estimates, which amount to a dual extrapolation of energy differences to the level of coupled cluster theory including single, double, and perturbative estimates of connected triple excitations [CCSD(T)] in the limit of an asymptotically complete basis set (cc-pVinfinityZ), the S(0)-T(1) vertical excitation energies of benzene (n=1) and n-acenes (n=2-7) amount to 100.79, 76.28, 56.97, 40.69, 31.51, 22.96, and 18.16 kcal/mol, respectively. Values of 87.02, 62.87, 46.22, 32.23, 24.19, 16.79, and 12.56 kcal/mol are correspondingly obtained at the CCSD(T)/cc-pVinfinityZ level for the S(0)-T(1) adiabatic excitation energies, upon including B3LYP/cc-PVTZ corrections for zero-point vibrational energies. In line with the absence of Peierls distortions, extrapolations of results indicate a vanishingly small S(0)-T(1) energy gap of 0 to approximately 4 kcal/mol (approximately 0.17 eV) in the limit of an infinitely large polyacene.
本文基于基于限制或非限制零阶波函数的单参考和多参考计算,对苯(n=1)和 n-并苯(C(4n+2)H(2n+4))的电子基态以及垂直和绝热单重态-三重态(ST)激发能进行了基准理论研究,范围从萘(n=2)到苝(n=7)。发现对于非限制单参考处理中存在的巨大但不合理的对称性破坏效应,必须进行高水平和大规模的电子相关处理。多组态波函数的组成、非限制 CASSCF 计算中自然轨道的拓扑结构、耦合簇理论的 T1 诊断以及进一步的基于能量的标准表明,所有研究的系统都表现出(1)A(g)单重态闭壳电子基态。因此,通过将焦点分析原理应用于一系列单点和对称限制计算的结果,可以非常准确地确定单重态-三重态(S(0)-T(1))能隙,这些计算使用相关一致的 cc-pVXZ 基组(X=D、T、Q、5)和单参考方法[HF、MP2、MP3、MP4SDQ、CCSD、CCSD(T)]进行质量改进。根据我们的最佳估计,这相当于耦合簇理论能量差异的双重外推,包括渐近完全基组(cc-pVinfinityZ)中单个、双个和连接三重激发的微扰估计[CCSD(T)],苯(n=1)和 n-并苯(n=2-7)的 S(0)-T(1)垂直激发能分别为 100.79、76.28、56.97、40.69、31.51、22.96 和 18.16 kcal/mol。在 CCSD(T)/cc-pVinfinityZ 水平上,对于 S(0)-T(1)绝热激发能,相应地得到 87.02、62.87、46.22、32.23、24.19、16.79 和 12.56 kcal/mol 的值,其中包括 B3LYP/cc-PVTZ 对零点振动能的校正。与不存在 Peierls 畸变一致,结果的外推表明,在无穷大聚并苯的极限中,S(0)-T(1)能隙几乎为零,约为 4 kcal/mol(约 0.17 eV)。