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.     

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.     

A new hybrid DFT approach to electronic excitation and first hyperpolarizabilities of transition metal complexes

The electronic excitations and the static first hyperpolarizability of three typical transition metal (TM) aromatic carbonyl complexes, two tungsten pentacarbonyl derivatives (W(CO)₅L, L = Py and PyCHO) and one chromium tricarbonyl arene derivative (Cr(CO)₃Bz, Bz = benzene), have been theoretically studied by a variety of density functional methods. The assessments reveal that most of the conventional DFT methods including local density approximation, generalized gradient approximation (GGA), and the various kinds of hybrid exchange‐correlation (xc) methods present the first hyperpolarizabilities of these TM‐containing molecules with large deviations from the experiments. A one‐parameter hybrid xc functional is introduced by using the Perdew‐Wang 1991 gradient‐corrected correlation functional (PW91) and the Barone's‐modified PW91 exchange functional (mPW). The ratio between the exact and the density functional exchange is determined to be 0.40 by the adiabatic connection method. The application of the new hybrid functional to the three organometallic carbonyl molecules results in the satisfactory agreement between the calculated first hyperpolarizabilities and the experimental ones. The second‐order nonlinear optical properties of the three organometallic complexes are addressed to the metal‐to‐ligand charge transfers, and the extended π‐delocalization ligands benefit the enhancement of the first hyperpolarizability. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Where is the spin? Understanding electronic structure and g-tensors for ruthenium complexes with redox-active quinonoid ligands

Understanding the bonding in transition metal complexes with redox-active ligands is a major challenge, for example in redox catalysis or in bioinorganic chemistry. In this work, electronic g-tensors, spin-density distributions, and electronic structure have been studied by different density functional methods for an extended series of complexes [Ru(acac)2(L)]n (n = -1, 0, +1; L = redox-active o-quinonoid ligand). Comparison is made with experimental g-tensors and g-tensor-based oxidation-state assignments for a number of experimentally studied examples, using both gradient-corrected (BP86) and hybrid functionals (B3LYP, BHLYP) representing a range of exact-exchange admixtures. Reasonable, albeit not perfect, agreement with experimental g-tensors is obtained in one-component DFT calculations with hybrid functionals. Analyses of spin densities confirm the assignment of the cationic complexes as predominantly d5-Ru(III) with a neutral quinonoid ligand. However, this conclusion is obtained only after inclusion of the appreciable spin polarization of the unrestricted determinant, while the singly occupied molecular orbital (SOMO) is localized more on the acac ligands. The anionic complexes turn out to be approximately halfway between a d6-Ru(II)/semiquinone and a d5-Ru(III)/catecholate formulation, but again only after taking into account the extensive spin polarization. Even the previous assignment of the neutral parent systems as d5-Ru(III)/semiquinone is not accurate, as a d6-Ru(II)/quinone resonance structure contributes to some extent. Very unusual trends in the spin contamination of the Kohn-Sham determinant with increasing exact-exchange admixture in some of the cationic complexes have been traced to an interplay between spin delocalization and spin polarization.     

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.     

Benchmarking approximate density functional theory for s/d excitation energies in 3d transition metal cations

Holthausen has recently provided a comprehensive study of density functional theory for calculating the s/d excitation energies of the 3d transition metal cations. This study did not include the effects of scalar relativistic effects, and we show here that the inclusion of scalar relativistic effects significantly alters the conclusions of the study. We find, contrary to the previous study, that local functionals are more accurate for the excitation energies of 3d transition method cations than hybrid functionals. The most accurate functionals, of the 38 tested, are SLYP, PBE, BP86, PBELYP, and PW91.     

Computational Study of Cage Like (ZnO)12 Cluster Using Hybrid and Hybrid Meta Functionals

