香豆素衍生物的荧光发射能计算及XC泛函的合理选择

采用含时密度泛函理论(TDDFT)与单激发组态相互作用(CIS)处理相结合的计算方案对香豆素系列15种已知荧光化合物的发射能进行了系统考察. 结果表明, 发射能与吸收能一样, 其计算值主要取决于交换-相关(XC)泛函的选择. 只要泛函选用得当, 在使用较小基组的TDDFT/6-31G(d)//CIS/3-21G(d)理论水平上即可使绝大部分化合物的实验发射能在精度达0.16 eV以内得以重现. 与吸收能计算不同的是, 无法选用单一的一种泛函来对全系列化合物的发射能作出满意的理论预测. 激发态无明显电荷转移的、7位上有氨(或胺)基取代或有氮原子相连的化合物, 其适用泛函为不含Hartree-Fock(HF)交换能的纯泛函OLYP和BLYP. 而激发态发生较大程度电荷转移的、3 位上有共轭取代基的衍生物, 其适用泛函则为含20%的HF交换成分的混合泛函B3LYP. 因此, 发射能计算中的XC泛函选择, 应同时考虑取代基团效应以及激发态的电子结构特征. 其中, 发射能计算值受XC泛函中HF交换能比例的影响十分敏感. 文中还对激发能计算中的溶剂效应校正方案和激发态几何优化精度的影响进行了讨论.     

Electronic Structure Characteristics of ESIPT and TICT Fluorescence Emissions and Calculations of Emitting Energies

A scheme of time-dependent density functional theory (TDDFT) combined with the single-excitation configuration interaction (CIS) approach was employed to investigate the first excited singlet state (S₁) for eight salicylanilide derivatives and analogues, which have similar structural formulas. The results showed that fluorescence-emitting mechanisms of these molecules were in two distinct manners (excited-state intramolecular proton transfer (ESIPT) and twisted intramolecular charge transfer (TICT)), which agreed with the well-known experiments. For ESIPT compounds with inconspicuous charge transfer (CT) during electron transition, pure functionals without Hartree-Fock (HF) exchange energy, such as OLYP and BLYP, were suitable to calculate emitting energies. For TICT compounds with large CT during electron transition, hybrid functionals with about 37% HF exchange energy, such as mPW1B95 and MPW1K, performed well. On condition that the exchange-correlation (XC) functionals were chosen properly according to the rules above, reliable emitting energies for salicylanilide derivatives and analogues could be obtained at the TDDFT/6-31G(d)//CIS/3-21G(d) level. The average accuracy reached about 0.2 eV. For the salicylanilides with both proton transfer (PT) and CT reaction channels, only one channel occurred actually according to the principle of energy minimum. This actual reaction decided proper XC functionals, whereas the reaction that did not occur actually was trivial. Eight appendent compounds were calculated to prove that this successful scheme is expected to be suitable for other ESIPT and TICT compounds.     

ESIPT和TICT荧光发射的电子结构特征及发射能计算

采用含时密度泛函理论(TDDFT)与单激发组态相互作用(CIS)相结合的计算方案对八种结构相似的水杨酰苯胺衍生物及其类似物第一激发单重态(S1)进行考察, 证实它们的荧光发射分属分子内质子转移(ESIPT)和分子扭转-电荷转移(TICT)两种不同机制且结论与已知实验事实相符. ESIPT发光的化合物在电子跃迁前后无明显的电荷转移发生, 发射能计算的适用泛函是OLYP和BLYP等无Hartree-Fock(HF)交换成分的纯泛函; TICT发光的化合物在电子跃迁前后发生明显的电荷转移, 其适用泛函为含约37% HF交换成分的混合型泛函(例如mPW1B95和MPW1K). 按上述原则来选择适用泛函, 即可在TDDFT/6-31G(d)//CIS/3-21G(d)理论水平上正确预测水杨酰苯胺衍生物和类似物的发射能, 平均精度可达0.2 eV. 兼具质子转移与电荷转移双反应通道的化合物, 两者的竞争遵从能量最小原理, 结果使荧光发射仅选择其中一个通道进行. 泛函的选择只与实际发生的反应有关, 与并未实际发生的反应通道无关. 附加的八个算例进一步表明, 此成功的计算方案可望推广应用于其它类型的ESIPT和TICT荧光有机物.     

