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 efficiency in the calculations of vibrational frequencies and molecular structures of β-diketones

Density functional theory (DFT) levels are employed to calculate the vibrational frequencies and geometrical data of β-diketones. We evaluate the relative performance of the different levels by comparing theoretical results to experimental values. The applied DFT levels in this work are B3LYP, BLYP, B3P86, B3PW91, BPW91, G96LYP, BP86, and G96PW91 with the standard 6-31G, 6-31G*, 6-31G**, 6-31+G**, 6-31++G**, 6-311G**, 6-311++G** basis sets. The best results are obtained at the B3LYP, B3PW91, and B3P86 levels.     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

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.     

Hydrogen and dihydrogen bonding of transition metal hydrides

Intermolecular interactions between a prototypical transition metal hydride WH(CO)₂NO(PH₃)₂ and a small proton donor H₂O have been studied using DFT methodology. The hydride, nitrosyl and carbonyl ligand have been considered as site of protonation. Further, DFT-D calculations in which empirical corrections for the dispersion energy are included, have been carried out. A variety of pure and hybrid density functionals (BP86, PW91, PBE, BLYP, OLYP, B3LYP, B1PW91, PBE0, X3LYP) have been considered, and our calculations indicate the PBE functional and its hybrid variation are well suited for the calculation of transition metal hydride hydrogen and dihydrogen bonding. Dispersive interactions make up for a sizeable portion of the intermolecular interaction, and amount to 20–30% of the bond energy and to 30–40% of the bond enthalpy. An energy decomposition analysis reveals that the H?H bond of transition metal hydrides contains both covalent and electrostatic contributions.     

Accurate prediction of vertical electronic transitions of Ni(II) coordination compounds via time dependent density functional theory

Time dependent density functional theory calculations are completed for five Ni(II) complexes formed by polydentate peptides to predict the electronic absorption spectrum. The ligands examined were glycyl‐glycyl‐glycine (GGG), glycyl‐glycyl‐glycyl‐glycine (GGGG), glycyl‐glycyl‐histidine (GGH), glycyl‐glycyl‐cysteine (GGC), and triethylenetetramine (trien). Fifteen functionals and two basis sets were tested. On the basis of the mean absolute percent deviation (MAPD), the ranking among the functionals is: HSE06 ∼ MPW1PW91 ∼ PBE0 > ω‐B97x‐D ∼ B3P86 ∼ B3LYP ∼ CAM‐B3LYP > PBE ∼ BLYP ∼ BP86 > TPSS > TPSSh > BHandHLYP > M06 ≫ M06‐2X. Concerning the basis sets, the triple‐ζ def2‐TZVP performs better than the double‐ζ LANL2DZ. With the functional HSE06 and basis set def2‐TZVP the MAPD with respect to the experimental λ_max is 1.65% with a standard deviation of 1.26%. The absorption electronic spectra were interpreted in terms of vertical excitations between occupied and virtual MOs based on Ni‐d atomic orbitals. The electronic structure of the Ni(II) species is also discussed.     

Electronic Structures of 5d Transition Metal Monoxides by Density Functional Theory

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the diatomic 5d transition metal (except La) monoxides and their positively and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, MPW1PW91, PBE1PBE, and SVWN. Our calculation shows that for each individual species, the calculated properties are quite sensitive to the method used. Compared with hybrid density functional method B3PW91 (B3P86), pure density functional method BPW91 (BP86) gives longer bond distance (lower vibrational frequency) from HfO to PtO for neutral species, HfO⁺ to IrO⁺ for cationic species, and HfO⁻ to AuO⁻ for anionic species. While for B3LYP and BLYP, the trend was observed for cationic species from HfO⁺ to IrO⁺ and anionic species from HfO⁻ to AuO⁻ (except TaO⁻), but not for neutrals. Pure density function methods BLYP, BPW91, and BP86 give larger dissociation energy compared with hybrid density functional methods B3LYP, B3PW91, and B3P86. SVWN in most cases gives the smallest bond distance, while BLYP gives the largest value. MPW1PW91 and PBE1PBE show the same performance in predicting the spectroscopic constants. In addition, useful empirical criteria that one has obtained the ground states of a species and its ions are the spin multiplicities of a neutral and its single charged ions which differs by ±1.     

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: benchmarks approaching the complete basis set limit

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Moller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size.     

Theoretical calculation of bond dissociation energies and heats of formation for nitromethane and polynitromethanes with density functional theory

