Bipolar membranes (BPMs), which have long seen usage in electrodialysis reactors for the generation of acid and base, have recently demonstrated potential to become critical components in electrochemical synthesis devices. 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 junction as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-devices for electrosynthesis applications. Prior work has modeled ion transport in bipolar membranes in neutral salt solutions for electrodialysis. However, no model exists for the BPM under the harsh applied pH gradients that would be present in electrosynthesis, and there is significant need to explore the effects of the internal hydration on the lifetime and performance of BPMs in such environments. Additionally, a mechanistic understanding of water dissociation catalysis will be required to develop interfacial catalysts that enable the high current density operation required for scalable electrosynthesis of fuels. In this talk, we discuss modeling methodologies and physics inherent in BPMs and present our recent model of ion transport and water dissociation catalysis in BPMs across the pH scale. Specifically, we simulate multi-ion transport for a BPM with various electrolyte combinations on each side of the membrane, demonstrating the significance of co- and counter-ion crossover in BPMs operating under harsh pH gradients. We then investigate effects caused by hydration gradients that occur due to internal ion-exchange and examine potential methods for improving performance and mitigating crossover. Finally, we examine the impact of the interfacial water dissociation catalyst and perform sensitivity analysis on the key properties (catalyst point of zero charge and pKa) that dictate catalyst performance. These results provide information that is critical to developing a comprehensive understanding of multi-component phenomena in BPMs and to informing the design and implementation of BPMs in next-generation devices for the numerous electrosynthesis chemistries that benefit from operation under an applied pH gradient.

The electrochemical reduction of CO2 (CO2R) to value-added products is an attractive technology for tackling the rising atmospheric CO2 levels and storing intermittent renewable energy into chemical bonds. Fundamental understanding of CO2R has progressed significantly in recent years and is critical in the development of CO2R to liquid-fuel electrolyzers, where gas-diffusion electrodes (GDEs) have been shown to be key enabling architectures. Various designs have been proposed and studied in the literature to enhance overall selectivity, rates, and maximize the conversion of CO2, the latter of which is only now being recognized as a critical issue. In this respect, there is a need to explore the governing phenomena inherent in these architectures to enable optimization. Mathematical modeling is ideally suited to tackle and explore these multiphysics interactions and provide virtual design analysis. In this talk, we discuss modeling methodologies and physics inherent in these devices and present our recent modeling of GDEs for CO2 reduction. We specifically examine the impacts of multiphase flow and related phenomena on overall cell performance. We then explore the performance and limitations of various cell designs guided by simulation results and examine potential methods for improving water management and tuning catalyst selectivity including the use of different anion-exchange and bipolar membranes. Finally, we discuss the disparities in local environments between aqueous and GDE devices and propose strategies to reduce the gap in knowledge between the two systems. Acknowledgements This material is based upon work performed by the Joint Center for Artificial Photosynthesis, a DOE Energy Innovation Hub, supported through the Office of Science of the U.S. Department of Energy under Award Number DE-SC0004993.

Bicarbonate electrolyzers are devices designed to convert CO2 captured from point sources or the atmosphere into chemicals and fuels without needing to first isolate pure CO2 gas. We report here an experimentally validated model that quantifies the reaction chemistry and mass transfer processes within the catalyst layer and cation exchange membrane layer of a bicarbonate electrolyzer. Our results demonstrate that two distinct chemical microenvironments are key to forming CO at high rates: an acidic membrane layer that promotes in situ CO2 formation and a basic catalyst layer that suppresses the hydrogen evolution reaction. We show that the rate of CO product formation can be increased by modulating the catalyst and membrane layer properties to increase the rate of in situ CO2 generation and transport to the cathode. These insights serve to inform the design of bicarbonate and BPM-based CO2 electrolyzers while demonstrating the value of modeling for resolving rate-determining processes in electrochemical systems.

The electrochemical reduction of CO2 (CO2R) to value-added products is an attractive technology for tackling the rising atmospheric CO2 levels and storing intermittent renewable energy into chemical bonds. Fundamental understanding of CO2R has progressed significantly in recent years and is critical in the development of CO2R to liquid-fuel electrolyzers, where gas-diffusion electrodes (GDEs) have been shown to be key enabling architectures. Various designs have been proposed and studied in the literature to enhance overall selectivity, rates, and maximize the conversion of CO2, the latter of which is only now being recognized as a critical issue. In this respect, there is a need to explore the governing phenomena inherent in these architectures to enable optimization. Mathematical modeling is ideally suited to tackle and explore these multiphysics interactions and provide virtual design analysis. In this talk, we discuss modeling methodologies and physics inherent in these devices and present our recent modeling of GDEs for CO2 reduction. We specifically examine the impacts of multiphase flow and related phenomena on overall cell performance. We then explore the performance and limitations of various cell designs guided by simulation results and examine potential methods for improving water management and tuning catalyst selectivity including the use of different anion-exchange and bipolar membranes. Finally, we discuss the disparities in local environments between aqueous and GDE devices and propose strategies to reduce the gap in knowledge between the two systems. Acknowledgements This material is based upon work performed by the Joint Center for Artificial Photosynthesis, a DOE Energy Innovation Hub, supported through the Office of Science of the U.S. Department of Energy under Award Number DE-SC0004993.

