Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2021

Influence of fine-scale B2 precipitation on dynamic compression and wear properties in hypo-eutectic Al0.5CoCrFeNi high-entropy alloy

The near equiatomic high-entropy alloys based on Al, Co, Cr, Fe, Ni have provided a novel microstructural template for alloy design. The face-centered cubic (FCC) matrix strengthened by body-centered cubic (BCC) based ordered B2 precipitates can be used to develop a new generation of precipitation strengthened alloys. The present study focuses on an Al0.5CoCrFeNi hypo-eutectic alloy consisting of an as-solidified FCC+B2 microstructure. Guided by solution thermodynamic modeling, this alloy has been isothermally aged for solid-state precipitation of a high fraction of fine-scale intragranular B2 particles. The dynamic compressive strength, hardness, and sliding wear response of the alloy as a function of B2 phase fraction is examined. The flow stress at 0.02% true strain increased from ~ 670 MPa to ~1350 (102%) MPa, the ultimate compressive strength increased from ~1160 MPa to 1500 MPa (20%), and the wear resistance increased more than 5 times with progressively increasing phase fraction of B2 precipitates.

Gwalani, Bharat↗
DOE OSTI2022

Spark plasma sintered, MoNbTi-based multi-principal element alloys with Cr, V, and Zr

MoNbTi, MoNbTiZr, CrMoNbTiZr, and MoNbTiVZr multi-principal element alloys (MPEAs) were fabricated via spark plasma sintering (SPS) and investigated for use in high-strength applications. The fabrication method by SPS and powder metallurgy differs from those presented in the prior literature, where most MPEAs are fabricated using arc melting (AM) methods. Here cryogenic milling was used to maximize potential defect sinks (grain boundaries) for radiation resistance and to increase ultimate tensile strength through the Hall-Petch effect. SPS was chosen for consolidation in order to maintain a fine-grained structure during densification. Each alloy was characterized using x-ray diffraction and scanning electron microscopy for phase identification and compositional homogeneity. The base ternary alloy MoNbTi presented a predominantly single BCC system with minor cubic phases. The introduction of additional alloying elements—Zr, V, and Cr—heightened the phase complexity and increased the fractions of a secondary BCC phase and an HCP phase from Zr. The addition of Cr induced a larger fraction of the Laves phase to form. The addition of V caused the precipitation of small Mo inclusions. Thermodynamic analysis was performed to understand the separation of phases in each alloy. Discrepancies among the phase predictions generated by thermodynamic models, phases previously presented in the literature, and the characterization results suggest that MPEA fabrication methods, especially solid-state methods, require significant investigation to ensure that alloys can remain stable throughout their anticipated service lifetimes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effect of Zr on precipitation in oxide dispersion strengthened FeCrAl alloys

Oxide dispersion strengthened (ODS) FeCrAl alloys are promising candidate materials for advanced nuclear reactor applications requiring high-temperature strength, corrosion resistance, and irradiation tolerance. As these alloys have increased in compositional complexity through attempts to use highly reactive elements such as Zr to refine particle sizes and optimize nanoprecipitate dispersion characteristics, much debate has ensued as to the effects of these alloying element additions on alloy properties. In an attempt to reconcile differences in nanoprecipitate distributions reported in the literature over the past decade, a detailed investigation of a recently developed ODS FeCrAl alloy with nominal composition Fe–10Cr-6.1Al-0.3Zr+0.3Y 2 O 3 is presented here using a combination of atom probe tomography (APT), scanning/transmission electron microscopy (S/TEM), and computational thermodynamics modeling. It is illustrated that based on the amount of Zr available in the lattice, Zr competes with Al and Cr to form carbides and nitrides as opposed to oxygen-rich precipitates. This alloy system has a high number density (>10 23 m -3 ) of ~2–4 nm diameter (Y,Al,O)-rich nanoprecipitates, but it is shown that due to the compositional spread and unknown partitioning of Al between the matrix and precipitates, significant challenges still exist for quantifying the exact compositions of these precipitates using APT. However, the noted compositional spread is supported by identified complex oxides yttrium aluminum monoclinic (YAM) and yttrium aluminum garnet (YAG) using S/TEM. As a result of these findings, researchers developing ODS FeCrAl with reactive element additions must pay careful attention to C and N impurities when optimizing reactive element additions.