Density Functional Theory employing hybrid and M06 functionals in combination with three different basis sets is used to calculate the ground state of a cage like (ZnO)₁₂ nanocluster which has been consistently reported as the more stable cluster for its particular size. B3LYP and B3PW91 hybrid functionals combined with 6‐31+G*, Lanl2dz and SDD basis sets are employed to treat the ZnO molecular system. Alternatively, three M06 functionals in combination with three basis sets are employed in the nanostructure calculations. Results obtained by treating ZnO sodalite cage nanocluster with M06 functionals demonstrated comparable quality to results obtained with hybrid functionals. Within this study, efficient theoretical DFT methods with the widely known hybrid and the recently created M06 meta‐hybrid functionals are employed to study nanostructured ZnO. Our resulting parameters provide a fresh approach performance wise on the different theoretical methods to treat transition metal nanostructures, particularly, ZnO nanoclusters geometry and electronic structure.     

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     

Performance assessment of density-functional methods for study of charge-transfer complexes

Various density functionals are applied to a number of weakly bound intermolecular pi-pi charge-transfer (CT) complexes. Most functionals, including the recently developed mPWPW91 and mPW1PW91, grossly underestimate experimental excitation energies; good agreement is obtained only with the half-and-half hybrid BH&HLYP functional. PW91PW91 provides the best agreement with intermolecular distances measured in crystal, while the BH&HLYP values are about 0.1 A too long. Various hybrid functionals with nonlocal exchange correction provide binding energies that compare favorably with the experimental heats of formation measured in solution.     

Computational studies of electron paramagnetic resonance parameters for paramagnetic molybdenum complexes. 1. Method validation on small and medium-sized systems

A variety of density functional methods have been evaluated in the computation of electronic g-tensors and molybdenum hyperfine couplings for systems ranging from the Mo atom through MoIIIN, [MoVOCl4]-, and [MoVOF5]2- to two larger MoV complexes MoXLCl2 (X=O, S; L=tris(3,5-dimethylpyrazolyl)hydroborate anion). In particular, the influence of the molybdenum basis set and of various exchange-correlation functionals with variable admixtures of Hartree-Fock exchange on the computed EPR parameters have been evaluated in detail. Careful basis-set studies have provided a moderate-sized 12s6p5d all-electron basis on molybdenum that gives hyperfine tensors in excellent agreement with much larger basis sets and that will be useful for calculations on larger systems. The best agreement with experimental data for both hyperfine and g-tensors is obtained with hybrid functionals containing approximately 30-40% Hartree-Fock exchange. Only for MoSLCl2 does increasing spin contamination with increasing exact-exchange admixture restrict the achievable computational accuracy. In all cases, spin-orbit corrections to the hyperfine tensors are sizable and have to be included in accurate calculations. Scalar relativistic effects enhance the isotropic Mo hyperfine coupling by approximately 15-20%. Two-component g-tensor calculations with variational inclusion of spin-orbit coupling show that the Deltag parallel components in [MoVOCl4]- and [MoVOF5]2- depend on higher-order spin-orbit contributions and are thus described insufficiently by the usual second-order perturbation approaches. Computed orientations of g- and hyperfine tensors relative to each other and to the molecular framework for the MoXLCl2 complexes provide good agreement between theory and single-crystal electron paramagnetic resonance experiments. In these cases, the hyperfine tensor orientations are influenced only slightly by spin-orbit effects.     

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.     

De novo prediction of the ground state structure of transition metal complexes using semiempirical and ab initio quantum mechanics. Coordination isomerism

