2,2'-Biphosphinines and 2,2'-bipyridines in homoleptic dianionic Group 4 complexes and neutral 2,2'-biphosphinine Group 6 d(6) metal complexes: octahedral versus trigonal-prismatic geometries

The geometric and electronic structure of formally d(6) tris-biphosphinine [M(bp)(3)](q) and tris-bipyridine [M(bpy)(3)](q) complexes were studied by means of DFT calculations with the B3LYP functional. In agreement with the available experimental data, Group 4 dianionic [M(bp)(3)](2-) complexes (1P-3P for M=Ti, Zr, and Hf, respectively) adopt a trigonal-prismatic (TP) structure, whereas the geometry of their nitrogen analogues [M(bpy)(3)](2-) (1N-3N) is nearly octahedral (OC), although a secondary minimum was found for the TP structures (1N'-3N'). The electronic factors at work in these systems are discussed by means of an MO analysis of the minima, MO correlation diagrams, and thermodynamic cycles connecting the octahedral and trigonal-prismatic limits. In all these complexes, pronounced electron transfer from the metal center to the lowest lying pi* ligand orbitals makes the d(6) electron count purely formal. However, it is shown that the bp and bpy ligands accommodate the release of electron density from the metal in different ways because of a change in the localization of the HOMO, which is a mainly metal-centered orbital in bp complexes and a pure pi* ligand orbital in bpy complexes. The energetic evolution of the HOMO allows a simple rationalization of the progressive change from the TP to the OC structure on successive oxidation of the [Zr(bp)(3)](2-) complex, a trend in agreement with the experimental structure of the monoanionic complex. The geometry of Group 6 neutral complexes [M(bp)(3)] (4P and 5P for M=Mo and W, respectively) is found to be intermediate between the TP and OC limits, as previously shown experimentally for the tungsten complex. The electron transfer from the metal center to the lowest lying pi* ligand orbitals is found to be significantly smaller than for the Group 4 dianionic analogues. The geometrical change between [Zr(bp)(3)](2-) and [W(bp)(3)] is analyzed by means of a thermodynamic cycle and it is shown that a larger ligand-ligand repulsion plays an important role in favoring the distortion of the tungsten complex away from the TP structure.     

Atom-molecule interactions on transition metal surfaces: a DFT Study of CO and several atoms on Rh(100), Pd(100) and Ir(100).

Density functional theory (DFT) calculations have been performed to determine the interaction energy between a CO probe molecule and all atoms from the first three rows of the periodic table coadsorbed on Rh(100), Pd(100) and Ir(100) metal surfaces. Varying the coverage of CO or the coadsorbed atom proved to have a profound effect on the strength of the interaction energy. The general trend, however, is the same in all cases: the interaction energy becomes more repulsive when moving towards the right along a row of elements, and reaches a maximum somewhere in the middle of a row of elements. The absolute value of the interaction energy between an atom-CO pair ranges from about -0.40 eV (39 kJ mol(-1)) attraction to +0.70 eV (68 kJ mol(-1)) repulsion, depending on the coadsorbate, the metal and the coverage. The general trend in interaction energies seems to be a common characteristic for several transition metals.     

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.     

Ambiphilic diphosphine-borane ligands: metal-->borane interactions within isoelectronic complexes of rhodium, platinum and palladium

Coordination of an ambiphilic diphosphine-borane (DPB) ligand to the RhCl(CO) fragment affords two isomeric complexes. According to X-ray diffraction analysis, each complex adopts a square-pyramidal geometry with trans coordination of the two phosphine buttresses and axial RhB contacts, but the two differ in the relative orientations around the rhodium and boron centres. DFT calculations on the actual complexes provide insight into the influence of the pi-accepting CO co-ligand, compared with previously reported complexes [Rh(mu-Cl)(dpb)]2 and [RhCl(dmap)(dpb)]. In addition, comparison of the nu(CO) frequency of [RhCl(CO)(dpb)] with that of the related borane-free complex [RhCl(CO)(iPr2PPh)2] substantiates the significant electron-withdrawing effect that the sigma-accepting borane moiety exerts on the metal. Valence isoelectronic [PtCl2(dpb)] and [PdCl2(dpb)] complexes have also been prepared and characterized spectroscopically and structurally. The pronounced influence of the transition metal on the magnitude of the M-->B interaction is highlighted by geometric considerations and NBO analyses.     

