Theoretical studies in palladium and platinum molecular chemistry

Bond activation by group-11 transition-metal cations

We have computationally explored C–X bond activation by the group-11 transition-metal cations Cu+, Ag+, and Au+, and, for comparison, Pd, using relativistic density functional theory (DFT) at ZORA-BLYP/TZ2P. Oxidative insertion of the second-row transition-metal species Ag+ and Pd leads, for a given bond, to the highest overall reaction barriers. On the other hand, if we compare the different bonds oxidative insertion into the C–F bond is associated with (one of the) highest overall barriers whereas insertion into the C–Cl bond leads to the lowest overall barrier for any transition metal. The main trends in reactivity are rationalized using the activation strain model of chemical reactivity, which is an extension of the fragment approach to reaction profiles. In this model, the shape of the reaction profile ΔE(ζ) and the height of the overall reaction barrier ΔE≠ = ΔE(ζ=ζTS) are interpreted in terms of the strain energy ΔEstrain(ζ) associated with deforming the reactants along the reaction coordinate ζ plus the interaction energy ΔEint(ζ) between these deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ).

Mechanisms of C-H bond activation: rich synergy between computation and experiment

Recent computational studies of C-H bond activation at late transition metal systems are discussed and processes where lone pair assistance via heteroatom co-ligands or carboxylates are highlighted as a particularly promising means of cleaving C-H bonds. The term 'ambiphilic metal ligand activation' (AMLA) is introduced to describe such reactions.

Detection of platinum dihydride bisphosphine complexes and studies of their reactivity through para-hydrogen-enhanced NMR methods