The C? NO₂ bond dissociation energies (BDEs) and the heats of formation (HOFs) of nitromethane and polynitromethanes (dinitromethane, trinitromethane, and tetranitromethane) system in gas phase at 298.15 K were calculated theoretically. Density functional theory (DFT) B3LYP, B3P86, B3PW91, and PBE0 methods in combination with different basis sets were employed. It was found that the C? NO₂ bond BDEs can be improved from B3LYP to B3PW91 to B3P86 or PBE0 functional. Levels of theory employing B3P86 and PBE0 functionals were found to be sufficiently reliable without the presence of diffusion functions. As the number of NO₂ groups on the same C atom increases, the PBE0 functional performs better than the B3P86 functional. Regarding the calculated HOFs, all four functionals can yield satisfactory results with deviations of <2 kcal mol^?1 from experimental ones for CH₂(NO₂)₂ and CH(NO₂)₃, when the diffusion functions are not augmented. For the C(NO₂)₄ molecule, the large basis sets augmented with polarization functions and diffusion functions are required to yield a good result. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density functional theory study of trans-dioxo complexes of iron, ruthenium, and osmium with saturated amine ligands, trans-[M(O)2(NH3)2(NMeH2)2]2+ (M=Fe, Ru, Os), and detection of [Fe(qpy)(O)2]n+ (n=1, 2) by high-resolution ESI mass spectrometry