The Challenge of Linear (E)‐Enones in the Rh‐Catalyzed,Asymmetric 1,4‐Addition Reaction of Phenylboronic Acid: A DFT Computational Analysis

Why are linear (E)‐enones such challenging substrates in the Rh‐catalyzed asymmetric arylation with boronic acids, which is one of the most important asymmetric catalysis methods? DFT computations show that these substrates adopt a specific conformation in which the largest substituent is antiperiplanar to Rh^I π‐complexed with the C?C bond within the enantioselectivity‐determining carborhodation transition state. Additionally, for such structures, there is a strong, but not exclusive, preference for s‐cis enone conformation. This folding minimizes steric interactions between the substrate and the ligand, and hence reduces the enantioselectivity. This idea is further confirmed by investigating three computation‐only substrate “probes”, one of which is capable of double asymmetric induction, and a recent computationally designed 1,5‐diene ligand. On average, excellent agreement between predicted and experimental enantioselectivity was attained by a three‐pronged approach: 1) thorough conformational search within ligand and substrate subunits to locate the most preferred carborhodation transition state; 2) including dispersion interaction and long‐range corrections by SMD/ωB97xD/DGDZVP level of theory; and 3) full substrate and ligand modeling. Based on the results, a theory‐enhanced enantioselectivity model that is applicable to both chiral diene and diphosphane ligands is proposed.     

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.     

Basis set error estimation for DFT calculations of electronic g‐tensors for transition metal complexes

We present a detailed study of the basis set dependence of electronic g‐tensors for transition metal complexes calculated using Kohn–Sham density functional theory. Focus is on the use of locally dense basis set schemes where the metal is treated using either the same or a more flexible basis set than used for the ligand sphere. The performance of all basis set schemes is compared to the extrapolated complete basis set limit results. Furthermore, we test the performance of the aug‐cc‐pVTZ‐J basis set developed for calculations of NMR spin‐spin and electron paramagnetic resonance hyperfine coupling constants. Our results show that reasonable results can be obtain when using small basis sets for the ligand sphere, and very accurate results are obtained when an aug‐cc‐pVTZ basis set or similar is used for all atoms in the complex. © 2014 Wiley Periodicals, Inc.     

Mechanism of the dehydrogenative silylation of alcohols catalyzed by cationic gold complexes: an experimental and theoretical study

The catalytic activity both of cationic [(XDPP)Au][X] (XDPP = bis-2,5-diphenylphosphole xantphos X = BF(4)) and of the isolated gold hydride complex [(XDPP)(2)Au(2)H][OTf] in the dehydrogenative silylation process is presented. A parallel theoretical study using density functional theory revealed a mechanism involving the counter anion as a co-catalyst, which was experimentally confirmed by testing various counterions (X = OTf, NTf(2), PF(6)). Finally, a "Au(2)H(+)" species was determined as the intermediate during the catalytic cycle, which correlates well with the experimental findings on the first example of catalytic activity of an isolated "Au-H" complex.     

Spectral, structural and DFT studies of platinum group metal 3,6-bis(2-pyridyl)-4-phenylpyridazine complexes and their ligand bonding modes

