Vibrational corrections to geometries of transition metal complexes from density functional theory

Zero-point vibrational corrections are computed at the BP86/AE1 level for the set of 50 transition-metal/ligand bonds that have recently been proposed as testing ground for DFT methods, because of the availability of precise experimental gas-phase geometries (Bühl and Kabrede, J Chem Theory Comput 2006, 2, 1282). These corrections are indicated to be transferable to a large extent between various density-functional/basis-set combinations, so that they can be used to estimate zero-point averaged r0g distances from re values optimized at other theoretical levels. Applying this approach to a number of popular DFT levels does not, in general, improve their overall accuracy in terms of mean and standard deviations from experiment. The hybrid variant of the meta-functional TPSS is confirmed as promising choice for computing structures of transition-metal complexes.     

Spin quenching of transition metals deposited on MgO insulator and CdO semiconductor density functional calculations

We have analyzed spin quenching of first row transition metals deposited on (001) defect‐free and defect‐containing surfaces of MgO insulator and CdO semiconductor by means of density functional calculations and embedded cluster model. Clusters of moderate sizes were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces, and relaxation of ions that surround the defect sites was taken into account. Spin states of metals deposited on the defect free surfaces were maintained as in the isolated metals except for Ti, V, and Co on MgO, and Ti, V, and Cr on CdO. On the defect containing surfaces, spin states were maintained too except for Fe on MgO, and V and Cr on CdO. The metal‐support interactions stabilize the low spin state of the adsorbed metal with respect to the isolated metal, but the effect was not in general enough to quench the spin. Spin polarization effects tend to preserve the spin states of the adsorbed metals relative to those of the isolated metals. Although charge transfer took place from the adsorbed metal to the insulator surface, it took place the other way round from the semiconductor surface to the adsorbed metal. The encountered variations in magnetic properties were attributed to the smaller band gap of the semiconductor, and the behavior of a single metal atom adsorbed on a particular surface was a result of a competition between Hund's rule for the adsorbed metal and the formation of a chemical bond at the interface. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Modeling a halogen dance reaction mechanism: A density functional theory study

Since the discovery of the halogen dance (HD) reaction more than 60 years ago, numerous insights into the mechanism have been unveiled. To date however, the reaction has not been investigated from a theoretical perspective. Density functional theory (DFT) was used to model the potential energy surface linking the starting reagents to the lithiated products for each step in the mechanism using a thiophene substrate. It was found that the lithium‐halogen exchange mechanism is critical to understand the HD mechanism in detail and yielded the knowledge that S_N2 transition states (TS) are favored over the four‐center type for the lithium‐bromine exchange steps. The overall driving force for the HD is thermodynamics, while the kinetic factors tightly control the reaction path through temperature. The S_N2 lithium‐bromide TS are barrierless, except the second, which is the limiting step. Finally, the model for the HD is discovered to be a pseudo‐clock type, due to a highly favorable bromide catalysis step and the reformation of 2‐bromothiophene. © 2016 Wiley Periodicals, Inc.     

Synthesis and characterization of new five-coordinate platinum nitrosyl derivatives: density functional theory study of their electronic structure

The neutral, five-coordinate platinum nitrosyl compounds [Pt(C(6)F(5))(3)(L)(NO)] (2) [L=CNtBu (2 a), NC(5)H(4)Me-4 (2 b), PPhMe(2) (2 c), PPh(3) (2 d) and tht (2 e)] have been prepared by the reaction of [NBu(4)][Pt(C(6)F(5))(3)(L)] (1) with NOClO(4) in CH(2)Cl(2). The ionic compound [N(PPh(3))(2)][Pt(C(6)F(5))(4)(NO)] (4) has been prepared in a similar way starting from the homoleptic species [N(PPh(3))(2)](2)[Pt(C(6)F(5))(4)] (3). Compounds 2 and 4 are all diamagnetic with [PtNO](8) electronic configuration and show nu(NO) stretching frequencies at around 1800 cm(-1). The crystal and molecular structures of 2 c and 4 have been established by X-ray diffraction methods. The coordination environment for the Pt center in both compounds can be described as square pyramidal (SPY-5). Bent nitrosyl coordination is observed in both cases with Pt-N-O angles of 120.1(6) and 130.2(7) degrees for 2 c and 4, respectively. The bonding mechanism of the nitrosyl ligand coordinated to various model [Pt(II)R(4)](2-) (R=H, Me, Cl, CN, C(6)F(5) or C(6)Cl(5)) and [Pt(C(6)F(5))(3)(L)](-) (L=CNMe, PH(3)) systems has been studied by density functional calculations at the B3LYP level of theory, using the SDD basis set. The R(4)Pt-NO and (C(6)F(5))(3)(L)Pt-NO interactions generally involve two components: i) a direct Pt-NO bonding interaction and ii) multicenter-bonding interactions between the N atom of the NO ligand and the donor atoms of the R and L ligands. Moreover, with the more complex R groups, C(6)F(5) or C(6)Cl(5), a third component has been found to arise, which involves multicenter electrostatic interactions between the positively charged NO ligand and the negatively charged halo-substituents in the ortho-position of the C(6)X(5) groups (X=F, Cl). The contribution of each component to the Pt-NO bonding in R(4)Pt-NO and (C(6)F(5))(3)(L)Pt-NO compounds seems to be modulated by the electronic and steric effects of the R and L ligands.     

