M-Ge-Si thermolytic molecular precursors and models for germanium-doped transition metal sites on silica.

The synthesis, thermolysis, and surface organometallic chemistry of thermolytic molecular precursors based on a new germanosilicate ligand platform, -OGe[OSi(OtBu)3]3, is described. Use of this ligand is demonstrated with preparation of complexes containing the first-row transition metals Cr, Mn, and Fe. The thermolysis and grafting behavior of the synthesized complexes, Fe{OGe[OSi(OtBu)3]3}2 (FeGe), Mn{OGe[OSi(OtBu)3]3}2(THF)2 (MnGe) and Cr{OGe[OSi(OtBu)3]3}2(THF)2 (CrGe), was evaluated using a combination of thermogravimetric analysis; nuclear magnetic resonance (NMR), ultraviolet-visible (UV-Vis), and electron paramagnetic resonance (EPR) spectroscopies; and single-crystal X-ray diffraction (XRD). Grafting of the precursors onto SBA-15 mesoporous silica and subsequent calcination in air led to substantial changes in transition metal coordination environments and oxidation states, the implications of which are discussed in the context of low-coordinate and low oxidation state thermolytic molecular precursors.

Elucidating the active phases of CoOx films on Au(111) in the CO oxidation reaction.

Noble metals supported on reducible oxides, like CoOx and TiOx, exhibit superior activity in many chemical reactions, but the origin of the increased activity is not well understood. To answer this question we studied thin films of CoOx supported on an Au(111) single crystal surface as a model for the CO oxidation reaction. We show that three reaction regimes exist in response to chemical and topographic restructuring of the CoOx catalyst as a function of reactant gas phase CO/O2 stoichiometry and temperature. Under oxygen-lean conditions and moderate temperatures (≤150 °C), partially oxidized films (CoOx<1) containing Co0 were found to be efficient catalysts. In contrast, stoichiometric CoO films containing only Co2+ form carbonates in the presence of CO that poison the reaction below 300 °C. Under oxygen-rich conditions a more oxidized catalyst phase (CoOx>1) forms containing Co3+ species that are effective in a wide temperature range. Resonant photoemission spectroscopy (ResPES) revealed the unique role of Co3+ sites in catalyzing the CO oxidation. Density function theory (DFT) calculations provided deeper insights into the pathway and free energy barriers for the reactions on these oxide phases. These findings in this work highlight the versatility of catalysts and their evolution to form different active phases, both topological and chemically, in response to reaction conditions exposing a new paradigm in the catalyst structure during operation.

Theory and Simulation of Metal-Insulator-Semiconductor (MIS) Photoelectrodes.

A metal-insulator-semiconductor (MIS) structure is an attractive photoelectrode-catalyst architecture for promoting photoelectrochemical reactions, such as the formation of H2 by proton reduction. The metal catalyzes the generation of H2 using electrons generated by photon absorption and charge separation in the semiconductor. The insulator layer between the metal and the semiconductor protects the latter element from photo-corrosion and, also, significantly impacts the photovoltage at the metal surface. Understanding how the insulator layer determines the photovoltage and what properties lead to high photovoltages is critical to the development of MIS structures for solar-to-chemical energy conversion. Herein, we present a continuum model for charge-carrier transport from the semiconductor to the metal with an emphasis on mechanisms of charge transport across the insulator. The polarization curves and photovoltages predicted by this model for a Pt/HfO2/p-Si MIS structure at different HfO2 thicknesses agree well with experimentally measured data. The simulations reveal how insulator properties (i.e., thickness and band structure) affect band bending near the semiconductor/insulator interface and how tuning them can lead to operation closer to the maximally attainable photovoltage, the flat-band potential. This phenomenon is understood by considering the change in tunneling resistance with insulator properties. The model shows that the best MIS performance is attained with highly symmetric semiconductor/insulator band offsets (e.g., BeO, MgO, SiO2, HfO2, or ZrO2 deposited on Si) and a low to moderate insulator thickness (e.g., between 0.8 and 1.5 nm). Beyond 1.5 nm, the density of filled interfacial trap sites is high and significantly limits the photovoltage and the solar-to-chemical conversion rate. These conclusions are true for photocathodes and photoanodes. This understanding provides critical insight into the phenomena enhancing and limiting photoelectrode performance and how this phenomenon is influenced by insulator properties. The study gives guidance toward the development of next-generation insulators for MIS structures that achieve high performance.

Ethene Hydroformylation Catalyzed by Rhodium Dispersed with Zinc or Cobalt in Silanol Nests of Dealuminated Zeolite Beta.