Assessment of the OLYP and O3LYP density functionals for first-row transition metals

We have investigated the performance of the OLYP and O3LYP density functionals for predicting atomic excitation energies and ionization potentials, and bond dissociation energies, geometries, and vibrational frequencies for selected first-row transition metal compounds, including hydrides (MH) and singly charged methylene and methyl cations. The OLYP and O3LYP functionals are similar to the well-known BLYP and B3LYP functionals, respectively, but use a new optimized exchange functional (OPTX) developed by Handy and Cohen (Mol Phys 2001, 99, 403) in place of the standard B88 exchange. A previous study by us on organic reactions (J Chem Phys 2002, 117, 1331) indicated that both OLYP and O3LYP gave results for heats of reaction and barrier heights that were overall superior to those using the popular B3LYP functional. For transition metals, however, although OLYP is overall superior to BLYP for molecular calculations, it is inferior to B3LYP. O3LYP provides results for molecules of about the same quality as B3LYP. For atomic excitation and 4s ionization energies, unless relativistic effects are included, OLYP and O3LYP are clearly worse than both BLYP and B3LYP. There is thus no real incentive to use either OLYP or O3LYP in place of B3LYP for calculations involving first-row transition metals.     

A TDDFT study of the low-lying excitation energies of polycyclic cinnolines and their carbocyclic analogues

The low-lying excited states of 21 compounds of polycyclic cinnoline, monoaza-hydrocarbon and their corresponding polycyclic aromatic hydrocarbon (PAH) analogues have been investigated in the framework of time-dependent density functional theory (TDDFT). The gradient corrected BLYP, hybrid B3LYP and B3P86 functionals together with 6-31G (d) basis set have been used. The hybrid-type B3LYP and B3P86 systematically overestimate the excitation energies for states with dominating ionic character (corresponding to group III band), with a mean absolute deviation (MAD) of 0.33 eV (B3LYP) and 0.34 eV (B3P86), respectively. However, they can accurately predict the excitation energies for states with covalent character (corresponding to group II and I bands). The MAD for group II (B3LYP: 0.05 eV; B3P86: 0.05 eV) and I bands (B3LYP: 0.12 eV; B3P86: 0.12 eV) are significantly smaller. The BLYP outperforms B3LYP and B3P86 for group III band (MAD: 0.09 eV), but has a worse performance for group II (MAD: 0.15 eV) and I (MAD: 0.13 eV) bands. Comparison of the lowest-lying excited states for polycyclic cinnolines with those of the corresponding PAH analogues, the first excited states of polycyclic cinnolines mainly result from n→π₀* transitions. Therefore, in non-polar solvents, the spectra of some polycyclic cinnolines exhibit an additional absorption band at longer wavelengths.     

Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes

It is well-known that time-dependent density functional theory (TDDFT) yields substantial errors for the excitation energies of charge-transfer (CT) excited states, when approximate standard exchange-correlation (xc) functionals are used, for example, SVWN, BLYP, or B3LYP. Also, the correct 1/R asymptotic behavior of CT states with respect to a distance coordinate R between the separated charges of the CT state is not reproduced by TDDFT employing these xc-functionals. Here, we demonstrate by analysis of the TDDFT equations that the first failure is due to the self-interaction error in the orbital energies from the ground-state DFT calculation, while the latter is a similar self-interaction error in TDDFT arising through the electron transfer in the CT state. Possible correction schemes, such as inclusion of exact Hartree-Fock or exact Kohn-Sham exchange, as well as aspects of the exact xc-functional are discussed in this context. Furthermore, a practical approach is proposed which combines the benefits of TDDFT and configuration interaction singles (CIS) and which does not suffer from electron-transfer self-interaction. The latter approach is applied to a (1,4)-phenylene-linked zincbacteriochlorin-bacteriochlorin complex and to a bacteriochlorophyll-spheroidene complex, in which CT states may play important roles in energy and electron-transfer processes. The errors of TDDFT alone for the CT states are demonstrated, and reasonable estimates for the true excitation energies of these states are given.     