A density functional theory study on the electrocyclization of 1,2,4,6-heptatetraene analogues: converting a pericyclic to a pseudopericyclic reaction

A comprehensive B3LYP/6-31+G* study on the electrocyclization of 1,2,4,6-heptatetraene analogues was conducted. Starting from the cyclization of (2Z)-2,4,5-hexatrienal, a pericyclic disrotatory process favored by the assistance of a electron lone pair, we incorporated small modifications in its molecular structure to obtain a truly pseudopericyclic process. To this purpose electronegative atoms (fluorine and nitrogen) were added to give a more electrophilic character on the carbon atom which is attacked by the electron lone pair of the oxygen atom. The complete pathway for each reaction was determined, and changes in magnetic properties were monitored with a view to estimating the aromatization associated with each process. This information, together with the energetic and structural results, allowed us to classify the reactions as pseudopericyclic or pericyclic. Among all studied reactions only one was a truly pseudopericyclic process and another was a borderline case. The features of this unequivocally pseudopericyclic case were analyzed in depth.     

Transition metal atoms on oxide supports density functional calculations

The adsorption properties of Cu, Ag, Ni, and Pd atoms on O^2?, F, and F⁺ sites of MgO, CaO, SrO, and BaO (001) surfaces have been studied by means of density functional calculations. The examined clusters were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces. The adsorption properties have been analyzed with reference to the basicity and energy gap of the oxide support in addition to orbital interactions. While the free Ni d⁹s¹ triplet ground state is preserved on adsorption on the O^2? sites of MgO, CaO, and SrO surfaces, it is no longer preserved on the O^2? site of BaO. For all adsorbates considered, adsorption is found to be stronger on F⁺ sites compared with regular O^2? sites. While on the O^2? site, Pd and Ni form the most stable complexes, on the F site, Pd and Cu form the most stable counterparts. On the F⁺ site, the single valence electron of Cu and Ag atoms couples with the unpaired electron of the vacancy forming a covalent bond. As a result, the adsorption energies of these atoms on the F⁺ site are stronger than those on the F and O^2? sites. The adsorption properties correlate linearly with the basicity and energy gap of the oxide support in addition to orbital interactions. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Actinium coordination chemistry: A density functional theory study with monodentate and bidentate ligands

In the current study, the coordination chemistry of nine-coordinate Ac(III) complexes with 35 monodentate and bidentate ligands was investigated using density functional theory (DFT) in terms of their geometries, charges, reaction energies, and bonding interactions. The energy decomposition analysis with naturals orbitals for chemical valence (EDA-NOCV) and the quantum theory of atoms in molecules (QTAIM) were employed as analysis methods. Trivalent Ac exhibits the highest affinities toward hard acids (such as charged oxophilic donors, fluoride), so its classification as a hard acid is justified. Natural population analysis quantified the involvement of 5f orbitals on Ac to be about 30% of total valence electron natural configuration indicating that Ac is a member of the actinide series. Pearson correlation coefficients were used to study the pairwise correlations among the bond lengths, ΔG reaction energies, charges on Ac and donor atoms, and data from EDA-NOCV and QTAIM. Strong correlations and anticorrelations were found between Voronoi charges on donor atoms with ΔG, EDA-NOCV interaction energies and QTAIM bond critical point densities.     

The challenge of the so-called electron configurations of the transition metals