Catalysts for hydroformylation of ethene were prepared by grafting Rh into nests of ≡SiOZn-OH or ≡SiOCo-OH species prepared in dealuminated BEA zeolite. X-ray absorption spectra and infrared spectra of adsorbed CO were used to characterize the dispersion of Rh. The Rh dispersion was found to increase markedly with increasing M/Rh (M = Zn or Co) ratio; further increases in Rh dispersion occurred upon use for ethene hydroformylation catalysis. The turnover frequency for ethene hydroformylation measured for a fixed set of reaction conditions increased with the fraction of atomically dispersed Rh. The ethene hydroformylation activity is 15.5-fold higher for M = Co than for M = Zn, whereas the propanal selectivity is slightly greater for the latter catalyst. The activity of the Co-containing catalyst exceeds that of all previously reported Rh-containing bimetallic catalysts. The rates of ethene hydroformylation and ethene hydrogenation exhibit positive reaction orders in ethene and hydrogen but negative orders in carbon monoxide. In situ IR spectroscopy and the kinetics of the catalytic reactions suggest that ethene hydroformylation is mainly catalyzed by atomically dispersed Rh that is influenced by Rh-M interactions, whereas ethene hydrogenation is mainly catalyzed by Rh nanoclusters. In situ IR spectroscopy also indicates that the ethene hydroformylation is rate limited by formation of propionyl groups and by their hydrogenation, a conclusion supported by the measured H/D kinetic isotope effect. This study presents a novel method for creating highly active Rh-containing bimetallic sites for ethene hydroformylation and provides new insights into the mechanism and kinetics of this process.

Catalytic deconstruction of waste polyethylene with ethylene to form propylene.

The conversion of polyolefins to monomers would create a valuable carbon feedstock from the largest fraction of waste plastic. However, breakdown of the main chains in these polymers requires the cleavage of carbon-carbon bonds that tend to resist selective chemical transformations. Here, we report the production of propylene by partial dehydrogenation of polyethylene and tandem isomerizing ethenolysis of the desaturated chain. Dehydrogenation of high-density polyethylene with either an iridium-pincer complex or platinum/zinc supported on silica as catalysts yielded dehydrogenated material containing up to 3.2% internal olefins; the combination of a second-generation Hoveyda-Grubbs metathesis catalyst and [PdP(tBu)3(µ-Br)]2 as an isomerization catalyst selectively degraded this unsaturated polymer to propylene in yields exceeding 80%. These results show promise for the application of mild catalysis to deconstruct otherwise stable polyolefins.

Continuum Modeling of Porous Electrodes for Electrochemical Synthesis.

Electrochemical synthesis possesses substantial promise to utilize renewable energy sources to power the conversion of abundant feedstocks to value-added commodity chemicals and fuels. Of the potential system architectures for these processes, only systems employing 3-D structured porous electrodes have the capacity to achieve the high rates of conversion necessary for industrial scale. However, the phenomena and environments in these systems are not well understood and are challenging to probe experimentally. Fortunately, continuum modeling is well-suited to rationalize the observed behavior in electrochemical synthesis, as well as to ultimately provide recommendations for guiding the design of next-generation devices and components. In this review, we begin by presenting an historical review of modeling of porous electrode systems, with the aim of showing how past knowledge of macroscale modeling can contribute to the rising challenge of electrochemical synthesis. We then present a detailed overview of the governing physics and assumptions required to simulate porous electrode systems for electrochemical synthesis. Leveraging the developed understanding of porous-electrode theory, we survey and discuss the present literature reports on simulating multiscale phenomena in porous electrodes in order to demonstrate their relevance to understanding and improving the performance of devices for electrochemical synthesis. Lastly, we provide our perspectives regarding future directions in the development of models that can most accurately describe and predict the performance of such devices and discuss the best potential applications of future models.

Assessing the stability of Pd-exchanged sites in zeolites with the aid of a high throughput quantum chemistry workflow.

Cation exchanged-zeolites are functional materials with a wide range of applications from catalysis to sorbents. They present a challenge for computational studies using density functional theory due to the numerous possible active sites. From Al configuration, to placement of extra framework cation(s), to potentially different oxidation states of the cation, accounting for all these possibilities is not trivial. To make the number of calculations more tractable, most studies focus on a few active sites. We attempt to go beyond these limitations by implementing a workflow for a high throughput screening, designed to systematize the problem and exhaustively search for feasible active sites. We use Pd-exchanged CHA and BEA to illustrate the approach. After conducting thousands of explicit DFT calculations, we identify the sites most favorable for the Pd cation and discuss the results in detail. The high throughput screening identifies many energetically favorable sites that are non-trivial. Lastly, we employ these results to examine NO adsorption in Pd-exchanged CHA, which is a promising passive NOx adsorbent (PNA) during the cold start of automobiles. The results shed light on critical active sites for NOx capture that were not previously studied.