Density functional theory (DFT) calculations on trans-dioxo metal complexes containing saturated amine ligands, trans-[M(O)2(NH3)2(NMeH2)2]2+ (M=Fe, Ru, Os), were performed with different types of density functionals (DFs): 1) pure generalized gradient approximations (pure GGAs): PW91, BP86, and OLYP; 2) meta-GGAs: VSXC and HCTH407; and 3) hybrid DFs: B3LYP and PBE1PBE. With pure GGAs and meta-GGAs, a singlet d2 ground state for trans-[Fe(O)2(NH3)2(NMeH2)2]2+ was obtained, but a quintet ground state was predicted by the hybrid DFs B3LYP and PBE1PBE. The lowest transition energies in water were calculated to be at lambda approximately 509 and 515 nm in the respective ground-state geometries from PW91 and B3LYP calculations. The nature of this transition is dependent on the DFs used: a ligand-to-metal charge-transfer (LMCT) transition with PW91, but a pi(Fe-O)-->pi*(Fe-O) transition with B3LYP, in which pi and pi* are the bonding and antibonding combinations between the dpi(Fe) and ppi(O(2-)) orbitals. The FeVI/V reduction potential of trans-[Fe(O)2(NH3)2NMeH2)2]2+ was estimated to be +1.30 V versus NHE based on PW91 results. The [Fe(qpy)(O)2](n+) (qpy=2,2':6',2':6',2':6',2'-quinquepyridine; n=1 and 2) ions, tentatively assigned to dioxo iron(V) and dioxo iron(VI), respectively, were detected in the gas phase by high-resolution ESI-MS spectroscopy.     

A DFT conformational analysis and VCD study on methyl tetrahydrofuran-2-carboxylate

DFT calculations were performed on (S)-methyl tetrahydrofuran-2-carboxylate to facilitate the interpretation of IR and VCD spectra. The potential energy surface could not be described unambiguously using the 6-31G* basis set in combination with different density functionals including B1LYP, B3LYP, B3P86, B3PW91, B98, BHandH, BHandHLYP, MPW1PW91 and PBE1PBE. In contrast, a uniform conformational picture could be found using the cc-pVTZ basis set. Using this large basis set and the collection of nine functionals from above, the dipole and rotational strengths were calculated, and compared to experimental values which were extracted from the experimental IR and VCD spectra for (+)-(S)-methyl tetrahydrofuran-2-carboxylate. A detailed analysis on the agreement between experiment and simulated spectra was performed by assigning the experimental bands based on the harmonic fundamentals obtained for all functionals except BHandH, which performs badly over the whole line. Assessing the dipole strengths, all tested functionals perform equally well. For the rotational strengths, differences can be observed: B3LYP, B1LYP and B98 give the highest correlation with experiment, while PBE1PBE gives the lowest correlation. Comparable conclusions are obtained using a neighborhood similarity measure.     

DFT study of the structural characteristics of the yttrium(3+) aqua ion

Density functional theory (DFT) using SVWN5, B3LYP, B3P86, O3LYP, B3PW91, B1LYP, B971, MPW1PW91, PBE1PBE, BHandH, and BHandHLYP density functionals was employed to study the structural characteristics of the Y(H₂O) ₈ ³⁺ yttrium aqua ion. The nonlocal hybrid GGA functionals show worse predictive ability in structural calculations of the Y(H₂O) ₈ ³⁺ aqua ion compared to the relatively simple combined functional BHandH and to the simplest SVWN5 functional in LSDA theory.     

铁磁性双核吲哚铜配合物的理论研究

采用密度泛函理论结合对称性破损态方法(DFT-BS),以双亚甲胺席夫碱为配体的双核吲哚铜配合物[Cu2L(AZ)(DMSO)](L=双甲亚胺三阴离子,由1-苯基-3-甲基-4-甲酰吡唑和1,3-二氨-2-丙醇衍生而成;AZ=7-氮杂吲哚阴离子)作为研究对象,通过与实验数据相比较,讨论了在不同密度泛函方法与基组对金属铜配物交换耦合常数的影响.结果表明,4种混合密度泛函方法(B3PW91,B3LYP,B3P86和PBE)及3种单一密度泛函(BPW91,BLYP和BP86)的计算结果都与实验值106cm-1有着相同的符号,但使用单一密度泛函所得到的结果与实验值比较接近.不论是混合密度泛函还是单一密度泛函计算所得到的耦合常数Jab对基组都有很大的依赖性,且BLYP/LANL2DZ水平下计算所得结果与实验数据吻合程度最好.研究表明,铁磁性配合物[Cu2L(AZ)(DMSO)]中存在着自旋离域和自旋极化效应.     

Benchmarking density functionals in conjunction with Grimme's dispersion correction for noble gas dimers (Ne2, Ar2, Kr2, Xe2, Rn2)

Eleven exchange‐correlational functionals of different types corrected for dispersion by Grimme's D3 correction in conjunction with the aug‐cc‐pVTZ basis set were tested on the following noble gas (Ng) dimers: Ne₂, Ar₂, Kr₂, Xe₂, and Rn₂. For comparison, the D2 and D3BJ corrections were probed with the B3LYP functional. From post‐HF wavefunction methods, CCSD(T) theory was also included. The investigated properties involved potential energy curves, equilibrium bond distances, and interaction energies. The B3LYP‐D3, B3LYP‐D3BJ, and PBE0‐D3 functionals performed overall best for bond distances, while B3LYP‐D3 and B97‐D3 performed best for interaction energies. The importance of fortunate error cancellations was seen in the often reduced agreement with reference data upon correction for BSSE. As several functionals performed well selectively for some noble gases (and poorly for others), we also analysed the performance on the Ng₂ dimers individually and recommended DFT‐D3 functionals for the calculation of large clusters of each Ng.     

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.     

Vibrational frequency scale factors for density functional theory and the polarization consistent basis sets

Calculated harmonic vibrational frequencies systematically deviate from experimental vibrational frequencies. The observed deviation can be corrected by applying a scale factor. Scale factors for: (i) harmonic vibrational frequencies [categorized into low (<1000 cm^?1) and high (>1000 cm^?1)], (ii) vibrational contributions to enthalpy and entropy, and (iii) zero‐point vibrational energies (ZPVEs) have been determined for widely used density functionals in combination with polarization consistent basis sets (pc‐n, n = 0,1,2,3,4). The density functionals include pure functionals (BP86, BPW91, BLYP, HCTH93, PBEPBE), hybrid functionals with Hartree‐Fock exchange (B3LYP, B3P86, B3PW91, PBE1PBE, mPW1K, BH&HLYP), hybrid meta functionals with the kinetic energy density gradient (M05, M06, M05‐2X, M06‐2X), a double hybrid functional with Møller‐Plesset correlation (B2GP‐PLYP), and a dispersion corrected functional (B97‐D). The experimental frequencies for calibration were from 41 organic molecules and the ZPVEs for comparison were from 24 small molecules (diatomics, triatomics). For this family of basis sets, the scale factors for each property are more dependent on the functional selection than on basis set level, and thus allow for a suggested scale factor for each density functional when employing polarization consistent basis sets (pc‐n, n = 1,2,3,4). A separate scale factor is recommended when the un‐polarized basis set, pc‐0, is used in combination with the density functionals. © 2012 Wiley Periodicals, Inc.     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

Performance of density functionals for calculation of reductive ring-opening reaction energies of Li<Superscript>+</Superscript>-EC and Li<Superscript>+</Superscript>-VC

Reaction energies were determined for reductive ring-opening reactions of Li⁺-coordinated ethylene carbonate (EC) and vinylene carbonate (VC) by using various density functional theory (DFT) and ab-initio methods applying the basis sets up to Dunnings aug-cc-pVQZ. The methods examined include the local density functional (SVWN), the pure gradient-corrected density functionals (BLYP and BPW91), and the hybrid density functionals (B3LYP, B1LYP, B3PW91, and mPW1PW91). Comparison of the DFT results with ab-initio results indicates that the mPW1PW91 approach introduced by Adamo and Barone, is superior to all the other DFT methods (including B3LYP). The performance of more cost-effective Pople-type basis sets ranging from 6-31G(d,p) to 6-311++G(3df,3pd) was assessed at DFT levels of theory by calibrating them with the aug-cc-pVQZ results