TD-CI simulation of the electronic optical response of molecules in intense fields II: comparison of DFT functionals and EOM-CCSD

Time-dependent configuration interaction (TD-CI) simulations can be used to simulate molecules in intense laser fields. TD-CI calculations use the excitation energies and transition dipoles calculated in the absence of a field. The EOM-CCSD method provides a good estimate of the field-free excited states but is rather expensive. Linear-response time-dependent density functional theory (TD-DFT) is an inexpensive alternative for computing the field-free excitation energies and transition dipoles needed for TD-CI simulations. Linear-response TD-DFT calculations were carried out with standard functionals (B3LYP, BH&HLYP, HSE2PBE (HSE03), BLYP, PBE, PW91, and TPSS) and long-range corrected functionals (LC-ωPBE, ωB97XD, CAM-B3LYP, LC-BLYP, LC-PBE, LC-PW91, and LC-TPSS). These calculations used the 6-31G(d,p) basis set augmented with three sets of diffuse sp functions on each heavy atom. Butadiene was employed as a test case, and 500 excited states were calculated with each functional. Standard functionals yield average excitation energies that are significantly lower than the EOM-CC, while long-range corrected functionals tend to produce average excitation energies slightly higher. Long-range corrected functionals also yield transition dipoles that are somewhat larger than EOM-CC on average. The TD-CI simulations were carried out with a three-cycle Gaussian pulse (ω = 0.06 au, 760 nm) with intensities up to 1.26 × 10(14) W cm(-2) directed along the vector connecting the end carbons. The nonlinear response as indicated by the residual populations of the excited states after the pulse is far too large with standard functionals, primarily because the excitation energies are too low. The LC-ωPBE, LC-PBE, LC-PW91, and LC-TPSS long-range corrected functionals produce responses comparable to EOM-CC.     

Performance of recently developed kinetic energy density functionals for the calculation of hydrogen binding strengths and hydrogen-bonded structures

 The accuracy of predicted hydrogen binding energies and equilibrium structures for a benchmark set of molecules is compared for some recently developed density functionals, Becke's three parameter hybrid method with the Lee, Yang, and Parr (LYP) correlation functional (B3LYP), Becke's half and half functional combined with the LYP correlation functional (BHLYP), Perdew, Burke and Ernzerhof functional (PBE), Van Voorhis, Scuseria exchange correlation functional (VSXC), the hybrid Perdew, Burke and Ernzerhof functional (PBE1PBE), and meta-generalized gradient approximation (meta-GGA). Overall, the hybrid functionals which contain a portion of Hartree–Fock exchange (B3LYP, BHLYP, and PBE1PBE) yield the most accurate results. The kinetic-energy-density-dependent functionals, VSXC and meta-GGA, are significantly less accurate. Received: 10 December 1999 / Accepted: 5 March 2000 / Published online: 21 June 2000     

DFT and TDDFT study related to electron transfer in nonbonded porphine...C60 complexes

Spectroscopic properties of a ground state nonbonded porphine-buckminsterfullerene (H2P...C60) complex are studied in several different relative orientations of C60 with respect to the porphine plane by using the density functional (DFT) and time-dependent density functional (TDDFT) theories. The geometries and electronic structures of the ground states are optimized with the B3LYP and PBE functionals and a SVP basis set. Excitation energies and oscillator strengths are obtained from the TDDFT calculations. The relative orientation of C60 is found to affect the equilibrium distance between H2P and C60 especially in the case of the PBE functional. The excitation energies of different H2P...C60 complexes are found to be practically the same for the same excitations when the B3LYP functional is used but to differ notably when PBE is used in calculations. Existence of the states related to a photoinduced electron transfer within a porphyrin-fullerene dyad is also studied. All calculations predict a formation of an excited charge-transfer complex state, a locally excited donor (porphine) state, as well as a locally excited acceptor (fullerene) state in the investigated H2P...C60 complexes.     

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.     

Density functional study of AuXq (X = O, S, Se, Te, q = +1, 0, -1) molecules