Reactions of 3,6-bis(2-pyridyl)-4-phenylpyridazine (L^ph) with [(η⁶-arene)Ru(μ-Cl)Cl]₂ (arene = C₆H₆, p-ⁱPrC₆H₄Me and C₆Me₆), [(η⁵-C₅Me₅)M(μ-Cl)Cl]₂, (M = Rh and Ir) and [(η⁵-Cp)Ru(PPh₃)₂Cl] (Cp = C₅H₅, C₅Me₅ and C₉H₇) afford mononuclear complexes of the type [(η⁶-arene)Ru(L^ph)Cl]PF₆, [(η⁵-C₅Me₅)M(L^ph)Cl]PF₆ and [(Cp)Ru(L^ph)(PPh₃)]PF₆ with different structural motifs depending on the π-acidity of the ligand, electronic properties of the central metal atom and nature of the co-ligands. Complexes [(η⁶-C₆H₆)Ru(L^ph)Cl]PF₆1, [(η⁶-p-ⁱPrC₆H₄Me)Ru(L^ph)Cl]PF₆2, [(η⁵-C₅Me₅)Ir(L^ph)Cl]PF₆5, [(η⁵-Cp)Ru(PPh₃)(L^ph)]PF₆, (Cp = C₅H₅, 6; C₅Me₅, 7; C₉H₇, 8) show the type-A binding mode (see text), while complexes [(η⁶-C₆Me₆)Ru(L^ph)Cl]PF₆3 and [(η⁵-C₅Me₅)Rh(L^ph)Cl]PF₆4 show the type-B binding mode (see text). These differences reflect the more electron-rich character of the [(η⁶-C₆Me₆)Ru(μ-Cl)Cl]₂ and [(η⁵-C₅Me₅)Rh(μ-Cl)Cl]₂ complexes compared to the other starting precursor complexes. Binding modes of the ligand L^ph are determined by ¹H NMR spectroscopy, single-crystal X-ray analysis as well as evidence obtained from the solid-state structures and corroborated by density functional theory calculations. From the systems studied here, it is concluded that the electron density on the central metal atom of these complexes plays an important role in deciding the ligand binding sites.     

Dissecting the intimate mechanism of cyanation of {2Fe3S} complexes related to the active site of all-iron hydrogenases by DFT analysis of energetics, transition states, intermediates and products in the carbonyl substitution pathway

A bridging carbonyl intermediate with key structural elements of the diiron sub-site of all-iron hydrogenase has been experimentally observed in the CN/CO substitution pathway of the {2Fe3S} carbonyl precursor, [Fe(2)(CO)(5){MeSCH(2)C(Me)(CH(2)S)(2)}]. Herein we have used density functional theory (DFT) to dissect the overall substitution pathway in terms of the energetics and the structures of transition states, intermediates and products. We show that the formation of bridging CO transitions states is explicitly involved in the intimate mechanism of dicyanation. The enhanced rate of monocyanation of {2Fe3S} over the {2Fe2S} species [Fe(2)(CO)(6){CH(2)(CH(2)S)(2)}] is found to rest with the ability of the thioether ligand to both stabilise a mu-CO transition state and act as a good leaving group. In contrast, the second cyanation step of the {2Fe3S} species is kinetically slower than for the {2Fe2S} monocyanide because the Fe2 atom is deactivated by coordination of the electron-donating thioether group. In addition, hindered rotation and the reaction coordinate of the approaching CN(-) group, are other factors which explain reactivity differences in {2Fe2S} and {2Fe3S} systems. The intermediate species formed in the second cyanation step of {2Fe3S} species is a mu-CO species, confirming the structural assignment made on the basis of FT-IR data (S. J. George, Z. Cui, M. Razavet, C. J. Pickett, Chem. Eur. J. 2002, 8, 4037-4046). In support of this we find that computed and experimental IR frequencies of structurally characterised {2Fe3S} species and those of the bridging carbonyl intermediate are in excellent agreement. In a wider context, the study may provide some insight into the reactivity of dinuclear systems in which neighbouring group on-off coordination plays a role in substitution pathways.     

Unexpected direct iron-fluorine bonds in trifluorophosphane iron complexes: an alternative to bridging trifluorophosphane and difluorophosphido groups

The iron trifluorophosphane complexes [Fe(PF(3))(n)] (n=4, 5), [Fe(2)(PF(3))(n)] (n=8, 9), [H(2)Fe(PF(3))(4)], and [Fe(2)(PF(2))(2)(PF(3))(6)] have been studied by density functional theory. The lowest energy structures of [Fe(PF(3))(4)] and [Fe(PF(3))(5)] are a triplet tetrahedron and a singlet trigonal bipyramid, respectively. Both cis and trans octahedral structures were found for [H(2)Fe(PF(3))(4)] with the cis isomer lying lower in energy by approximately 10 kcal mol(-1). The lowest energy structure for [Fe(2)(PF(3))(8)] has two [Fe(PF(3))(4)] units linked only by an iron-iron bond of length 2.505 A consistent with the formal Fe=Fe double bond required to give both iron atoms the favored 18-electron configuration. In the lowest energy structure for [Fe(2)(PF(3))(9)] one of the iron atoms has inserted into a P-F bond of one of the PF(3) ligands to give a structure [(F(3)P)(4)Fe<--PF(2)Fe(F)(PF(3))(4)] with a bridging PF(2) group and a direct Fe-F bond. A bridging PF(3) group is found in a considerably higher energy [Fe(2)(PF(3))(9)] structure at approximately 30 kcal mol(-1) above the global minimum. However, this bridging PF(3) group keeps the two iron atoms too far apart (approximately 4 A) for the direct iron-iron bond required to give the iron atoms the favored 18-electron configuration. The preferred structure for [Fe(2)(PF(2))(2)(PF(3))(6)] has a bridging PF(2) group, as expected. However, this bridging PF(2) group bonds to one of the iron atoms through an P-Fe covalent bond and to the other iron through an F-->Fe dative bond, leaving an uncomplexed phosphorus lone pair.     