electron microscopy↗
DOE OSTI2022

On the high-temperature stability of the Al 8 Cu 3 Ce intermetallic in an additively manufactured Al-Cu-Ce-Zr alloy

High-temperature resistant eutectic Al alloys are crucial materials for lightweight and energy efficient design in the automotive and aviation industries. Additive manufacturing offers a pathway to refine eutectic microstructures and develop novel alloys with superior high-temperature strength. High-volume fraction intermetallic Al-Cu-Ce alloys have been developed to deliver high-temperature strength in combination with reduced hot-tearing susceptibility. Zr is added to provide additional strengthening via nanoscale Al 3 Zr precipitation, and to stabilize and avoid coarsening of the Al 8 Cu 3 Ce phase. However, the detailed interaction between Zr and Al 8 Cu 3 Ce remains unexplored. In this work, we show with synchrotron X-ray diffraction that laser powder bed fusion fabricated Al-Cu-Ce and Al-Cu-Ce-Zr alloys contain predominantly the Al 8 Cu 3 Ce intermetallic in the as-fabricated condition. Heat treatment of the Al-Cu-Ce alloy results in the Al 8 Cu 3 Ce → Al 8 Cu 4 Ce phase transformation. In the Al-Cu-Ce-Zr alloy, minor fractions of (Al,Cu,Si) 4 Ce and Al 2 Cu-θ are found in the as-fabricated condition, while Al 8 Cu 3 Ce remains stable during heat treatment. Atom probe microscopy quantifies intermetallic stoichiometries and reveals how Zr is enriched at the Al-matrix/Al 8 Cu 3 Ce interface acting as a diffusion barrier against solute exchange. Calibrated thermodynamic modeling underpins this as a kinetic effect. Here, a qualitative microstructural model summarizes, how Zr stabilizes Al 8 Cu 3 Ce against phase transformations and coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗
DOE OSTI2021

Thermodynamic analysis of hydrogel swelling in aqueous sodium chloride solutions

Here, a comprehensive thermodynamic model is presented for swelling and salt partitioning of poly (n-isopropyl acrylamide) hydrogels in aqueous sodium chloride (NaCl) solutions at 298 K. Together with a Helmholtz energy expression for network elastic energy, a modified electrolyte Nonradom Two-Liquid activity coefficient model is used to represent both the long-range electrostatic interactions and the short-range van der Waals interactions present in the gel phase. The model parameters include one network parameter per hydrogel system plus two binary parameters per interaction pair in the polymer–solvent-salt systems. With semiquantitave agreement with experimental data for both hydrogel swelling and salt partitioning up to salt saturation, the analysis suggests the ion hydration between NaCl and water is the root cause for hydrogel deswelling in aqueous NaCl solution at high NaCl concentrations, i.e., NaCl weight fraction >0.03. On the other hand, the strong repulsive interaction between NaCl and hydrogel polymer is the dominant factor for salt partitioning.

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

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis of pumped thermal energy storage using particle media integrated with concentrating solar power

Pumped Thermal Energy Storage (PTES) is an electricity storage system that converts electricity into thermal energy which is stored and later transformed back into electricity. Previous work has illustrated that particles have low capital costs and can be operated over a wide range of temperatures. PTES with particle storage achieves higher round-trip efficiency and specific power output than when molten salt thermal energy storage is used. This article explores hybrid systems that combine PTES with Concentrating Solar Power (CSP). Hybrid systems share the majority of components thereby reducing costs compared to two stand-alone devices. In addition, hybrid systems can provide multiple services (such as renewable power generation and electricity storage services). Using particle thermal storage in these hybrid concepts provides freedom in choosing the design conditions, since a wide range of operating temperatures is allowable, therefore making it possible to identify hybrid system designs that have good performance. In this article, two concepts for hybrid PTES-CSP are introduced. Thermodynamic models are developed and these are used to evaluate the performance of two hybrid systems. These models account for turbomachinery efficiency, and approach temperature and pressure loss in heat exchangers, as well as other sources of inefficiency, such as motor-generator losses, and air fan power. The “Solar Top-Up” Concept uses CSP to increase the temperature delivered by the charging heat pump. The discharging system uses a topping gas cycle and a bottoming steam cycle to fully exploit the available energy. Using solar heat to increase the maximum temperature from 750 K to 1100 K increases the round-trip efficiency from 41 % to 62 % and the specific work output from 88 kJ/kg to 301 kJ/kg. The second concept is a “Dual-Mode” device which provides both electricity storage and solar electricity generation with the same set of components. The PTES-mode and CSP-mode performance are optimized at different design values but careful parameter selection leads to good performance of both modes: one design produces PTES round-trip efficiency > 60 %, CSP heat engine efficiency > 40 %, and specific work outputs > 150 kJ/kg (for both cycles), when the particle receiver temperature is > 1200 K and maximum heat pump temperature is 1100 K.

