Miura Masanori, Aoki Yuriko, Champagne Benoît
Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, Kasuga, Fukuoka 816-8580, Japan.
J Chem Phys. 2007 Aug 28;127(8):084103. doi: 10.1063/1.2761886.
In present study the relevance of using the time-dependent density functional theory (DFT) within the adiabatic approximation for computing oscillator strengths (f) is assessed using different LDA, GGA, and hybrid exchange-correlation (XC) functionals. In particular, we focus on the lowest-energy valence excitations, dominating the UV/visible absorption spectra and originating from benzenelike HOMO(pi)-->LUMO(pi(*)) transitions, of several aromatic molecules: benzene, phenol, aniline, and fluorobenzene. The TDDFT values are compared to both experimental results obtained from gas phase measurements and to results determined using several ab initio schemes: random phase approximation (RPA), configuration interaction single (CIS), and a series of linear response coupled-cluster calculations, CCS, CC2, and CCSD. In particular, the effect of the amount of Hartree-Fock (HF) exchange in the functional is highlighted, whereas a basis set investigation demonstrates the need of including diffuse functions. So, the hybrid XC functionals--and particularly BHandHLYP--provide f values in good agreement with the highly correlated CCSD scheme while these can be strongly underestimated using pure DFT functionals. These results also display systematic behaviors: (i) larger f and squares of the transition dipole moments (mid R:mumid R:(2)) are associated with larger excitation energies (DeltaE); (ii) these relationships present generally a linear character with R>0.9 in least-squares fit procedures; (iii) larger amounts of HF exchange in the XC functional lead to larger f, R:mumid R:(2), as well as DeltaE values; (iv) these increases in f, mid R:mumid R:(2), and DeltaE are related to increased HOMO-LUMO character; and (v) these relationships are, however, not universal since the linear regression parameters (the slopes and the intercepts at the origin) depend on the system under investigation as well as on the nature of the excited state.
在本研究中,我们使用不同的局域密度近似(LDA)、广义梯度近似(GGA)和杂化交换相关(XC)泛函,评估了在绝热近似下使用含时密度泛函理论(DFT)计算振子强度(f)的相关性。特别地,我们关注了几种芳香族分子(苯、苯酚、苯胺和氟苯)的最低能量价激发,这些激发主导了紫外/可见吸收光谱,并且源自类似苯的最高占据分子轨道(HOMO(π))到最低未占据分子轨道(LUMO(π*))的跃迁。将含时密度泛函理论(TDDFT)值与从气相测量获得的实验结果以及使用几种从头算方法(随机相位近似(RPA)、单组态相互作用(CIS)以及一系列线性响应耦合簇计算方法,如CCSD、CC2和CC)确定的结果进行了比较。特别强调了泛函中哈特里 - 福克(HF)交换量的影响,而基组研究表明需要包含弥散函数。因此,杂化XC泛函——特别是BHandHLYP——提供的f值与高度相关的CCSD方法得到的值吻合良好,而使用纯DFT泛函时这些值可能会被严重低估。这些结果还呈现出系统的行为:(i)较大的f和跃迁偶极矩的平方(μ²)与较大的激发能(ΔE)相关;(ii)在最小二乘拟合过程中,这些关系通常呈现线性特征,相关系数R>0.9;(iii)XC泛函中HF交换量的增加会导致f、μ²以及ΔE值增大;(iv)f、μ²和ΔE的这些增加与HOMO - LUMO特征的增强有关;(v)然而,这些关系并非普遍适用,因为线性回归参数(斜率和原点处的截距)取决于所研究的体系以及激发态的性质。