Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Oxygen evolution”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
DOE OSTI2021

Chromium Isotopic Evidence for Mixing of NC and CC Reservoirs in Polymict Ureilites: Implications for Dynamical Models of the Early Solar System

Nucleosynthetic isotope anomalies show that the first few million years of solar system history were characterized by two distinct cosmochemical reservoirs, CC (carbonaceous chondrites and related differentiated meteorites) and NC (the terrestrial planets and all other groups of chondrites and differentiated meteorites), widely interpreted to correspond to the outer and inner solar system, respectively. At some point, however, bulk CC and NC materials became mixed, and several dynamical models offer explanations for how and when this occurred. We use xenoliths of CC materials in polymict ureilite (NC) breccias to test the applicability of such models. Polymict ureilites represent regolith on ureilitic asteroids but contain carbonaceous chondrite-like xenoliths. We present the first 54 Cr isotope data for such clasts, which, combined with oxygen and hydrogen isotopes, show that they are unique CC materials that became mixed with NC materials in these breccias. It has been suggested that such xenoliths were implanted into ureilites by outer solar system bodies migrating into the inner solar system during the gaseous disk phase ~3–5 Myr after CAI, as in the "Grand Tack" model. However, combined textural, petrologic, and spectroscopic observations suggest that they were added to ureilitic regolith at ~50–60 Myr after CAI, along with ordinary, enstatite, and Rumuruti-type chondrites, as a result of the breakup of multiple parent bodies in the asteroid belt at this time. This is consistent with models for an early instability of the giant planets. The C-type asteroids from which the xenoliths were derived were already present in inner solar system orbits.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2017

Application of New Partition Coefficients to Modeling Plagioclase

Previously, studies that determined the partition coefficient for an element, i, between plagioclase and the residual basaltic melt (Di plag) have been conducted using experimental conditions dissimilar from the Moon, and thus these values are not ideal for modeling plagioclase fractionation in a lunar system. However, recent work [1] has determined partition coefficients for plagioclase at lunar oxygen fugacities, and resulted in plagioclase with Anorthite contents ≥An90; these are significantly more calcic than plagioclase in previous studies, and the An content has a profound effect on partition coefficient values [2,3]. Plagioclase D-values, which are dependent on the An content of the crystal [e.g., 2-6], can be determined using published experimental data and the correlative An contents. Here, we examine new experimental data from [1] to ascertain their effect on the calculation of equilibrium liquids from Apollo 16 sample 60635,2. This sample is a coarse grained, subophitic impact melt composed of ~55% plagioclase laths with An94.4-98.7 [7,8], distinctly more calcic than of previous partition coefficient studies (e.g., [3-6, 9-10]). Sample 60635,2 is notable as having several plagioclase trace element analyses containing a negative Europium anomaly (-Eu) in the rare-earth element (REE) profile, rather than the typical positive Eu anomaly (+Eu) [7-8] (Fig. 1). The expected +Eu is due to the similarity in size and charge with Ca2+, thereby allowing Eu2+ to be easily taken up by the plagioclase crystal structure, in contrast to the remaining REE3+. Some 60635,2 plagioclase crystals only have +Eu REE profiles, some only have -Eu REE profiles, and some +Eu and -Eu analyses in different areas on a single crystal [7, 8]. Moreover, there does not seem to be any core-rim association with the +Eu or -Eu analyses, nor does there appear to be a correlation between the size, shape, or location of a particular crystal within the sample and the sign of its Eu anomaly, which suggests a complex evolution. In order to investigate this sample further, we can calculate the equilibrium liquids, but with An contents distinct from previous experimental studies, we must calculate the appropriate partition coefficients for each trace element analysis.

Fagan, A. L.↗
DOE OSTI2023

Laisk measurements in the nonsteady state: Tests in plants exposed to warming and variable CO2 concentrations