14 SOLAR ENERGY↗
DOE OSTI2025

Facet-Dependent Doping and Dopant-Dependent Faceting in Si-Doped GaAsSb Nanowires

Here, we correlate the spatial distributions of Si, Sb, and rotational twins in Si-doped GaAs 1-x Sb x nanowires. GaAs 1-x Sb x nanowires were grown epitaxially on Si(111) substrates by tuning process conditions to achieve repeated nucleation of rotational twins and growth along the [111]B direction; dilute Sb and Si fluxes were chosen to create a sufficient twin density to achieve high yield while avoiding growth of the wurtzite phase. While the impact of Si and Sb on twin density and nanowire growth rate has been previously reported, the facet-dependent incorporation of these species has not been established. Scanning transmission electron microscopy was used to confirm that Sb incorporates preferentially on the (111)B facets relative to {1 ̅1 ̅0} facets prior to nucleation of a rotational twin. With periodic twinning, this facet dependence leads to alternating regions of enriched and depleted Sb concentration attributed to a growth rate-dependent Sb-As-exchange mechanism. Atom probe tomography measurements establish that while Si doping is not perturbed by twinning on (111)B facets, Si and Sb concentrations are anti-correlated for growth on non-(111)B facets. Density functional theory calculations underpin a thermodynamic model that explains the observed anisotropies in dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

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

Shape Control of Colloidal Semiconductor Nanocrystals through Thermodynamically Driven Aggregative Growth

The optoelectronic properties of colloidal semiconductor nanocrystals (NCs) can be manipulated by changing their geometric shapes. The precise synthetic control over particle morphologies, however, has remained elusive. Conventional growth techniques rely on the kinetic assembly of atomic units, where supersaturation and precipitation processes can lead to a broad distribution of particle shapes. In this paper, we demonstrate that replacing atomic precursors with small-size nanocrystals as building blocks for larger colloids offers an easier, more predictive control over nanoparticle shape evolution. The reported growth strategy is illustrated via shape-selective syntheses of CdSe and CdS NC cubes, spheres, rods, as well as unprecedented “donut” and ring-like structures. Growth of desired particle morphologies was achieved by choosing an appropriate reaction solvent that minimizes the surface free energy through a unique shape formation. This so-called "aggregative growth mechanism" is explained using a thermodynamic model for interacting viscous colloids. We expect that this novel growth strategy will offer a viable platform for shape–selective synthesis of many inorganic colloids, providing pathways to optimizable optoelectronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nitrilotriacetamides: Nitric Acid and Water Extraction and Physicochemical Properties

Nitrilotriacetamides are an emerging class of extractants that have been identified as promising candidates for a variety of separation schemes. Systems incorporating these ligands are prone to third-phase formation at high loading, but this can be mitigated by increasing the alkyl chain length or adding a phase modifier (e.g., 1-octanol). In this work, the impacts of these modifications on HNO 3 extraction, H 2 O extraction, and physiochemical properties (interfacial tension, viscosity, and density) are examined. The studies show that ligand alkyl chain length has very little impact on HNO 3 and H 2 O extraction, whereas the addition of 10% octanol significantly increases extraction of both HNO 3 and H 2 O. Viscosity and interfacial ten-sion are impacted by HNO 3 concentration, alkyl chain length, and the presence of phase modifier but remain in a range that should be compatible with commonly used solvent contacting equipment. Additional studies of physiochemical properties at variable ligand concentration suggest that a change in aggregation occurs above approximately 10-50 mM ligand. Thermodynamic modeling of distribution ratios and spectroscopic studies aimed at understanding speciation in the presence of HNO 3 are also reported. Further, these experiments suggest that the first equivalent of extracted HNO 3 interacts with both the central amine and carbonyl groups of the ligand. The spectroscopic changes at higher HNO 3 concentrations are relatively minor, suggesting additional HNO 3 is associated in an outer-sphere manner

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

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-exchange Resin at Initial pH ~3 – 7 and Temperature 25 – 55 °C