Engineering Catalyst-Electrolyte Microenvironments to Optimize the Activity and Selectivity for the Electrochemical Reduction of CO2 on Cu and Ag.

The electrochemical reduction of carbon dioxide (CO2R) driven by renewably generated electricity (e.g., solar and wind) offers a promising means for reusing the CO2 released during the production of cement, steel, and aluminum as well as the production of ammonia and methanol. If CO2 could be removed from the atmosphere at acceptable costs (i.e., <$100/t of CO2), then CO2R could be used to produce carbon-containing chemicals and fuels in a fully sustainable manner. Economic considerations dictate that CO2R current densities must be in the range of 0.1 to 1 A/cm2 and selectivity toward the targeted product must be high in order to minimize separation costs. Industrially relevant operating conditions can be achieved by using gas diffusion electrodes (GDEs) to maximize the transport of species to and from the cathode and combining such electrodes with a solid-electrolyte membrane by eliminating the ohmic losses associated with liquid electrolytes. Additionally, high product selectivity can be attained by careful tuning of the microenvironment near the catalyst surface (e.g., the pH, the concentrations of CO2 and H2O, and the identities of the cations in the double layer adjacent to the catalyst surface).We begin this Account with a discussion of our experimental and theoretical work aimed at optimizing catalyst microenvironments for CO2R. We first examine the effects of catalyst morphology on the production of multicarbon (C2+) products over Cu-based catalysts and then explore the role of mass transfer combined with the kinetics of buffer reactions in the local concentration of CO2 and pH at the catalyst surface. This is followed by a discussion of the dependence of the local CO2 concentration and pH on the dynamics of CO2R and the formation of specific products over both Cu and Ag catalysts. Next, we explore the impact of electrolyte cation identity on the rate of CO2R and the distribution of products. Subsequently, we look at utilizing pulsed electrolysis to tune the local pH and CO2 concentration at the catalyst surface. The last part of the discussion demonstrates that ionomer-coated catalysts in combination with pulsed electrolysis can enable the attainment of very high (>90%) selectivity to C2+ products over Cu in an aqueous electrolyte. This part of the Account is then extended to consider the difference in the catalyst-nanoparticle microenvironment, present in the catalyst layer of a membrane electrode assembly (MEA), with respect to that of a planar electrode immersed in an aqueous electrolyte.

Propane Dehydrogenation Catalyzed by Isolated Pt Atoms in ≡SiOZn-OH Nests in Dealuminated Zeolite Beta.

Atomically dispersed noble metal catalysts have drawn wide attention as candidates to replace supported metal clusters and metal nanoparticles. Atomic dispersion can offer unique chemical properties as well as maximum utilization of the expensive metals. Addition of a second metal has been found to help reduce the size of Pt ensembles in bimetallic clusters; however, the stabilization of isolated Pt atoms in small nests of nonprecious metal atoms remains challenging. We now report a novel strategy for the design, synthesis, and characterization of a zeolite-supported propane dehydrogenation catalyst that incorporates predominantly isolated Pt atoms stably bonded within nests of Zn atoms located within the nanoscale pores of dealuminated zeolite Beta. The catalyst is stable in long-term operation and exhibits high activity and high selectivity to propene. Atomic resolution images, bolstered by X-ray absorption spectra, demonstrate predominantly atomic dispersion of the Pt in the nests and, with complementary infrared and nuclear magnetic resonance spectra, determine a structural model of the nested Pt.

Critical Role of Thermal Fluctuations for CO Binding on Electrocatalytic Metal Surfaces.

This work considers the evaluation of density functional theory (DFT) when comparing against experimental observations of CO binding trends on the strong binding Pt(111) and intermediate binding Cu(111) and for weak binding Ag(111) and Au(111) surfaces important in electrocatalysis. By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we find that the RPBE and B97M-rV DFT functionals yield qualitatively better metal surface strain trends and CO enthalpies of binding for Cu(111) and Pt(111) than found at 0 K, thereby correcting the overbinding by 0.2 to 0.3 eV to yield better agreement with the enthalpies determined from experiment. The importance of dispersion effects are manifest for the weak CO binding Ag(111) and Au(111) surfaces at finite temperatures in which the RPBE functional does not bind CO at all, while the B97M-rV functional shows that the CO-metal interactions are a mixture of chemisorbed and physisorbed species with binding enthalpies that are within ∼0.05 eV of experiment. Across all M(111) surfaces, we show that the B97M-rV functional consistently predicts the correct atop site preference for all metals due to thermally induced surface distortions that preferentially favor the undercoordinated site. This study demonstrates the need to fully account for finite temperature fluctuations to make contact with the binding enthalpies from surface science experiments and electrocatalysis applications.