Abstract Light respiration (RL) is an important component of plant carbon balance and a key parameter in photosynthesis models. RL is often measured using the Laisk method, a gas exchange technique that is traditionally employed under steady-state conditions. However, a nonsteady-state dynamic assimilation technique (DAT) may allow for more rapid Laisk measurements. In 2 studies, we examined the efficacy of DAT for estimating RL and the parameter Ci* (the intercellular CO2 concentration where Rubisco's oxygenation velocity is twice its carboxylation velocity), which is also derived from the Laisk technique. In the first study, we compared DAT and steady-state RL and Ci* estimates in paper birch (Betula papyrifera) growing under control and elevated temperature and CO2 concentrations. In the second, we compared DAT-estimated RL and Ci* in hybrid poplar (Populus nigra L. × P. maximowiczii A. Henry “NM6”) exposed to high or low CO2 concentration pre-treatments. The DAT and steady-state methods provided similar RL estimates in B. papyrifera, and we found little acclimation of RL to temperature or CO2; however, Ci* was higher when measured with DAT compared to steady-state methods. These Ci* differences were amplified by the high or low CO2 pre-treatments. We propose that changes in the export of glycine from photorespiration may explain these apparent differences in Ci*.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Performance of Base and Noble Metals for Electrocatalytic Hydrogenation of Bio-Oil-Derived Oxygenated Compounds

Electrocatalytic hydrogenation is a particularly attractive approach for converting the most unstable compounds in biogenic feedstocks at ambient conditions without external hydrogen. Here, we synthesized a variety of carbon-supported transition metal catalysts and characterized their activity for the electrocatalytic hydrogenation of a series of model compounds and pyrolysis oil. Carbonyl compounds, especially aromatic aldehydes, such as furfural and benzaldehyde, are particularly inclined to hydrogenation driven by an applied current. This was verified with pure solutions of the model compounds and with pyrolysis oil, where we achieved stable and steady continuous operation on Pd. By choosing the optimal catalyst composition, the conversion of benzaldehyde shifted from alcohol production (e.g., on Pd and Cu) to dimerization (e.g., on Co, Ni, and Zn). Pd and Cu were shown to offer the best compromise between reaction rates and efficiency although, in general, base metals offer similar conversions but better efficiencies that noble metals. Thus, the present work offers foundational results and guidelines for choosing the optimal metal catalyst and the applied potential for processing organic feedstocks as a function of its composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗
NASA NTRSDate not supplied

First N 2 Profile for Venus’ Deep Lower Atmosphere

We present the first N 2 profile for Venus’ deep lower atmosphere (<15 km) and a constrained isotopic composition for cloud N 2 [1]. This work directly addresses unresolved questions for Venus using legacy observations. Prior to this work, there were no reported measurements for N 2 abundances at <15 km and the isotopic composition remained unconstrained [2]. The N 2 parameters are critical to understanding the evolution and thermal properties of the atmosphere [3-7]. Our N 2 results were obtained by re-analysis of data acquired in 1978 by the Pioneer Venus Large Probe Neutral Mass Spectrometer (LNMS) [8]. The archived mass spectra from 64.2 to 0.2 km were treated using the analytical procedures specifically developed for the LNMS [9-11]. Judicious peak fitting permitted disambiguation of (A) N 2 + and CO + at 28 u and (B) 14 N 15 N+, 13 CO + , and C 2 H 5 + at 29 u. Quality controls included comparing the (A) LNMS CO + /CO 2 + ratios to the NIST database and literature and (B) fitted counts for CO + to the expected counts of CO + calculated from C 18 O + and 13 CO + using the LNMS 16 O/ 18 O and 12 C/ 13 C ratios (obtained from CO 2 ). Our results show that N 2 is uniformly mixed across the deep lower atmosphere between ~0.2 and 15 km (2.49 ± 0.10 v%). In contrast, N 2 is non-uniformly mixed across the sub-cloud atmosphere and clouds (~15–59 km), where N 2 abundances increase by ~ 2-fold between ~15 km (2.45 ± 0.32 v%) and ~59-51 km (5.21 ± 0.18 v%). Using the cloud data, we also obtained a constrained 15 N/ 14 N ratio (2.93×10 -3 ± 0.13×10 -3 ) and δ 15 N value (-204 ± 35‰). Thus, the LNMS results [1] suggest that (A) the atmosphere is unstable at <15 km, (B) N 2 is not well-mixed >15 km, and (C) the cloud δ 15 N falls between Earth and the solar wind [12, 13]. Comparisons of the N 2 abundances and isotopic compositions for nitrogen, carbon, and oxygen to other Venus measurements will be discussed.

deep lower atmosphere↗
NASA NTRS2020

Cycling phosphorus on the Archean Earth: Part I. Continental weathering and riverine transport of phosphorus

Phosphorus (P) is the key nutrient thought to limit primary productivity on geological timescales. Phosphate levels in Archean marine sediments are low, but quantification of the P cycle and how it changed through a billion years of recorded Archean history remain a challenge, hindering our understanding of the role played by P in biosphere/geosphere co-evolution on the early Earth. Here, we design kinetic and thermodynamic models to quantitatively assess one key component of the early P cycle – continental weathering – by considering the emergence and elevation of continents, as well as the evolution of climate, the atmosphere, and the absence of macroscopic vegetation during the Archean Eon. Our results suggest that the weathering rate of apatite, the major P-hosting mineral in the rocks, was at least five times higher in the Archean Eon than today, attributable to high levels of pCO2,g. Despite this, the weathering flux of P to the oceans was negligible in the early Archean Eon, increasing to a level comparable to or greater than the modern by the end of eon, a consequence of accelerating continental emergence. Furthermore, our thermodynamic calculations indicate high solubilities of primary and secondary P-hosting minerals in the acidic weathering fluids on land, linking to high Archean pCO2,g. Thus, weathering of P was both kinetically and thermodynamically favorable on the Archean Earth, and river water could transport high levels of dissolved P to the oceans, as also supported by the observed P-depletion in our new compilation of Archean paleosols. Lastly, we evaluated the relative rates of physical erosion and chemical weathering of silicates during the Archean Eon. The results suggest that continental weathering on the early and middle Archean Earth might have been transport-limited due to low erosion rates associated with limited subaerial emergence and low plateau elevations; by the late Archean, however, continental weathering would have transited to kinetically-limited state because of continental emergence, increased plateau elevation, and weakening weathering rates. Overall, our weathering calculations together with paleosol evidence indicate an increasing flux of bioavailable P to the oceans through time, associated with late Archean continental emergence, reaching levels comparable to or higher than modern values by the end of the eon. Increased P fluxes could have fueled increasing rates of primary production, including oxygenic photosynthesis, through time, facilitating the irreversible oxidation of the Earth’s atmosphere early in the Proterozoic Eon.

Phosphorus↗
NASA NTRS1994

Paleolimnology of the McMurdo Dry Valleys, Antarctica

The McMurdo Dry Valleys presently contain more than 20 permanent lakes and ponds, which vary markedly in character. All, with the exception of a hypersaline pond, have a perennial ice-cover. The dry valley lakes, and lakes in other ice-free regions of continental Antarctica, are unique on this planet in that they consistently maintain a thick year-round ice cover (2.8-6.0 m) over liquid water. The persistent ice covers minimize wind-generated currents and reduce light penetration, as well as restricting sediment deposition into a lake and the exchange of atmospheric gases between the water column and the atmosphere. From a paleolimnological perspective, the dry valley lakes offer an important record of catchment and environmental changes. These lakes are also modern-day equivalents of periglacial lakes that were common during glacial periods at temperate latitudes. The present lakes are mostly remnants of larger glacial lakes that occupied the valleys in the past, perhaps up to 4.6 Ma ago. Two of the valleys contain evidence of being filled with large glacial lakes within the last 10000 years. Repeated drying and filling events since then have left a characteristic impression on the salt profiles of some lakes creating a unique paleo-indicator within the water column. These events are also marked in the sediments by the concentration and dilution of certain chemical constituents, particularly salts, and are also corroborated by carbonate speciation and oxygen isotope analysis. Stratigraphic analysis of dry valley lake sediments is made difficult by the occurrence of an 'old carbon' reservoir creating spurious radiocarbon dates, and by the high degree of spatial variability in lake sedimentation. From a biological perspective, the lakes are relatively simple, containing various taxa of planktonic and benthic microorganisms, but no higher forms of life, which is an advantage to paleolimnologists because there is no bioturbation in the sediments. Useful biological paleo-indicators found in the sediments include cyanobacterial filament sheaths, diatom frustules and other eukaryotic algal cells, protozoan cysts, photosynthetic pigments, and minerals (e.g. carbonates) associated with microbial activity. Future work will benefit from fully characterizing the connection between the ice covers, environmental conditions, and paleo-indicators, thereby allowing refinement of inferences made concerning the paleoenvironment. New dating techniques need to be tested in this environment to overcome the problems associated with radiocarbon dating. The establishment of a detailed and focused paleolimnological campaign is proposed.

Non-NASA Center↗
DOE OSTI2021

Electrocatalysts Development for Hydrogen Oxidation Reaction in Alkaline Media: from Mechanism Understanding to Materials Design

Anion exchange member (AEM) fuel cells have gained great attention recently mainly owing to the advantages of non-precious metal catalysts application. Despite of the similar oxygen reduction reaction (ORR) kinetics in acid and base, the reaction kinetics of hydrogen oxidation reaction (HOR) is two orders of magnitude slower in alkaline systems than in acid. To understand the slower kinetics of HOR in base, two major interpretations have been proposed, including 1) hydrogen binding energy as a major descriptor for HOR in both acid and base; and 2) bifunctional theory based on the contributions from both hydrogen and hydroxide adsorption for HOR in alkaline electrolyte. Here, we discuss the possible HOR mechanisms in alkaline electrolytes with the corresponding Tafel behavior change. Besides the discussion of traditional Tafel-Volmer and Heyrovsky-Volmer HOR mechanisms, the recently proposed hydroxide adsorption step is discussed to illustrate the HOR mechanism difference in acid and base. We further correlate the hydroxide adsorption step involved mechanism with HOR catalyst design, and summarize its role in HOR and HER activities improvement in alkaline electrolyte over the same catalyst. The representative works of alkaline HOR catalyst development were briefly introduced, and the strategy of strengthening interaction between H2 desorption and H2O recombination was finally proposed for future HOR catalyst design in alkaline environment.

