skip to main content

Title: Selective electrochemical oxidative coupling of methane mediated by Sr2Fe1.5Mo0.5O6-δ and its chemical stability

Efficient conversion of methane to value-added products such as olefins and aromatics has been in pursuit for the past few decades. The demand has increased further due to the recent discoveries of shale gas reserves. Oxidative and non-oxidative coupling of methane (OCM and NOCM) have been actively researched, although catalysts with commercially viable conversion rates are not yet available. Recently,$${{{{{{{\mathrm{Sr}}}}}}}}_2Fe_{1.5 + 0.075}Mo_{0.5}O_{6 - \delta }$$Sr2Fe1.5+0.075Mo0.5O6δ(SFMO-075Fe) has been reported to activate methane in an electrochemical OCM (EC-OCM) set up with a C2 selectivity of 82.2%1. However, alkaline earth metal-based materials are known to suffer chemical instability in carbon-rich environments. Hence, here we evaluated the chemical stability of SFMO in carbon-rich conditions with varying oxygen concentrations at temperatures relevant for EC-OCM. SFMO-075Fe showed good methane activation properties especially at low overpotentials but suffered poor chemical stability as observed via thermogravimetric, powder XRD, and XPS measurements where SrCO3was observed to be a major decomposition product along with SrMoO3and MoC. Nevertheless, our study demonstrates that electrochemical methods could be used to selectively activate methane towards partial oxidation products such as ethylene at low overpotentials while higher applied biases result in the complete oxidation of methane to carbon dioxide and water.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Communications Chemistry
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Perovskite oxides (ternary chemical formula ABO3) are a diverse class of materials with applications including heterogeneous catalysis, solid-oxide fuel cells, thermochemical conversion, and oxygen transport membranes. However, their multicomponent (chemical formula$${A}_{x}{A}_{1-x}^{\text{'}}{B}_{y}{B}_{1-y}^{\text{'}}{O}_{3}$$AxA1x'ByB1y'O3) chemical space is underexplored due to the immense number of possible compositions. To expand the number of computed$${A}_{x}{A}_{1-x}^{{\prime} }{B}_{y}{B}_{1-y}^{{\prime} }{O}_{3}$$AxA1xByB1yO3compounds we report a dataset of 66,516 theoretical multinary oxides, 59,708 of which are perovskites. First, 69,407$${A}_{0.5}{A}_{0.5}^{{\prime} }{B}_{0.5}{B}_{0.5}^{{\prime} }{O}_{3}$$A0.5A0.5B0.5B0.5O3compositions were generated in theab+aGlazer tilting mode using the computationally-inexpensive Structure Prediction and Diagnostic Software (SPuDS) program. Next, we optimized these structures with density functional theory (DFT) using parameters compatible with the Materials Project (MP) database. Our dataset contains these optimized structures and their formation (ΔHf) and decomposition enthalpies (ΔHd) computed relative to MP tabulated elemental references and competing phases, respectively. This dataset can be mined, used to train machine learning models, and rapidly and systematically expanded by optimizing more SPuDS-generated$${A}_{0.5}{A}_{0.5}^{{\prime} }{B}_{0.5}{B}_{0.5}^{{\prime} }{O}_{3}$$A0.5A0.5B0.5B0.5O3perovskite structures using MP-compatible DFT calculations.

    more » « less
  2. Abstract

    Massive gully land consolidation projects, launched in China’s Loess Plateau, aim to restore 2667$$\mathrm{km}^2$$km2agricultural lands in total by consolidating 2026 highly eroded gullies. This effort represents a social engineering project where the economic development and livelihood of the farming families are closely tied to the ability of these emergent landscapes to provide agricultural services. Whether these ‘time zero’ landscapes have the resilience to provide a sustainable soil condition such as soil organic carbon (SOC) content remains unknown. By studying two watersheds, one of which is a control site, we show that the consolidated gully serves as an enhanced carbon sink, where the magnitude of SOC increase rate (1.0$$\mathrm{g\,C}/\mathrm{m}^2/\mathrm{year}$$gC/m2/year) is about twice that of the SOC decrease rate (− 0.5$$\mathrm{g\,C}/\mathrm{m}^2/\mathrm{year}$$gC/m2/year) in the surrounding natural watershed. Over a 50-year co-evolution of landscape and SOC turnover, we find that the dominant mechanisms that determine the carbon cycling are different between the consolidated gully and natural watersheds. In natural watersheds, the flux of SOC transformation is mainly driven by the flux of SOC transport; but in the consolidated gully, the transport has little impact on the transformation. Furthermore, we find that extending the surface carbon residence time has the potential to efficiently enhance carbon sequestration from the atmosphere with a rate as high as 8$$\mathrm{g\,C}/\mathrm{m}^2/\mathrm{year}$$gC/m2/yearcompared to the current 0.4$$\mathrm{g\,C}/\mathrm{m}^2/\mathrm{year}$$gC/m2/year. The success for the completion of all gully consolidation would lead to as high as 26.67$$\mathrm{Gg\,C}/\mathrm{year}$$GgC/yearsequestrated into soils. This work, therefore, not only provides an assessment and guidance of the long-term sustainability of the ‘time zero’ landscapes but also a solution for sequestration$$\hbox {CO}_2$$CO2into soils.

    more » « less
  3. Abstract

    In this work, we develop a deep neural network model for the reaction rate of oxidative coupling of methane from published high-throughput experimental catalysis data. A neural network is formulated so that the rate model satisfies the plug flow reactor design equation. The model is then employed to understand the variation of reactant and product composition within the reactor for the reference catalyst Mn–Na2WO4/SiO2at different temperatures and to identify new catalysts and combinations of known catalysts that would increase yield and selectivity relative to the reference catalyst. The model revealed that methane is converted in the first half of the catalyst bed, while the second part largely consolidates the products (i.e. increases ethylene to ethane ratio). A screening study of3400combinations of pairs of previously studied catalysts of the form M1(M2)12M3Ox/support (where M1, M2 and M3 are metals) revealed that a reactor configuration comprising two sequential catalyst beds leads to synergistic effects resulting in increased yield of C2compared to the reference catalyst at identical conditions and contact time. Finally, an expanded screening study of 7400 combinations (comprising previously studied metals but with several new permutations) revealed multiple catalyst choices with enhanced yields of C2products. This study demonstrates the value of learning a deep neural network model for the instantaneous reaction rate directly from high-throughput data and represents a first step in constraining a data-driven reaction model to satisfy domain information.

    more » « less
  4. Abstract

    Thin film evaporation is a widely-used thermal management solution for micro/nano-devices with high energy densities. Local measurements of the evaporation rate at a liquid-vapor interface, however, are limited. We present a continuous profile of the evaporation heat transfer coefficient ($$h_{\text {evap}}$$hevap) in the submicron thin film region of a water meniscus obtained through local measurements interpreted by a machine learned surrogate of the physical system. Frequency domain thermoreflectance (FDTR), a non-contact laser-based method with micrometer lateral resolution, is used to induce and measure the meniscus evaporation. A neural network is then trained using finite element simulations to extract the$$h_{\text {evap}}$$hevapprofile from the FDTR data. For a substrate superheat of 20 K, the maximum$$h_{\text {evap}}$$hevapis$$1.0_{-0.3}^{+0.5}$$1.0-0.3+0.5 MW/$$\text {m}^2$$m2-K at a film thickness of$$15_{-3}^{+29}$$15-3+29 nm. This ultrahigh$$h_{\text {evap}}$$hevapvalue is two orders of magnitude larger than the heat transfer coefficient for single-phase forced convection or evaporation from a bulk liquid. Under the assumption of constant wall temperature, our profiles of$$h_{\text {evap}}$$hevapand meniscus thickness suggest that 62% of the heat transfer comes from the region lying 0.1–1 μm from the meniscus edge, whereas just 29% comes from the next 100 μm.

    more » « less
  5. A theoretical analysis on crack formation and propagation was performed based on the coupling between the electrochemical process, classical elasticity, and fracture mechanics. The chemical potential of oxygen, thus oxygen partial pressure, at the oxygen electrode-electrolyte interface (μO2OE∣El) was investigated as a function of transport properties, electrolyte thickness and operating conditions (e.g., steam concentration, constant current, and constant voltage). Our analysis shows that: a lower ionic area specific resistance (ASR),riOE,and a higher electronic ASR (reOE) of the oxygen electrode/electrolyte interface are in favor of suppressing crack formation. TheμO2OEEl,thus local pO2, are sensitive towards the operating parameters under galvanostatic or potentiostatic electrolysis. Constant current density electrolysis provides better robustness, especially at a high current density with a high steam content. While constant voltage electrolysis leads to greater variations ofμO2OEEl.Constant current electrolysis, however, is not suitable for an unstable oxygen electrode becauseμO2OEElcan reach a very high value with a gradually increasedriOE.A crack may only occur under certain conditions whenpO2TPB>pcr.

    more » « less