Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interfacial force modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Designing Anion-Exchange Ionomers with Oriented Nanoscale Phase Separation at a Silver Interface

We report performance of polymer electrolyte-based energy systems is significantly impacted by transport within the electrode catalyst layer, where ionomer thin films coat catalyst particles. Proton exchange ionomer thin films have been thoroughly characterized, but few studies have critically examined anion-exchange ionomer (AEI) thin films. Further, none have reported nanoscale phase separation for hydrocarbon AEIs, which is critical to mitigate transport resistances. In this work, a set of hydrocarbon-based AEIs with nanoscale phase separation are developed from tunable block copolymer systems composed of polyisoprene (PIp) and polychloromethylstyrene (PCMS). The effect of the PIp/PCMS ratio, architecture, and thickness on the thin-film morphology of the neutral block copolymer precursors on silicon and silver substrates is investigated using grazing-incidence small-angle X-ray scattering (GISAXS) and atomic force microscopy (AFM). AEIs are prepared by quaternizing with trimethylamine or methylpiperidine and their cation-dependent morphology is characterized at 60 °C and 95% RH. A perpendicularly aligned morphology is observed on silver, while no phase separation is observed on silicon, indicating that silver-polymer interfacial interactions drive phase separation. After quaternization, dipole-dipole interactions induce some disorder, but nanoscale phase separation is still maintained. GISAXS patterns are modeled using a Unified Fit approach to understand water uptake and swelling, and recommendations for AEI design are presented.

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

Control of Two Solid Electrolyte Interphases at the Negative Electrode of an Anode‐Free All Solid‐State Battery based on Argyrodite Electrolyte

Abstract Anode‐free all solid‐state batteries (AF‐ASSBs) employ “empty” current collector with three active interfaces that determine electrochemical stability; lithium metal – Solid electrolyte (SE) interphase (SEI‐1), lithium – current collector interface, and collector – SE interphase (SEI‐2). Argyrodite Li 6 PS 5 Cl (LPSCl) solid electrolyte (SE) displays SEI‐2 containing copper sulfides, formed even at open circuit. Bilayer of 140 nm magnesium/30 nm tungsten (Mg/W‐Cu) controls the three interfaces and allows for state‐of‐the‐art electrochemical performance in half‐cells and fullcells. AF‐ASSB with NMC811 cathode achieves 150 cycles with Coulombic efficiency (CE) above 99.8%. With high mass‐loading cathode (8.6 mAh cm −2 ), AF‐ASSB retains 86.5% capacity after 45 cycles at 0.2C. During electrodeposition of Li, gradient Li‐Mg solid solution is formed, which reverses upon electrodissolution. This promotes conformal wetting/dewetting by Li and stabilizes SEI‐1 by lowering thermodynamic driving force for SE reduction. Inert refractory W underlayer is required to prevent ongoing formation of SEI‐2 that also drives electrochemical degradation. Inert Mo and Nb layers likewise protect Cu from corroding, while Li‐alloying layers (Mg, Sn) are less effective due to ongoing volume changes and associated pulverization. Mechanistic explanation for observed Li segregation within alloying Li x Mg layer is provided through mesoscale modelling, considering opposing roles of diffusivity differences and interfacial stresses.