Bond distances, vibrational frequencies, electron affinity, ionization potential, and dissociation energies of the title molecules were studied by use of density functional methods B3LYP, B3P86, B3PW91, BHLYP, BLYP, BP86, mPW1PW91, and PBE1PBE. It was found that the ground electronic state is doublet for neutral species, singlet for the anion, and triplet for the cation, in agreement with experiments and previous theoretical studies. The calculated properties are highly dependent on the functionals employed, in particular for the dissociation energy. The predicted bond distances and vibrational frequencies are in agreement with experiments and previous theoretical results. BP86 and BLYP have relatively good performance in reproducing the experimental results, while BHLYP is the worst functional method compared with the other density functional methods used for the title molecules.     

The electronic structure and energetics of V(+)-benzene half-sandwiches of different multiplicities: Comparative multireference and single-reference theoretical study

Multireference [complete active space self-consistent field (CASSCF) and multiconfigurational quasidegenerate perturbation theory (MCQDPT)] and single-reference ab initio (Moller-Plesset second order perturbation theory (MP2) and coupled clusters with singles, doubles and noniterative triples [CCSD(T)]) and density functional theory (PBE and B3LYP) electronic structure calculations of V(C(6)H(6))(+) half-sandwich in the states of different multiplicities are described and compared. Detailed analyses of the geometries and electronic structures of the all found states are given; adiabatic and diabatic dissociation energies are estimated. The lowest electronic state of V(C(6)H(6))(+) half-sandwich was found to be the quintet (5)B(2) state with a slightly deformed upside-down-boat-shaped benzene ring and d(4) configuration of V atom, followed by a triplet (3)A(2) state lying about 4 kcal/mol above. The lowest singlet state (1)A(1)(d(4)) lies much ( approximately 28 kcal/mol) higher. MCQDPT calculated adiabatic dissociation energy (53.6 kcal/mol) for the lowest (5)B(2)(d(4)) state agrees well with the current 56.4 (54.4) kcal/mol experimental estimate, giving a preference to the lower one. Compared to MCQDPT, B3LYP hybrid exchange-correlation functional provides the best results, while CCSD(T) performs usually worse. Gradient-corrected PBE calculations tend to systematically overestimate metal-benzene binding in the row quintet相似文献   

Designing new free-radical reducing reagents: Theoretical study on Si–H bond dissociation energies of organic silanes

Bond dissociation energies of a series of substituted silanes were studied with the density functional theory methods. The performances of six different density functional methods including B3LYP, B3P86, BH&HLYP, B1LYP, PBE1KCIS, and TPSSLYP1W were examined for the prediction of Si–H bond dissociation energies. The results showed that B3P86 was the most accurate theoretical procedure among these six DFT methods. Using the B3P86 method, we then carried out a systematic study about the substituent effects on Si–H bond dissociation energies, with a focus to identify the possible approaches to weaken the Si–H bond strength. On the basis of the knowledge learned from the systematic study on model systems, we proposed some new silicon-based radical reducing reagents which may be used to replace toxic tin hydride reagents.     

Characterizing the potential energy surface of the water dimer with DFT: failures of some popular functionals for hydrogen bonding

Ten stationary points on the water dimer potential energy surface have been examined with ten density functional methods (X3LYP, B3LYP, B971, B98, MPWLYP, PBE1PBE, PBE, MPW1K, B3P86, and BHandHLYP). Geometry optimizations and vibrational frequency calculations were carried out with the TZ2P(f,d)+dif basis set. All ten of the density functionals correctly describe the relative energies of the ten stationary points. However, correctly describing the curvature of the potential energy surface is far more difficult. Only one functional (BHandHLYP) reproduces the number of imaginary frequencies from CCSD(T) calculations. The other nine density functionals fail to correctly characterize the nature of at least one of the ten (H(2)O)(2) stationary points studied here.     

Structure and stability of the (001) alpha-quartz surface

The structure and surface energies of the cleaved, reconstructed, and fully hydroxylated (001) alpha-quartz surface of various thicknesses are investigated with periodic density functional theory (DFT). The properties of the cleaved and hydroxylated surface are reproduced with a slab thickness of 18 atomic layers, while a thicker 27-layer slab is necessary for the reconstructed surface. The performance of the hybrid DFT functional B3LYP, using an atomic basis set, is compared with the generalised gradient approximation, PBE, employing plane waves. Both methodologies give similar structures and surface energies for the cleaved and reconstructed surfaces, which validates studying these surfaces with hybrid DFT. However, there is a slight difference between the PBE and B3LYP approach for the geometry of the hydrogen bonded network on the hydroxylated surface. The PBE adsorption energy of CO on a surface silanol site is in good agreement with experimental values, suggesting that this method is more accurate for hydrogen bonded structures than B3LYP. New hybrid functionals, however, yield improved weak interactions. Since these functionals also give superior activation energies, we recommend applying the new functionals to contemporary issues involving the silica surface and adsorbates on this surface.     

