Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes’ innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Science & Technology - Other Topics↗
DOE OSTI2024

Polymer Layer-Accelerated CO 2 Absorption in Aqueous Amino Acid Solutions

Direct air capture (DAC) of CO 2 via solvent-based absorption is considered a promising negative-emission technology. However, the low concentration of CO 2 in the air and slow transport into the solvent make DAC notoriously challenging to implement without costly investments. In this study, we explore the fundamental role that the bulk and surface properties of CO 2 -permeable polymer membranes play in enhancing the efficiency of the solution sorption process in passive DAC of CO 2 . This work leverages various spectroscopic and computational studies to demonstrate that a hybrid system, comprising a reusable CO 2 -permeable polymer layer placed atop an aqueous amino acid (AA) solution, can outperform a pure aqueous AA system by 2-fold. Here we show how the enhanced solubility of CO 2 in the polymer layer can improve the transport of CO 2 into the aqueous phase, while the chemistry of the polymer can control the interfacial barrier for CO 2 permeation and the interfacial concentration of reactive AAs. The derived knowledge of the material properties achieved here can aid in the design of DAC systems with improved performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Analysis of Flow Area in the Axial Cross-Vessel Annular Duct Following a Breach of the Primary Helium Boundary of an HTGR

The lock-exchange flow phenomenon has been heavily studied for decades, including the lock-exchange phenomenon that takes place during an air ingress accident in a High-Temperature Gas-cooled Reactor (HTGR) during a break on the axial cross-vessel duct. The present study aims to characterize the gas dynamics in the cavity region of an HTGR, particularly the axial cross-vessel outer duct region. This study computes the velocity and flow area of the heavy and light currents when a break on the primary helium boundary occurs. This analysis is supported using experimental data from experiments conducted in a 1/20th length-scale HTGR. The following study summarizes the conditions of the experiment, the approach adopted to predict the flow area, and future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multimetallic Metal-Organic Frameworks as Heterogeneous Catalysts for Gas Phase Hydroformylation and Hydrogenation Reactions

This project focused on the development bimetallic metal-organic frameworks (MOFs) as gas phase heterogeneous catalysts for hydrogenation and hydroformylation reactions. MOFs are a new class of hybrid inorganic/organic materials that are highly crystalline, with structures consisting of metal nodes of specific geometries connected by organic linkers. Although there have been a number of studies of catalysis at MOF nodes in solution, there is little experimental data in the literature for gas phase reactions despite the fact that industrial heterogeneous catalysis on MOFs is more economically viable than homogeneous catalysis. The use of MOFs as heterogeneous catalysts presents the unique opportunity to carefully control the composition, geometry and ensemble sizes of the active sites, which are all critical factors for the rational design of new catalysts. Specific objectives of the project are as follows: (1) to tailor the geometry, composition and ensemble size of the active sites; (2) to understand how the adsorption of molecules at metal sites can be modified by interactions with a neighboring metal site; (3) to determine how oxidation states of the metal change during reaction and how these states can be modified by metal-metal electronic interactions; and (4) to elucidate reaction mechanisms and intermediates. For these studies, we have chosen to investigate selective hydrogenation of hydrocarbons and hydroformylation, which are industrially relevant reactions. Our interdisciplinary team of Chen, Shustova and Vogiatzis/Henkelman provides critical expertise in novel MOF synthesis (Shustova), atomic-scale surface science and catalysis (Chen) and computational studies of reaction mechanisms (Vogiatzis/Henkelman) that are necessary for the development of these catalysts. We believe that this work will lead directly to the rational design of new catalysts that are versatile, highly active/selective and suitable for industrial processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Deciphering the Olefin Isomerization-Polymerization Paradox of Palladium(II) Diimine Catalysts: Discovery of Simultaneous and Independent Pathways of Olefin Isomerization and Living Polymerization