Ab initio and DFT benchmark study for nucleophilic substitution at carbon (SN2@C) and silicon (SN2@Si)

To obtain a set of consistent benchmark potential energy surfaces (PES) for the two archetypal nucleophilic substitution reactions of the chloride anion at carbon in chloromethane (S(N)2@C) and at silicon in chlorosilane (S(N)2@Si), we have explored these PESes using a hierarchical series of ab initio methods [HF, MP2, MP4SDQ, CCSD, CCSD(T)] in combination with a hierarchical series of six Gaussian-type basis sets, up to g polarization. Relative energies of stationary points are converged to within 0.01 to 0.56 kcal/mol as a function of the basis-set size. Our best estimate, at CCSD(T)/aug-cc-pVQZ, for the relative energies of the [Cl(-), CH(3)Cl] reactant complex, the [Cl-CH(3)-Cl](-) transition state and the stable [Cl-SiH(3)-Cl](-) transition complex is -10.42, +2.52, and -27.10 kcal/mol, respectively. Furthermore, we have investigated the performance for these reactions of four popular density functionals, namely, BP86, BLYP, B3LYP, and OLYP, in combination with a large doubly polarized Slater-type basis set of triple-zeta quality (TZ2P). Best overall agreement with our CCSD(T)/aug-cc-pVQZ benchmark is obtained with OLYP and B3LYP. However, OLYP performs better for the S(N)2@C overall and central barriers, which it underestimates by 2.65 and 4.05 kcal/mol, respectively. The other DFT approaches underestimate these barriers by some 4.8 (B3LYP) to 9.0 kcal/mol (BLYP).     

Low-lying excited states of C120 and C151: a multireference perturbation theory study

Excited states of two 7-aminocoumarin derivatives, coumarin 120 (7-amino-4-methylcoumarin) and coumarin 151 (7-amino-4-trifluoromethylcoumarin), were investigated using generalized multiconfigurational quasidegenerate perturbation theory (GMC-QDPT), multiconfigurational quasidegenerate perturbation theory (MC-QDPT) and time-dependent density functional theory (TDDFT) with the B3LYP and CAM-B3LYP functionals. The absorption and fluorescence spectra of C120 and C151 were calculated. We elucidated the characters of the low-lying states of C120 and C151. The absorption spectra calculated with GMC-QDPT and TDDFT B3LYP agreed well with the experimental data, while for the fluorescence spectra, the TDDFT calculations overestimated the fluorescence spectra compared to GMC-QDPT calculations. Utilizing active spaces with large numbers of electrons and orbitals for reference functions, GMC-QDPT showed a better performance than MC-QDPT with a complete active space self-consistent field (CASSCF) reference of active space with smaller number of electrons and orbitals. In our gas phase calculation, we found that the optimized structures for the first excited states have a planar amino group with a CN single bond, while the amino group is pyramidal in the ground state.     

Electronic structures of 4d transition metal monoxides by density functional theory

Bond distances, dissociation energies, ionization potentials and electron affinities of 4d transition metal monoxides from YO to CdO and their positive and negative ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, SVWN, MPW1PW91 and PBE1PBE. It was found that calculated properties are highly dependent on the functionals employed, especially for dissociation energy. For most neutral species, pure density functionals BLYP, BPW91 and BP86 have good performance in predicting dissociation energy than hybrid density functionals B3LYP, B3PW91 and B3P86. In addition, BLYP gives the largest bond distance compared with other density functional methods, while SVWN gives shortest bond distance, largest dissociation energy and electron affinity. For the ground state, the spin multiplicity of the charged species can be obtained by ± 1 of their corresponding neutral species.     

Scaling dynamical correlation energy from density functional theory correlation functionals

