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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Morphological Evolution and Dealloying During Corrosion of Ni20Cr (wt.%) in Molten FLiNaK Salts

The dealloying corrosion behavior of the FCC Ni20Cr (wt%) in molten LiF-NaF-KF (FLiNaK) salts at 600 °C under varying applied potentials was investigated. Using in-operando electrochemical techniques and a multi-modal suite of characterization methods, we connect electrochemical potential, thermodynamic stability, and electro-dissolution kinetics to the corrosion morphologies. Notably, under certain potential regimes, a micron-scale bicontinuous structure, characterized by a network of interconnected pores and ligaments riched with the composition of the more noble (MN) element, becomes prominent. At other potentials both MN and less noble (LN) elements dealloy but at different rates. The dealloying process consists of lattice and grain boundary diffusion of Cr to the metal/salt interface, interphase Cr oxidation, accompanied by surface diffusion of Ni to form interconnected ligaments. At higher potentials, the bicontinuous porous structure undergoes further surface coarsening. Concurrently, Cr(II), Cr(III), and Ni(II) begin to dissolve, with the dissolution of Ni occurring at a significantly slower rate. When solid-state transport of Cr is exceeded by the interfacial rates, dealloying depths are limited.

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

Primary biomolecular adsorption energetics of core–shell nanocomplexes: Implications for biological interactions

Molecular organization at the nano-bio interface governing the colloidal stability, reactivity, immune recognition, and drug delivery performance of nanoparticles remains difficult to predict. Quantifying the primary hydration energetics of biomolecule-coated nanomaterials can determine those interactions and provide a basis for engineered nanocarriers with tailored behavior in biological systems. Here, we measured the thermodynamics of water adsorption on patchy dry magnetite (Fe 3 O 4 ) nanoparticles coated with three model biomolecules, bovine serum albumin, potato starch, and lauric acid and compared these properties to the hydration energetics of the corresponding free dry biomolecules. The results demonstrate how the surface functionalization alters the hydrophilicity, the accessible hydrophilic surface, and the interaction potential of the nanocomplex surface with biological media. The protein coating increases the interaction potential of the surface of the nanocomplex. The weaker interaction potential of the polysaccharide coating and the relatively large hydrophilic surface area allow dynamic and reversible binding, while the fatty acid rearranges into a partial bilayer with very strong hydrophilicity. The findings establish the hydration enthalpy as a quantitative basis to determine and interpret nanoparticle interactions with proteins, membranes, and biological fluids, and provide a thermodynamic foundation for designing nanocarriers with predictable biological reactivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

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

Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite

Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.

Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗
DOE OSTI2025

Bulk Anion Recognition Kinetically Holds Back Interfacial Adsorption

The competition between bulk and interfacial phenomena underlies many key processes in complex chemical phenomena and transport. While competitive processes are often framed in a thermodynamic context, opportunities to leverage transient species found away from equilibrium can provide a kinetic handle to achieve unconventional reaction outcomes. In this work, we outfit an iminoguanidinium headgroup capable of selective SO 4 2– complexation with alkyl tails of varying complexity to probe competitive bulk and interfacial reaction pathways and tune kinetic pathways for selective chemical separations. Using sum frequency generation (SFG) vibrational spectroscopy we unexpectedly find that adsorption of ligands to the air–aqueous interface was dramatically slowed down for species with increasingly hydrophobic tails. Underlying this phenomenon, we show that the formation of bulk colloidal species with differing propensities for SO 4 2– inhibited surface adsorption via a kinetic bottleneck in the exchange of molecular extractants with colloidal aggregates. This kinetic effect could open up avenues to access unconventional selectivity via complexation of strongly coordinating species in the bulk phase, allowing for more weakly coordinating species to transport via interfacial mechanisms. Furthermore, this work broadly probes nonequilibrium phenomena in chemical separations that arise through unexpected interfacial events that are neglected in traditional equilibrium descriptions.

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

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

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

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al 3+ ) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Furthermore, our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

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

Epitaxial 2D Magnet and Topological Insulator Heterostructures