This work elucidates a long-standing unexplained paradox commonly observed within the polymerization of α-olefin using palladium (Pd)(II)–diimine catalysts, in which isomerization and living polymerization of α-olefins are both observed. With a classical mechanistic understanding of these complexes, this behavior is often dismissed and interpreted as experimental error. Herein, we present a comprehensive mechanistic investigation into this phenomenon that supports the existence of a novel mechanistic pathway for Pd(II)–diimine complexes. Part one of the mechanistic study lays the foundation of the proposed mechanism, in which neutral Pd(II)–diimine complexes were found to exhibit a moderate to good catalytic activity for olefin isomerization of α-olefins despite the established notion that catalyst activation is required. Extensive experimental and computational studies reveal the possibility of a partial dissociation of the diimine ligand, which frees up one coordination site and enables coordination–insertion. This finding is significant as the coexistence of two reactive coordination sites at the palladium center becomes a valid proposal for the activated cationic Pd(II)–diimine complexes. In part two, we examined and validated the simultaneously observed α-olefin isomerization and living polymerization using the cationic Pd(II)–diimine catalyst, which supports the presence of two independent reaction pathways of isomerization and polymerization, respectively. Furthermore, the addition of a strong Lewis acid, such as AlCl 3 , accelerates the ligand dissociation and the consequential isomerization as it weakens the palladium–nitrogen bond through competitive binding. In part three, Lewis acid-triggered olefin isomerization-polymerization is employed to prepare living olefinic block copolymers and further synthesize novel polyolefin-polar block copolymers with unique architectures, distinct levels of branching, crystallinity, and polar functionality in a one-pot manner.

Catalysts↗
DOE OSTI2026

Computational Analysis of Hydraulic Efficiency of Michigan DOT Covers J and K

Drainage structures are used in urban street and highway systems to capture stormwater runoff. These structures typically consist of catch basins fitted with grates, inlets, or combination grate/inlet configurations that collect runoff and convey it through buried drainage systems. They are strategically placed within curb-and-gutter systems to enhance public safety by efficiently removing water from roadways and thereby reducing the risk of hydroplaning. The performance of drainage structures is commonly evaluated in terms of hydraulic efficiency, defined as the percentage of flow captured by the basin relative to the total flow reaching the structure. Understanding the hydraulic performance of these structures allows designers to properly space inlets, resulting in cost-effective designs that also help ensure the safety of the traveling public. MDOT uses a variety of drainage structures for runoff capture, as documented in Michigan Department of Transportation (MDOT) Drainage Manual. Many of these structures incorporate sinusoidal-type grates that are not addressed in HEC-22. Because physical modeling of these structures has been limited, further analysis is needed to verify their capture efficiency. Under current MDOT practice, the capture efficiency of these grates is estimated by assuming performance similar to that of a comparably sized reticuline grate described in HEC22. The first phase of this effort, titled Computational Analysis of Hydraulic Efficiency of Michigan DOT Cover C, focused on evaluating the hydraulic performance of MDOT’s Cover C grate. Cover C was selected as the initial test candidate because its sinusoidal pattern is representative of other MDOT grates, while it is typically used in high-volume, higher-speed applications. A similar version, Cover CX, is used on interstate highways but does not include transverse bars for bicycle safety. The current second phase of the study expands this work to evaluate MDOT’s Covers J and K. These grates were selected for additional analysis to further assess the hydraulic performance of MDOT drainage structures that are not directly represented by grate configurations in HEC-22. The results of this phase will build on the findings from the Cover C analysis and support improved understanding of the capture efficiency of MDOT’s standard drainage grates. This report is intended to serve as a companion document to the earlier study, Computational Analysis of Hydraulic Efficiency of Michigan DOT Cover C [3]. The present work applies the same overall CFD-based evaluation approach to MDOT Covers J and K and compares the resulting performance trends with those previously identified for Cover C. In particular, both studies assess on-grade interception efficiency, sag-location hydraulic capacity, and the effects of partial obstruction relative to HEC-22-based design estimates.

42 ENGINEERING↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Observation of single-defect memristor in an MoS2 atomic sheet

