Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface dynamics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Many-Body Perturbation Theory Approach to Energy Band Alignment at the Crystalline Tetracene–Silicon Interface

Hybrid inorganic–organic semiconductor (HIOS) interfaces are of interest for new photovoltaic devices operating above the Shockley–Queisser limit. Predicting energy band alignment at the interfaces is crucial for their design, but represents a challenging problem due to the large scales of the system, the energy precision required and a wide range of physical phenomena that occur at the interface. To tackle this problem, many-body perturbation theory in the non-self-consistent GW approximation, orbital relaxation corrections for organic semiconductors, and line-up potential method for inorganic semiconductors which allows for tractable and accurate computing of energy band alignment in crystalline van-der-Waals hybrid inorganic–organic semiconductor interfaces are used. In this work, crystalline tetracene physisorbed on the clean hydrogen-passivated 1 × 2 reconstructed (100) silicon surface is studied. Using this computational approach, it is found that the energy band alignment is determined by an interplay of the mutual dynamic dielectric screening of two materials and the formation of a dipole layer due to a weak hybridization of atomic/molecular orbitals at the interface. The significant role of the exchange-correlation effects in predicting band offsets for the hybrid inorganic–organic semiconductor interfaces is also emphasized.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Direct Visualization of Defect‐Controlled Diffusion in van der Waals Gaps

Abstract Diffusion processes govern fundamental phenomena such as phase transformations, doping, and intercalation in van der Waals (vdW) bonded materials. Here, the diffusion dynamics of W atoms by visualizing the motion of individual atoms at three different vdW interfaces: hexagonal boron nitride (BN)/vacuum, BN/BN, and BN/WSe 2 , by recording scanning transmission electron microscopy movies is quantified. Supported by density functional theory (DFT) calculations, it is inferred that in all cases diffusion is governed by intermittent trapping at electron beam‐generated defect sites. This leads to diffusion properties that depend strongly on the number of defects. These results suggest that diffusion and intercalation processes in vdW materials are highly tunable and sensitive to crystal quality. The demonstration of imaging, with high spatial and temporal resolution, of layers and individual atoms inside vdW heterostructures offers possibilities for direct visualization of diffusion and atomic interactions, as well as for experiments exploring atomic structures, their in situ modification, and electrical property measurements of active devices combined with atomic resolution imaging.

Chemistry↗
DOE OSTI2024

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine‐Residue Low Complexity Peptide, as Probed by 2 H Solid‐State NMR

Abstract Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid‐state NMR to assess rotameric inter‐conversions and other dynamic modes of the methionine in the context of a nine‐residue random‐coil peptide (RC9) with the low‐complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine‐containing systems, characterized previously. Namely, the M35 of hydrated amyloid‐β 1–40 in the three‐fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)‐methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Vugmeyster, Liliya↗
DOE OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

Li-Metal Anode in Dilute Electrolyte LiFSI/TMP: Electrochemical Stability Using Ab Initio Molecular Dynamics