Quite different meanings are attached by chemists to the words element, atom, orbital, order of orbitals or configurations. This causes conceptual inconsistencies, in particular with respect to the transition-metal elements and their atoms or ions. The different meanings will here be distinguished carefully. They are analyzed on the basis of empirical atomic spectral data and quasi-relativistic density functional calculations. The latter are quite reliable for different average configuration energies of transition-metal atoms. The so-called "configurations of the chemical elements", traditionally displayed in periodic tables, are the dominant configurations of the lowest spin-orbit levels of the free atoms. They are chemically rather irrelevant. In many-electron systems the ns and np AOs are significantly below the more hydrogen-like nd ones. Even (n+1)s is below nd for all light neutral atoms from C onwards, but only up to the first elements of the respective long rows! The most common orbital order in transition-metal atoms is 3p < 3d < 4s etc. The chemically relevant configuration in group g is always d(g) instead of d(g-2) s(2). Conceptually clear reasoning eliminates apparent textbook inconsistencies between simple quantum-chemical models and the empirical facts. The empirically and theoretically well-founded Rydberg (n-deltal) rule is to be preferred instead of the historical Madelung (n+l) rule with its large number of exceptions.     

Theoretical study of rhodium(I) carbene complexes: the structural versatility of phosphino- compared with aminocarbenes

Density functional calculations are reported for complexes of general formula [(carbene)RhClL(2)] featuring model phosphino- and aminocarbenes. Both the cis and trans isomers of the rhodium(I) eta(1)-complexes (1-9) were investigated, and the influence of the rhodium co-ligands (L=ethylene, phosphine, or carbon monoxide) was evaluated. In the case of phosphinocarbenes and carbon monoxide as a ligand, a somewhat unusual coordination mode was observed, in which a significant intramolecular Cl-->C(carbene) interaction is present. The propensity of phosphino- and aminocarbenes to behave as four electron donors was also investigated both structurally and energetically on the related eta(2)-complexes 10-18. These results as a whole emphasize the structural versatility of phosphino- compared with aminocarbene complexes.     

Hydroformylation and isomerization of allene and propyne: a density functional theory study

The [HCo(CO)3]-catalyzed hydroformylation of allene and propyne has been investigated at the B3LYP level of density functional theory. It is found that hydroformylation of allene favors the linear anti-Markovnikov product in high regioselectivity both kinetically and thermodynamically. The origin of this regioselectivity comes from the enhanced stability of the eta3-allylic intermediate [(eta3-CH2CHCH2)Co(CO)3]. By contrast, propyne does not show any regioselectivity. The possible interconversion between allene and propyne mediated by [HCo(CO)3] has been explored.     

Thermal effects and vibrational corrections to transition metal NMR chemical shifts

Both zero-point and classical thermal effects on the chemical shift of transition metals have been calculated at appropriate levels of density functional theory for a number of complexes of titanium, vanadium, manganese and iron. The zero-point effects were computed by applying a perturbational approach, whereas classical thermal effects were probed by Car-Parrinello molecular dynamics simulations. The systematic investigation shows that both procedures lead to a deshielding of the magnetic shielding constants evaluated at the GIAO-B3 LYP level, which in general also leads to a downfield shift in the relative chemical shifts, delta. The effect is small for the titanium and vanadium complexes, where it is typically on the order of a few dozen ppm, and is larger for the manganese and iron complexes, where it can amount to several hundred ppm. Zero-point corrections are usually smaller than the classical thermal effect. The pronounced downfield shift is due to the sensitivity of the shielding of the metal centre with regard to the metal-ligand bond length, which increase upon vibrational averaging. Both applied methods improve the accuracy of the chemical shifts in some cases, but not in general.     

Aggregation behavior and non-covalent functionalization of borofullerenes B28, B38, and B40: A density functional theory investigation

Fullerene-like boron clusters (borofullerenes) are rising stars in the field of cluster chemistry. In this work, density functional theory calculations revealed that the recently reported small borofullerenes B_n (n = 28, 38, and 40) are all highly reactive and tend to form dimers and even trimers spontaneously. In addition, the non-covalent modification of these borofullerenes by various cycloparaphenylene nanorings can form stable host-guest systems with substantial intermolecular charge transfer at both ground and excited states. Our results demonstrate that the borofullerenes are versatile platform for exohedral functionalization, and are very promising candidates for the design of novel nanomaterials with desirable properties.     

Formaldehyde encapsulated in lithium-decorated metal-organic frameworks: a density functional theory study