The scaling of dynamical correlation energy in molecules obtained by the correlation functionals of density functional theory (DFT) is examined. The approach taken is very similar to the scaled external correlation method of Brown and Truhlar but is based on the observation that DFT correlation functionals, especially the LYP, appear to represent the dynamical portion of the correlation energy in molecules. We examine whether higher accuracy in atomization energies can be gained by scaling without significant deterioration of the structural and spectroscopic properties of the molecules using four DFT functionals (BLYP, OLYP, B3LYP, and O3LYP) on 19 molecules including the six molecule AE6 database, the latter being representative of a much larger, 109 molecule training set. We show that, with molecule specific scale factors, nearly perfect agreement with experiment can be achieved in atomization energies without increasing the average errors in other molecular properties relative to the DFT calculation. We further show that it is possible to find optimal scale factors which reduce the mean unsigned error per bond to levels comparable to those of some multilevel multicoefficient methods.     

硅杂环戊二烯衍生物电子结构与光谱性质的理论研究

采用密度泛函理论(DFT)B3LYP/6-31G*方法, 对2,5位取代的硅杂环戊二烯(Silole)系列衍生物进行几何结构优化, 通过计算得到电离能、电子亲和势、空穴抽取能及电子抽取能等相关能量, 并使用TD-DFT方法研究其吸收光谱, 分析相关能量及光谱的变化规律. 采用单组态相互作用(CIS/6-31G*)方法优化得到它们的最低单重激发态(S1)结构, 在此基础上, 使用TD-DFT方法计算对应的发射光谱. 分析2,5位芳基取代硅杂环戊二烯衍生物(DADPS)激发态与基态的结构差异及原因, 研究前线分子轨道的分布情况, 并讨论发光特征及载流子传输性能. 研究结果表明, 激发态结构弛豫主要发生在Silole环和直接与2,5位芳基相连的部位; 前线轨道主要分布在Silole环和2,5位芳基上; 二吡咯取代物有望成为空穴传输材料, 二噻吩取代物和二呋喃取代物有望在发光器件中表现出较高的发光效率.     

DFT calculations on the spin-crossover complex Fe(salen)(NO): a quest for the best functional

DFT calculations on the spin-crossover complex Fe(salen)(NO) provide a striking illustration of the comparative performance of different exchange-correlation functionals vis-à-vis the issue of transition metal spin state energetics. Thus, although the "classic" pure functionals PW91 and BLYP favor the S = 1/2 state by about 10 kcal/mol, relative to the S = 3/2 state, the hybrid functional B3LYP favors the latter state by nearly the same margin. In contrast, the newer pure functionals OLYP and OPBE, based on the OPTX exchange functional, as well as the B3LYP* hybrid functional (which has 15% Hartree-Fock exchange, compared with 20% for B3LYP) predict nearly isoenergetic S = 1/2 and 3/2 states, as required for a spin-crossover complex. Intriguingly, the OLYP and B3LYP* spin density profiles for the S = 1/2 state of Fe(salen)(NO) are substantially dissimilar.     

Theoretical study on infrared vibrational spectra of boric-acid in gas-phase using density functional methods

Density functional theory (DFT) methods with various exchange-correlation functionals such as SVWN, BVWN, BVWN5, BLYP, B1LYP, B3LYP, B3PW91, and BH and H are employed in a theoretical study of molecular boric-acid in gas-phase. In the calculations, the split valence 6-311++G** and 6-31G* basis sets were used. The geometry, zero-point vibrational energies (ZPVEs), and harmonic infrared vibrational (IR) frequencies are predicted. The calculated C_3h-symmetry geometrical parameters are compared with Hartree–Fock (HF) calculation results and experimental data. IR frequencies predicted by the BLYP, B3LYP, and B3PW91 calculations are in good agreement with experimental data. The frequency calculations presented here also suggest that the C_3h-symmetrical structure corresponds to a minimum in the potential energy surface, but neither is D_3h- or C₃-symmetrical structure.     

Assessing alkyl-, silyl-, and halo-substituent effects on the electron affinities of silyl radicals