In-situ NMR studies on the reactions of Pt{CH2 = CHSi(Me)2}2O)(PCy3) with phosphines, HSiEt3 and--hydrogen or Pt(L)(L')(Me)(2) alone enable the detection of cis-Pt(L)(L')(H)2 [L = PCy3 and L' = PCy2H, PPh3 or PCy3] which then undergo hydride site interchange and H2 reductive elimination on the NMR timescale.

The platinum-olefin binding energy in series of (PH3)2Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R1R2C = CR3R4)Pt(PH3)2 [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents.

Catalytic carbophilic activation: catalysis by platinum and gold pi acids

The ability of platinum and gold catalysts to effect powerful atom-economic transformations has led to a marked increase in their utilization. The quite remarkable correlation of their catalytic behavior with the available structural data, coordination chemistry, and organometallic reactivity patterns, including relativistic effects, allows the underlying principles of catalytic carbophilic activation by pi acids to be formulated. The spectrum of reactivity extends beyond their utility as catalytic and benign alternatives to conventional stoichiometric pi acids. The resulting reactivity profile allows this entire field of catalysis to be rationalized, and brings together the apparently disparate electrophilic metal carbene and nonclassical carbocation explanations. The advances in coupling, cycloisomerization, and structural reorganization--from the design of new transformations to the improvement to known reactions--are highlighted in this Review. The application of platinum- and gold-catalyzed transformations in natural product synthesis is also discussed.

Transition-State Energy and Position along the Reaction Coordinate in an Extended Activation Strain Model

We investigate palladium-induced activation of the C-H, C-C, C-F, and C-Cl bonds in methane, ethane, cyclopropane, fluoromethane, and chloromethane, using relativistic density functional theory (DFT) at ZORA-BLYP/TZ2P. Our purpose is to arrive at a qualitative understanding, based on accurate calculations, of the trends in activation barriers and transition state (TS) geometries (e.g. early or late along the reaction coordinate) in terms of the reactants' properties. To this end, we extend the activation strain model (in which the activation energy Delta E(not equal) is decomposed into the activation strain Delta E(not equal)(strain) of the reactants and the stabilizing TS interaction Delta E(not equal)(int) between the reactants) from a single-point analysis of the TS to an analysis along the reaction coordinate zeta, that is, Delta E(zeta)=Delta E(strain)(zeta)+Delta E(int)(zeta). This extension enables us to understand qualitatively, trends in the position of the TS along zeta and, therefore, the values of the activation strain Delta E(not equal)(strain)=Delta E(strain)(zeta(TS)) and TS interaction Delta E(not equal)(int)=Delta E(int)(zeta(TS)) and trends therein. An interesting insight that emerges is that the much higher barrier of metal-mediated C-C versus C-H activation originates from steric shielding of the C-C bond in ethane by C-H bonds. Thus, before a favorable stabilizing interaction with the C-C bond can occur, the C-H bonds must be bent away, which causes the metal-substrate interaction Delta E(int)(zeta) in C-C activation to lag behind. Such steric shielding is not present in the metal-mediated activation of the C-H bond, which is always accessible from the hydrogen side. Other phenomena that are addressed are anion assistance, competition between direct oxidative insertion (OxIn) versus the alternative S(N)2 pathway, and the effect of ring strain.

Ester versus Polyketone Formation in the Palladium−Diphosphine Catalyzed Carbonylation of Ethene

The origin of the chemoselectivity of palladium catalysts containing bidentate phosphine ligands toward either methoxycarbonylation of ethene or the copolymerization of ethene and carbon monoxide was investigated using density functional theory based calculations. For a palladium catalyst containing the electron-donating bis(dimethylphosphino)ethane (dmpe) ligand, the rate determining step for chain propagation is shown to be the insertion of ethene into the metal-acyl bond. The high barrier for chain propagation is attributed to the low stability of the ethene intermediate, (dmpe)Pd(ethene)(C(O)CH3). For the competing methanolysis process, the most likely pathway involves the formation of (dmpe)Pd(CH3OH)(C(O)CH3) via dissociative ligand exchange, followed by a solvent mediated proton-transfer/reductive- elimination process. The overall barrier for this process is higher than the barrier for ethene insertion into the palladium-acetyl bond, in line with the experimentally observed preference of this type of catalyst toward the formation of polyketone. Electronic bite angle effects on the rates of ethene insertion and ethanoyl methanolysis were evaluated using four electronically and sterically related ligands (Me)2P(CH2)nP(Me)2 (n = 1-4). Steric effects were studied for larger tert-butyl substituted ligands using a QM/MM methodology. The results show that ethene coordination to the metal center and subsequent insertion into the palladium-ethanoyl bond are disfavored by the addition of steric bulk around the metal center. Key intermediates in the methanolysis mechanism, on the other hand, are stabilized because of electronic effects caused by increasing the bite angle of the diphosphine ligand. The combined effects explain successfully which ligands give polymer and which ones give methyl propionate as the major products of the reaction.

Pd–H elimination reactions in palladium(ii) allylic complexes

The ease of H elimination from the 4- (beta-) position in a series of allylic complexes [Pd(5-C(6)F(5)-1-3-eta(3)-cyclohexenyl)XL](n+) (X, L = Br, N-, P-donor, C(6)Cl(2)F(3); n = -1, 0, +1) was compared by analyzing their decomposition products at 50 degrees C. Pd-H elimination does not occur for cationic complexes, whereas it is the main decomposition pathway for neutral and anionic complexes. In addition to the charge of the complex, the ease of this Pd-H elimination is determined by the trans influence of the ligands (aryl > PMe(3) > Br, N-donor).

The relativistic complete active-space second-order perturbation theory with the four-component Dirac Hamiltonian

The relativistic complete active-space second-order perturbation theory (CASPT2) is developed for the four-component relativistic Hamiltonian. The present method can describe the near-degenerated and dissociated electronic states of molecules involving atoms of heavy elements. The present approach is less expensive than the relativistic multireference configuration interaction method. The ground and low-lying excited states of TlH, Tl(2), and PtH molecules are calculated with the Dirac-Coulomb (DC) CASPT2 method and their spectroscopic constants are obtained. These spectroscopic constants are compared with experimental findings and previous theoretical work. For all the molecules, the spectroscopic constants of DC-CASPT2 show good agreement with the experimental or previous theoretical spectroscopic constants. The present theory provides accurate descriptions of bonding or dissociation states and of ground and excited states in a well-balanced way.

DFT Study on the Palladium-Catalyzed Allylation of Primary Amines by Allylic Alcohol

The palladium-catalyzed allylation of primary amines has been investigated by DFT calculations (B3PW91, PCM method), and two potential mechanisms were studied. The first mechanism relies on the formation of cationic hydridopalladium complexes. Their formation involves a metal-assisted formal (1,3) shift of a proton from the nitrogen atom of an ammonium to the Cbeta carbon atom. The second part of the cycle relies on a ligand exchange through a pentacoordinated 18VE hydridopalladium complex. The last step likely proceeds through a bimolecular pathway and formally consists of a proton transfer from the allylammonium to the alcohol group of the complex. The second mechanism, which is closer to that currently admitted for nucleophilic allylic substitutions, relies on the decomplexation of the coordinated allylammonium and appears to be favored. This catalytic cycle was recomputed on model complexes varying the ligands, and a charge decomposition analysis was carried out to assess the influence of the electronic properties of the ligands. To compare our results with competitive experiments, CDA calculations were also performed on real ligands. In agreement with experimental observations, this process was found to be strongly ligand dependent, decomplexation being favored by strong pi-acceptor ligands. These calculations led us to show experimentally that complex [Pd(P(OPh)(3))(2)(eta(3)-C(3)H(5))][OTf] is an efficient catalyst for this allylation. Finally, this catalytic process proved to be sensitive to the nature of the amine, with poorly basic amines favoring the re-formation of the catalytic precursor.

Natural orbitals for chemical valence as descriptors of chemical bonding in transition metal complexes

Natural orbitals for chemical valence (NOCV) are defined as the eigenvectors of the chemical valence operator defined by Nalewajski et al.; they decompose the deformation density (differential density, Deltarho) into diagonal contributions. NOCV were used in a description of the chemical bond between the organometallic fragment and the ligand in example transition-metal complexes: heme-CO ([FeN(5)C(20)H(15)]-CO), [Ni-diimine hydride]-ethylene ([N;N-Ni-H]-C(2)H(4), N;N = -NH-CH-CH-NH-), and [Ni(NH(3))(3)]-CO. DFT calculations were performed using gradient-corrected density functional theory (DFT) in the fragments resolution, using the fragment/ligand Kohn-Sham orbitals as a basis set in calculations for the whole fragment-ligand complex. It has been found that NOCV lead to a very compact description of the fragment-ligand bond, with only a few orbitals exhibiting non-zero eigenvalues. Results of NOCV analysis, compared with Mulliken populations analysis and Zigler-Rauk interaction-energy decomposition, demonstrate that the use of the natural valence orbitals allows for a separation of the sigma-donation and pi-back-donation contributions to the ligand-fragment bond. They can be also useful in comparison of these contributions in different complexes.

Acetylenic Carbon-13 Chemical Shift Tensors for Diphenylacetylene and (η2-Diphenylacetylene)Pt(PPh3)2: A Solid-State NMR and Theoretical Study

The structure of (eta2-diphenylacetylene)Pt(PPh3) (2), as well as those of its dichloromethane and benzene solvates, is determined via X-ray crystallography. An investigation of the chemical shift (CS) tensors of the 13C-labeled carbons in Ph13C13CPh and (eta2-Ph13C13CPh)Pt(PPh3)2.(C6H6) is carried out via analysis of 13C NMR spectra from stationary solid samples. The principal components of the CS tensors as well as their orientations with respect to the 13C,13C internuclear vector are determined. DFT calculations of these CS tensors are in close agreement with the experimental values. For diphenylacetylene (tolane), the orientations and principal-component magnitudes of the alkynyl carbon CS tensors are comparable to those for other alkynyl carbons, although the CS tensor is not axially symmetric in this case. Coordination to platinum causes a change in the CS tensor orientation and a net increase in the isotropic chemical shift, resulting from a significant increase in two principal components (delta11 and delta33) while the third (delta22) decreases only slightly. The measured carbon CS tensors in the platinum complex bear a striking similarity to those of the alkenyl carbons in trans-Ph(H)C=C(H)Ph, and a short theoretical discussion of these observations is presented.

Synthesis and Theoretical Study of a Series of Dipalladium(I) Complexes Containing the Pd2(μ-CO)2 Core

The complex [PBu4]2[Pd2(mu-CO)2Cl4] has been prepared in high yields by carbonylation of [PBu4]2[Pd2Cl6]. Methanol, potassium acetate, or CO readily reacted under ambient conditions to quantitatively afford a series of dipalladium(I) complexes, namely [Pd2(mu-CO)2Cl3(OCH3)]2-, [Pd2(mu-CO)2Cl3(OC(O)CH3)]2-, [Pd2(mu-CO)2Cl3(CO)]-, and [Pd2(mu-CO)2Cl2(OCH3)(CO)]-, all of which have the Pd2(mu-CO)2 core preserved. All these complexes have been characterized by infrared and NMR spectroscopies; the high nu(CO) stretching wavenumbers observed and the diamagnetic character of these complexes prompted us to perform theoretical calculations to describe the electronic structure of the Pd2(mu-CO)2 core and to gain an intimate description of the Pd-CO bonds. The pairing of the two lonely electrons of the Pd d9 atoms is due to the delocalization along the CO bridging ligands.

A Two-State Computational Investigation of Methane CH and Ethane CC Oxidative Addition to [CpM(PH3)]n+ (M=Co, Rh, Ir;n=0, 1)

Reductive elimination of methane from methyl hydride half-sandwich phosphane complexes of the Group 9 metals has been investigated by DFT calculations on the model system [CpM(PH(3))(CH(3))(H)] (M = Co, Rh, Ir). For each metal, the unsaturated product has a triplet ground state; thus, spin crossover occurs during the reaction. All relevant stationary points on the two potential energy surfaces (PES) and the minimum energy crossing point (MECP) were optimized. Spin crossover occurs very near the sigma-CH(4) complex local minimum for the Co system, whereas the heavier Rh and Ir systems remain in the singlet state until the CH(4) molecule is almost completely expelled from the metal coordination sphere. No local sigma-CH(4) minimum was found for the Ir system. The energetic profiles agree with the nonexistence of the Co(III) methyl hydride complex and with the greater thermal stability of the Ir complex relative to the Rh complex. Reductive elimination of methane from the related oxidized complexes [CpM(PH(3))(CH(3))(H)](+) (M = Rh, Ir) proceeds entirely on the spin doublet PES, because the 15-electron [CpM(PH(3))](+) products have a doublet ground state. This process is thermodynamically favored by about 25 kcal mol(-1) relative to the corresponding neutral system. It is essentially barrierless for the Rh system and has a relatively small barrier (ca. 7.5 kcal mol(-1)) for the Ir system. In both cases, the reaction involves a sigma-CH(4) intermediate. Reductive elimination of ethane from [CpM(PH(3))(CH(3))(2)](+) (M = Rh, Ir) shows a similar thermodynamic profile, but is kinetically quite different from methane elimination from [CpM(PH(3))(CH(3))(H)](+): the reductive elimination barrier is much greater and does not involve a sigma-complex intermediate. The large difference in the calculated activation barriers (ca. 12.0 and ca. 30.5 kcal mol(-1) for the Rh and Ir systems, respectively) agrees with the experimental observation, for related systems, of oxidatively induced ethane elimination when M = Rh, whereas the related Ir systems prefer to decompose by alternative pathways.

DFT studies of the methyl exchange reaction between Cp2M–CH3or Cp*2M–CH3(Cp = C5H5, Cp* = C5Me5, M = Y, Sc, Ln) and CH4. Does M ionic radius control the reaction?

The activation energies for the methyl exchange reactions between Cp2M-CH3 and H-CH3 have been calculated for M = Sc, Y and representative metals of the lanthanide family (La, Ce, Sm, Ho, Yb and Lu) with DFT(B3PW91) calculations with large-core pseudopotentials for M. The sigma-bond metathesis reactions are calculated to have lower activation energies for early lanthanides than for late lanthanides and any of group 3 metals. The relative activation barriers are analyzed using the NBO charge distributions in the reactant and in the transition states. It is shown that the methane needs to be polarized in the transition state as H((+delta))-CH3((-delta)) by the reactant, because this sigma-bond metathesis is best viewed as heterolytic cleavage of methane, leading to a proton transfer between two methyl groups in the field of an electropositive M metal. Early lanthanides, which are involved in strongly ionic metal-ligands bonds are thus associated with the lowest activation energies. The ionic radius and the steric effects influence the relative rates of reaction for the complexes of Sc, Y and Lu. In agreement with earlier works of Sherer et al., the experimental reactivity trends found by Tilley are reproduced best with Cp*2M-CH3 (Cp* = C5Me5) rather than Cp2M-CH3 (Cp = C5H5) because the steric bulk of C5Me5 deactivates most the complex where the metal has the smallest ionic radius (Sc). While the steric effects and the influence of the metal ionic radius cannot be neglected, these factors are not the only ones involved in determining the activation barriers of the sigma-bond metathesis reaction.

Theoretical Study of C−H and C−F Activation in CH4-nFn (n =1−4) Molecules by Platinum

CASSCF followed by MRMP2 calculations have been carried out to analyze the reactions of a naked platinum atom with the fluorocarbon compounds CH(4-n)Fn (n = 1-4). For each of these interactions the potential-energy surfaces which correlate with the triplet ground state and the first excited singlet state of the free fragments were investigated for representative states evolving from different approaching modes of the reactants. For all the fluorinated fragments activation of the C-H and C-F bonds by the metal is strongly determined by the low-multiplicity channels arising from the first excited asymptote. Although stable products are predicted for insertion of the metallic atom into both the C-H and the C-F bonds of the different fluorocarbon compounds, comparison between the calculated energy barriers for reactions taking place in the same fluorinated molecule suggests in all cases a kinetic preference for the C-H bond oxidative addition to the platinum atom.

Oxidative Addition of the Fluoromethane C−F Bond to Pd. An ab Initio Benchmark and DFT Validation Study

We have computed a state-of-the-art benchmark potential energy surface (PES) for two reaction pathways (oxidative insertion, OxIn, and S(N)2) for oxidative addition of the fluoromethane C-F bond to the palladium atom and have used this to evaluate the performance of 26 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing these reactions. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) in combination with a hierarchical series of seven Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -5.3 (-6.1) kcal/mol for the formation of the reactant complex, 27.8 (25.4) kcal/mol for the activation energy for oxidative insertion (OxIn) relative to the separate reactants, 37.5 (31.8) kcal/mol for the activation energy for the alternative S(N)2 pathway, and -6.4 (-7.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.4-2.7 kcal/mol and errors in activation energies ranging from 0.3 to 2.8 kcal/mol. The B3LYP functional compares very well with a slight underestimation of the overall barrier for OxIn by -0.9 kcal/mol. For comparison, the well-known BLYP functional underestimates the overall barrier by -10.1 kcal/mol. The relative performance of these two functionals is inverted with respect to previous findings for the insertion of Pd into the C-H and C-C bonds. However, all major functionals yield correct trends and qualitative features of the PES, in particular, a clear preference for the OxIn over the alternative S(N)2 pathway.