I will discuss our latest advances on the epitaxial growth of 2D van der Waals (vdW) magnets and their integration with topological insulators (TI). This work is motivated by the realization of topological phases such as the quantum anomalous Hall effect and highly efficient spin-orbit torque produced by TIs. Our initial studies of MnSe2 growth on Bi2Se3 showed a tendency for the interdiffusion of Mn into the Bi2Se3. This ultimately led to the synthesis of MnBi2Se4 (MBS), a new magnetic TI. Interestingly, the vdW phase is not the thermodynamically stable phase and bulk crystals do not exist, so the epitaxial stabilization of MBS creates the opportunity to explore the magnetic and topological properties of this material. We find that MBS is a layered antiferromagnet, similar to MnBi2Te4, but a difference is that the magnetic moments lie in the plane of the film. Angle resolved photoemission experiments show the presence of a topological surface state with Dirac dispersion. For bilayers of 2D magnets and TIs, we have developed FGT films on Bi2Te3. We first optimized FGT by studying its growth on Ge(111) substrates, where we find that kinetic considerations play a major role. Using cross-sectional scanning transmission electron microscopy and scanning tunneling microscopy, we optimize the FGT films to have atomically smooth surfaces and abrupt interfaces with the Ge(111). Subsequently, we have developed the growth of FGT on Bi2Te3 for the integration of 2D magnets with Tis. Interestingly, we observe room temperature ferromagnetism in FGT/Bi2Te3 heterostructures by varying the growth conditions.

Kawakami, Roland↗
DOE OSTI2026

Modulating Operational Conditions to Mitigate Deactivation in Formate Dehydrogenation on Pd Phases

In heterogeneous catalysis, poisoning by surface-bound intermediates poses a major barrier to sustained catalyst performance in (de)hydrogenation reactions. Formate/bicarbonate systems, as liquid organic hydrogen carriers (LOHCs), offer a CO 2 -integrated, low-temperature pathway for hydrogen storage and release, making them attractive for circular energy applications. However, their lower hydrogen density and susceptibility to catalyst deactivation limit their competitiveness compared to conventional LOHCs like methylcyclohexane. Here, this study investigates the mechanistic origins of formate (HCOO – ) dehydrogenation and associated deactivation on Pd interfaces. Using density functional theory (DFT) simulations, we show that under thermocatalytic conditions, strongly bound formate accumulates on the catalyst surface (Pd(111)), blocking active sites, raising activation barriers, and leading to progressive performance loss. Because formate adsorption involves charge transfer, we exploit its sensitivity to electronic structure by modulating the electrochemical potential of the catalyst. Our results reveal that hydrogen transfer from water and formate exhibits opposing potential dependencies, providing insights into the opposing driving forces behind both catalytic activity and poisoning. To further probe this phenomenon, we examine electrochemically induced phase transitions in Pd, focusing on PdO(100) and PdH(110), which are stable under oxidizing and reducing potentials, respectively, and demonstrate enhanced dehydrogenation activity between −0.4 and 0.2 V vs standard hydrogen electrode (SHE). Complementary thermal treatments help decouple kinetic and thermodynamic contributions to intermediate binding. These findings underscore the critical role of the catalyst phase and external stimuli in dictating poison-active site interactions and highlight phase engineering as a promising strategy to mitigate deactivation. This work offers mechanistic insights and design principles for developing more resilient and efficient catalysts for LOHC applications under realistic operating conditions.

Pd phase↗
DOE OSTI2025

Stimulus-Responsive Modulation of Solvation Environments in Solid Catalysts