Hydrogen oxidation reaction, Hydrogen evolution re↗
DOE OSTI2022

Gas-phase metallicity determinations in nearby AGNs with SDSS-IV MaNGA: evidence of metal-poor accretion

ABSTRACT We derive the metallicity (traced by the O/H abundance) of the narrow-line region (NLR) of 108 Seyfert galaxies as well as radial metallicity gradients along their galaxy discs and of these of a matched control sample of no active galaxies. In view of that, observational data from the SDSS-IV MaNGA survey and strong emission-line calibrations taken from the literature were considered. The metallicity obtained for the NLRs was compared to the value derived from the extrapolation of the radial oxygen abundance gradient, obtained from H ii region estimates along the galaxy disc, to the central part of the host galaxies. We find that, for most of the objects ($\sim 80\, {{\ \rm per\ cent}}$), the NLR metallicity is lower than the extrapolated value, with the average difference ($\langle$D$\rangle$) between these estimates ranging from 0.16 to 0.30 dex. We suggest that $\langle$D$\rangle$ is due to the accretion of metal-poor gas to the AGN that feeds the nuclear supermassive black hole (SMBH), which is drawn from a reservoir molecular and/or neutral hydrogen around the SMBH. Additionally, we look for correlations between D and the electron density (Ne), [O iii]λ5007, and H α luminosities, extinction coefficient (AV) of the NLRs, as well as the stellar mass (M*) of the host galaxies. Evidence of an inverse correlation between the D and the parameters Ne, M*, and Av was found.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Temperature Effect on Silicate Melt-Sulfide-Metal Trace Element Partitioning in the Presence of Sulfur Under Reduced Conditions