Neutral anion energy differences for a large class of alpha-substituted silyl radicals have been computed to determine the effect of alkyl, silyl, and halo substituents on their electron affinities. In particular, we report theoretical predictions of the adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) for a series of methyl-, silyl-, and halo-substituted silyl radical compounds. This work utilizes the carefully calibrated DZP++ basis set, in conjunction with the pure BLYP and OLYP functionals, as well as with the hybrid B3LYP, BHLYP, PBE1PBE, MPW1K, and O3LYP functionals. Bromine has the largest effect in stabilizing the anions, and the BLYP/DZP++ AEA for SiBr(3) is 3.29 eV. The other predicted electron affinities are for SiH(3) (1.37 eV), SiH(2)CH(3) (1.09 eV), SiH(2)F (1.54 eV), SiH(2)Cl (1.94 eV), SiH(2)Br (2.05 eV), SiH(2)(SiH(3)) (1.77 eV), SiH(CH(3))(2) (0.92 eV), SiHF(2) (1.86 eV), SiHCl(2) (2.53 eV), SiHBr(2) (2.67 eV), Si(CH(3))(3) (0.86 eV), SiF(3) (2.66 eV), SiCl(3) (3.21 eV), Si(SiH(3))(3) (2.25 eV), and SiFClBr (3.13 eV). For the five silyl radicals where experimental data are available, the BLYP functional gives the most accurate determination of AEAs; the average absolute error is 0.04(1) eV, whereas the corresponding errors for the O3LYP, MPW1K, PBE1PBE, B3LYP, OLYP, and BHLYP functionals are 0.05(8), 0.06(0), 0.06(3), 0.08(5), 0.11(5), and 0.15(3) eV, respectively.     

Optical absorption and emission properties of fluoranthene, benzo[k]fluoranthene, and their derivatives. A DFT study

Fluoranthene and benzo[k]fluoranthene-based oligoarenes are good candidates for organic light-emitting diodes (OLEDs). In this work, the electronic structure and optical properties of fluoranthene, benzo[k]fluoranthene, and their derivatives have been studied using quantum chemical methods. The ground-state structures were optimized using the density functional theory (DFT) methods. The lowest singlet excited state was optimized using time-dependent density functional theory (TD-B3LYP) and configuration interaction singles (CIS) methods. On the basis of ground- and excited-state geometries, the absorption and emission spectra have been calculated using the TD-DFT method with a variety of exchange-correlation functionals. All the calculations were carried out in chloroform medium. The results show that the absorption and emission spectra calculated using the B3LYP functional is in good agreement with the available experimental results. Unlikely, the meta hybrid functionals such as M06HF and M062X underestimate the absorption and emission spectra of all the studied molecules. The calculated absorption and emission wavelength are more or less basis set independent. It has been observed that the substitution of an aromatic ring significantly alters the absorption and emission spectra.     

Computational Assessment of MLCT versus MC Stabilities in First-to-Third-Row d6 Pseudo-Octahedral Transition Metal Complexes

An accurate modeling of metal-to-ligand-charge-transfer (MLCT) and metal-centered (MC) excited state energies is key to predict the photoinduced response in transition metal complexes (TMCs). Herein, the importance of the ground state and excited state reference geometries is addressed for three-prototype d⁶ pseudo-octahedral TMCs, each displaying a different potential energy landscape of MLCT versus MC relative stabilities. Several functionals are used within the time-dependent density functional theory (TDDFT), as well as multireference wave-function theory (MS-CASPT2), applied to [Mn(im)(CO)₃(phen)]⁺, [Ru(im)₂(bpy)₂]²⁺, and [Re(im)(CO)₃(phen)]⁺, (im: imidazole, bpy: bypiridine, phen: phenantroline). The results revel that TDDFT is robust except when using B3LYP functional for first-row d⁶ TMCs. In contrast, MS-CASPT2 calculations are strongly biased in those cases with competitive MLCT/MC states. The results reinforce the reliability of B3LYP to describe the excited states in d⁶ TMCs, but question the validity of assessing the density functional theory (DFT)/TDDFT performance via direct comparison with MS-CASPT2 performed at the same DFT reference geometry as a standard strategy. © 2019 Wiley Periodicals, Inc.     