Ab initio molecular dynamics simulations were performed for Li + conducting electrolytes based on trimethyl phosphates (TMP) and lithium bis(fluorosulfonyl)imide (Li + FSI – ) salt in contact with a Li-metal electrode. We focused on the transient-state behavior at the electrolyte, interfacial electrolyte–Li-metal electrode, and lithium reference electrode–electrolyte–Li-metal electrode to study dynamics and activation energy barriers of the Li + ion, electrochemical and thermal stability of the interface electrode-electrolyte, and potential behavior of the Li-metal electrode, respectively. Our results show that in the most stable state, Li + ions are tetrahedrally coordinated to three TMPs and one FSI – . The inner solvation shell of a Li-ion is composed of three TMP and one FSI – in one contact ion-pair and four TMPs in a solvent-separated ion-pair. On the other hand, Li-ions transport through electrolyte cages takes place when they are coordinated with three or less molecules that could be a combination of TMPs and FSI – . The decomposition pathway of the LiFSI salt when in direct contact with the Li-metal anode starts with defluorination of FSI – , rapidly losing F– to the lithium surface, forming LiF species. The remaining FSO 2 NSO 2 –2 with the addition of 2e – from the Li-metal decomposes into SO 2 –2 and NFSO 2 –2 . SO 2 –2 deposits on the Li-surface and decomposes into Li 2 O and Li 2 S. The remaining NFSO 2 –2 defluorinates, losing F – ion to the lithium surface, resulting in LiF and the remaining NSO 2 –1 deposits on the lithium surface and decomposes in the following picoseconds, forming several binary compounds such as Li 3 N, Li 2 S and Li 2 O. In contrast, when the salt is solvated by the TMP molecules, avoiding a direct contact with the Li metal electrode, only one defluorination occurs, decomposing the FSI – into FSO 2 NSO 2 –2 and a F – . The two anions remain stable as they are solvated by the TMP molecules. Furthermore, we also analyzed the open circuit potential energy (OCPE) of the Li-metal electrode during the SEI formation. OCPE is calculated from the average local potential profile difference within the Li-metal electrode and a pristine Li-crystal reference electrode (LRE). When no SEI is formed, the Li-metal electrode has an average OCPE of +0.36 eV vs LRE. Due to the formation of a SEI, the Li-metal electrode has an average OCPE between -0.07 and -0.21 eV vs LRE. The OCPE of the Li-metal electrode decreases by ~0.42 eV when a SEI is formed.

25 ENERGY STORAGE↗
DOE OSTI2022

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target–apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix–helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

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

Electro-chemo-mechanical evolution of sulfide solid electrolyte/Li metal interfaces: operando analysis and ALD interlayer effects

Sulfide solid electrolytes (SE) show promise for high-performance solid-state batteries because of their high ionic conductivity and ease of processing. However, sulfide electrolytes have suffered from chemical and electrochemical instability against Li metal anodes. In this paper, we use a suite of in situ/operando microscopy and spectroscopy techniques to demonstrate that Al 2 O 3 interlayers deposited using atomic layer deposition (ALD) delay degradation at the Li/LGPS interface by modifying the solid electrolyte interphase (SEI) chemistry and morphology. Optical and scanning electron microscopy are used to rationalize the electrochemical response of the system, which is attributed to a delayed onset of mechanical degradation at the interface when ALD interlayers are used. Operando X-ray photoelectron spectroscopy demonstrates that the dynamic evolution of SEI chemistry is impacted by the presence of the ALD interlayer. In situ Auger spectroscopy and operando optical microscopy provide visual evidence of spatial heterogeneity of Li plating, which is attributed to the chemo-mechanical degradation of the ALD interphase. The implications of these observations provide valuable insights toward the development of robust interlayers for solid-state batteries.

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

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

OpenFerro v0.1.0

OpenFerro is a Python package for on-lattice atomistic dynamics simulation of ferroic materials. OpenFerro is based on JAX, a high-performance linear algebra package supporting auto-differentiation and GPU acceleration. OpenFerro is designed to minimize the effort required to build on-lattice Hamiltonian models, and to perform molecular dynamics (MD) and Landau-Lifshitz-Gilbert simulations. Unlike existing codes, OpenFerro provides a unified interface to model different types of local order parameters.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2025

Achieving the hydrogen shot: Interrogating ionomer interfaces

