Kinetics of the Trans Effect in Ruthenium Complexes Provide Insight into the Factors That Control Activity and Stability in CO2 Electroreduction.

Comparative kinetic studies of a series of new ruthenium complexes provide a platform for understanding how strong trans effect ligands and redox-active ligands work together to enable rapid electrochemical CO2 reduction at moderate overpotential. After synthesizing isomeric pairs of ruthenium complexes featuring 2'-picolinyl-methyl-benzimidazol-2-ylidene (Mebim-pic) as a strong trans effect ligand and 2,2':6',2″-terpyridine (tpy) as a redox-active ligand, chemical and electrochemical kinetic studies examined how complex geometry and charge affect the individual steps and overall catalysis of CO2 reduction. The relative trans effect of picoline vs the N-heterocyclic carbene (NHC) was quantified through a kinetic analysis of reductively triggered chloride dissociation, revealing that chloride loss is 1000 times faster in the isomer with the NHC trans to chloride. The kinetics of CO dissociation from a site trans to the NHC were examined in a systematic study of isostructural carbonyl complexes across four different overall charges. The rate constants for CO loss span 12 orders of magnitude and are fastest upon two-electron reduction, leading to a hypothesis that redox-active ligands play a key role in promoting reductive CO dissociation during catalysis. Analogous studies of complexes featuring the picoline ligand trans to the carbonyl reveal the importance of the trans effect of the CO ligand itself, with picoline ligand dissociation observed upon reduction. The complexes with NHC trans to the active site proved to be active electrocatalysts capable of selective CO2 electroreduction to CO. In acidic solutions under a N2 atmosphere, on the other hand, H2 evolution proceeds via an intermediate that positions a hydride ligand trans to picoline. The mechanistic insight and quantitative kinetic parameters that arise from these studies help establish general principles for molecular electrocatalyst design.

The Trans Effect in Electrocatalytic CO2 Reduction: Mechanistic Studies of Asymmetric Ruthenium Pyridyl-Carbene Catalysts.