Liquid environments play a crucial role in the biological processes occurring in living organisms as well as in many human-made processes involving electrochemistry, photo-, and thermocatalysis. In the majority of these systems, aqueous phases are ubiquitous due to water’s natural abundance. Water molecules, however, can exert large changes in the chemical environment of catalytically active sites, altering the reaction rates, selectivity, and catalyst stability. These solvation effects induced by water molecules near catalytic sites can drastically change the energy landscape and unlock unique reaction pathways with far more favorable kinetics. In nature, living organisms couple these complex interactions with detection, communication, and actuation mechanisms to induce self-regulatory behavior, ensuring stability of the system and thus long-term durability. Extrapolating this behavior to heterogeneous catalysis is desirable because the resulting “smart materials” can potentially unlock new chemical conversion processes with higher atom efficiency, rates, and stability. The combination of polymer chemistry and heterogeneous catalysis has introduced versatile approaches for creating materials that can respond to cues in the reaction medium that alter the accessibility, intrinsic activity, and selectivity of the catalyst. To achieve this, one could combine stimulus-responsive polymers, which undergo a large volumetric phase transition in response to an external stimulus, with a solid catalyst. This chemo-mechanical response has been employed to create a variety of nanoreactor vessels with stimulus-responsive character that turn on- and off- depending on the reaction conditions. In this Account, we focus on the impact of these polymer coatings on the solvation environment around the active site and the implications of these effects on the reaction energy landscape, molecular arrangement of the solvent, electric fields at the catalyst–liquid interface, binding energy, and mobility of surface reaction intermediates. These seemingly subtle changes in solvent molecules induced by the presence of polymers can have a tremendous impact on the development of bioinspired heterogeneous catalysts, reliable chemical clocks, micro/nanoreactors, and robots. The large library of polymer chemistries offers a plethora of combinations of stimulus-responsive mechanisms (e.g., temperature, pH, light, magnetic field, solvent composition), providing the possibility of creating homeostatic catalysts à la carte.

catalysts↗
DOE OSTI2025

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

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

The mechanics of nucleation and growth and the surface tensions of active matter

Homogeneous nucleation, a textbook transition path for phase transitions, is typically understood on thermodynamic grounds through the prism of classical nucleation theory. However, recent studies have suggested the applicability of classical nucleation theory to systems far from equilibrium. Here, in this article, we formulate a purely mechanical perspective of homogeneous nucleation and growth, elucidating the criteria for the properties of a critical nucleus without appealing to equilibrium notions. Applying this theory to active fluids undergoing motility-induced phase separation, we find that nucleation proceeds in a qualitatively similar fashion to equilibrium systems, with concepts such as the Gibbs–Thomson effect and nucleation barriers remaining valid. We further demonstrate that the recovery of such concepts allows us to extend classical theories of nucleation rates and coarsening dynamics to active systems upon using the mechanically derived definitions of the nucleation barrier and surface tensions. Three distinct surface tensions—the mechanical, capillary, and Ostwald tensions—play a central role in our theory. While these three surface tensions are identical in equilibrium, our work highlights the distinctive role of each tension in the stability of active interfaces and the nucleation and growth of motility-induced phases.

Capillary waves↗
DOE OSTI2025

Distinct Fracture Mineralogy That is Out of Equilibrium With Modern Groundwaters Provides Important Context for Subsurface Life

Rock fracture surfaces in the crust are essential habitat for microorganisms. Fracture‐groundwater interfaces provide physical substrates for biofilm growth and are sources of carbon, nutrients, and electron donors and acceptors. To better understand geochemical processes impacting fracture surfaces and the subsurface microbiome, we identified fractures in archived rock cores from the Soudan formation, which is known to host saline groundwaters and isolated microbial communities dependent on rock‐water interactions. Cores with open fractures were thin sectioned and studied via electron microprobe and synchrotron X‐ray fluorescence microprobe. Most fracture surfaces had mineralogy distinct from that of the bulk rock. Chlorite minerals were abundant on fracture surfaces and had elemental compositions suggesting deposition during late‐stage hydrothermal alteration. Fracture‐lining chlorites likely limit access to iron oxide and sulfide minerals that are active in subsurface biogeochemical cycles. Calcium‐rich rinds were also observed along fracture edges. These rinds were too thin and poorly ordered to be identified via light microscopy or X‐ray diffraction; however, Ca K‐edge micro‐X‐ray absorption near‐edge structure spectroscopy identified them as carbonates, minerals not observed in the bulk rock. Thermodynamic modeling shows that carbonate precipitation is largely unfavorable in Soudan groundwaters, indicating that fracture edge conditions differed from those in modern water samples. Because of the low carbon concentrations in Soudan groundwaters, carbonate rinds likely play an important role in subsurface carbon cycling and may mark fracture surfaces that once hosted biofilms. Overall, this study suggests that fracture alteration can both play an active role in and suppress rock‐water interactions essential to subsurface life.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nickel hydroxide–nickel carbonate competitive growth on carbonate surfaces