The aim of this study is to enable the hydrogen economy and decarbonize various sectors in our environment that requires less expensive and more durable water electrolyzers, which can meet the Hydrogen-Shot target. The key is to improve the ionomer interfaces in low-temperature water electrolyzers as rapidly as possible, but to do so, it requires a systematic and holistic campaign combining both experiments and theory. In this perspective, we discuss the issues of electrolyzers and needs for translational science. We then present the approach that the Energy EarthShot Research Center: Center for Ionomer-based Water Electrolysis is taking in hopes of inspiring the community with this approach that can be leveraged to multiple problems and technologies.Graphical abstractHighlightsOne way to achieve the Hydrogen-Shot goal of low-cost, clean hydrogen, is advancing research and development on the interfaces of water electrolyzers for both performance and lifetime. The Center for Ionomer-based Water Electrolysis is exploring new techniques and strategies to not only interrogate interfacial phenomena in water electrolyzers to increase efficiency and durability, but also a new paradigm related to synergistic, cojoined experimental and theoretical research.DiscussionCatalyst\ionomer interfaces are complex and not fully understood, but through investigating different interfaces and utilizing digital and physical twins, we can elucidate key mechanisms and understanding.Understanding the dynamic double layer in electrochemical systems that use solid electrolytes is crucial to identifying and mitigating the controlling phenomena to enable increased performance and durability at the technology level.Studying the time and length scales of interfacial changes can be a powerful tool to understand reaction mechanisms and changes in the electrolyzer performance and durability.

Fornaciari, Julie C↗
DOE OSTI2020

Deciphering the PMI Surface Chemistry of Lithium-Based PFCs and its Effects on High Performance Plasmas in NSTX

The primary aims or objectives of the proposed work here focus on the dynamic measurement (Illinois PI Allain) and atomistic-based multi-scale computational atomistic simulations (UTK PI: Wirth) of: 1) lithium coatings with re-deposited hydrogenated carbon surfaces and high-Z substrates, 2) lithium coatings on boronized ATJ graphite surfaces and high-Z substrates and 3) lithium coatings on ATJ graphite and high-Z substrates with variation in temperature. Three main primary tasks are proposed to address each of the primary objectives listed above: 1) dynamic in-situ irradiation with D+, He+ species of lithium coatings with variation in interface substrate morphology (e.g. smooth vs rough, fuzz vs nanostructured); 2) validation and coupling with multi-scale computational simulations connecting the irradiated surface to the plasma edge; 3) extrapolation to long-pulse conditions (e.g. function of flux vs fluence studies). In primary objective #3, temperature-based effects will be studied to evaluate three primary mechanisms: a) D uptake and recycling, b) erosion, and c) surface melting and evaporation as a function of substrate morphology and chemistry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗
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

Nanoscale Tracking of the High-Temperature Spin-State Transition in LaCoO 3

The high-temperature spin and electronic transitions in LaCoO 3 have recently been leveraged to create neuromorphic (brain-inspired) devices. While these devices have shown the potential for impactful functionality in next-generation computing systems, the nanoscale dynamics of the spin and electronic transitions that underlie their operation are not well understood. Inhomogeneities related to interfaces, electrode contacts, strain, and crystal defects can all affect device performance, making nanoscale characterization of the transitions essential for producing consistent and reliable devices. Here, we demonstrate the first nanoscale in situ measurement of the spin transition in LaCoO 3 at device-relevant temperatures (25–325 °C) over length scales of tens of nanometers using STEM-EELS. This measurement is enabled by an Al 2 O 3 coating, which prevents unwanted reduction of the LaCoO 3 specimen at high temperature and vacuum. The detailed understanding of LaCoO 3 transition dynamics enabled by such measurements will be crucial for optimizing LaCoO 3 -based neuromorphic devices and increasing reliability for real-world application.

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

pH-Driven Restructuring of Hydration Layers and Cation Ad-sorption at the Alumina-Water Interface

Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.

Electrical double layer (EDL), Surface protonation↗
DOE OSTI2022

Using Ensemble Data Assimilation to Estimate Transient Hydrologic Exchange Flow Under Highly Dynamic Flow Conditions

Abstract Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high‐frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady‐state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub‐daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation‐based approach to estimate the sub‐daily flux under highly dynamic flow conditions using multi‐depth temperature observations at a 5‐min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi‐depth temperature data at 1‐ or 2‐hr time intervals assuming one‐dimensional flow and heat transport governing equation. By comparing estimated fluxes with model‐generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation‐based flux estimation method was able to capture the vertical sub‐daily fluxes using multi‐depth high‐resolution temperature data alone, even in the presence of multi‐dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1‐D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

54 ENVIRONMENTAL SCIENCES↗