This work presents systematic aqueous-phase adsorption equilibrium, kinetics, and column breakthrough measurements with three key biointermediates that are common compounds in many bioprocesses. Adsorption equilibrium experiments were carried out with acetic acid, butyric acid, and lactic acid on a commercial ion-exchange resin, Amberlite IRN-78, at wide ranges of acid concentration (8-500 mmol/L), initial pH (similar to 3-7), and temperature (25-55 degrees C), simulating the effluent characteristics from different fermenter operations. The kinetics and column breakthrough experiments were conducted at an initial pH of 6 and a concentration of 200 mmol/L. The equilibrium study shows a higher loading at the initial pH < pK(a) and a lower loading at the initial pH > pK(a). Overall removal varies between 16 and 99% depending on the initial pH, temperature, and organic acid concentration and type. The study further indicates monolayer adsorption at the equilibrium pH > 10 and multilayer adsorption at the equilibrium pH < 6. The thermodynamic modeling of adsorption isotherm data was carried out using Langmuir and Freundlich isotherms. IRN-78 presents fast adsorption kinetics as the maximum loading was attained in <= 10 min and nearly the same breakthrough time for all three organic acids involved in this study.

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

Complex Fluorine Chemical Potential Effects on the Shape and Compositional Heterogeneity of KTa 1– x Nb x O 3 Nanoparticles

When kinetic conditions dominate and dictate the growth in a nanoparticle synthesis, properties of the synthesis environment can have considerable effects on the properties of the products. Such effects were studied here, where the solution environment was changed via the addition of KF to the hydrothermal syntheses of KTa 1– x Nb x O 3 and KTaO 3 . One result demonstrated the straightforward cause-and-effect relationship between the solution and reaction kinetics: KF directly increased the solution stability of the Ta species and therefore decreased its reaction rate, resulting in a change in composition heterogeneity of Ta and Nb in the KTa 1– x Nb x O 3 particles. However, not all effects are so simple; changing the chemical potential of the solution with KF can also promote the formation of particles with anisotropic defect enhanced kinetic Wulff shapes instead of cuboidal shapes. The increased F chemical potential in the solution enabled the formation of planar defects in the bulk, which accelerated growth in-plane to form particles characterized by flat rectangular flake geometries. In this work, thermodynamic modeling with density functional theory calculations confirmed that sufficient KF concentrations can drive the formation of a defect phase K n +1 Ta n O 3 n F.

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

Simultaneous Absorption and Desorption Isotherms of Various Hydrogen and Deuterium Mixtures between 20 and 120 °C Are Used to Determine the Activity Coefficients for Palladium-Hydride Solutions

A test bed was constructed to measure the absorption and desorption isotherms for palladium hydride using H 2 , D 2 , and various H 2 /D 2 mixtures for temperatures in the range 20 °C ≤ T ≤ 120 °C. The pressure–composition–temperature isotherms were measured. The pressures obtained with mixtures between each pure isotope were monotonic, yet nonlinear. This nonlinear dependence of total pressure with feed gas protium concentration reveals the mixed isotope hydride system behaves nonideally. A thermodynamic model was adapted from the literature, which accounts for the nonideal nature of the mixed-isotope system. The measured data were used as constraints in this model in order to calculate the protium mole fractions in the hydride phase and the protium activity coefficients for palladium hydride at various temperatures and protium concentrations in the system. Knowing the protium mole fraction and activity coefficients allows for a priori calculation of the isotopologue distribution in the gas phase and the isotope distribution in the hydride phase, given a palladium temperature and equilibrium pressure.

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

Ion Pairing and Redissociaton in Low-Permittivity Electrolytes for Multivalent Battery Applications

Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.

25 ENERGY STORAGE↗
DOE OSTI2025

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