Syntheses, reactivity and DFT studies of group 2 and group 12 metal complexes of tris(pyrazolyl)methanides featuring "free" pyramidal carbanions

Reactions of HC(Me2pz)3 with Grignard reagents, dialkyl magnesium compounds and dimethylzinc are reported, together with a DFT study on some of the aspects of this chemistry. Reactions of HC(Me2pz)3 with MeMgX (X=Cl or Br) gave the half-sandwich zwitterionic compounds [Mg((Me)Tpmd)X] (X=Cl (2) or Br (3); (Me)Tpmd(-)=[C(Me2pz)3](-)). Addition of HCl to 2 gave the structurally characterised half-sandwich compound [Mg{HC(Me2pz)3}Cl2(thf)] (4). The zwitterionic sandwich compound [Mg(MeTpmd)2] (5) formed in low yields in the reaction of MeMgX with HC(Me2pz)3 but was readily prepared from HC(Me2pz)3 and either MgnBu2 or MgPh2. The structurally characterised compound 5 contains two "naked" sp3-hybridised carbanions fully separated from the dicationic metal centre. Only by using MgPh2 as starting material could the half-sandwich compound [Mg(MeTpmd)Ph(thf)] (6) be isolated. The zwitterionic sandwich compound 5 reacted with HOTf (OTf(-)=[O3SCF3](-)) to form the dication [Mg{HC(Me2pz)3}2]2+ (7(2+)), which was structurally characterised. Pulsed field gradient spin-echo (PGSE) diffusion NMR spectroscopy revealed both compounds to be intact in solution. In contrast to the magnesium counterparts, HC(Me2pz)3 reacted only slowly with ZnMe2 (and not at all with ZnPh2) to form the half-sandwich zwitterion [Zn(MeTpmd)Me] (8), which contains a cationic methylzinc moiety separated from a single sp3-hybridised carbanion. Density functional calculations on the zwitterions [M(MeTpmd)Me] and [M(MeTpmd)2] (M=Mg, Zn) revealed that the HOMO in each case is a (Me)Tpmd-based carbanion lone pair. The kappa 1C isomers of [M(MeTpmd)Me] were calculated to be considerably less stable than their kappa 3N-bound counterparts, with the largest gain in energy for Mg due to the greater ease of electron transfer from metal to the (Me)Tpmd apical carbon atom on formation of the zwitterion. Moreover, the computed M-C bond dissociation enthalpies of the kappa 1C isomers of [M(MeTpmd)Me] are considerably higher than expected by simple extrapolation from the corresponding computed H-C bond dissociation enthalpy.     

Extension of the PS-GVB electronic structure code to transition metal complexes

We have developed a parameterization enabling ab initio electronic structure calculation via the PS-GVB program on transition-metal-containing systems using two standard effective core potential basis sets. Results are compared with Gaussian-92 for a wide range of complexes, and superior performance is demonstrated with regard to computational efficiency for single-point energies and geometry optimization. Additionally, the initial guess strategy in PS-GVB is shown to provide considerably more reliable convergence to the ground state. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1863–1874, 1997     

Electronic structure and coordination chemistry of phenanthridine ligand in first‐row transition metal complexes: A DFT study

