Direct Observation of Micelle Fragmentation via In Situ Liquid-Phase Transmission Electron Microscopy
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.
The production of acrylic acid by coupling CO 2 and ethylene is highly desirable, as it uses a C 1 molecule to introduce the carboxylic group into a large commodity chemical, forming an unsaturated carboxylic acid. Molecular complexes have been used as catalysts in the literature for this reaction; however, these homogeneous studies were all limited by low activity and stability, and efficient heterogeneous catalysts remain to be explored. Here, we report density functional theory calculations, through which we investigated metal-exchanged MFI zeolites with metals from 3d, 4d, and 5d series for this reaction. Screening of reaction thermodynamics has shown that Sc, Y, Cr, Mo, W, Mn, Tc, and Re dispersed into MFI zeolites are potential catalysts. We further calculated the activation barriers for the β–H transfer, a key elementary step that normally determines the overall rates, and found the lowest barrier over Sc- and Y-MFI; such a feasible kinetics is attributed to the low electronegativity of Sc and Y that strongly affects the nature of the transition states and the charge redistribution during the β–H transfer. In conclusion, this work thus provides guidelines for designing a zeolite-based catalyst for this valuable reaction and identifies promising metal centers to be used.
Electrostatic interactions play a pivotal role in enzymatic catalysis and are increasingly modeled explicitly in computational enzyme design; nevertheless, they are challenging to measure experimentally. Using vibrational Stark effect (VSE) spectroscopy, we have measured electric fields inside the active site of the enzyme ketosteroid isomerase (KSI). These studies have shown that these fields can be unusually large, but it has been unclear to what extent they specifically stabilize the transition state (TS) relative to a ground state (GS). In the following, we use crystallography and computational modeling to show that KSI’s intrinsic electric field is nearly perfectly oriented to stabilize the geometry of its reaction’s TS. Moreover, we find that this electric field adjusts the orientation of its substrate in the ground state so that the substrate needs to only undergo minimal structural changes upon activation to its TS. This work provides evidence that the active site electric field in KSI is preorganized to facilitate catalysis and provides a template for how electrostatic preorganization can be measured in enzymatic systems.
The trityl-substituted bisoxazoline ( TrH BOX) was prepared as a chiral analogue to a previously reported nickel dipyrrin system capable of ring-closing amination catalysis. Ligand metalation with divalent NiI 2 (py) 4 followed by potassium graphite reduction afforded the monovalent (TrH BOX)Ni(py) (4). Slow addition of 1.4 equiv of a benzene solution of 1-adamantylazide to 4 generated the tetrazido ( TrH BOX)Ni(κ 2 -N 4 Ad 2 ) (5) and terminal iminyl adduct ( TrH BOX)Ni(NAd) (6). Investigation of 6 via single-crystal X-ray crystallography, NMR and EPR spectroscopies, and computations revealed a Ni(II)-iminyl radical formulation, similar to its dipyrrinato congener. Complex 4 exhibits enantioselective intramolecular C–H bond amination to afford N-heterocyclic products from 4-aryl-2-methyl-2-azidopentanes. Catalytic C–H amination occurs under mild conditions (5 mol % catalyst, 60 °C) and provides pyrrolidine products in decent yield (29%–87%) with moderate ee (up to 73%). Substrates with a 3,5-dialkyl substitution on the 4-aryl position maximized the observed enantioselectivity. Kinetic studies to probe the reaction mechanism were conducted using 1H and 19F NMR spectroscopies. A small, intermolecular kinetic isotope effect (1.35 ± 0.03) suggests an H-atom abstraction step with an asymmetric transition state while the reaction rate is measured to be first order in catalyst and zeroth order in substrate concentrations. Enantiospecific deuterium labeling studies show that the enantioselectivity is dictated by both the H-atom abstraction and radical recombination steps due to the comparable rate between radical rotation and C–N bond formation. Furthermore, the competing elements of the two-step reaction where H-removal from the pro-R configuration is preferred while the preferential radical capture occurs with the Si face of the carboradical likely lead to the diminished ee observed, as corroborated by theoretical calculations. Furthermore, based on these enantio-determining steps, catalytic enantioselective synthesis of 2,5-bis-tertiary pyrrolidines is demonstrated with good yield (50–78%) and moderate ee (up to 79%).
Ethylene-forming enzyme (EFE) is an ambifunctional iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase. In its major (EF) reaction, it converts carbons 1, 2, and 5 of 2OG to CO 2 and carbons 3 and 4 to ethylene, a four-electron oxidation drastically different from the simpler decarboxylation of 2OG to succinate mediated by all other Fe/2OG enzymes. EFE also catalyzes a minor reaction, in which the normal decarboxylation is coupled to oxidation of L-arginine (a required activator for the EF pathway), resulting in its conversion to L-glutamate semialdehyde and guanidine. Here we show that, consistent with precedent, the L-Arg-oxidation (RO) pathway proceeds via an iron(IV)-oxo (ferryl) intermediate. Use of 5,5-[ 2 H 2 ]-L-Arg slows decay of the ferryl complex by >16-fold, implying that RO is initiated by hydrogen-atom transfer (HAT) from C5. That this large substrate deuterium kinetic isotope effect has no impact on the EF:RO partition ratio implies that the same ferryl intermediate cannot be on the EF pathway; the pathways must diverge earlier. Consistent with this conclusion, the variant enzyme bearing the Asp191Glu ligand substitution accumulates ~ four times as much of the ferryl complex as the wild-type enzyme and exhibits a ~ 40-fold diminished EF:RO partition ratio. Here, the selective detriment of this nearly conservative substitution to the EF pathway implies that it has unusually stringent stereoelectronic requirements. An active-site, like-charge guanidium pair, which involves the L-Arg substrate/activator and is unique to EFE among four crystallographically characterized L-Arg modifying Fe/2OG oxygenases, may serve to selectively stabilize the transition state leading to the unique EF branch.
The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.
Ga-chabazite zeolites (Ga-CHA) have been found to efficiently catalyze propane dehydrogenation with high propylene selectivity (96%). In situ FTIR spectroscopy and pulse titrations are employed to determine that upon reduction, surface Ga 2 O 3 is reduced and diffuses into the zeolite pores, displacing the Brønsted acid sites (BAS) and forming extra-framework Ga + sites. This isolated Ga + site reacts reversibly with H 2 to form GaHx (2034 cm -1 ) with an enthalpy of formation of ~ -51.2 kJ·mol -1 , a result supported by Density functional theory (DFT) calculations. The initial C 3 H 6 dehydrogenation rates decrease rapidly (40%) during the first 100 min and then decline slowly afterward, while the C 3 H 6 selectivity is stable at ~ 96%. The reduction in the reaction rate is correlated with the formation of polycyclic aromatics inside the zeolite (using UV-vis spectroscopy) indicating that the accumulation of polycyclic aromatics is the main cause of the deactivation. The carbon species formed can be easily oxidized at 600 °C with complete recovery of the PDH catalytic properties. The correlations between GaH x vs. Ga/Al ratio, and PDH rates vs. Ga/Al ratio show that extra-framework Ga + is the active center catalyzing propane dehydrogenation. The higher reaction rate on Ga + than In + in CHA zeolites, by a factor of 43, is the result of differences in the stabilization of the transition state due to the higher stability of Ga 3+ vs. In 3+ . The uniformity of the Ga + sites in this material makes it an excellent model for the molecular understanding of metal cation exchanged hydrocarbon interactions in zeolites.
We recently reported a reaction sequence that activates C–H bonds in simple arenes as well as the N–N triple bond in N 2 , delivering the aryl group to N 2 to form a new N–C bond ( Nature 2020 , 584, 221). This enables the transformation of abundant feedstocks (arenes and N 2 ) into N-containing organic compounds. The key N–C bond forming step occurs upon partial silylation of N 2 . However, the pathway through which reduction, silylation, and migration occurred was unknown. Here, in this work, we describe synthetic, structural, magnetic, spectroscopic, kinetic, and computational studies that elucidate the steps of this transformation. N 2 must be silylated twice at the distal N atom before aryl migration can occur, and sequential silyl radical and silyl cation addition is a kinetically competent pathway to a formally iron(IV)–NN(SiMe 3 ) 2 intermediate that can be isolated at low temperature. Kinetic studies show its first-order conversion to the migrated product, and DFT calculations indicate a concerted transition state for migration. The electronic structure of the formally iron(IV) intermediate is examined using DFT and CASSCF calculations, which reveal contributions from iron(II) and iron(III) resonance forms with oxidized NNSi 2 ligands. The depletion of electron density from the Fe-coordinated N atom makes it electrophilic enough to accept the incoming aryl group. This new pathway for the N–C bond formation offers a method for functionalizing N 2 using organometallic chemistry.
Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.
Second sphere coordination effects ubiquitous in enzymatic catalysis occur through direct interactions, either covalent or non-covalent, with reaction intermediates and transition states. Here, we present herein evidence of indirect second sphere coordination effects in which ligation of water/alkanols far removed from the primary coordination sphere of the active site nevertheless alter energetic landscapes within catalytic redox cycles in the absence of direct physicochemical interactions with surface species mediating catalytic turnovers. Density functional theory, in situ X-ray absorption and infrared spectroscopy, and a wide array of steady-state and transient CO oxidation rate data suggest that the presence of peripheral water renders oxidation half-cycles within two-electron redox cycles over μ m -oxo-bridged trimers in MIL-100(M) more kinetically demanding. Communication between ligated water and the active site appears to occur through the Fe-O-Fe backbone, as inferred from spin density variations on the central μ m -oxygen 'junction'. Evidence is provided for the generality of these second sphere effects in that they influence different types of redox half-cycles or metals, and can be amplified or attenuated through choice of coordinating ligand. Specifically in the case of MIL-100(M) materials, the Cr isostructure can be made to kinetically mimic the Fe variant by disproportionately hindering oxidation half-cycles relative to the reduction half-cycles. Kinetic and spectroscopic inferences presented here significantly expand both the conceptual definition of second sphere effects as well as the palette of synthetic levers available for tuning catalytic redox performance through chemical ligation.
Single atoms of uranium supported on molybdenum sulfide surfaces (U@MoS 2 ) have been recently demonstrated to facilitate the hydrogen evolution reaction (HER) through electrocatalysis. Theoretical calculations have predicted uranium hydroxide moieties bound to edge-sulfur atoms of MoS 2 as a proposed transition state involved in the HER process. However, the isolation of relevant intermediates involved in this process remains a challenge, rendering mechanistic hypotheses unverified. The present work describes the isolation and characterization of a uranium-hydroxide intermediate on molybdenum sulfide surfaces using [(Cp* 3 Mo 3 S 4 )UCp*], a molecular model of a reduced uranium center supported at MoS 2 . Mechanistic investigations highlight the metalloligand cooperativity with uranium involved in the water activation pathway. The corresponding uranium-oxo analogue, [(Cp* 3 Mo 3 S 4 )Cp*U(═O)], was also accessed from the hydroxide cluster via hydrogen atom transfer and from [(Cp* 3 Mo 3 S 4 )UCp*] through an alternative direct oxygen atom transfer. These results provide an atomistic perspective on the reactivity of low-valent uranium at molybdenum sulfide surfaces toward water, modeling key intermediates associated with the HER of U@MoS 2 catalysts.
Here we report a series of isomeric, dicationic Re(bpy 2+ )(CO) 3 I com-plexes with bpy (2,2'-bipyridine) modified by two phenyl-CH 2 -(NMe 3 )+ pendants with cations located at variable distances from the active site for electrocatalytic CO 2 reduction in CH 3 CN/2.8 M H 2 O. The position of the cationic groups dramatically increases the rate of catalysis by 800-fold, from 1.2 to 950 s -1 , with a minor increase in overpotential. Acceleration is due to stabilization of the initial CO 2 adduct and lowering of ΔG ‡ for C-OH bond cleavage by Coulombic stabilization of anionic charges. Performance may be enhanced by accumulation in the electrochemical double layer. Transition state stabilization in the optimized isomer unlocks the low overpotential "protonation-first" pathway, highlighting the sizable effects of subtle structural optimization.
The significant increase in CO 2 emissions from heavy fossil fuel utilization has raised serious concerns, highlighting the need for effective methods to convert CO 2 into value-added chemicals. Here, in this work, we report a computational investigation on the catalytic activity of ZrO 2 -on-Cu inverse catalysts for CO 2 hydrogenation to methanol, considering highly dispersed ZrO 2 trimers on Cu (111). Such clusters present a large ensemble of formate-containing configurations, Zr 3 O n (OH) m (OCHO) l , making the evaluation of the catalytic activity very challenging. We found that the sites on the various catalyst configurations exhibit markedly different activities for formate hydrogenation, despite their similar free energy and composition. To understand these differences in reactivity, we examined the structural and electronic nature of the low free-energy catalyst configurations and identified that the energy of the lowest unoccupied orbital of the reacting formate, modified by its binding with the catalytic site, is a descriptor for the reaction energy of the formate hydrogenation step. From there, we screened an ensemble of catalyst structures using this descriptor to predict highly active metastable catalyst configurations and computed the reaction pathways and transition states for formate hydrogenation. From this investigation, we distinguished reactive from nonreactive sites and formate species on the ZrO 2 /Cu inverse catalyst based on structural and electronic features. We showed that rare metastable configurations control the activity. Additionally, an efficient method for examining the reactivity of a large number of coexisting catalyst structures was developed.
Isotactic poly(propylene oxide) (iPPO) is a semicrystalline polyether that has emerged as a high-strength, photodegradable material for marine applications. To improve the accessibility of i PPO, catalysts with higher activity and selectivity are required. Using rational catalyst design informed by computational insights, we developed a flexibly tethered, bimetallic chromium catalyst exhibiting high enantioselectivity ( k rel ∼ 100) and unprecedented activity (TOF ∼ 50,000 h −1 ) for propylene oxide (PO) polymerization. Mechanistic studies reveal that high enantioselectivity originates from increased steric bulk at the ortho position of the salicylimine moiety, which increases steric repulsion between the alkoxide chain end and the ligand in the disfavored transition state. Furthermore, introducing geminal dimethyl groups that rigidify the flexible tether between the two ligand moieties significantly enhances catalyst activity by destabilizing the resting state during polymerization. The catalyst remains active at loadings as low as 0.5 ppm, enabling the synthesis of colorless, tough i PPO.
A reaction usually involves a few active modes while the other modes are largely preserved throughout the reaction as spectators. Excitation of an active mode is expected to promote the reaction, but depositing energy in a spectator mode typically has no effect, because of the differing ability for energy flow to the reaction coordinate. Here, we report a surprising case of mode specificity in a key radical–radical reaction OH + HO 2 → H 2 O + O 2 , where such canonical expectations fail to hold. Despite its spectator nature, the vibrational excitation of the OH reactant is shown at low collision energies to enhance the reactivity significantly. This unique effect can be attributed to the increased attraction with HO 2 due to the larger dipole of the stretched OH. At low collision energies, the stronger attraction increases the chance of capturing the reactants to form a hydrogen-bonded complex, thus of passing through the submerged barrier. The novel mechanism differs from the conventional vibrational enhancement via coupling to the reaction coordinate at the transition state, enriching our understanding of mode specificity in chemistry.
Aldol condensations of carbonyl compounds for C–C bond formation are a very important class of reactions in organic synthesis and upgrading of biomass-derived feedstocks. However, the atomic level understanding of reaction mechanisms and structure-activity correlation on widely used transition metal oxide catalysts are limited due to the high degree of struc-tural heterogeneity of catalysts such as commercial TiO 2 powders. Here, we provide a deep understanding of the reaction mechanisms, kinetics, and structure–function relationships for vapor phase acetone aldol condensation through the con-trolled synthesis of two catalysts with high surface areas and clean, dominant facets, coupled with detailed characterization and kinetic studies that are further assisted by density functional theory (DFT) calculations. Temperature-dependent diffuse reflectance infrared Fourier transform spectroscopy showed the existence of abundant acetone bonded to surface hydroxyl groups (acetone-O s H) and acetone bonded to Lewis acid sites (acetone-Ti 5c ) on surface of both {101} and {001} facet dom-inant TiO 2 . Intermolecular C–C coupling of enolate intermediate from acetone-Ti 5c and a vicinal acetone-O s H is a kinetically relevant step, which is consistent with kinetic and isotopic studies and DFT calculations. The {001} facet showed a lower apparent activation energy (or higher reactivity) than the {101} facet. This is likely caused by the weaker Lewis acid and Brønsted base strengths of the {001} facet which favors the reprotonation–desorption of the coupled intermediate, making the C–C coupling step more exothermic on the {001} facet and resulting in an earlier transition state with a lower activation barrier. It is also possible that the {001} facet has a smoother surface configuration and less steric hindrance during intermo-lecular C–C bond formation than the {101} facet.
Carbon dioxide is a potentially valuable feedstock for carbon-based fuels or commodities but is only available in dilute streams. Many studies have focused on either the capture and concentration of CO 2 or the reduction of pure CO 2 streams. The direct reduction of sorbent-captured CO 2 in an integrated process would skip the energy-intensive CO 2 concentration and sorbent regeneration step. Herein, we report the electrocatalytic reduction of 1,3-bis(2,6-diisopropylphenyl)imidazolium-2-carboxylate (IPr·CO 2 ), which forms quantitatively from the reaction of sorbent 1,3-bis(2,6-diisopropylphenyl)imidazol-2-ylidene (IPr) with 10% and 0.04% CO 2 streams, by catalyst [Pt(dmpe) 2 ](PF 6 ) 2 (dmpe = 1,2-bis(dimethylphosphino)ethane) to formate with >70% Faradaic efficiencies. Unexpectedly, experimental studies indicate that the proton source phenol facilitates rapid decarboxylation of IPr·CO 2 to release CO 2 , which is the substrate for reduction. Kinetic studies determined the rate of hydride transfer from a catalytic intermediate [HPt(dmpe) 2 ](PF 6 ) to form the C–H bond in formate to be 0.22 M –1 s –1 . Further details on the mechanism, transition state energy, and structure for hydride transfer to CO 2 , a common step in CO 2 reduction, were explored using computational methods.