Scaling up CO2 electrolysis is a vital aspect in the transition to manufacturing sustainable fuels and chemical compounds, satisfying the demand for chemicals and demand for storing renewable electricity. Depending on the employed catalyst, different products can be produced, such as carbon monoxide and ethylene, by applying a voltage on a CO2 and H2O fed electrolyzer. CO is a desirable product according to techno-economic analysis, because it can be produced selectively using a silver catalyst and is a precursor for hydrocarbons in the Fischer–Tropsch process. In a state-of-the-art CO2 electrolyzer, two electrodes are directly pressed against an ion-exchange membrane – this is called a zero-gap configuration. Therefore, the membrane is a crucial component for the system, since it has the role of providing a conductive medium between the electrodes. One of the challenges in CO2 electrolysis is that water is consumed in the reaction. At high current density, this may cause the membrane's surface near the cathode to dry out, lowering efficiency or perhaps stopping the process entirely, since there is no longer a conductive medium. As a result, water management is critical for this process. In this work, we approach the drying-out challenge by studying the novel concept of a membrane with internal microchannels. These channels allow the circulation of water or an electrolyte inside the membrane, which reduces the water diffusion path and affects the membrane’s conductivity. The effects of channel geometry, location, and concentrations of electrolyte inside on water content, conductivity and overall performance are studied in a 2D COMSOL model. In addition, the effect of internal concentration of electrolyte on the membrane’s resistance, on the process performance and the K+ cross-over to cathode side were investigated experimentally. Our modeling results prove that the presence of the channels can keep the membrane hydrated. The highest current densities are observed when the channel is closest to the cathode, and with smaller pores. Smaller pores are advantageous due to the trade-off between enhanced membrane conductivity and the lower conductivity of the liquid itself. If water is circulated in a large channel, it increases the membrane’s ionomer conductivity due to hydration but the overall conductivity is decreased since water is not highly conductive. Nonetheless, the results also show that a higher concentration of electrolyte inside the microchannels can significantly increase the total conductivity of the membrane, and therefore the energy efficiency of the process. These effects are most significant at higher current densities. In the experimental results, we’ve observed similar effects in terms of membrane conductivity and current density of the process – the higher the electrolyte concentration the higher the current density. Furthermore, a low concentration decreases the amount of potassium which crosses over to the cathode side, inhibiting salt deposition. We’ve concluded that a small channel, up to 90 µm wide, close to the pore with an electrolyte with a concentration of up to 10 mM could be very beneficial for the water management and energy efficiency of the process. This helps to keep the membrane hydrated at higher current densities, improves the conductivity of the membrane, and it doesn’t have significant impacts on the salt deposition. Figure 1

Highlights Photoelectrochemical CO2 reduction to multicarbon products with high efficiency. Co-design of photoelectrochemical CO2 reduction devices for multicarbon product generation. Optimal semiconductor bandgaps for high multicarbon product selectivity. Limitations in multicarbon product generation from photoelectrochemical CO2 reduction. Optimizing photoelectrochemical CO2 reduction systems for multicarbon product generation.

Bipolar membranes (BPMs) possess the potential to optimize pH environments for electrochemical synthesis applications when employed in reverse bias. Unfortunately, the performance of BPMs in reverse bias has long been limited by the rate of water dissociation (WD) occurring at the interface of the BPM. Herein, we develop a continuum model of the BPM that agrees with experiment to understand and enhance WD catalyst performance by considering multiple kinetic pathways for WD in the BPM junction catalyst layer. The model reveals that WD catalysts with a more highly alkaline or acidic pH at the point of zero charge (pHPZC) exhibit accelerated WD kinetics because the more acidic or alkaline pHPZC catalysts possess greater surface charge, enhancing the local electric field and rate of WD. The model is then employed to explore the sensitivity of the BPM performance to various BPM physical parameters. Finally, the model is used to simulate the operation of bimetallic WD catalysts, demonstrating that an optimal bimetallic catalyst has an acidic pHPZC catalyst matched with the cation-exchange layer and an alkaline pHPZC catalyst matched with the anion-exchange layer. The study provides insight into the operation of BPM WD catalysts and gives direction toward the development of next-generation WD catalysts for optimal BPM performance under water-splitting and related conditions.

Pulsed electrolysis has been demonstrated to improve the faradaic efficiency (FE) to C2+ products during the electrochemical reduction of CO2 over a Cu catalyst, but the nature of this enhancement is poorly understood. Herein, we developed a time-dependent continuum model of pulsed CO2 electrolysis on Cu in 0.1 M CsHCO3 that faithfully represents the experimentally observed effects of pulsed electrolysis. This work shows that pulsing results in dynamic changes in the pH and CO2 concentration near the Cu surface, which lead to an enhanced C2+ FE as a consequence of repeatedly accessing a transient state of heightened pH and CO2 concentration at high cathodic overpotential. Using these insights, a variety of pulse shapes were explored to establish operating conditions that maximize the rate of C2+ product formation and minimize the rates of H2 and C1 product formation.

There is increasing interest in the possibility of photoelectrochemical (PEC) reduction of CO2to C2+products; however, the criteria for maximizing PEC solar-to-C2+(STC2+) rates are not well understood. We report here a continuum-scale model of PEC CO2reduction (CO2R) on Cu in 0.1 M CsHCO3and use it to optimize the design and operating conditions for generating C2+products. We demonstrate that the potential-dependent product distribution of CO2R on Cu requires operating near the potential that maximizes C2+generation rates (Vid), unlike PEC water splitting, which desires operation at the maximum photocurrent density. Because of this requirement, the criterion for a high STC2+rate includes high-photocurrent semiconductors with photovoltages near Vidand low series resistance. The STC2+rate in these systems is enhanced by optimal CO2transport and exhibits low sensitivity to dirunal solar irradiance variations.