Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interface & surface thermodynamics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nanoionics Drastically Accelerating Mass Transfer at Elevated Temperatures over 750 °C

Nanoionics were previously considered thermally unstable and infeasible for devices operating above 500 °C. Here, we elucidate the design principle for establishing stable nanoionics from various oxides. We utilized reversible solid oxide cells (SOCs) as the test bed and implemented nanoionics using atomic layer deposition (ALD). We demonstrate a straightforward, interface-controlled, practical approach to render a conformal, ∼15 nm thick ALD film, which initially thermodynamically favors the formation of a solid solution with the substrate into surface nanoionics with single or double layers of nanograins with random crystal orientations. The nanoionics exhibited conductivity estimated to be 7 orders of magnitude higher than that of their bulk-scale counterpart. They demonstrated conformability with uniform grain sizes of ∼15 nm, even after electrochemical operation for ∼500 h at 750 °C and 1000 h at 850 °C. The thermal stability and conductivity of such nanoionics represent a conceptual and technological framework in nanoionics.

atomic layer deposition↗
DOE OSTI2025

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗
DOE OSTI2025

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2025

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2025

Facet-Dependent Cold Welding of Au Nanorods Revealed by Liquid Cell Transmission Electron Microscopy

Cold welding of metals at the nanoscale has been demonstrated to play a significant role in bottom-up manufacturing and self-healing processes of nanostructures and nanodevices. However, the welding mechanism at the nanoscale is not well understood. In this study, a comprehensive demonstration of the cold welding process of gold nanorods with different modes is presented through in situ liquid cell transmission electron microscopy. The experimental results and molecular dynamics simulations reveal that the nanorods are welded through the facet-dependent atomic surface diffusion and rearrangement along {100} facets. The density functional theory calculations indicate that the preferred coalescence of two {100} surfaces is thermodynamically favorable. Unlike the prevalent “oriented attachment” in the nanoparticle coalescence, the misalignment of nanorod orientations and local stresses can induce grain boundaries and stacking faults in the welded interface.

42 ENGINEERING↗
DOE OSTI2025

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗
DOE OSTI2026

The Role of Cooperative Interactions Among Surfaces, Solvents, and Reactive Intermediates on Catalysis at Liquid–Solid Interfaces (Final Report DE-SC0020224)

This project established quantitative links between inner‑sphere chemistry (active metal identity, coordination, and zeolite topology) and outer‑sphere organization (solvent identity, hydrogen‑bond networks, and pore condensation) that govern rates, activation barriers, and selectivities for alkene epoxidation and epoxide ring‑opening at solid–liquid and quasi‑liquid–solid interfaces. We deconvoluted contributions from covalent interactions at active sites and noncovalent, solvent‑mediated interactions within pores by pairing well‑defined metal substituted zeolites with controlled solvent environments. We then mapped those contributions onto measurable kinetics (ΔH‡, ΔS‡), adsorption thermodynamics (ITC), and in situ spectroscopy. The transferrable outcomes include a set of design rules that include the following understandings. First, tune silanol ((SiOH)x) density and pore topology to organize solvent networks that selectively stabilize transition states. Second, exploit activity‑coefficient‑normalized rates and adsorption– barrier correlations to diagnose when solvent reorganization rather than surface chemistry limits performance. Third, use partial pore condensation (e.g., acetonitrile, water but also generalizable to other solvents) to elicit liquid‑like stabilization effects even in nominally vapor‑phase reactors. Collectively, these results provide strategies to increase epoxidation rates, improve oxidant utilization (i.e., selectivities), and steer regioselectivity in zeolite‑based catalytic processes relevant to sustainable oxidation chemistry. These outcomes should be transferable to other classes of reactions that proceed in microporous materials and under confinement provided by organized solvents (e.g., electrochemical double layers).

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

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2025

Phase Separation, Capillarity, and Odd-Surface Flows in Chiral Active Matter

Active phase separations evade canonical thermodynamic descriptions and have thus challenged our understanding of coexistence and interfacial phenomena. Considerable progress has been made towards a nonequilibrium theoretical description of these traditionally thermodynamic concepts. Spatial parity symmetry is conspicuously assumed in much of this progress, despite the ubiquity of chirality in experimentally realized systems. Here, in this Letter, we derive a theory for the phase coexistence and interfacial fluctuations of a system that microscopically violates spatial parity. We find suppression of the phase separation as chirality is increased as well as the development of steady-state currents tangential to the interface dividing the phases. These odd flows are irrelevant to stationary interfacial properties, with stability, capillary fluctuations, and surface area minimization determined entirely by the capillary surface tension. Using large-scale Brownian dynamics simulations, we find excellent agreement with our theoretical scaling predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗
DOE OSTI2025

Perovskite-driven solar C 2 hydrocarbon synthesis from CO 2