Wang, Yixian [Materials Science and Engineering Pr↗
DOE OSTI2026

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;↗
DOE OSTI2024

Electron/hole selectivity in organic semiconductor contacts for solar energy conversion (Final Report)

Four major and accomplishments for the project "Electron/hole selectivity in organic semiconductor contacts for solar energy conversion.'' are reported. The first is a development of a model describing the contact-determined behavior of a solar cell. The most basic solar cell consists of an intrinsic (pure) semiconductor that acts as an absorber to which contacts are made that provide the asymmetry needed to create a driving force (voltage) for the directional flow of electrons (current), in other words, to convert sunlight into electrical energy. The asymmetry at the contacts comes from the different rate at which electrons and holes, the charge carriers generated by illumination, are collected. We developed a model that describes the ideal current-voltage, and hence energy converting properties, of an ideal photovoltaic that is entirely determined by the kinetics for these charge collection processes. Second, we measured how organic semiconductor interfacial layers impact contact electron and hole transfer rates at semiconductor interfaces and described how these concepts do or do not determine the open-circuit voltage of a solar cell. The work clearly controverts common general misconceptions about the action of contact interfacial layers, such as the idea that improved selectivity for one carrier over the other results from decreased recombination, and quantitatively demonstrates how specific interfacial layer materials act to improve the efficiency of a solar cell. Third, we measured the sub-band gap external quantum efficiency (EQE) of a series of organometal halide perovskite (“perovskite”) solar cells. These measurements quantified band-tailing and revealed defect states, which can cause recombination, that correlate with composition, performance, and hysteresis. Finally, we developed a semiconductor bipolar membrane that uses light to pump ions. This structure is unique among systems that drive ion gradients using light in that it is designed to pump salt rather than one sign of ion; it is a photochemical salt pump.

14 SOLAR ENERGY↗
DOE OSTI2022

A generalized 3D elastic model for nanoscale, self-assembled oxide-metal thin films with pillar-in-matrix configurations

In recent years, functional oxide-metal based vertically aligned nanocomposite (VAN) thin films have gained interest due to their intriguing physical properties and multifunctionalities stemming from the complex interactions between the two phases in the film and the substrate. In this work, we develop a model for studying the energetics of these thin film systems, including the effects of both lattice mismatch and capillary forces due to interface curvature. Each phase is incorporated into the model using a phase indicator function, and we introduce the capillary forces as body forces using a vector density representation of the interface. The model is implemented using the finite element method to study the deformation of the thin film which is composed of Au nanopillars embedded in a La 0.7 Sr 0.3 MnO 3 (LSMO) matrix on an SrTiO 3 (STO) substrate. Further, the results suggest that the total energy is lowest for random configurations of pillars compared to ordered square and hexagonal lattice configurations, consistent with the random distribution of pillars found in experiments. Furthermore, we find that the interfacial energy dominates the total energy of each configuration, suggesting that interfacial energy in the system is an important design parameter for nanocomposite growth, along with the lattice mismatch.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Mechanism-Resolved PFM of Ferroionic and Ferroelectric Responses in Thickness-Gradient Hf 0.5 Zr 0.5 O 2 Libraries

Resolving growth mechanisms and thickness evolution of functional properties is one of the key tasks in materials discovery and optimization involving thin-film materials, traditionally requiring significant experimental budgets. Here we introduce the combination of thickness-gradient libraries and automated scanning probe microscopy as a systematic pathway to elucidate growth modes and disentangle ferroelectric and electrochemical contributions in ferroelectric thin films. As a model system, we explore the Hf 0.5 Zr 0.5 O 2 (HZO) gradient thin films grown on La x Sr 1-x MnO 3 (LSMO) bottom electrode thin films. Automated piezoresponse force microscopy, spectroscopy, and lithography reveals that irreversible topographic deformation arises from electrochemical activity at the LSMO surface, whereas reversible phase inversion in HZO reflects ferroelectric switching. Automated topography height-map scans are used to further quantify nucleation density, particle-size evolution, and roughness correlations across the thickness-gradient, demonstrating that improved plume stabilization during growth suppresses interfacial reactions and promotes dense, fine-grained HZO conducive to ferroelectric phase formation. This combined materials-engineering and automated-SPM framework establishes a platform for high-throughput, mechanism-resolved characterization of ferroionic and ferroelectric responses in complex oxide films.

FOS: Physical sciences↗
DOE OSTI2026

Heteroepitaxial control of thickness, strain, and domain architecture in few-layer ferroelectric tin monochalcogenides

Thin-film epitaxy and epitaxial strain have been widely exploited to tune domain configurations, switching behavior, and ferroic properties in conventional three-dimensional ferroelectric thin films; however, its application to controlling the properties of two-dimensional (2D) ferroelectrics has remained largely unexplored. Here, using SnX (X = Se, S) as a model system, we demonstrate heteroepitaxial control of thickness, strain state, and domain architecture in few-layer ferroelectric SnX via growth on monolayer MoS2 van der Waals (vdW) templates. Compared with conventional growth, MoS2-templated heteroepitaxy promotes epitaxial alignment, yielding ultrathin SnSe films with improved crystalline quality, full areal coverage, and enlarged lateral dimensions. Strong interfacial epitaxial coupling induces pronounced in-plane strain and stabilizes a hierarchical ferroelastic domain architecture, in which long-range 90° stripe domains are further subdivided into nanoscale rotational variants, as revealed by scanning transmission electron microscopy and synchrotron X-ray microscopy. Piezoresponse force microscopy and second-harmonic polarimetry confirm robust in-plane polarization, while ferroelectricity in SnSe is established through polarization-electric field hysteresis and nonvolatile ferroelectric resistive switching with on/off ratios approaching 1000. A nonvolatile, switchable ferroelectric diode effect further evidences direct coupling between polarization and charge transport. Notably, ferroelectric switchability exhibits a strong thickness dependence and is preserved only below ∼10 layers. This vdW heteroepitaxial strategy is further extended to ferroelectric SnS. The seamless heteroepitaxial integration of 2D ferroelectrics with CMOS-compatible, wafer-scale MoS2 templates provides a general and scalable route for strain-enabled structural and ferroelectric engineering in emerging memory and low-power optoelectronic applications.

Wang, Yueyin↗
DOE OSTI2021

Learning intermolecular forces at liquid–vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

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

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2023

Quantifying interfacial energetics of 2D semiconductor electrodes using in situ spectroelectrochemistry and many-body theory

Hot carrier extraction occurs in 2D semiconductor photoelectrochemical cells. Boosting the energy efficiency of hot carrier-based photoelectrochemical cells requires maximizing the hot carrier extraction rate relative to the cooling rate. One could expect to tune the hot carrier extraction rate constant (k ET ) via a Marcus–Gerischer relationship, where k ET depends exponentially on ΔG°' (the standard driving force for interfacial electron transfer). ΔG°' is defined as the energy level difference between a semiconductor's conduction/valence band (CB/VB) minima/maxima and the redox potential of reactant molecules in solution. A major challenge in the electrochemistry community is that conventional approaches to quantify ΔG°' for bulk semiconductors (e.g., Mott–Schottky measurements) cannot be directly applied to ultrathin 2D electrodes. The specific problem is that enormous electronic bandgap changes (>0.5 eV) and CB/VB edge movement take place upon illuminating or applying a potential to a 2D semiconductor electrode. Here, we develop an in situ absorbance spectroscopy approach to quantify interfacial energetics of 2D semiconductor/electrolyte interfaces using a minimal many-body model. Our results show that band edge movement in monolayer MoS 2 is significant (0.2–0.5 eV) over a narrow range of applied potentials (0.2–0.3 V). Such large band edge shifts could change k ET by a factor of 10–100, which has important consequences for practical solar energy conversion applications. We discuss the current experimental and theoretical knowledge gaps that must be addressed to minimize the error in the proposed optical approach.

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

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2022

Electrochemically-Induced Phase Transformations in Battery Storage Compounds (Final Technical Report)

Compounds of interest for ion storage in advanced batteries frequently exhibit phase transformations as the working ion concentration varies. Under large electrochemical driving forces inherent to practical use, systems are often driven far from equilibrium. This program combines experiments and theory to understand the phase transition behavior of ion insertion compounds when electrochemically driven far from equilibrium. As model systems, we focus on alkaline metal phosphates AMPO4 (A = alkali; M = first row transition metal) of olivine structure, which are both technologically interesting and ideally suited for fundamental study due to the ability to systematically tune transformation strain, and along with it, the phase transformation pathway. Behavior in compositions having large transformation strains (~15 vol%) requiring plasticity for strain accommodation is emphasized. Experimental techniques include operando characterization of structure while simultaneously varying electrokinetic parameters, and high resolution microscopy of nanoscale and interfacial phenomena. Phase-field modeling is used to model the thermodynamics and kinetics of competing transformation pathways, extended to include the effects of plasticity, and integrated with porous electrode kinetic theory to treat multi-particle effects. Success in this project will lead to an ability to design ion storage compounds with predictable transformation pathways, electrochemical kinetics, capacity utilization, and durability. New technologically important compounds may also be discovered.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Deterministic Fabrication of Large-Area, High-Crystallinity Oxide Moiré Superlattices

Oxide twistronics extends moiré engineering beyond van der Waals materials, offering a promising platform for accessing emergent interfacial phenomena arising from the strong coupling of lattice, charge, and orbital degrees of freedom in complex oxides. However, deterministic fabrication of high-crystallinity oxide moiré superlattices over large lateral dimensions remains challenging due to the three-dimensional bonding network of oxides. Here, we demonstrate a scalable, generalized fabrication strategy that enables the formation of high-crystallinity oxide moiré superlattices with clean, chemically bonded interfaces and precisely controlled twist angles down to nominal values of 0.1°, achieving subdegree twist-angle accuracy across large contiguous lateral dimensions approaching the millimeter scale. Using NaNbO3 as a model system, we show that the resulting interlayer coupling drives pronounced structural reconstruction that modifies both the phase structure and ferroelectric domain configuration. Synchrotron-based X-ray 3D reciprocal space mapping reveals the emergence of a single-phase state in twisted bilayers, in contrast to the mixed-phase structure observed in single-layer membranes prior to twist assembly. The structural signatures are further consistent with gradual lattice rotation distributed along the thickness direction that may accommodate interfacial shear strain, distinct from reconstruction observed in van der Waals moiré systems which primarily occurs through in-plane stacking rearrangement. This collective lattice response is correlated with twist-dependent nanoscale electromechanical modulations observed by piezoresponse force microscopy. These results establish a scalable materials platform for oxide twistronics and support the implementation of twist-engineered functionalities in practical, macroscale device architectures.

Ghanbari, Reza [North Carolina State University (N↗
DOE OSTI2022

Nucleation of Grain Boundary Phases

Here, we derive a theory that describes homogeneous nucleation of grain boundary (GB) phases. Our analysis takes account of the energy resulting from the GB phase junction, the line defect separating two different GB structures, which is necessarily a dislocation as well as an elastic line force due to the jump in GB stresses. The theory provides analytic forms for the elastic interactions and the core energy of the GB phase junction that, along with the change in GB energy, determines the nucleation barrier. We apply the resulting nucleation model to simulations of GB phase transformations in tungsten. Our theory explains why under certain conditions GBs cannot spontaneously change their structure even to a lower energy state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating Abnormal Grain Growth in Thermomagnetic Processed Materials with Transfer Learning and Reinforcement Learning

The goal of this research program is to establish the mechanism governing local grain boundary motion, which is needed to design and process desirable microstructures for better performance, by identifying the relative contributions of grain boundary (GB) energy and mobility to grain growth. Classical models for grain growth assume that the primary mechanism for reducing the total interfacial energy is area reduction and that GB restructuring is not significant. This assumption implies that grain growth is locally driven by curvature. However, recent experimental observations using new non-destructive 3D x-ray diffraction microscopy techniques (3D-XRM) reveal that classic descriptors (i.e., curvature, number of neighbors, grain size) do not predict real grain growth. Instead, local GB motion appears to be governed by its energy relative to its neighbors such that low-energy boundaries replace those of higher energy. However, simulations that incorporate GB energy anisotropy still fail to reproduce these observations. These discrepancies suggest that the common assumption for grain growth theory must be re-examined to predict and, thus, control microstructure evolution in real polycrystals. A significant challenge to testing this assumption is due to anisotropic GB mobility. Mobility may cause abnormal grain growth or affect the final grain shapes or growth rate but its true contributions are unknown because it is difficult to measure. For example, observations in Fe have found that grains associated with high energy and high mobility boundaries tend to experience abnormal grain growth, whereas abnormal grain growth is associated with low energy and high mobility boundaries in alumina. As mobility and energy both control GB motion, it is challenging to isolate the local driving forces necessary to test the common assumption that the primary mechanism is area reduction. The novelty of this work is the use of machine learning tools to capture GB mobility and energy from 3D-XRM measurements in polycrystals to test the common assumption used in grain growth models. Machine learning can capture high-order correlations in dynamic systems like those found in the evolving GB topology. The PIs have developed a physics-regularized interpretable machine learning microstructure evolution (PRIMME) model that accurately replicates the grain growth behavior of its trained data set.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Heterogeneous catalysis: Optimal performance at a phase boundary?

Most of the industrially used heterogeneous catalysts have been discovered by trial and error, and despite decades of experience, the discovery of new catalysts continues to be extremely challenging. The drive to uncover guiding principles in catalyst design is more present than ever. We share a series of observations indicating that optimal catalysts typically function at characteristic phase boundaries (e.g., abrupt changes in adsorbate coverage, catalyst structure, etc.) accessed in the reaction conditions. The catalyst exploits the associated instability—the desire to exist in multiple states simultaneously—as a driving force for chemical transformations. In other words, phase boundaries are good places to start the catalyst search, and indeed, we should focus on at least two phases at once rather than just one. Here, we substantiate this claim with several studies that combine statistical operando modeling and experiments. Transpiring from these observations is a hitherto unrecognized vector in catalyst discovery.

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

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