The reduced nature of Mercury, enstatite chondrites, and the aubrite parent bodies (APB) have raised many questions regarding the geochemical behavior of typically lithophile, heat-producing, and rare-earth elements (REE) in magmas at low oxygen fugacity (fO2). Due to decreasing O availability at these low fO2, and an abundance of S(sup 2(-)), sulfur (S) acts as an important anion that changes the partitioning behavior of many elements and modifies the physical properties of silicate melts. Preliminary observations suggest that major and minor elements exhibit different geochemical affinities in highly reduced, S-rich systems compared to terrestrial rocks. The speciation and bonding environment of S, dictated by P/T/fO2 conditions, may strongly influence the degree to which S affects partitioning behavior. Here we investigate the partitioning behavior of major, minor, and trace elements between silicate melt, sulfide melt, and metal as well as the coordination chemistry of S in highly reduced silicate melts. Our work is focused on investigating solely the entropy-dependent temperature effect on partitioning of elements for which we currently have MESSENGER data (K, Na, Th, U, Si, Mg, Fe, Ti, Ca, Al, Cr, Mn, S, Cl) as well as a host of geochemically relevant trace elements such as REEs (P, Co, Ni, Mo, Ce, Nd, Sm, Eu, Gd, Dy, Yb). Previous studies in which temperature, pressure, and fO2 were co-varied found that as fO2 decreases, heat-producing elements U and Th become more chalcophile, while K becomes less chalcophile. Concurrently, nominally lithophile elements Mg and Ca become more chalcophile and were observed as minor elements in exsolved sulfides and bonded with S species in silicate melt. These studies, however, could not disentangle entropic effects from changes in the fO2. New temperature-dependent partitioning data from our work will be used to determine which elements are most likely to retain their lithophile character and hence be incorporated into silicates, and which elements are likely contained within the sulfide (chalcophile) and metal core (siderophile), setting the stage for the thermal and magmatic evolution of reduced planetary bodies.