The ground state coordination isomers for 30 different trigonal bipyramidal transition metal complexes have been predicted using different levels of quantum mechanics: semiempirical (PM3(tm)), ab initio (MP2//HF), pure (BPW91) and hybrid (B3PW91) density functional theory (DFT) methods. For species where these methods failed to reproduce crystallographic data, hybrid quantum mechanics/molecular mechanics (QM/MM) methods were used to study more exact experimental models. Literature deficiencies regarding ground state multiplicity of these species were supplemented by spin predictions using previously tested PM3(tm) methods. Geometry optimization calculations were performed for each possible coordination isomer. The predicted ground state minima provided by the different methods are compared to each other and with crystallographic data. Pure DFT functionals outperformed hybrid functionals and MP2//HF. The very rapid PM3(tm) parameterization method provided accurate predictions in comparison to other levels of theory. An integrated MM/PM3(tm)/DFT de novo scheme accurately reproduced crystallographic data for species where the individual methods failed.     

Understanding the EPR parameters of glycine-derived radicals: the case of N-acetylglycyl in the N-acetylglycine single-crystal environment

As a step toward an in-depth understanding of the electron paramagnetic resonance parameters of glycyl radicals in proteins, the hyperfine tensors and, particularly, the g-tensor of N-acetylglcyl in the environment of a single crystal of N-acetylglycine have been studied by systematic state-of-the-art quantum chemical calculations on various suitable model systems. The quantitative computation of the g-tensors for such glycyl-derived radicals is a veritable challenge, mainly because of the very small g-anisotropy combined with a nonsymmetrical, delocalized spin-density distribution and several atoms with comparable spin-orbit contributions to the g-tensors. The choice of gauge origin of the magnetic vector potential, and of approximate spin-orbit operators, both turn out to be more critical than found in previous studies of g-tensors for organic radicals. Environmental effects, included by supermolecular hydrogen-bonded models, were found to be moderate, because of a partial compensation between the influences from intramolecular and intermolecular hydrogen bonds. The largest effects on the g-tensor are caused by the conformation of the radical. The density functional theory methods employed systematically overestimate both the Delta gx and Delta gy components of the g-tensor. This is important for parallel investigations on the protein-glycyl radicals. The 1H alpha and 13C alpha hyperfine couplings depend only slightly on the supermolecular model chosen and appear less sensitive probes of detailed structure and environment.     

Relativistic two-component calculations of electronic g-tensors that include spin polarization

The first two-component relativistic density-functional approach for the calculation of electronic g-tensors is reported that includes spin polarization using noncollinear spin-density functionals. The method is based on the relativistic Douglas-Kroll-Hess Hamiltonian and has been implemented into the ReSpect program package. Using three self-consistent-field calculations with orthogonal orientations of total magnetization J, the full g-matrix may be obtained. In contrast to previous spin-restricted two-component treatments, results with the new approach agree excellently with spin-polarized one-component calculations for light-atom radicals. Additionally, unlike one-component approaches, the method also reproduces successfully the negative deltag(parallel)-values of heavy-atom 2sigma radicals and the negative deltag(perpendicular) components in cysteinyl. The new method removes effectively the dilemma existing up to now regarding the simultaneous inclusion of spin polarization and higher-order spin-orbit effects in g-tensor calculations. It is straightforwardly applicable to higher than doublet spin multiplicities and has been implemented with hybrid functionals.     

Comparison of DFT methods for molecular orbital eigenvalue calculations

We report how closely the Kohn-Sham highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenvalues of 11 density functional theory (DFT) functionals, respectively, correspond to the negative ionization potentials (-IPs) and electron affinities (EAs) of a test set of molecules. We also report how accurately the HOMO-LUMO gaps of these methods predict the lowest excitation energies using both time-independent and time-dependent DFT (TD-DFT). The 11 DFT functionals include the local spin density approximation (LSDA), five generalized gradient approximation (GGA) functionals, three hybrid GGA functionals, one hybrid functional, and one hybrid meta GGA functional. We find that the HOMO eigenvalues predicted by KMLYP, BH&HLYP, B3LYP, PW91, PBE, and BLYP predict the -IPs with average absolute errors of 0.73, 1.48, 3.10, 4.27, 4.33, and 4.41 eV, respectively. The LUMOs of all functionals fail to accurately predict the EAs. Although the GGA functionals inaccurately predict both the HOMO and LUMO eigenvalues, they predict the HOMO-LUMO gap relatively accurately (approximately 0.73 eV). On the other hand, the LUMO eigenvalues of the hybrid functionals fail to predict the EA to the extent that they include HF exchange, although increasing HF exchange improves the correspondence between the HOMO eigenvalue and -IP so that the HOMO-LUMO gaps are inaccurately predicted by hybrid DFT functionals. We find that TD-DFT with all functionals accurately predicts the HOMO-LUMO gaps. A linear correlation between the calculated HOMO eigenvalue and the experimental -IP and calculated HOMO-LUMO gap and experimental lowest excitation energy enables us to derive a simple correction formula.     

The confirmation of accurate combination of functional and basis set for transition‐metal dimers: Fe2, Co2, Ni2, Ru2, Rh2, Pd2, Os2, Ir2, and Pt2

Eight kinds of density functionals named B3LYP, PBE1PBE, B1B95, BLYP, BP86, G96PW91, mPWPW91, and SVWN along with two different valence basis sets (LANL2DZ and CEP‐121g) are employed to study the transition‐metal dimers for the elements of group VIII. By comparing the equilibrium bond distances, vibrational frequencies, and dissociation energies of the ground state of these dimers with the available experimental values and theoretical data, we show that the “pure” DFT methods (G96PW91, BLYP, and BP86) with great‐gradient approximation always give better results relative to the hybrid HF/DFT schemes (B3LYP, PBE1PBE, and B1B95). The striking case found by us is that the G96PW91 functional, which is not tested in previous systemic studies, always predicts the dissociation energy to be well. The Ru₂ and Os₂ dimers are sensitive to not only the functionals employed but also the valence basis sets adopted. The natural bond orbital population is analyzed, and the molecular orbitals of the unpaired electrons are determined. Furthermore, our results indicate that the s and d orbitals of these dimers always hybridize with each other except for Rh₂ and Pt₂ molecules. And by analyzing the electron configuration of the bonding atom, the dissociation limit of the ground state is obtained. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Does the gradient-regulated connection improve the description of correlated metal bond properties?

Gradient-regulated connection (GRAC) is a generalized gradient approximation exchange density functional designed by combining the revPBE and PW91 exchange functionals to impose their behaviors in the slowly- and fast-varying density regions, respectively. Such a construction allows one single density functional to accurately estimate both covalent and weak interactions occurring in main-group-based molecular systems. For the first time, the assessment of the performance of the GRAC exchange functional is extended to the modeling of various metal bond energy and structure properties. This assessment shows that when GRAC is coupled with the Perdew, Burke, Ernzerhof (PBE) correlation, the resulting exchange-correlation density functional is an excellent alternative to global hybrids to model bond dissociation energy, atomic electronic excitation energy, and bond length structure properties of single-reference metal bonds. It also shows that coupling with the Tognetti, Cortona, Adamo (TCA) correlation constitutes a robust approach to tackle energy bond properties of organometallic complexes with multi-reference character.     

B2-PPW91: a promising double-hybrid density functional for the electric response properties

A new double-hybrid density functional, termed B2-PPW91, is presented which includes the Becke88 (B88) exchange in conjunction with Perdew-Wang91 (PW91) gradient-corrected correlation functional. The fitting parameters are obtained by minimization of mean absolute error of the static dipole polarizability of 4d transition metal monohalides against the CCSD(T)∕aug-cc-pVTZ∕SDD results. The performance of proposed functional has been assessed for estimation of other response properties, such as dipole moment and excitation energy, for the same species. We then proceed to explore the validity of B2-PPW91 method for calculation of the dipole polarizability of some 5d transition metal monofluorides. In all cases, the improvement compared to common density functional methods and even previously reported double-hybrid functionals such as B2-PLYP and mPW2-PLYP has been observed. This indicates that the utility of double-hybrid density functional methods can be further extended to study linear and non-linear optical properties of transition metal containing molecules.