Non-volatile resistive switching, also known as memristor effect, where an electric field switches the resistance states of a two-terminal device, has emerged as an important concept in the development of high-density information storage, computing and reconfigurable systems. The past decade has witnessed substantial advances in non-volatile resistive switching materials such as metal oxides and solid electrolytes. It was long believed that leakage currents would prevent the observation of this phenomenon for nanometre-thin insulating layers. However, the recent discovery of non-volatile resistive switching in two-dimensional monolayers of transition metal dichalcogenide and hexagonal boron nitride sandwich structures (also known as atomristors) has refuted this belief and added a new materials dimension owing to the benefits of size scaling. Here we elucidate the origin of the switching mechanism in atomic sheets using monolayer MoS 2 as a model system. Atomistic imaging and spectroscopy reveal that metal substitution into a sulfur vacancy results in a non-volatile change in the resistance, which is corroborated by computational studies of defect structures and electronic states. These findings provide an atomistic understanding of non-volatile switching and open a new direction in precision defect engineering, down to a single defect, towards achieving the smallest memristor for applications in ultra-dense memory, neuromorphic computing and radio-frequency communication systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Singlet Fission

This report summarizes key findings from DOE Award DE‑SC0007004 – Singlet Fission, documenting more than a decade of research led by Professor Josef Michl. The project advanced foundational understanding of singlet fission (SF) through integrated synthetic, spectroscopic, crystallographic, and computational studies. Across a diverse set of chromophores—including tetracenes, cibalackrot derivatives, BODIPY dimers, TDPPs, and fluorinated isobenzofurans—the team identified structural and packing features that enable or inhibit SF, clarified the roles of charge‑transfer intermediates and excimer formation, and developed simplified theoretical models and computational screening tools for predicting optimal geometries and excited‑state-energetics. The work significantly broadened design principles for SF‑active materials and provided insights relevant to next‑generation solar energy conversion and molecular photophysics.

14 SOLAR ENERGY↗
DOE OSTI2020

A missense variant remote from the active site impairs stability of human phosphoglucomutase 1

Abstract Missense variants of human phosphoglucomutase 1 (PGM1) cause the inherited metabolic disease known as PGM1 deficiency. This condition is categorised as both a glycogen storage disease and a congenital disorder of glycosylation. Approximately 20 missense variants of PGM1 are linked to PGM1 deficiency, and biochemical studies have suggested that they fall into two general categories: those affecting the active site and catalytic efficiency, and those that appear to impair protein folding and/or stability. In this study, we characterise a novel variant of Arg422, a residue distal from the active site of PGM1 and the site of a previously identified disease‐related variant (Arg422Trp). In prior studies, the R422W variant was found to produce insoluble protein in a recombinant expression system, precluding further in vitro characterisation. Here we investigate an alternative variant of this residue, Arg422Gln, which is amenable to experimental characterisation presumably due to its more conservative physicochemical substitution. Biochemical, crystallographic, and computational studies of R422Q establish that this variant causes only minor changes in catalytic efficiency and 3D structure, but is nonetheless dramatically reduced in stability. Unexpectedly, binding of a substrate analog is found to further destabilise the protein, in contrast to its stabilising effect on wild‐type PGM1 and several other missense variants. This work establishes Arg422 as a lynchpin residue for the stability of PGM1 and supports the impairment of protein stability as a pathomechanism for variants that cause PGM1 deficiency. Synopsis Biochemical and structural studies of a missense variant far from the active site of human PGM1 identify a residue with a key role in enzyme stability.

Stiers, Kyle M.↗
DOE OSTI2023

Sulfur Lone Pairs Control Topology in Heterotrimetallic Complexes: An Experimental and Theoretical Study