Few-Unit-Cell MFI Zeolite Synthesized using a Simple Di-quaternary Ammonium Structure-Directing Agent.

Synthesis of a pentasil-type zeolite with ultra-small few-unit-cell crystalline domains, which we call FDP (few-unit-cell crystalline domain pentasil), is reported. FDP is made using bis-1,5(tributyl ammonium) pentamethylene cations as structure directing agent (SDA). This di-quaternary ammonium SDA combines butyl ammonium, in place of the one commonly used for MFI synthesis, propyl ammonium, and a five-carbon nitrogen-connecting chain, in place of the six-carbon connecting chain SDAs that are known to fit well within the MFI pores. X-ray diffraction analysis and electron microscopy imaging of FDP indicate ca. 10 nm crystalline domains organized in hierarchical micro-/meso-porous aggregates exhibiting mesoscopic order with an aggregate particle size up to ca. 5 µm. Al and Sn can be incorporated into the FDP zeolite framework to produce active and selective methanol-to-hydrocarbon and glucose isomerization catalysts, respectively.

Mechanism and Kinetics of Acetone Conversion to Isobutene over Isolated Hf Sites Grafted to Silicalite-1 and SiO2.

Isolated hafnium (Hf) sites were prepared on Silicalite-1 and SiO2 and investigated for acetone conversion to isobutene. Characterization by IR, 1H MAS NMR, and UV-vis spectroscopy suggests that Hf atoms are bonded to the support via three O atoms and have one hydroxyl group, i.e, (≡SiO)3Hf-OH. In the case of Hf/Silicalite-1, Hf-OH groups hydrogen bond with adjacent Si-OH to form (≡SiO)3Hf-OH···HO-Si≡ complexes. The turnover frequency for isobutene formation from acetone is 4.5 times faster over Hf/Silicalite-1 than Hf/SiO2. Lewis acidic Hf sites promote the aldol condensation of acetone to produce mesityl oxide (MO), which is the precursor to isobutene. For Hf/SiO2, both Hf sites and Si-OH groups are responsible for the decomposition of MO to isobutene and acetic acid, whereas for Hf/Silicalite-1, the (≡SiO)3Hf-OH···HO-Si≡ complex is the active site. Measured reaction kinetics show that the rate of isobutene formation over Hf/SiO2 and Hf/Silicalite-1 is nearly second order in acetone partial pressure, suggesting that the rate-limiting step involves formation of the C-C bond between two acetone molecules. The rate expression for isobutene formation predicts a second order dependence in acetone partial pressure at low partial pressures and a decrease in order with increasing acetone partial pressure, in good agreement with experimental observation. The apparent activation energy for isobutene formation from acetone over Hf/SiO2 is 116.3 kJ/mol, while that for Hf/Silicalite-1 is 79.5 kJ/mol. The lower activation energy for Hf/Silicalite-1 is attributed to enhanced adsorption of acetone and formation of a C-C bond favored by the H-bonding interaction between Hf-OH and an adjacent Si-OH group.

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals.

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO2 or CO reduction. Here, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO-water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

Computational Modeling Predicts the Stability of Both Pd+ and Pd2+ Ion-Exchanged into H-CHA.

Passive NOx adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NOx adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework, before the introduction of NOx. Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next-next-nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd+ instead of Pd2+. The stability of Pd+ in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-VIS, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

Siloxyaluminate and Siloxygallate Complexes as Models for Framework and Partially Hydrolyzed Framework Sites in Zeolites and Zeotypes.

Anionic molecular models for nonhydrolyzed and partially hydrolyzed aluminum and gallium framework sites on silica, M[OSi(OtBu)3 ]4 - and HOM[OSi(OtBu)3 ]3 - (where M=Al or Ga), were synthesized from anionic chlorides Li{M[OSi(OtBu)3 ]3 Cl} in salt metathesis reactions. Sequestration of lithium cations with [12]crown-4 afforded charge-separated ion pairs composed of monomeric anions M[OSi(OtBu)3 ]4 - with outer-sphere [([12]crown-4)2 Li]+ cations, and hydroxides {HOM[OSi(OtBu)3 ]3 } with pendant [([12]crown-4)Li]+ cations. These molecular models were characterized by single-crystal X-ray diffraction, vibrational spectroscopy, mass spectrometry and NMR spectroscopy. Upon treatment of monomeric [([12]crown-4)Li]{HOM[OSi(OtBu)3 ]3 } complexes with benzyl alcohol, benzyloxide complexes were formed, modeling a possible pathway for the formation of active sites for Meerwin-Ponndorf-Verley (MPV) transfer hydrogenations with Al/Ga-doped silica catalysts.