B A Anzures↗
NASA NTRS2017

The Role of Carbon in Core Formation Under Highly Reducing Conditions With Implications for the Planet Mercury

Results from the MErcury Surface, Space ENvironment, GEochemistry and Ranging (MESSENGER) spacecraft have shown elevated abundances of carbon on the surface of Mercury. Furthermore, the X-Ray Spectrometer on board MESSENGER measured elevated abundances of sulfur and low abundances of iron, suggesting the planet's oxygen fugacity (fO2) is several log10 units below the Iron-Wüstite (IW) buffer. Similar to the role of other volatiles (e.g. sulfur) on highly reducing planetary bodies, carbon is expected to behave differently than it would under higher fO2. As discussed by Nittler et al. and Hauck et al., under such highly reducing conditions, the majority of the iron partitions into the core. On Mercury, this resulted in a relatively large core and a thin mantle. Using a composition similar to the largest volcanic field on the planet (the northern volcanic plains), Vander Kaaden and McCubbin conducted sink-float experiments to determine the density of melts and minerals on Mercury. They showed that graphite would be the only buoyant mineral in a mercurian magma ocean. Therefore, Vander Kaaden and McCubbin proposed a possible primary flotation crust on the planet composed of graphite. Concurrently, Peplowski et al. used GRS data from MESSENGER to show an average northern hemisphere abundance of C on the planet of 1.4 +/- 0.9 wt%. However, as this result was only at the one-sigma detection limit, possible carbon abundances at the three-sigma detection limit for Mercury range from 0 to 4.1 wt% carbon. Additionally, Murchie et al. investigated the possible darkening agent on Mercury and concluded that coarse-grained graphite could darken high reflectance plains to the low reflectance material. To further test the possibility of elevated abundances of carbon in Mercury's crust, Peplowski et al. used the low-altitude MESSENGER data to show that carbon is the only material consistent with both the visible to near-infrared spectra and the neutron measurements of low reflectance material on Mercury, confirming that C is the primary darkening agent on Mercury. Confirmation of carbon on the planet prompts many questions regarding the role of carbon during the differentiation and evolution of Mercury. Given the elevated abundances of both S and C on Mercury's surface, it begs the question, what is the core composition of the planet? This study seeks to understand the impact of C as a light element on potential core compositions on Mercury.

Vander Kaaden, Kathleen E..↗
DOE OSTI2023

Computationally Accelerated Discovery and Experimental Demonstration of High-Performance Materials for Advanced Solar Thermochemical Hydrogen Production