Electronic spectroscopy of UO2(2+), NUO(+) and NUN: an evaluation of time-dependent density functional theory for actinides

The performance of the time-dependent density functional theory (TDDFT) approach has been evaluated for the electronic spectrum of the UO(2)(2+), NUO(+) and NUN molecules. Different exchange-correlation functionals (LDA, PBE, BLYP, B3LYP, PBE0, M06, M06-L, M06-2X, CAM-B3LYP) and the SAOP model potential have been investigated, as has the relative importance of the adiabatic local density approximation (ALDA) to the exchange-correlation kernel. The vertical excitation energies have been compared with reference data obtained using accurate wave-function theory (WFT) methods.     

Density functional theory studies on the electronic and vibrational spectra of octaethylporphyrin diacid

The ground-state structure and electronic and vibrational spectra of octaethylporphyrin diacid (H4OEP2+) have been studied with the density functional theory. The geometrical parameters computed with B3LYP, PBE1PBE and mPW1PW91 functionals and 6-31G* basis sets are well consistent with the experimental values. Electronic absorption spectrum of H4OEP2+ has been studied with the time-dependent DFT method, and the calculated excitation energies and oscillator strengths are compared with the experimental results. The Raman and IR spectra of H4OEP2+ and the Raman spectrum of its N-deuterated analogue (D4OEP2+) were measured. The observed Raman and IR bands have been assigned based on the frequency calculations at the B3LYP/6-31G* level of theory.     

Absorption and Fluorescence Emission Spectra of Molecules Containing Azo and/or Oxadiazole Units

The ground state structures of a series of organic molecules containing azo and/or oxadiazole units were obtained by means of density functional theory B3LYP/6-31G(d) method. The first singlet excited state structures were optimized by virtue of singlet-excitation configuration interaction CIS/6-31G(d) method. The absorption and fluorescence emission spectra were then evaluated via the time-dependent density functional theory B3LYP and PBE1PBE methods with 6-311++G(3df,2p) basis set. The calculation results show that compared with those of their parent molecules A-H, B-H, C-H, D-H, the absorption and emission wavelength values of all the derivatives show red shifts. The derivatives containing both the oxadiazole and methoxyl units are good candidates for longer absorption wavelength materials. The effects of azo, oxadiazole, and methoxyl units on the absorption and emission wavelength were discussed.     

Calculations of ionization energies and electron affinities for atoms and molecules: A comparative study with different methods

In the present work, we examined the performance of 36 density functionals, including the newly developed doubly hybrid density functional XYG3 (Y. Zhang, X. Xu, and W. A. Goddard III, Proc. Natl. Acad. Sci, USA, 2009, 106, 4963), to calculate ionization energies (IEs) and electron affinities (EAs). We used the well-established G2-1 set as reference, which contains 14 atoms and 24 molecules for IE, along with 7 atoms and 18 molecules for EA. XYG3 leads to mean absolute deviations (MADs) of 0.057 and 0.080 eV for IEs and EAs, respectively, using the basis set of 6-311 + G (3df,2p). In comparison with some other functionals, MADs for IEs are 0.109 (B2PLYP), 0.119 (M06-2X), 0.159 (X3LYP), 0.161 (PBE), 0.162 (B3LYP), 0.165 (PBE0), 0.173 (TPSS), 0.200 (BLYP), and 0.215 eV (LC-BLYP). MADs for EAs are 0.090 (X3LYP), 0.090 (B2PLYP), 0.102 (PBE), 0.103 (M06-2X), 0.104 (TPSS), 0.105 (BLYP), 0.106 (B3LYP), 0.126 (LC-BLYP), and 0.128 eV (PBE0).     

XYG3型双杂化密度泛函方法新进展：从能量到能量的解析梯度