The thermodynamic and kinetic factors controlling the competitive heterogeneous nucleation and growth of ubiquitous metal carbonate and hydroxide phases are poorly understood. In this work, calcite (CaCO 3 ) and magnesite (MgCO 3 ) powders were reacted with NiCl 2 (0–600 μM) for 7 days at 22 °C and 5 °C. The reacted powders were analyzed with X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and energy-dispersive X-ray spectroscopy to characterize the Ni surface precipitates formed. Evidence from these techniques pointed to the formation of mixed Ni carbonate-Ni hydroxide amorphous surface precipitates. On calcite, XPS detected primarily Ni(OH) 2 despite the initial solutions being more supersaturated with respect to gaspéite (NiCO 3 ) than to theophrastite (Ni(OH) 2 ) by a factor of 17–18. In contrast, NiCO 3 was the dominant component detected by XPS on magnesite in the same conditions. Decreasing the temperature had the effect of increasing the proportion of NiCO 3 to the detriment of Ni(OH) 2 . The experimental observations were consistent with low lattice/cation size mismatch favoring NiCO 3 nucleation and temperature most influencing Ni(OH) 2 nucleation. Comparison to previous work on Co-reacted powders indicated the differences in lattice/cation size mismatch and/or water exchange rate impacted the composition of the surface precipitates more than the relative thermodynamic stabilities of the competing minerals. This work explored the heterogeneous growth regime of Ni carbonate and hydroxide phases on carbonate surfaces and shed light on the factors that control the competition between surface precipitates when mineral surfaces are in contact with aqueous solutions supersaturated with respect to multiple mineral phases. In conclusion, these results contribute to geochemists’ efforts toward interpreting data from geochemical systems with elevated Ni concentrations, improving Ni environmental remediation and recovery strategies, and predicting the fate and transport of Ni in geochemical systems.

Heterogeneous nucleation↗
DOE OSTI2024

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2025

Redox Couples Control Band Bending, Photovoltage, and Quasi-Fermi Levels in Tungsten Oxide (WO 3 ) Photoanodes

Tungsten oxide (WO 3 ) is a well-known photoanode and photocatalyst for photoelectrochemical (PEC) water oxidation. Because the compound has a deep valence band, it can facilitate the oxygen evolution reaction without added cocatalysts, and it can drive the oxidation of species with much higher electrochemical potentials, including the conversion of water to hydrogen peroxide, sulfate to persulfate, and iodate to meta-periodate. Here, we use the liquid vibrating Kelvin probe surface photovoltage (liquid VK-SPV) technique in combination with open circuit potential (OCP) and photoelectrochemical (PEC) scans to assess the possibility of reaching such oxidizing potentials in aqueous electrolytes and at open circuit. Here, this is done by mapping the quasi-Fermi levels of electrons and holes at the interfaces as a function of the light intensity. Nanostructured WO 3 photoelectrodes for this purpose were fabricated by thermal annealing of a tungstic acid solution on fluorine-doped tin oxide. Electrochemical measurements are conducted at open circuit and 400 nm LED light illumination in electrolytes containing fast (O 2 /H 2 O 2 ), slow (O 2 /H 2 O), and very oxidizing (NaIO 4 /NaIO 3 ) redox couples. Photovoltage values scale with the light intensity and with the built-in potential for each redox couple and reach values up to 0.61 V under 20 mW cm –2 illumination for the NaIO 4 electrolyte. This shows that the photoelectrodes behave like Schottky-type diodes whose maximum possible energy output is determined mainly by the built-in voltage of each junction. For slow redox couples, the quasi-Fermi level of the holes increases with light intensity due to hole accumulation at the WO 3 –liquid interface. For example, for the O 2 /H 2 O electrolyte, interfacial hole accumulation and removal occur on the 90–300 s time scale. For the fast hole acceptor H 2 O 2 , on the other hand, the quasi-Fermi level of the photoholes is pinned to the electrochemical potential of the O 2 /H 2 O 2 couple. This limits the energy conversion efficiency of the electrode. Overall, these results reveal the influence of charge transfer thermodynamics and kinetics on the photovoltage of WO 3 . Furthermore, the work further establishes VK-SPV as a contactless method to observe the photovoltage, carrier dynamics, and quasi-Fermi levels of semiconductor-liquid junctions.

Electrodes↗