The geometric parameters, electronic structures, and haptotropic migration of a series of hypothetical compounds of general formula CpM(C₁₃H₉N) and (CO)₃M(C₁₃H₉N) (M = fist row transition metal, Cp = C₅H₅, and C₁₃H₉N = phenanthridine ligand) are investigated by means of the density functional theory. The phenanthridine ligand can bind to the metal through η¹ to η⁶ coordination mode, in agreement with the electron count and the nature of the metal, showing its capability to adapt itself to the electronic demand of the metal as well as to the polycyclic aromatic hydrocarbons. In the investigated species, the most favored closed‐shell count is 18‐electron except for the Ti and V models which are deficient open‐shell 16‐electron configuration. This study has shown the difference in coordination ability of this heteropolycyclic ligand: the coordination of the central C₅N ring is less favored than the terminal C₆ rings, in agreement with the π‐electron density localization. Most of the investigated complexes are expected to exhibit a rich fluxional behavior. This flexibility favors the possibility for the existence of several isomers as well as their interconversion through haptotropic shifts. © 2012 Wiley Periodicals, Inc.     

Description of the ground-state covalencies of the bis(dithiolato) transition-metal complexes from X-ray absorption spectroscopy and time-dependent density-functional calculations

The electronic structures of [M(L(Bu))(2)](-) (L(Bu)=3,5-di-tert-butyl-1,2-benzenedithiol; M=Ni, Pd, Pt, Cu, Co, Au) complexes and their electrochemically generated oxidized and reduced forms have been investigated by using sulfur K-edge as well as metal K- and L-edge X-ray absorption spectroscopy. The electronic structure content of the sulfur K-edge spectra was determined through detailed comparison of experimental and theoretically calculated spectra. The calculations were based on a new simplified scheme based on quasi-relativistic time-dependent density functional theory (TD-DFT) and proved to be successful in the interpretation of the experimental data. It is shown that dithiolene ligands act as noninnocent ligands that are readily oxidized to the dithiosemiquinonate(-) forms. The extent of electron transfer strongly depends on the effective nuclear charge of the central metal, which in turn is influenced by its formal oxidation state, its position in the periodic table, and scalar relativistic effects for the heavier metals. Thus, the complexes [M(L(Bu))(2)](-) (M=Ni, Pd, Pt) and [Au(L(Bu))(2)] are best described as delocalized class III mixed-valence ligand radicals bound to low-spin d(8) central metal ions while [M(L(Bu))(2)](-) (M=Cu, Au) and [M(L(Bu))(2)](2-) (M=Ni, Pd, Pt) contain completely reduced dithiolato(2-) ligands. The case of [Co(L(Bu))(2)](-) remains ambiguous. On the methodological side, the calculation led to the new result that the transition dipole moment integral is noticeably different for S(1s)-->valence-pi versus S(1s)-->valence-sigma transitions, which is explained on the basis of the differences in radial distortion that accompany chemical bond formation. This is of importance in determining experimental covalencies for complexes with highly covalent metal-sulfur bonds from ligand K-edge absorption spectroscopy.     

Transition metal hydrazide-based hydrogen-storage materials: the first atoms-in-molecules analysis of the Kubas interaction

Molecular models of the M-H(2) binding sites of experimentally characterised amorphous vanadium hydrazide gels are studied computationally using gradient corrected density functional theory, to probe the coordination number of the vanadium in the material and the nature of the interaction between the metal and the H(2) molecules. The H(2) is found to bind to the vanadium through the Kubas interaction, and the first quantum theory of atoms-in-molecules analysis of this type of interaction is reported. Strong correlation is observed between the electron density at the H-H bond critical point and the M-H(2) interaction energy. Four coordinate models give the best reproduction of the experimental data, suggesting that the experimental sites are four coordinate. The V-H(2) interaction is shown to be greater when the non-hydrazine based ligand, THF, of the experimental system is altered to a poorer π-acceptor ligand. Upon altering the metal to Ti or Cr the M-H(2) interaction energy changes little but the number of H(2) which may be bound decreases from four (Ti) to two (Cr). It is proposed that changing the metal from V to Ti may increase the hydrogen storage capacity of the experimental system. A 9.9 wt% maximum storage capacity at the ideal binding enthalpy for room temperature performance is predicted when the Ti metal is combined with a coordination sphere containing 2 hydride ligands.