This project achieved its overarching goal of accelerating the discovery and validation of solar thermochemical hydrogen (STCH) materials through a tightly integrated approach that combined high-throughput computational screening, advanced machine learning (ML), and experimental testing. Guided by the objectives outlined in the Statement of Project Objectives (SOPO), our work fulfilled all major milestones across four technical tasks and delivered scientific breakthroughs and practical tools that significantly exceeded the original scope of the project. We began by addressing the challenge of predicting material phase stability through machine learning. A novel Python module was developed to generate thousands of meaningful features from composition, structure, and electronic properties, enabling rapid and reproducible ML model development. Using these tools, we trained a model to predict temperature-dependent Gibbs energies (G(T)) for inorganic crystalline materials with near-chemical accuracy—roughly 40 meV/atom—marking the first such descriptor of its kind. We also introduced a new machine-learned tolerance factor, τ, that accurately predicted perovskite formability with over 90% success, outperforming traditional heuristic models, such as the Goldschmidt tolerance factor. These capabilities allowed for rapid and accurate predictions of phase stability across a vast oxide composition space, setting the stage for high-throughput thermodynamic screening. Building on this foundation, we conducted an extensive computational screening of candidate STCH oxide materials. Over 1.1 million perovskite compositions were evaluated using the τ descriptor, leading to the identification of more than 27,000 predicted stable structures. Using density functional theory (DFT), we refined over 68,000 multinary perovskite structures and computed oxygen vacancy formation energies for over 1,300 ternary and double perovskites. These calculations enabled us to isolate compounds with redox behavior consistent with STCH requirements and resulted in a public dataset now hosted on the Materials Project. Recognizing that thermodynamic screening alone is insufficient, we addressed kinetic limitations by developing a suite of tools to estimate transition state (TS) energies for key redox reactions. We implemented a novel bounding approach that provides lower and upper estimates of TS energies with dramatically reduced computational cost, requiring less than 10% of the CPU time of a full nudged elastic band (NEB) calculation while maintaining high accuracy. This enabled rapid evaluation of over 200 reaction pathways across 90 materials. To further accelerate screening, we developed a SISSO-based ML model to predict diffusion barriers with a 96.7% success rate in classifying fast vs. slow materials, supporting a robust, data-driven framework for assessing redox kinetics. Experimental validation was critical to confirming the predictive power of our models. We synthesized and tested a wide array of candidate materials, including Mn-doped hercynite and several Gd- and La-based perovskites. Notably, Sr 0.4 Gd 0.6 Mn 0.6 Al 0.4 O 3 (SGMA) and Gd 0.5 La 0.5 Co 0.5 Fe 0.5 O 3 (GLCF) emerged as leading STCH materials, exhibiting robust redox cycling and high hydrogen yields exceeding 150 µmol H 2 /g per cycle. These materials also retained over 50% of their hydrogen productivity under high-conversion conditions (H 2 O:H 2 = 1333:1), demonstrating strong thermodynamic favorability and promising performance under industrially relevant scenarios. Additional candidates, such as La 2 MnNiO 6 (L2MN), were found to produce even higher yields than ceria under standard STCH conditions. Our collaborators at Sandia National Laboratories confirmed these findings using high-temperature X-ray diffraction and thermogravimetric analysis, observing stable phase evolution and reversible redox activity. In several respects, the project went beyond the goals initially outlined in the SOPO. We published 17 peer-reviewed articles, including a large dataset of over 66,000 theoretical perovskites and a new structure prediction method (SPuDS-DFT) that accurately identifies ground-state structures at a fraction of the cost of traditional DFT. We demonstrated that our machine-learned G(T) model offers accuracy rivaling quasiharmonic calculations while being orders of magnitude faster. In partnership with the Materials Project, we made our datasets openly available, providing a powerful new resource for the broader materials science community. The combined computational and experimental advances of this project represent a significant advance in STCH materials discovery. By creating a robust, generalizable, and open workflow for thermodynamic and kinetic screening, and validating key findings through synthesis and reactor testing, we have provided a practical and scalable pathway for the rapid identification of new redox-active materials. The tools, data, and materials developed under this project are already supporting ongoing research and have laid the groundwork for the next generation of solar fuel technologies.

08 HYDROGEN↗
DOE OSTI2021

Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus-like Atmosphere

Exoplanets orbiting M dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our solar system, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO 2 -dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O 2 + ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. Furthermore, the extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Effects of NaCl on metabolic heat evolution rates by barley roots