Charge-transfer excitations in uranyl tetrachloride ([UO2Cl4]2-): how reliable are electronic spectra from relativistic time-dependent density functional theory?

Four-component relativistic time-dependent density functional theory (TD-DFT) is used to study charge-transfer (CT) excitation energies of the uranyl molecule as well as the uranyl tetrachloride complex. Adiabatic excitation energies and vibrational frequencies of the excited states are calculated for the lower energy range of the spectrum. The results for TD-DFT with the CAM-B3LYP exchange-correlation functional for the [UO(2)Cl(4)](2-) system are in good agreement with the experimentally observed spectrum of this species and agree also rather well with other theoretical data. Use of the global hybrid B3LYP gives qualitatively correct results, while use of the BLYP functional yields results that are qualitatively wrong due to the too low CT states calculated with this functional. The applicability of the overlap diagnostic of Peach et al. (J. Chem. Phys.2008, 128, 044118) to identify such CT excitations is investigated for a wide range of vertical transitions using results obtained with three different approximate exchange-correlation functionals: BLYP, B3LYP, and CAM-B3LYP.     

DFT/TDDFT studies of the geometry, electronic structure and spectra of (12S)-1,4,7,10-tetraazadicyclo[10,3,0]-pentadecane-3,11-dione and its derivatives

The FT-Raman and UV-visible spectra of (12S)-1,4,7,10-tetraazadicyclo[10,3,0]-pentadecane-3,11-dione and its derivatives were obtained and discussed. The harmonic vibrational wavenumbers and the corresponding Raman scattering activities in their electronic ground-states were calculated at the DFT-B3LYP/6-31G(d) level of theory. The calculated wavenumbers were then scaled and compared with the experimental values. The 7-(2,4-dinitrophenyl)-(12S)-1,4,7,10-tetrazadicyclo[10,3,0]-pentadecane-3,11-dione derivative has mainly an amide (II) character, while the others have an amide (I) character. Moreover, the different substituents do not cause a significant shift of the vibrational mode of the macrocyclic plane. The electronic vertical excitation energy and the oscillator strength were determined with the help of TDDFT calculations and by employing pure (BLYP) and hybrid (B3LYP, B3P86, and mPW1PW91) functionals together with the 6-31G(d) basis set. The BLYP functional reproduces the UV-vis absorption spectra better than the B3LYP, B3P86, or mPW1PW91 hybrid functionals. A dimolecular model, which considers hydrogen-bonded structures, proved that strong inter- and intramolecular hydrogen bonds are present in these compounds. Due to the transannular effect, the UV-vis absorption spectrum of macrocyclic dioxotetraamines is completely different from that of single amide compounds.     

Ground and excited states of the monomer and dimer of certain carboxylic acids

The ground-state properties of the monomer and the dimer of formic acid, acetic acid, and benzoic acid have been investigated using Hartree-Fock (HF) and density functional theory (DFT) methods using the 6-311++G(d,p) basis set. Some of the low-lying excited states have been studied using the time-dependent density functional theory (TDDFT) with LDA and B3LYP functionals and also employing complete-active-space-self-consistent-field (CASSCF) and multireference configuration interaction (MRCI) methodologies. DFT calculations predict the ground-state geometries in quantitative agreement with the available experimental results. The computed binding energies for the three carboxylic acid dimers are also in accord with the known thermodynamic data. The TDDFT predicted wavelengths corresponding to the lowest energy n-pi* transition in formic acid (214 nm) and acetic acid (214 nm) and the pi-pi* transition in benzoic acid (255 nm) are comparable to the experimentally observed absorption maxima. In addition, TDDFT calculations predict qualitatively correctly the blue shift (4-5 nm) in the excitation energy for the pi-pi* transition in going from the monomer to the dimer of formic acid and acetic acid and the red shift (approximately 19 nm) in pi-pi* transition in going from benzoic acid monomer to dimer. This also indicates that the electronic interaction arising from the hydrogen bonds between the monomers is marginal in all three carboxylic acids investigated.