本文回顾了现代密度泛函理论的基础,着重评述了XYG3型双杂化（XYG3typeofdoublyhybrid,xDH）泛函的最新进展,解析能量梯度的实现．XYG3是首个依照绝热途径理论建立的双杂化泛函,在具体实现上具有独特的构架．该类型泛函利用常用泛函（如B3LYP或PBE0等）作母泛函来进行自洽计算,以期获得更好的密度和轨道,然后将所得到的轨道和密度信息带入到xDH泛函中以得到最终能量．由于自洽泛函和最终能量泛函不同,因而在计算解析能量梯度时需要求解耦合微扰Kohn—Sham方程．在此基础上,还评述了xDH泛函在能量,尤其是构型优化方面的具体表现．测试的构型集包括以共价键键合的分子和非键相互作用体系的平衡结构,以及反应过渡态结构．结果表明,xDH双杂化泛函总体上给出了比母泛函更好的能量和几何构型．     

Potential Function and Thermodynamic Property of UO (cited: 2)

Potential energy scan for uranium oxide (UO) was performed by ab initio configuration inter-action (CI) method and density functional theory methods at the PBE1 and the B3LYP levels in combination with the (ECP80MWB_AVQZ+2f) basis set for uranium and 6-311+G* foroxygen. The dissociation energies of UO, after being corrected for the zero-point vibrational energy, are 2.38, 3.76, and 3.31 eV at the CI, PBE1, and B3LYP levels, respectively. The calculated energy was fitted to potential functions of Morse, Lennard-Jones, and Rydberg. Only the Morse function is eligible for the potential. The anharmonicity constant is 0.00425. The anharmonic frequency is 540.95 cm^-1 deduced from the PBE1 results. Thermodynamic properties of entropy and heat capacity at 298.2-1500 K were calculated using DFT-UPBE1 results and Morse parameters. The relationship between entropy and temperature was es-tablished.     

Nickel-catalyzed alkyl coupling reactions: evaluation of computational methods

The B3LYP, M06, M06L, M062X, MPW1K, and PBE1PBE DFT methods were evaluated for modeling nickel-catalyzed coupling reactions. The reaction consists of a nucleophilic attack by a carbanion equivalent on the nickel complex, S(N)2 attack by the anionic nickel complex on an alkyl halide, and reductive elimination of the coupled alkane product, regenerating the nickel catalyst. On the basis of CCSD(T)//DFT single-point energies, the B3LYP, M06, and PBE1PBE functionals were judged to generate the best ground state geometries. M06 energies are generally comparable or superior to B3LYP and PBE1PBE energies for transition state calculations. The MP2 and CCSD methods were also evaluated for single-point energies at the M06 geometries. The rate-determining step of this reaction was found to be nucleophilic attack of a L(2)NiR anion on the alkyl halide.     

Benchmark calculations on the adiabatic ionization potentials of M-NH(3) (M=Na,Al,Ga,In,Cu,Ag)

The ground states of the M-NH(3) (M=Na,Al,Ga,In,Cu,Ag) complexes and their cations have been studied with density functional theory and coupled cluster [CCSD(T)] methods. The adiabatic ionization potentials (AIPs) of these complexes are calculated, and these are compared to results from high-resolution zero-electron kinetic energy photoelectron spectroscopy. By extrapolating the CCSD(T) energies to the complete basis set (CBS) limit and including the core-valence, scalar relativistic, spin-orbit, and zero-point corrections, the CCSD(T) method is shown to be able to predict the AIPs of these complexes to better than 6 meV or 0.15 kcal/mol. 27 exchange-correlation functionals, including one in the local density approximation, 13 in the generalized gradient approximation (GGA), and 13 with hybrid GGAs, were benchmarked in the calculations of the AIPs. The B1B95, mPW1PW91, B98, B97-1, PBE1PBE, O3LYP, TPSSh, and HCTH93 functionals give an average error of 0.1 eV for all the complexes studied, with the B98 functional alone yielding a maximum error of 0.1 eV. In addition, the calculated metal-ammonia harmonic stretching frequencies with the CCSD(T) method are in excellent agreement with their experimental values, whereas the B3LYP method tends to underestimate these stretching frequencies. The metal-ammonia binding energies were also calculated at the CCSD(T)/CBS level, and are in excellent agreement with the available experimental values considering the error limits, except for Ag-NH(3) and Ag(+)-NH(3), where the calculations predict stronger bond energies than measured by about 4 kcal/mol, just outside the experimental error bars of +/-3 kcal/mol.     

X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。