The stability of monomeric formaldehyde encapsulated in the lithium-decorated metal-organic framework Li-MOF-5 was investigated by means of density functional calculations with the M06-L functional and the 6-31G(d,p) basis set. To assess the efficiency of Li-MOF-5 for formaldehyde preservation, we consider the reaction kinetics and the thermodynamic equilibrium between formaldehyde and its trimerized product, 1,3,5-trioxane. We propose that trimerization of encapsulated formaldehyde takes place in a single reaction step with an activation energy of 34.5 kcal mol(-1). This is 17.2 kcal mol(-1) higher than the corresponding activation energy in the bare system. In addition, the reaction energy of the system studied herein is endothermic by 6.1 kcal mol(-1) and the Gibbs free energy (ΔG) of the reaction becomes positive (11.0 kcal mol(-1)). Consequently, the predicted reverse rate for the trimerization reaction in the Li-MOF-5 is significantly faster than the forward rate. The calculations show that the oligomerization of formaldehyde in Li-MOF-5 is a reversible reaction, suggesting that such a zeolite might be a good candidate material for preserving formaldehyde in its monomeric form.     

The 1deltag dioxygen ene reaction with propene: a density functional and multireference perturbation theory mechanistic study

This study aims to determine whether a balance between concerted and non-concerted pathways exists, and in particular to ascertain the possible role of diradical/zwitterion or peroxirane intermediates. Three non-concerted pathways, via 1) diradical or 2) peroxirane intermediates, and 3) by means of hydrogen-abstraction/radical recoupling, plus one concerted pathway (4), are explored. The intermediates and transition structures (TS) are optimized at the DFT(MPW1K), DFT(B3LYP) and CASSCF levels of theory. The latter optimizations are followed by multireference perturbative CASPT2 energy calculations. (1) The polar diradical forms from the separate reactants by surmounting a barrier (deltaE(++)(MPW1K)=12, deltaE++(B3LYP)=14, and deltaE(++)(CASPT2)=16 kcal mol(-1) and can back-dissociate through the same TS, with barriers of 11 (MPW1K) and 8 kcal mol(-1) (B3LYP and CASPT2). The diradical to hydroperoxide transformation is easy at all levels (deltaE(++)(MPW1K)<4, deltaE(++)(B3LYP)=1 and deltaE(++)(CASPT2)=1 kcal mol(-1)). (2) Peroxirane is attainable only by passing through the diradical intermediate, and not directly, due to the nature of the critical points involved. It is located higher in energy than the diradical by 12 kcal mol(-1), at all theory levels. The energy barrier for the diradical to cis-peroxirane transformation (deltaE(++)=14-16 kcal mol(-1)) is much higher than that for the diradical transformation to the hydroperoxide. In addition, peroxirane can very easily back-transform to the diradical (deltaE(++)<3 kcal mol(-1)). Not only the energetics, but also the qualitative features of the energy hypersurface, prevent a pathway connecting the peroxirane to the hydroperoxide at all levels of theory. (3) The last two-step pathway (hydrogen-abstraction by (1)O(2), followed by HOO-allyl radical coupling) is not competitive with the diradical mechanism. (4) A concerted pathway is carefully investigated, and deemed an artifact of restricted DFT calculations. Finally, the possible ene/[pi2+pi2] competition is discussed.     

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

The performance of a number of different implementations of density functional theory (DFT) for predicting the s/d interconfigurational energies of the 3d transition metal cations is investigated. Systematic comparisons of computed results with experimental data indicate that gradient corrected correlation functionals, like the LYP GGA, efficiently correct the flaws of the LDA, but reveal shortcomings in the treatment of exchange by currently available GGAs. The admixture of exact exchange in hybrid functionals eventually leads to largely reduced errors. Several basis sets available for the 3d elements are tested in combination with the B3LYP functional. Finally, the influence of variations of the admixture of exact exchange is systematically tested. The results reveal that computed s/d excitation energies obtained for the individual ions depend in markedly different ways on the amount of exact exchange admixture and that there is no single optimal and transferable exchange parameter to create a hybrid functional that yields improved results for all ions alike. Several of the recently devised functionals perform as good as or slightly better than the B3LYP functional in the present study. But given the fact that the B3LYP functional has been identified as the most successful DFT method in an overwhelming number of systematic investigations in very many areas of chemical research, there is no persuasive motivation to recommend its replacement by one of the other functionals, as much less is known about their robustness.