A comprehensive mechanistic study of electrocatalytic CO2 reduction by ruthenium 2,2':6',2″-terpyridine (tpy) pyridyl-carbene catalysts reveals the importance of stereochemical control to locate the strongly donating N-heterocyclic carbene ligand trans to the site of CO2 activation. Computational studies were undertaken to predict the most stable isomer for a range of reasonable intermediates in CO2 reduction, suggesting that the ligand trans to the reaction site plays a key role in dictating the energetic profile of the catalytic reaction. A new isomer of [Ru(tpy)(Mebim-py)(NCCH3)]2+ (Mebim-py is 1-methylbenzimidazol-2-ylidene-3-(2'-pyridine)) and both isomers of the catalytic intermediate [Ru(tpy)(Mebim-py)(CO)]2+ were synthesized and characterized. Experimental studies demonstrate that both isomeric precatalysts facilitate electroreduction of CO2 to CO in 95/5 MeCN/H2O with high activity and high selectivity. Cyclic voltammetry, infrared spectroelectrochemistry, and NMR spectroscopy studies provide a detailed mechanistic picture demonstrating an essential isomerization step in which the N-trans catalyst converts in situ to the C-trans variant. Insight into molecular electrocatalyst design principles emerge from this study. First, the use of an asymmetric ligand that places a strongly electron-donating ligand trans to the site of CO2 binding and activation is critical to high activity. Second, stereochemical control to maintain the desired isomer structure during catalysis is critical to performance. Finally, pairing the strongly donating pyridyl-carbene ligand with the redox-active tpy ligand proves to be useful in boosting activity without sacrificing overpotential. These design principles are considered in the context of surface-immobilized electrocatalysis.

Origin of molecular conformational stability: perspectives from molecular orbital interactions and density functional reactivity theory.

To have a quantitative understanding about the origin of conformation stability for molecular systems is still an unaccomplished task. Frontier orbital interactions from molecular orbital theory and energy partition schemes from density functional reactivity theory are the two approaches available in the literature that can be used for this purpose. In this work, we compare the performance of these approaches for a total of 48 simple molecules. We also conduct studies to flexibly bend bond angles for water, carbon dioxide, borane, and ammonia molecules to obtain energy profiles for these systems over a wide range of conformations. We find that results from molecular orbital interactions using frontier occupied orbitals such as the highest occupied molecular orbital and its neighbors are only qualitatively, at most semi-qualitatively, trustworthy. To obtain quantitative insights into relative stability of different conformations, the energy partition approach from density functional reactivity theory is much more reliable. We also find that the electrostatic interaction is the dominant descriptor for conformational stability, and steric and quantum effects are smaller in contribution but their contributions are indispensable. Stable molecular conformations prefer to have a strong electrostatic interaction, small molecular size, and large exchange-correlation effect. This work should shed new light towards establishing a general theoretical framework for molecular stability.

Stability and dynamic processes in 16VE iridium(III) ethyl hydride and rhodium(I) σ-ethane complexes: experimental and computational studies.

Iridium(I) and rhodium(I) ethyl complexes, (PONOP)M(C2H5) (M = Ir (1-Et), Rh (2-Et)) and the iridium(I) propyl complex (PONOP)Ir(C3H7) (1-Pr), where PONOP is 2,6-(tBu2PO)2C5H3N, have been prepared. Low-temperature protonation of the Ir complexes yields the alkyl hydrides, (PONOP)Ir(H)(R) (1-(H)(Et)(+) and 1-(H)(Pr)(+)), respectively. Dynamic (1)H NMR characterization of 1-(H)(Et)(+) establishes site exchange between the Ir-H and Ir-CH2 protons (ΔG(exH)(‡)(-110 °C) = 7.2(1) kcal/mol), pointing to a σ-ethane intermediate. By dynamic (13)C NMR spectroscopy, the exchange barrier between the α and ß carbons ("chain-walking") was measured (ΔG(exC)(‡)(-110 °C) = 8.1(1) kcal/mol). The barrier for ethane loss is 17.4(1) kcal/mol (-40 °C), to be compared with the reported barrier to methane loss in 1-(H)(Me)(+) of 22.4 kcal/mol (22 °C). A rhodium σ-ethane complex, (PONOP)Rh(EtH) (2-(EtH)(+)), was prepared by protonation of 2-Et at -150 °C. The barrier for ethane loss (ΔG(dec)(‡)(-132 °C) = 10.9(2) kcal/mol) is lower than for the methane complex, 2-(MeH)(+), (ΔG(dec)(‡)(-87 °C) = 14.5(4) kcal/mol). Full spectroscopic characterization of 2-(EtH)(+) is reported, a key feature of which is the upfield signal at -31.2 ppm for the coordinated CH3 group in the (13)C NMR spectrum. The exchange barrier of the hydrogens of the coordinated methyl group is too low to be measured, but the chain-walking barrier of 7.2(1) kcal/mol (-132 °C) is observable by (13)C NMR. The coordination mode of the alkane ligand and the exchange pathways for the Rh and Ir complexes are evaluated by DFT studies. On the basis of the computational studies, it is proposed that chain-walking occurs by different mechanisms: for Rh, the lowest energy path involves a η(2)-ethane transition state, while for Ir, the lowest energy exchange pathway proceeds through the symmetrical ethylene dihydride complex.

Selective electrocatalytic reduction of CO2 to formate by water-stable iridium dihydride pincer complexes.

Iridium dihydride complexes supported by PCP-type pincer ligands rapidly insert CO(2) to yield κ(2)-formate monohydride products in THF. In acetonitrile/water mixtures, these complexes become efficient and selective catalysts for electrocatalytic reduction of CO(2) to formate. Electrochemical and NMR spectroscopic studies have provided mechanistic details and structures of key intermediates.

Ring-slippage and multielectron redox properties of Fe/Ru/Os-bis(arene) complexes: does hapticity change really cause potential inversion?

Bis(hexamethylbenzene) complexes of the group 8 metals (Fe, Ru, Os) show surprising diversity in their electron-transfer mechanisms and associated thermodynamics for the M(II) → M(I) → M(0) redox series. In electrochemical experiments, the Fe complex exhibits normally ordered potentials separated by ∼1 V, the Ru system shows nearly overlapping one-electron redox events, and Os demonstrates a one-step, two-electron transfer with a peak potential separation suggestive of highly inverted potentials. It has been conjectured that the sequential one-electron transfers observed for Fe are due to the lack of an accessible η(4):η(6) Fe(0) state, destabilizing the fully reduced species. Using an established model chemistry based on DFT, we demonstrate that the hapticity change is a consequence of the bonding throughout this transition metal triad and that apparent multielectron behavior is controlled by the vertical electron attachment component of the M(II) → M(I) redox event. Furthermore, the η(6):η(6) Fe(0) triplet state is more favorable than the hypothetical η(4):η(6) singlet state, emphasizing that the hapticity change is not sufficient for multielectron behavior. Despite both displaying two-electron redox responses, Ru and Os traverse fundamentally different mechanisms based on whether the first (Os) or second (Ru) electron transfer induces the hapticity change. While the electronic structure analysis is limited to the Fe triad here, the conceptual model that we developed provides a general understanding of the redox behavior exhibited by d(6) bis(arene) compounds.

Role of coordination geometry in dictating the barrier to hydride migration in d6 square-pyramidal iridium and rhodium pincer complexes.

Syntheses of the olefin hydride complexes [(POCOP)M(H)(olefin)][BAr(f)(4)] (6a-M, M = Ir or Rh, olefin = C(2)H(4); 6b-M, M = Ir or Rh, olefin = C(3)H(6); POCOP = 2,6-bis(di-tert-butylphosphinito)benzene; BAr(f) = tetrakis(3,5-trifluoromethylphenyl)borate) are reported. A single-crystal X-ray structure determination of 6b-Ir shows a square-pyramidal coordination geometry for Ir, with the hydride ligand occupying the apical position. Dynamic NMR techniques were used to characterize these complexes. The rates of site exchange between the hydride and the olefinic hydrogens yielded ΔG(++) = 15.6 (6a-Ir), 16.8 (6b-Ir), 12.0 (6a-Rh), and 13.7 (6b-Rh) kcal/mol. The NMR exchange data also established that hydride migration in the propylene complexes yields exclusively the primary alkyl intermediate arising from 1,2-insertion. Unexpectedly, no averaging of the top and bottom faces of the square-pyramidal complexes is observed in the NMR spectra at high temperatures, indicating that the barrier for facial equilibration is >20 kcal/mol for both the Ir and Rh complexes. A DFT computational study was used to characterize the free energy surface for the hydride migration reactions. The classical terminal hydride complexes, [M(POCOP)(olefin)H](+), are calculated to be the global minima for both Rh and Ir, in accord with experimental results. In both the Rh ethylene and propylene complexes, the transition state for hydride migration (TS1) to form the agostic species is higher on the energy surface than the transition state for in-place rotation of the coordinated C-H bond (TS2), while for Ir, TS2 is the high point on the energy surface. Therefore, only for the case of the Rh complexes is the NMR exchange rate a direct measure of the hydride migration barrier. The trends in the experimental barriers as a function of M and olefin are in good agreement with the trends in the calculated exchange barriers. The calculated barriers for the hydride migration reaction in the Rh complexes are ∼2 kcal/mol higher than for the Ir complexes, despite the fact that the energy difference between the olefin hydride ground state and the agostic alkyl structure is ∼4 kcal/mol larger for Ir than for Rh. This feature, together with the high barrier for interchange of the top and bottom faces of the complexes, is proposed to arise from the unique coordination geometry of the agostic complexes and the strong preference for a cis-divacant octahedral geometry in four-coordinate intermediates.

Density functional reactivity theory characterizes charge separation propensity in proton-coupled electron transfer reactions.

Proton-coupled electron transfer (PCET) reactions occur in many biological and artificial solar energy conversion processes. In these reactions the electron is often transferred to a site distant to the proton acceptor site. In this work, we employ the dual descriptor and the electrophilic Fukui function from density functional reactivity theory (DFRT) to characterize the propensity for an electron to be transferred to a site other than the proton acceptor site. The electrophilic regions of hydrogen bond or van der Waal reactant complexes were examined using these DFRT descriptors to determine the region of space to which the electron is most likely to be transferred. This analysis shows that in PCET reactions the electrophilic region of the reactant complex does not include the proton acceptor site.

Theory of hydride-proton transfer (HPT) carbonyl reduction by [Os(III)(tpy)(Cl)(NH═CHCH(3))(NSAr)].

Quantum mechanical analysis reveals that carbonyl reduction of aldehydes and ketones by the imine-based reductant cis-[Os(III)(tpy)(Cl)(NH═CHCH(3))(NSAr)] (2), which is accessible by reduction of the analogous nitrile, occurs by hydride-proton transfer (HPT) involving both the imine and sulfilimido ligands. In carbonyl reduction, water or alcohol is necessary to significantly lower the barrier for proton shuttling between ligands. The -N(H)SAr group activates the carbonyl group through hydrogen bonding while the -NC(H)CH(3) ligand delivers the hydride.

Characterization of a rhodium(I) sigma-methane complex in solution.

Numerous transition metal-mediated reactions, including hydrogenations, hydrosilations, and alkane functionalizations, result in the cleavage of strong sigma bonds. Key intermediates in these reactions often involve coordination of the sigma bond of dihydrogen, silanes (Si-H), or alkanes (C-H) to the metal center without full scission of the bond. These sigma complexes have been characterized to varying degrees in solid state and solution. However, a sigma complex of the simplest hydrocarbon, methane, has eluded full solution characterization. Here, we report nuclear magnetic resonance spectra of a rhodium(I) sigma-methane complex obtained by protonation of a rhodium-methyl precursor in CDCl2F solvent at -110 degrees C. The sigma-methane complex is shown to be more stable than the corresponding rhodium(III) methyl hydride complex. Even at -110 degrees C, methane rapidly tumbles in the coordination sphere of rhodium, exchanging free and bound hydrogens. Kinetic studies reveal a half-life of about 83 minutes at -87 degrees C for dissociation of methane (free energy of activation is 14.5 kilocalories per mole).

Molecular acidity: A quantitative conceptual density functional theory description.

Accurate predictions of molecular acidity using ab initio and density functional approaches are still a daunting task. Using electronic and reactivity properties, one can quantitatively estimate pKa values of acids. In a recent paper [S. B. Liu and L. G. Pedersen, J. Phys. Chem. A 113, 3648 (2009)], we employed the molecular electrostatic potential (MEP) on the nucleus and the sum of valence natural atomic orbital (NAO) energies for the purpose. In this work, we reformulate these relationships on the basis of conceptual density functional theory and compare the results with those from the thermodynamic cycle method. We show that MEP and NAO properties of the dissociating proton of an acid should satisfy the same relationships with experimental pKa data. We employ 27 main groups and first to third row transition metal-water complexes as illustrative examples to numerically verify the validity of these strong linear correlations. Results also show that the accuracy of our approach and that of the conventional method through the thermodynamic cycle are statistically similar.

Chemical functionalization of silica and alumina particles for dispersion in carbon dioxide.

The steric stabilization and flocculation of modified silica and alumina particle suspensions in condensed CO(2) were studied. Silica particles (average diameters of 7 and 12 nm) were functionalized using chlorosilanes of the form C(n)F(2n+1)CH(2)CH(2)Si(CH(3))(2)Cl (n = 8, 4, or 1) to give C(n)F(2n+1)-silica. Alumina particles (diameter of 8-14 nm) were grafted with C(8)F(17)CH(2)CH(2)Si(OEt)(3) and chemically modified with perfluorononanoic acid to yield C(8)F(17)-alumina and C(8)F(17)COOH-alumina, respectively. Elemental analysis and thermogravimetric analysis on the derivatized particles were carried out, and surface coverage was calculated. The stabilization of these modified particles in condensed CO(2) was quantified using turbidimetry. Particle stability was found to increase with increasing fluorinated tail length, temperature, and CO(2) density. Unmodified particles and those modified with only -CF(3) tails were unstable in condensed CO(2). Stabilization in supercritical CO(2) is continuous up to 24 h for the C(n)F(2n+1)-silica (n >/= 4) particles and 96 h for the C(8)F(17)-alumina particles. The C(8)F(17)COOH-alumina particles gave a significantly higher graft density than the C(8)F(17)-alumina particles but are not as stable in CO(2). The C(8)F(17)-alumina particles were stable at lower CO(2) densities than the modified silica particles. This stability difference may be attributed to the precursor organosilanes being monofunctional (modified silica) versus trifunctional (modified alumina), producing different structures on the surface.

Oxidative dissolution of copper and zinc metal in carbon dioxide with tert-butyl peracetate and a beta-diketone chelating agent.

A series of beta-diketone ligands, R(1)COCH(2)COR(2) [tmhdH (R(1) = R(2) = C(CH(3))(3)); tfacH (R(1) = CF(3); R(2) = CH(3)); hfacH (R(1) = R(2) = CF(3))], in combination with tert-butyl peracetate (t-BuPA), have been investigated as etchant solutions for dissolution of copper metal into carbon dioxide solvent. Copper removal in CO(2) increases in the order tfacH < tmhdH < hfacH. A study of the reactions of the hfacH/t-BuPA etchant solution with metallic copper and zinc was conducted in three solvents: scCO(2) (supercrical CO(2)); hexanes; CD(2)Cl(2). The etchant solution/metallic zinc reaction produced a diamagnetic Zn(II) complex, which allowed NMR identification of the t-BuPA decomposition products as tert-butyl alcohol and acetic acid. Gravimetric analysis of the amount of zinc consumed, together with NMR studies, confirmed the 1:1:2 Zn:t-BuPA:hfacH reaction stoichiometry, showing t-BuPA to be an overall two-electron oxidant for Zn(0). The metal-containing products of the copper and zinc reactions were characterized by elemental analysis, IR spectroscopy, and, as appropriate, NMR spectroscopy and single-crystal X-ray diffraction [trans-M(hfac)(2)(H(2)O)(CH(3)CO(2)H) (1, M = Cu; 2, M = Zn)]. On the basis of the experimental results, a working model of the oxidative dissolution reaction is proposed, which delineates the key chemical variables in the etching reaction. These t-BuPA/hfacH etchant solutions may find application in a CO(2)-based chemical mechanical planarization (CMP) process.

Etchant solutions for the removal of Cu(0) in a supercritical CO2-based "dry" chemical mechanical planarization process for device fabrication.

The microelectronics industry is focused on increasing chip complexity, improving the density of electron carriers, and decreasing the dimensions of the interconnects into the sub-0.25 mum regime while maintaining high aspect ratios. Water-based chemical mechanical planarization or polishing (CMP) faces several technical and environmental challenges. Condensed CO2 has significant potential for replacing current CMP solvents as a "dry" etching medium because of its unique properties. In working toward a condensed CO2-based CMP process, we have successfully investigated the oxidation and chelation of solid copper metal in liquid and supercritical CO2 using ethyl peroxydicarbonate and a beta-diketone chelating agent.

Density functional theory study of redox pairs: 2. Influence of solvation and ion-pair formation on the redox behavior of cyclooctatetraene and nitrobenzene.

A study of the electrochemical behavior of cyclooctatetraene (COT) and nitrobenzene with Density Functional Theory and the conductor like solvation model (COSMO) is reported. The two-electron reduction of the tub-shaped COT molecule is accompanied by a structural change to a planar structure of D(4)(h)() symmetry in the first electron addition step, and to a fully aromatic structure of D(8)(h)() symmetry in the second electron addition step. Theoretical models are examined that are aimed at understanding the electrolyte- and solvent-dependent redox behavior of COT, in which a single 2e(-) redox wave is observed with KI electrolyte in liquid ammonia solution (DeltaDeltaE(disp) = [E(-2) - E(-1)] - [E(-1) - E(0)] < 0, inverted potential), while two 1e(-) redox waves are observed (DeltaDeltaE(disp) > 0) with NR(4)(+)X(-) (R = butyl, propyl; X(-) = perchlorate) electrolyte in dimethylformamide solution. In all cases, the computed reaction energy profiles are in fair agreement with the experimental reduction potentials. A chemically intuitive theoretical square scheme method of energy partitioning is introduced to analyze in detail the effects of structural changes and ion-pair formation on the relative energies of the redox species. The structural relaxation energy for conversion of tub-COT to planar-COT is mainly apportioned to the first reduction step, and is therefore a positive contribution to DeltaDeltaE(disp). The effect of the structural change on the disproportionation energy for COT is counteracted by the substantially more positive reduction potential for planar-(COT)(-1) in comparison to tub-(COT)(-1). Ion pairing of alkali metal counterions with the anionic reduction products gives rise to a negative contribution to DeltaDeltaE(disp) because the second ion-pairing step is more exothermic than the first, and the reduction of [KA] (A = COT, NB) is more exothermic than the reduction of A(-1). For COT, this negative energy differential term as a result of ion pairing predicts the experimentally observed inversion in the two 1e(-) potentials (DeltaDeltaE(disp) < 0). Nitrobenzene is treated with the same computational protocol to provide a system for comparison that is not complicated by the major structural change that influences the COT energy profile.