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.     

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.     

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.     

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.     

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.     

On the theory of ligand substitution in transition metal complexes

Expressions are given for the calculation of activation energies for ligand substitution in transition metal complexes. A comparison between calculated and experimental activation energies for the substitution of water molecules shows reasonably good agreement.     

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.     

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.     

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     

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.     

Fluorine... and pi...alkali metal interactions control in the stereoselective amide enolate alkylation with fluorinated oxazolidines (Fox) as a chiral auxiliary: an experimental and theoretical study

The alpha-alkylation of amide enolates by using a pseudo-C(2) symmetry trans 4-phenyl-2-trifluoromethyloxazolidine (trans-Fox) as a chiral auxiliary occurs with an extremely high diastereoselectivity (>99 % de). The origin of this excellent stereocontrol was investigated by an experimental and theoretical (DFT) study. With this trans chiral auxiliary, both F...metal and pi...metal interactions compete to give the same diastereomer through Re face alkylation of the enolate. A 5.5 kcal mol(-1) energy difference found between the Re face and the Si face attack transition states is consistent with the complete diastereoselectivity that has been experimentally achieved. On the other hand, in the case of the cis chiral auxiliary (cis-Fox) the competition between the F...metal and pi...metal interactions is unfavourable to the diastereoselectivity. In this case, the Re face and the Si face attack transition states were found to be nearly isoenergetic (0.3 kcal mol(-1) difference), which is in good agreement with the very low diastereoselectivity observed.     

Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni,Pd)-Based Transition-Metal Catalysts

The performance of 23 density functionals, including one LDA, four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid meta-GGAs, and ten double-hybrid functionals, was investigated for the computation of activation energies of various covalent main-group single bonds by four catalysts: Pd, PdCl⁻, PdCl₂, and Ni (all in the singlet state). A reactant complex, the barrier, and reaction energy were considered, leading to 164 energy data points for statistical analysis. Extended Gaussian AO basis sets were used in all calculations. The best functional for the complete benchmark set relative to estimated CCSD(T)/CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal mol⁻¹ followed by PW6B95-D3, the double hybrid PWPB95-D3, and B3LYP-D3 (1.9 kcal mol⁻¹ each). The other tested hybrid meta-GGAs perform less well (M06-HF: 7.0 kcal mol⁻¹; M06-2X: 6.3 kcal mol⁻¹; M06: 4.9 kcal mol⁻¹) for the investigated reactions. In the Ni case, some double hybrids show larger errors due to partial breakdown of the perturbative treatment for the correlation energy in cases with difficult electronic structures (partial multi-reference character). Only double hybrids either with very low amounts of perturbative correlation (e.g., PBE0-DH) or that use the opposite-spin correlation component only (e.g., PWPB95) seem to be more robust. We also investigated the effect of the D3 dispersion correction. While the barriers are not affected by this correction, significant and mostly positive results were observed for reaction energies. Furthermore, six very recently proposed double-hybrid functionals were analyzed regarding the influence of the amount of Fock exchange as well as the type of perturbative correlation treatment. According to these results, double hybrids with <50–60 % of exact exchange and ∼30 % perturbative correlation perform best.     

Origin of the regio- and stereoselectivity in palladium-catalyzed electrophilic substitution via bis(allyl)palladium complexes

Palladium-catalyzed allylic substitution of aryl allyl chlorides with aromatic and heteroaromatic aldehydes was performed in the presence of hexamethylditin. This procedure involves palladium-catalyzed formation of transient allylstannanes followed by generation of a bis(allyl)palladium intermediate, which subsequently reacts with the aldehyde electrophile. The catalytic substitution reaction proceeds with high regio- and stereoselectivity. The stereoselectivity is affected by the steric and electronic properties of the allylic substituents. Various functionalities including NO(2), COCH(3), Br, and F groups are tolerated under the applied catalytic conditions. Density functional calculations at the B3PW91/DZ+P level of theory were applied to study the steric and electronic effects controlling the regio- and stereoselectivity of the electrophilic addition. The development of the selectivity was studied by modeling the various bis(allyl)palladium species occurring in the palladium-catalyzed substitution of cinnamyl chloride with benzaldehyde. It was found that the electrophilic attack proceeds via a six-membered cyclic transition state, which has a pronounced chair conformation. The regioselectivity of the reaction is controlled by the location of the phenyl group on the eta(1)-allyl moiety of the complex. The stereoselectivity of the addition process is determined by the relative configuration of the phenyl substituents across the developing carbon-carbon bond. The lowest energy path corresponds to the formation of the branched allylic isomer with the phenyl groups in anti configuration, which is in excellent agreement with the experimental findings.     

The rate-determining step in the rhodium-xantphos-catalysed hydroformylation of 1-octene

The rate-determining step in the hydroformylation of 1-octene, catalysed by the rhodium-Xantphos catalyst system, was determined by using a combination of experimentally determined (1)H/(2)H and (12)C/(13)C kinetic isotope effects and a theoretical approach. From the rates of hydroformylation and deuterioformylation, a small (1)H/(2)H isotope effect of 1.2 was determined for the hydride moiety of the rhodium catalyst. (12)C/(13)C isotope effects of 1.012(1) and 1.012(3) for the alpha-carbon and beta-carbon atoms of 1-octene were determined, respectively. Both quantum mechanics/molecular mechanics (QM/MM) and full quantum mechanics calculations were carried out on the key catalytic steps, for "real-world" ligand systems, to clarify whether alkene coordination or hydride migration is the rate-determining step. Our calculations (21.4 kcal mol(-1)) quantitatively reproduce the experimental energy barrier for CO dissociation (20.1 kcal mol(-1)) starting at the (bisphosphane)RhH(CO)(2) resting state. The barrier for hydride migration lies 3.8 kcal mol(-1) higher than the barrier for CO dissociation (experimentally determined trend approximately 3 kcal mol(-1)). The computed (1)H/(2)H and (12)C/(13)C kinetic isotope effects corroborate the results of the energy analysis.