Understanding Multi-Ion Transport Mechanisms in Bipolar Membranes.

Bipolar membranes (BPMs) have the potential to become critical components in electrochemical devices for a variety of electrolysis and electrosynthesis applications. Because they can operate under large pH gradients, BPMs enable favorable environments for electrocatalysis at the individual electrodes. Critical to the implementation of BPMs in these devices is understanding the kinetics of water dissociation that occurs within the BPM as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-electrolyzers. In this study, a continuum model of multi-ion transport in a BPM is developed and fit to experimental data. Specifically, concentration profiles are determined for all ionic species, and the importance of a water-dissociation catalyst is demonstrated. The model describes internal concentration polarization and co- and counter-ion crossover in BPMs, determining the mode of transport for ions within the BPM and revealing the significance of salt-ion crossover when operated with pH gradients relevant to electrolysis and electrosynthesis. Finally, a sensitivity analysis reveals that the performance and lifetime of BPMs can be improved substantially by using of thinner dissociation catalysts, managing water transport, modulating the thickness of the individual layers in the BPM to control salt-ion crossover, and increasing the ion-exchange capacity of the ion-exchange layers in order to amplify the water-dissociation kinetics at the interface.

Facing the Challenges of Borderline Oxidation State Assignments Using State-of-the-Art Computational Methods.

The oxidation state (OS) of metals and ligands in inorganic complexes may be defined by carefully curated rules, such as from IUPAC, or by computational procedures such as the effective oxidation state (EOS) or localized orbital bonding analysis (LOBA). Such definitions typically agree for systems with simple ionic bonding and innocent ligands but may disagree as the boundary between ionic and covalent bonds is approached, or as the role of ligand noninnocence becomes nontrivial, or high oxidation states of metals are supported by heavy dative bonding, and so on. This work systematically compares IUPAC, EOS, and LOBA across a series of complexes where OS assignment is challenging. These systems include high-valent transition metal oxides, transition metal complexes with noninnocent ligands such as dithiolate and nitrosyl, metal sulfur dioxide adducts, and two transition metal carbene complexes. The differences in OS assignment by the three methods are carefully discussed, demonstrating the synergy between EOS and LOBA. In addition, a clarity index for LOBA OS assignments is introduced that provides an indication of whether or not its predictions are close to the ionic-covalent boundary.

Ethanol Conversion to Butadiene over Isolated Zinc and Yttrium Sites Grafted onto Dealuminated Beta Zeolite.

Zinc and Yttrium single sites were introduced into the silanol nests of dealuminated BEA zeolite to produce Zn-DeAlBEA and Y-DeAlBEA. These materials were then investigated for the conversion of ethanol to 1,3-butadiene. Zn-DeAlBEA was found to be highly active for ethanol dehydrogenation to acetaldehyde and exhibited low activity for 1,3-butadiene generation. By contrast, Y-DeAlBEA was highly active for 1,3-butadiene formation but exhibited no activity for ethanol dehydrogenation. The formation of 1,3-butadine over Y-DeAlBEA and Zn-DeAlBEA does not occur via aldol condensation of acetaldehyde but, rather, by concerted reaction of coadsorbed acetaldehyde and ethanol. The active centers for this process are ≡Si-O-Y(OH)-O-Si≡ or ≡Si-O-Zn-O-Si-O≡ groups closely associated with adjacent silanol groups. The active sites in Y-DeAlBEA are 70 times more active than the Y sites supported on silica, for which the Y site is similar to that in Y-SiO2 but which lacks adjacent hydroxyl groups, and are 7 times more active than the active sites in Zn-DeAlBEA. We propose that C-C bond coupling in Y-DeAlBEA proceeds via the reaction of coadsorbed acetaldehyde and ethanol to form crotyl alcohol and water. The dehydration of crotyl alcohol to 1,3-butadiene is facile and occurs over the mildly Brønsted acidic ≡Si-OH groups present in the silanol nest of DeAlBEA. The catalysts reported here are notably more active than those previously reported for both the direct conversion of ethanol to 1,3-butadiene or the formation of this product by the reaction of ethanol and acetaldehyde.