The effect of salinity stress on metabolic heat output of barley (Hordeum vulgare L.) root tips was measured by isothermal microcalorimetry. Several varieties differing in tolerance to salinity were compared and differences quantified. Two levels of inhibition by increasing salt were found. Following the transition from the initial rate of the first level, inhibition remained at about 50% with further increases in salt concentration up to 150 millimolar. The concentration of salt required to inhibit to this level was cultivar dependent. At highter concentrations (>150 millimolar) of salt, metabolism was further decreased. This decrease was not cultivar dependent. The decreased rate of metabolic heat output at the first transition could be correlated with decreases in uptake of NO3-, NH4+, and Pi that occurred as the salt concentration was increased. The high degree of dependence of the inhibition of metabolic heat output on NaCl concentration points to a highly cooperative reaction responsible for the general inhibition of metabolism and nutrient uptake. The time required to attain the first level of salt inhibition is less than 20 minutes. Inhibition of root tips was not reversible by washing with salt free solutions. In addition to revealing these features of salt inhibition, isothermal microcalorimetry is a promising method for convenient and rapid determination of varietal differences in response to increasing salinity.

NASA Discipline Life Support Systems↗
DOE OSTI2023

Topological Bismuth (1 1̅ 0) Facet for Efficient Oxygen Reduction Cathode in Fuel Cells

The performance of proton-exchange membrane fuel cells is strongly dependent on the efficiency of the electrochemical oxygen reduction reaction (ORR). Here, we report both two-dimensional (4 monolayer, ML) and three-dimensional Bi (16 ML) (11̅0) facets deliver a strong correlation between nontrivial topological surface states and low electrochemical overpotential toward ORR. Both topological Bi (11̅0) film and slab have two slightly gapped (9–45 meV) Dirac cones with the invariant Z 2 =1, derived from the parity products for the occupied states at four time-reversal invariant momentum points. By contrast, the Bi (110) film is stabilized as an ordinary insulator due to atomic buckling. The much lower thermodynamic ORR overpotential (η = 0.58 V) and higher H 2 O selectivity by four steps of H + /e – transfer are obtained on the topological Bi (11̅0) surface compared with the nontopological (110) facet (η = 0.84 V). This promotional effect can be attributed to the accumulation of topological surface bands (Bi 6p orbitals) near the Fermi energy, which provide stronger electronic transport channels along the surface Bi–Bi and Bi–O bonds. Finally, a 16 ML Bi (11̅0) slab is also constructed to simulate the topological (11̅0) surface of three-dimensional bulk Bi as a topological crystalline insulator (TCI). Furthermore, the consistent topological properties and slightly lower thermodynamic ORR overpotential (η = 0.53 V) are also observed on the thicker slab.

25 ENERGY STORAGE↗
DOE OSTI2020

Oxidation of the Interiors of Carbide Exoplanets

Astrophysical measurements have shown that some stars have sufficiently high carbon-to-oxygen ratios such that the planets they host would be mainly composed of carbides instead of silicates. We studied the behavior of silicon carbide in the presence of water under the high pressure–temperature conditions relevant to planetary interiors in the laser-heated diamond-anvil cell. When reacting with water, silicon carbide converts to silica (stishovite) and diamond at pressures up to 50 GPa and temperatures up to 2500 K: $\mathrm{SiC}+2{{\rm{H}}}_{2}{\rm{O}}\to {\mathrm{SiO}}_{2}+{\rm{C}}+2{{\rm{H}}}_{2}$. Therefore, if water can be incorporated into carbide planets during their formation or through later delivery, they could be oxidized and have mineralogy dominated by silicates and diamond in their interiors. The reaction could produce CH 4 at shallower depths and H 2 at greater depths that could be degassed from the interior, causing the atmospheres of the converted carbon planets to be rich in reducing gases. Excess water after the reaction can be stored in dense silica polymorphs in the interiors of the converted carbon planets. Such conversion of mineralogy to diamond and silicates would decrease the density of carbon-rich planet, making the converted planets distinct from silicate planets in mass–radius relations for the 2–8 Earth mass range.

79 ASTRONOMY AND ASTROPHYSICS↗