Photoelectrochemistry (PEC) presents a direct pathway to solar fuel synthesis by integrating light absorption and catalysis into compact electrodes. Yet, PEC hydrocarbon production remains elusive due to high catalytic overpotentials and insufficient semiconductor photovoltage. Here we demonstrate ethane and ethylene synthesis by interfacing lead halide perovskite photoabsorbers with suitable copper nanoflower electrocatalysts. The resulting perovskite photocathodes attain a 9.8% Faradaic yield towards C 2 hydrocarbon production at 0 V against the reversible hydrogen electrode. The catalyst and perovskite geometric surface areas strongly influence C 2 photocathode selectivity, which indicates a role of local current density in product distribution. The thermodynamic limitations of water oxidation are overcome by coupling the photocathodes to Si nanowire photoanodes for glycerol oxidation. These unassisted perovskite–silicon PEC devices attain partial C 2 hydrocarbon photocurrent densities of 155 µA cm -2 , 200-fold higher than conventional perovskite–BiVO 4 artificial leaves for water and CO 2 splitting. These insights establish perovskite semiconductors as a versatile platform towards PEC multicarbon synthesis.

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

Activation of propane on Ag–PdO(101) model surfaces

Oxidation of alkanes remains a central challenge in catalysis due to the high activation barriers of C–H bonds and the thermodynamic favorability of complete oxidation. Palladium oxide (PdO), particularly its (101) facet, is known for its high reactivity in alkane oxidation, which is attributed to its coordinatively unsaturated palladium (Pd) and O atoms. In this study, we investigate the effect of silver (Ag) incorporation on the oxidation behavior of propane over PdO(101) using temperature-programmed reaction spectroscopy (TPRS) under controlled conditions. While pristine PdO(101) exhibits complete oxidation of propane with CO₂ and H₂O desorption at high temperatures (approximately 475 K), Ag incorporation induces a new CO₂ desorption peak at significantly lower temperatures (approximately 330 K). This shift is attributed to the formation of new active sites at the Ag–PdO(101) interface. Quantitative analysis reveals that low-temperature activity correlates with Ag coverage, while overall CO₂ production decreases, suggesting a redistribution of reactivity rather than an increase in active surface area. Activation energy estimations using the Redhead method confirm that C–H bond activation becomes more facile at the interface, with a 46 kJ/mol reduction compared to pristine PdO(101). These findings demonstrate that incorporating a less reactive metal such as Ag into PdO surfaces not only modifies the reaction energetics but also enables the design of bimetallic catalysts with improved selectivity for partial oxidation reactions.

Chemistry↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

Phase selecting nearly degenerate crystalline polymorphs during epitaxial growth can be challenging yet critical to targeting physical properties for specific applications. Here, we establish how phase selectivity of altermagnetic and noncentrosymmetric polymorphs of MnTe can be programmed by subtle changes to the surface of lattice-matched InP substrates in molecular beam epitaxy growth. Bulk altermagnetic MnTe is thermodynamically stable in the hexagonal NiAs-structure and is synthesized here on the polar (111)A surface (In-terminated) of InP, while the noncentrosymmetric, cubic ZnS-structure with wide band gap (>3 eV), which epitaxially matches III–V materials, is stabilized on the (111)B surface (P-terminated). Electron microscopy, X-ray photoemission spectroscopy, and reflection high-energy electron diffraction indicate that phase selection is triggered at the interface and proceeds along the growing surface. First-principles calculations suggest that interfacial termination and strain have a significant effect on the interfacial energy; stabilizing the NiAs polymorph on the In-terminated surface and the ZnS structure on the P-terminated surface. Here, selectively grown, high-quality, phase pure films of both MnTe polymorphs will enable our understanding of the novel properties of these materials, thereby facilitating their use in new applications ranging from spintronics to microelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The Role of Interfaces and Charge for Chemical Reactivity in Microdroplets

A wide variety of reactions are reported to be dramatically accelerated in aqueous microdroplets, making them a promising platform for environmentally clean chemical synthesis. However, to fully utilize the microdroplets for accelerating chemical reactions requires a fundamental understanding of how microdroplet chemistry differs from that of a homogeneous phase. Here, in this work, we provide our perspective on recent progress to this end, both experimentally and theoretically. We begin by reviewing the many ways in which microdroplets can be prepared, creating water/hydrophobic interfaces that have been frequently implicated in microdroplet reactivity due to preferential surface adsorption of solutes, persistent electric fields, and their acidity or basicity. These features of the interface interplay with specific mechanisms proposed for microdroplet reactivity, including partial solvation, possible gas phase channels, and the presence of highly reactive intermediates. We especially highlight the role of droplet charge and associated electric fields, which appears to be key to understanding how certain reactions, like the formation of hydrogen peroxide and reduced transition metal complexes, are thermodynamically possible in microdroplets. Lastly, we emphasize opportunities for theoretical advances and suggest experiments that would greatly enhance our understanding of this fascinating subject.

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

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