Heterotrimetallic complexes with (N 2 S 2 )M metallodithiolates, M = Ni 2+ , [Fe(NO)] 2+ , and [Co(NO)] 2+ , as bidentate chelating ligands to a central trans-Cr(NO)(MeCN) unit were characterized as the first members of a new class, NiCrNi, FeCrFe, CoCrCo. The complexes exhibit a cisoid structural topology, ascribed to the stereoactivity of the available lone pair(s) on the sulfur donors, resulting in a dispersed, electropositive pocket from the N/N and N/S hydrocarbon linkers wherein the Cr-NO site is housed. Computational studies explored alternative isomers (transoid and inverted cisoid) that suggest a combination of electronic and steric effects govern the geometrical selectivity. Electrostatic potential maps readily display the dominant electronegative potential from the sulfurs which force the NO to the electropositive pocket. The available S lone pairs work in synergy with the π-withdrawing ability of NO to lift Cr out of the S 4 plane toward the NO and stabilize the geometry. The metallodithiolate ligands bound to Cr(NO) thus find structural consistency across the three congeners. Although the dinitrosyl [(bme-dach)Co(NO)-Mo(NO)(MeCN)-(bme-dach)Co(MeCN)][PF 6 ] 2 (CoMoCo') analogue displays chemical noninnocence and a partial Mo–Co bond toward (N 2 S 2 )Co'(NCCH 3 ) in an “asymmetric butterfly” topology, the stability of the {Cr(NO)} 5 unit prohibits such bond rearrangement. Magnetism and EPR studies illustrate spin coupling across the sulfur thiolate sulfur bridges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Solution Chemistry to Control Boron-Containing Monolayers on Silicon: Reactions of Boric Acid and 4-Fluorophenylboronic Acid with H- and Cl-terminated Si(100)

In this work, the reactions of boric acid and 4-fluorophenylboronic acid with H- and Cl-terminated Si(100) surfaces in solution were investigated. X-ray photoelectron spectroscopy (XPS) studies reveal that both molecules react preferentially with Cl–Si(100) and not with H–Si(100) at identical conditions. On Cl–Si(100), the reactions introduce boron onto the surface, forming a Si–O–B structure. The quantification of boron surface coverage demonstrates that the 4-fluorophenylboronic acid leads to ~2.8 times higher boron coverage compared to that of boric acid on Cl–Si(100). Consistent with these observations, density functional theory studies show that the reaction of boric acid and 4-fluorophenylboronic acid is more favorable with the Cl- versus H-terminated surface and that on Cl–Si(100) the reaction with 4-fluorophenylboronic acid is ~55.3 kJ/mol more thermodynamically favorable than the reaction with boric acid. The computational studies were also used to demonstrate the propensity of the overall approach to form high-coverage monolayers on these surfaces, with implications for selective-area boron-based monolayer doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗
DOE OSTI2022

Visible-Light-Driven, Iridium-Catalyzed Hydrogen Atom Transfer: Mechanistic Studies, Identification of Intermediates, and Catalyst Improvements

The harvesting of visible light is a powerful strategy for the synthesis of weak chemical bonds involving hydrogen that are below the thermodynamic threshold for spontaneous H 2 evolution. Piano-stool iridium hydride complexes are effective for the bluelight- driven hydrogenation of organic substrates and contrathermodynamic dearomative isomerization. In this work, a combination of spectroscopic measurements, isotopic labeling, structure–reactivity relationships, and computational studies has been used to explore the mechanism of these stoichiometric and catalytic reactions. Photophysical measurements on the iridium hydride catalysts demonstrated the generation of long-lived excited states with principally metal-to-ligand charge transfer (MLCT) character. Transient absorption spectroscopic studies with a representative substrate, anthracene revealed a diffusion-controlled dynamic quenching of the MLCT state. The triplet state of anthracene was detected immediately after the quenching events, suggesting that triplet–triplet energy transfer initiated the photocatalytic process. The key role of triplet anthracene on the post-energy transfer step was further demonstrated by employing photocatalytic hydrogenation with a triplet photosensitizer and a HAT agent, hydroquinone. DFT calculations support a concerted hydrogen atom transfer mechanism in lieu of stepwise electron/proton or proton/electron transfer pathways. Kinetic monitoring of the deactivation channel established an inverse kinetic isotope effect, supporting reversible C(sp 2 )–H reductive coupling followed by rate-limiting ligand dissociation. Mechanistic insights enabled design of a piano-stool iridium hydride catalyst with a rationally modified supporting ligand that exhibited improved photostability under blue light irradiation. The complex also provided improved catalytic performance toward photoinduced hydrogenation with H 2 and contra-thermodynamic isomerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Development of Martini 2.2 parameters for N -glycans: a case study of the HIV-1 Env glycoprotein dynamics

N-linked glycans are ubiquitous in nature and play key roles in biology. For example, glycosylation of pathogenic proteins is a common immune evasive mechanism, hampering the development of successful vaccines. Due to their chemical variability and complex dynamics, an accurate molecular understanding of glycans is still limited by the lack of effective resolution of current experimental approaches. Here, we have developed and implemented a reductive model based on the popular Martini 2.2 coarse-grained force field for the computational study of N-glycosylation. We used the HIV-1 Env as a direct applied example of a highly glycosylated protein. Finally, our results indicate that the model not only reproduces many observables in very good agreement with a fully atomistic force field but also can be extended to study large amount of glycosylation variants, a fundamental property that can aid in the development of drugs and vaccines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Carbonate Management to Enable Energy- and Carbon-Efficient CO 2 Electrolysis (Final Technical Report)

The rapid growth and plummeting cost of solar energy have spurred growing interest in using CO 2 electrolysis to produce chemicals and fuels as an alternative to conventional petrochemical processes. High-temperature (>800 °C) solid oxide electrolyzers that convert CO 2 into CO and O 2 have recently become commercially available. Low-temperature electrolysis cells offer the prospect of more convenient and flexible operation, which is critical for utilizing intermittent solar energy, and provide access to more valuable C 2+ products such as ethylene, ethanol, and propanol. Over the past 10 years, research in this area has yielded substantial progress in both fundamental understanding of the requisite electrocatalytic reactions and design of prototype devices. Leveraging insights from fuel cells and membrane water electrolyzers, researchers have developed electrolysis cells with gas diffusion electrodes (GDE) that have demonstrated high CO 2 electrolysis current densities (>100 mA cm –2 ) as well as promising selectivity and stability. Despite these advances, the energy efficiency (electrical energy-to-product) and carbon efficiency (CO 2 -to-product) of low-temperature CO 2 electrolysis remain far too low for large-scale deployment. A preponderance of evidence indicates that the principal source of efficiency losses is the rapid and thermodynamically favorable reaction of CO 2 with hydroxide (OH – ) to form carbonate (CO 3 2– ). Carbonate formation imposes steady-state electrolysis conditions that result in large voltage and CO 2 losses for all known (photo)electrochemical CO 2 cells. While much current research remains focused on CO 2 reduction catalyst design, this largely overlooked CO 3 2– problem presents a fundamental scientific barrier to creating a viable electrochemical option for converting solar energy into chemicals and fuels. The project pursues an integrated, multi-PI research effort that establishes a fundamental science of CO 3 2– management. PI Kanan and Co-PI Mani’s contribution to the project is to evaluate strategies to mitigate the CO 3 2– problem by changing the properties of the electrolyte and the environment in which CO 2 reduction catalysis takes place. Experimental studies showed that electrolytes composed of a high concentration of both CO 3 2– and HCO 3 – , which serve as moderately alkaline buffers, improved the cell voltage by compared to all-HCO 3 – electrolytes, but these buffered systems still show substantial CO 2 uptake that reduces pH over time. Computational studies developed a homogenized model of a CO 2 reduction catalyst layer that permits a low-cost exploration of the high-dimensional parameter space associated with catalyst layers on gas diffusion electrodes. The model was validated by accurately reproducing experimental data for the related but simpler reaction of CO reduction and then used to probe the effects of catalyst layer architecture on CO 2 reduction. Minimizing the size of catalyst and hydrophobic domains in the catalyst layer is predicted to mitigate CO 3 2– formation and thereby enable prolonged operation at elevated pH. In support of the CO 2 electrolysis studies, a new method for rapidly prototyping electrochemical cells was developed and validated. The method uses a combination of 3D printing and electroless plating to generate conductive cell components for evaluating new cell designs. The carbonate problem encompasses mass transport processes and acid-base reactions that are relevant to many other electrochemical systems. Investigation of strategies to address the carbonate problem led to an additional line of inquiry into the physicochemical phenomena that determine the efficiency of electrochemical acid-base production, which has numerous applications in the broader field of carbon management. New strategies for using the supporting electrolyte to inhibit H + /OH – recombination in electrochemical acid-base production were evaluated, leading to the development of a novel acid-base producing system that eliminates the need for ion exchange membranes and exhibits promising efficiency and current densities for scalable applications.

25 ENERGY STORAGE↗