Search NASA⌕ Search

SEARCH · Search NASA

Results for “evolvable”

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

Quantifying structural errors in cloud condensation nuclei activity from reduced representation of aerosol size distributions

Aerosol effects on clouds and radiation are the dominant contribution to uncertainty in radiative forcing relative to the pre-industrial atmosphere. While previous studies have assessed the impact of parametric uncertainty on modeled forcing, structural errors from the numerical representation of particle distributions have not been well quantified. Here we present a framework for quantifying error in aerosol size distributions and cloud condensation nuclei activity, which we apply to the widely used 4-mode version of the Modal Aerosol Module (MAM4). Box model predictions from the MAM4 are evaluated against the Particle Monte Carlo Model for Simulating Aerosol Interactions and Chemistry (PartMC-MOSAIC), a benchmark model that tracks the evolution of individual particles. We show that size distributions simulated by MAM4 diverge from those simulated by PartMC-MOSAIC after only a few hours of aging by condensation and coagulation in polluted conditions, which leads to large errors in modeled cloud condensation nuclei concentrations. We find that differences between MAM4 and PartMC-MOSAIC are largest under polluted conditions, where the size distribution evolves rapidly though aging by condensation of semi-volatile substances and coagulation among particles. These findings suggest that structural error in modeled aerosol properties contributes to the large inter-model variability in aerosol radiative forcing.

Fierce, Laura M.↗
DOE OSTI2026

Engineering an aldoxime dehydratase with high activity and isomer tolerance for biosynthesis of an O -protected primary cyanohydrin

O-protected primary cyanohydrins (glycolonitriles) are important building blocks for many difunctionalized compounds and precursors to known bioactive molecules. Their synthesis, however, utilizes toxic cyanide, which raises significant safety concerns for industrial synthesis. Here, in this study, we present a cyanide-free enzymatic synthesis of an o-benzyl protected primary cyanohydrin from an (E)- or (Z)-α-oxygen protected aldoxime using an engineered aldoxime dehydratase enzyme from Bacillus sp. OxB-1 (OxdB). In contrast to many evolved enzymes that tend to “specialize” as their activity increases, we used directed evolution to engineer OxdB for efficient dehydration of both isomers in a mixture of (E)- or (Z)-α-oxygen aldoximes with high activity and substrate loading to achieve near quantitative yield. Using this enzyme, we further demonstrate a cyanide-free chemoenzymatic pathway to an o-protected primary cyanohydrin starting from a readily available aldehyde, where the aldehyde is first condensed with hydroxylamine, followed by dehydration using our evolved enzyme. This pathway was readily scaled up to 1 g scale with high substrate loading, demonstrating its utility in industrial synthesis of these important building block functional groups.

Aldoxime dehydratase↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2026

Stoichiometric effects on grain growth in zirconium carbide coatings for high-temperature nuclear fuel

Interest in coated particle fuel for space nuclear propulsion (SNP) has expanded in recent years due to successful demonstrations of the resiliency of the coatings to extreme environments. For SNP applications, the coating layer for the particle design needs to be able to withstand exposure to high temperature hydrogen during operating conditions. ZrC has been proposed as a protective layer, however, it is important to understand the high temperature behavior to ensure adequate protection to this fuel. In this study, surrogate ZrC coated particles were heat treated at 1900 °C up to 300 min, to examine how the microstructure evolves when exposed to high temperature. Scanning electron microscopy and electron backscatter diffraction (EBSD) were conducted to determine grain size and grain boundary character and orientation to determine the degree of change in the ZrC layer post heat treatment. Raman spectroscopy provided insight to understand how the as-fabricated stoichiometry of each sample contributed to the differences in grain growth behavior. Despite the as-fabricated samples showing a similar initial grain size and grain boundary character, the samples with a higher amount of excess carbon exhibited smaller grain areas and slower growth rates when exposed to 1900 °C. In conclusion, this investigation details the as-fabricated microstructure of the ZrC layer, specifically grain size, evolved under high temperature as this can impact the performance of the fuel under operating conditions.

EBSD↗
DOE OSTI2025

Scalability analysis of heavy-duty gas turbines using data-driven machine learning

With the increasing integration of variable renewable energy sources into power systems, the role of flexible power generation technologies like gas turbines (GT) in rapid grid balancing remains crucial. This sustained importance underscores the need for scaled and precise modeling of GT to ensure effective integration within evolving energy frameworks. While physics-driven GT models integrate thermodynamics, fluid dynamics, and combustion principles, they often rely on approximate mathematical representations to accommodate scaling that may not capture the actual complex dynamics for GTs and inertial effects associated to GTs with different ratings. In this study, a data-driven model is proposed using machine learning (ML) techniques to conduct GT scalability analysis and performance evaluation with high accuracy. The ML model, trained on data from various operating conditions and performance parameters, aims to uncover intricate relationships and patterns, resembling GT characteristics at different scales (ratings). The model is developed to capture complex system interaction and to adapt to changing operational scenarios at different capacities, providing valuable insights of power system dynamics. In this study, the real-time digital simulator platform was employed to generate training data for the ML model and assess its dynamic characteristics. The ultimate objective was to develop a detailed modeling framework based on governing equations and data-driven ML capable of predicting key performance indicators, in thermal systems such as GTs, including power output, speed, fuel consumption, and exhaust temperature under diverse operating conditions at different scales. The developed ML framework demonstrated high accuracy, with mean relative errors for GT power prediction, reference speed, exhaust temperature, and compressor pressure ratio (CPR) parameters consistently below 0.1% across typical load fluctuation scenarios. Maximum deviations were limited to approximately 0.5 K for exhaust temperature and 0.009 for CPR, underscoring the model’s ability to replicating dynamic GT behavior with high precision. The adaptability of the ML model enables its application across diverse operational conditions and its extension to other thermal systems. By leveraging advanced ML techniques, this study presents a robust and scalable modeling framework that enhances GT simulation precision, facilitating improved integration into evolving power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Characterization of the procoagulant phenotype of amniotic fluid across gestation in rhesus macaques and humans

Background: Amniotic fluid (AF) plays a key role in fetal development, yet the evolving composition of AF and its effects of hemostasis and thrombosis are poorly understood. Objectives: Here, we aim to determine how the evolving molecular composition of AF relates to its procoagulant properties. Methods: We analyzed the proteomes, lipidomes and procoagulant properties of AF obtained by amniocentesis from rhesus macaque and human pregnancies at gestational-age matched timepoints. Results When added to human plasma, both rhesus and human AF accelerated clotting time and fibrin generation. We identified proteomic modules associated with clotting time and enriched for coagulation-related pathways. Proteins known to be involved in hemostasis were highly correlated with each other and their intensity of expression varied across gestation in both rhesus and humans. Inhibition of contact pathway did not affect the procoagulant effect of AF. Blocking tissue factor pathway inhibitor reversed the ability of AF to block the generation of FXa. The prothrombinase activity of AF was inhibited by phospholipid inhibitors. The levels of phosphatidylserine in AF were inversely correlated with clotting time. AF promoted platelet activation and secretion in plasma. Conclusions: The addition of AF to plasma enhances coagulation in a manner dependent on phospholipids as well as the presence of proteases and other proteins that directly regulate coagulation. We describe a correlation between clotting time and expression of coagulation proteins and phosphatidylserine in both rhesus and human AF, supporting the use rhesus models for future studies of AF biology.

amniotic fluid↗
DOE OSTI2024

Evolution and engineering of pathways for aromatic O -demethylation in Pseudomonas putida KT2440

In this study, biological conversion of lignin from biomass offers a promising strategy for sustainable production of fuels and chemicals. However, aromatic compounds derived from lignin commonly contain methoxy groups, and O-demethylation of these substrates is often a rate-limiting reaction that influences catabolic efficiency. Several enzyme families catalyze aromatic O-demethylation, but they are rarely compared in vivo to determine an optimal biocatalytic strategy. Here, two pathways for aromatic O-demethylation were compared in Pseudomonas putida KT2440. The native Rieske non-heme iron monooxygenase (VanAB) and, separately, a heterologous tetrahydrofolate-dependent demethylase (LigM) were constitutively expressed in P. putida, and the strains were optimized via adaptive laboratory evolution (ALE) with vanillate as a model substrate. All evolved strains displayed improved growth phenotypes, with the evolved strains harboring the native VanAB pathway exhibiting growth rates ~1.8x faster than those harboring the heterologous LigM pathway. Enzyme kinetics and transcriptomics studies investigated the contribution of selected mutations toward enhanced utilization of vanillate. The VanAB-overexpressing strains contained the most impactful mutations, including those in VanB, the reductase for vanillate O-demethylase, PP_3494, a global regulator of vanillate catabolism, and fghA, involved in formaldehyde detoxification. These three mutations were combined into a single strain, which exhibited approximately 5x faster vanillate consumption than the wild-type strain in the first 8 h of cultivation. Overall, this study illuminates the details of vanillate catabolism in the context of two distinct enzymatic mechanisms, yielding a platform strain for efficient O-demethylation of lignin-related aromatic compounds to value-added products.

09 BIOMASS FUELS↗
DOE OSTI2024

Toward Enhancing Wastewater Treatment with Resource Recovery in Integrated Assessment and Computable General Equilibrium Models

Sustainable water management is essential to increasing water availability and decreasing water pollution. The wastewater sector is expanding globally and beginning to incorporate technologies that recover nutrients from wastewater. Nutrient recovery increases energy consumption but may reduce the demand for nutrients from virgin sources. We estimate the increase in annual global energy consumption (1,100 million GJ) and greenhouse gas emissions (84 million t CO2e) for wastewater treatment in the year 2030 compared to today’s levels to meet sustainable development goals. To capture these trends, integrated assessment and computable general equilibrium models that address the energy-water nexus must evolve. We reviewed 16 of these models to assess how well they capture wastewater treatment plant energy consumption and GHG emissions. Only three models include biogas production from the wastewater organic content. Four explicitly represent energy demand for wastewater treatment, and eight include explicit representation of wastewater treatment plant greenhouse gas emissions. Of those eight models, six models quantify methane emissions from treatment, five include representation of emissions of nitrous oxide, and two include representation of emissions of carbon dioxide. Our review concludes with proposals to improve these models to better capture the energy-water nexus associated with the evolving wastewater treatment sector.

42 ENGINEERING↗
DOE OSTI2026

Modulating Mid-Gap Electronic States Through Site-Selective Modification in B-Pbx/B'-CuyV2O5/CdS Heterostructures for Photocatalytic Hydrogen Evolution

We interfaced ..beta..-Pbx/..beta..'-CuyV2O5 compounds, with varying stoichiometries of precisely positioned Pb-ions (x) and Cu-ions (y) in interstitial sites along a tunnel-structured ?-V2O5 framework, with cysteine-capped CdS (cysCdS) quantum dots (QDs) to yield heterostructured photocatalysts. ..beta..-Pbx/..beta..'-CuyV2O5 compounds exhibit midgap electronic states with orbital contributions from both Cu 3d and stereochemically active Pb 6s states that show distinctive light-initiated reactivity with photoexcited QDs. ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures were prepared by linker-assisted assembly (LAA). Scanning and transmission electron microscopy, energy-dispersive X-ray spectroscopy, and Raman spectroscopy revealed that cysCdS QDs were deposited onto surfaces of ..beta..-Pbx/..beta..'-CuyV2O5 via LAA. HAXPES revealed that the site-selective positioning of Pb-ions and Cu-ions promoted close energetic alignment of the midgap states of ..beta..-Pbx/..beta..'-CuyV2O5 compounds with the valence-band maximum of cysCdS QDs. Transient absorption spectroscopy revealed that photogenerated holes were transferred from CdS QDs to midgap states of ..beta..-Pbx/..beta..'-CuyV2O5 compounds on time scales <50 ps. Finally, photoelectrochemical and photochemical experiments revealed that ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures promoted the photocatalytic reduction of H+ to H2. In photoelectrochemical experiments, under oxidative conditions, for all ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures, H2 was evolved at a Pt counter electrode while a sacrificial donor was oxidized at the heterostructure-functionalized working electrode. In contrast, under reductive conditions, for ..beta..-Pb0.152V2O5/CdS and ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures, H2 was evolved at the working electrode. In photochemical experiments, dispersed ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures promoted the reduction of H+ to H2 under white-light illumination; ..beta..'-Cu0.55V2O5/CdS and ..beta..-Pbx/..beta..'-CuyV2O5/CdS heterostructures, for which midgap states have Cu 3d orbital character, generated 2-fold more H2 than ..beta..-Pb0.152V2O5/CdS heterostructures. Cu-ion insertion thus appends additional acceptor surface states that improve ligand-mediated hole transfer from photoexcited QDs, but such states are intrinsically limited in mediating hole transport to the substrate as a result of the low mobility of holes in narrow Cu 3d-states. Our results reveal that the density and orbital character of midgap states of ..beta..-Pbx/..beta..'-CuyV2O5 compounds, tunable through recently developed site-selective ion insertion strategies, determine efficiencies of charge-transfer and charge-transport mechanisms that underpin photocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Extreme Thermal Stabilization of Carborane–Cyanate Composites via B–N Dative-Bonded Boroxine Barriers

Enhancing the thermal stability of cyanate ester (CE) resins is crucial for high-temperature aerospace applications. This study elucidates the degradation mechanisms of CE composites reinforced with carborane additives, focusing on the formation of boron–nitrogen (B–N) dative bonds and their role in improving thermal resistance. Using thermogravimetric analysis, we examined char formation and evolved gases at multiple degradation stages. Reaction sequences and char products were characterized via Fourier transform infrared spectroscopy (FTIR), solid-state 13C and 11B NMR, and evolved gas analysis using mass spectrometry and FTIR. The solid-state 11B NMR was instrumental in detecting various boron oxidation states and the formation of B–N dative bonds during decomposition. These bonds facilitate cross-linking in the char phase, enhancing material integrity at elevated temperatures. However, in inert environments, the volatility of boron additives like carborane limits their effectiveness. These findings advance the understanding of CE composite degradation mechanisms and offer critical insights for developing high-temperature-resistant materials for aerospace applications.

Fourier transform infrared spectroscopy↗
DOE OSTI2025

Unraveling the Surface Termination and Evolution of Surface States for Electrocatalyst PtSn 4 in Alkaline HER

Semimetal PtSn 4 has been experimentally demonstrated as a promising topological electrocatalyst for the hydrogen evolution reaction (HER) under both acidic and alkaline conditions. While two possible mechanisms have been proposed to explain its activity, the role of its surface states in HER remains unclear. It is indeed in question how the surface states of this alloy evolve as HER proceeds. In this study, we investigate the surface termination that sustains conducting surface states on PtSn 4 , and we track their evolution during HER catalysis. We show that a reconstructed surface with a Sn-poor termination reproduces the scanning tunneling microscopy pattern observed in experiments and sustains a conducting surface. Through phase diagram and geometric structure analysis, we outline the HER profile following the Volmer–Heyrovsky mechanism. As hydrogen atoms adsorb onto the surface, the structure undergoes further reconstruction to an equilibrium phase with a coverage of two hydrides per unit cell. Meanwhile, the surface electronic bands evolve in response to interactions with the adsorbed hydrogen atoms. A hybridization diagram is further proposed for understanding the surface state evolution based on wave function and chemical bonding analyses. While the Pt atoms serve as conventional sites for hydrogen binding, the surface states of PtSn 4 are essential for stabilizing the hydrogen antibonding states via in-phase electronic interactions with the Sn components. This stabilization results in frontier surface bands that are responsible for driving the HER catalysis. Here, our findings provide a detailed description for the direct involvement of surface states on PtSn4 when employed as a catalyst for HER.

catalysts↗
DOE OSTI2026

Theoretical Insights into Reaction-Induced Transformation and Tuning of Catalytic Behavior in Heterogenous Catalysis

Reaction-induced transformations in heterogenous catalysis represent diverse phenomena that challenge traditional views of static catalyst surfaces. From surface adsorbate dynamics, atomic rearrangements, to composition and phase transitions, these processes reveal the profound differences between idealized model systems under ultrahigh vacuum and the complex, evolving interfaces that govern real catalytic behaviors under reaction conditions. Here, this perspective reviews recent theoretical efforts to provide atomic-level mechanistic insights into significant reaction-induced transformations and their impact on catalytic activity and selectivity. It underscores the need for an integrated framework that combines predictive simulations with operando characterization to uncover active sites and mechanisms under realistic operating conditions. Achieving this requires accelerating existing simulations to fully capture diverse reaction-induced surface dynamics, enabling scalable and accurate modeling of catalysts as condition-dependent, dynamically evolving systems. Such approaches are critical to bridge the gap between theory and practice, offering a pathway to more impactful and predictive catalyst design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Impact of Interposer Microstructure on Ionic Transport in Liquid-Phase Bicarbonate Electrolysis

The electrochemical reduction of CO 2 (CO 2 RR) is a potentially scalable approach for converting captured carbon dioxide into value-added products. Conventional gas-phase electrolysis systems can suffer from carbonate crossover, which limits the efficiency of the system. Liquid-phase (bi)carbonate electrolysis using bipolar membrane electrode assemblies (BPMMEA) has emerged as a promising alternative. The interposer layer, a porous mass-transport material between the BPM and the catalyst, is an essential component of the MEA, as it allows evolved CO 2 to reach the catalyst surface for reaction. In the absence of this layer, evolved CO 2 generated by the pH swing process at the BPM can be converted back into (bi)carbonate (CO 2 recapture) due to the high bulk pH. Thus, clear design guidelines are needed to maximize CO 2 conversion, minimize CO 2 recapture in the catholyte, and improve energy efficiency. Here, the transport properties of the interposer are systematically characterized by X-ray tomography and symmetric-cell impedance spectroscopy to quantify porosity, tortuosity, and the resulting MacMullin number. We then examine the correlation between these material properties and the electrolyzer performance. We focus on characterizing two commercial porous membrane filters, mixed cellulose ester (MCE) and poly(ether sulfone) (PES).

CO2 electrolysis↗
DOE OSTI2024

Organosilica Nanoparticles

The dynamic field of nanotechnology is continually evolving, with organosilica nanoparticles (OSNPs) standing out as a significant and versatile class of materials. Combining the advantageous properties of organic and inorganic components, OSNPs offer unique capabilities that drive innovation across multiple scientific disciplines. Here, this primer is designed to be an accessible and informative guide for researchers embarking on their journey into the world of organosilica nanoparticles. Whether you want to incorporate OSNPs into your research or seek to understand the latest developments, this primer will provide the foundational knowledge and context needed to navigate this exciting and rapidly evolving field.

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

A Lens into the Cu Nanograin by In Situ Vibrational Spectroscopy

Cu-based catalysts are uniquely capable of C-C coupling during electrochemical CO 2 reduction (CO 2 R), yet further mechanistic understanding remains hampered by the lack of spectroscopically resolved descriptors that demonstrate how surface adsorbates emerge and evolve within their catalytic environment. Here, in this study, we correlate in situ surface-enhanced Raman spectroscopy (SERS) and surface-enhanced infrared absorption spectroscopy (SEIRAS) to resolve the potential-dependent dynamics during CO 2 R on Cu nanograin catalysts. By building on previous benchmarking of low overpotential performance and nanograin structural evolution, we offer a diagnostic framework linking vibrational signatures to catalytic function, unveiling which species appear, persist, and turnover as the electrified surface and interfacial environment evolve under bias. The onset of linear CO is marked below -0.45 V, coincident with persistent adsorbed *OH/*O domains beyond the CO 2 R onset. In this context, Cu nanograins serve as a platform to dissect contributions of adsorbate coverage. By carefully dissecting the potential dependence of emergent twin-defect/step CO stretch bands (P1-P2), alongside the prototypical terrace-site CO stretch band (P3), we provide important context for interpreting coupled spectroscopic trends driven by coverage effects and resolve this for the evidently complex nanograin morphology. Together, these observations highlight the intertwined roles of surface stabilization and interfacial flux in steering multicarbon product formation. By directly linking vibrational signatures to catalytic behavior, this work aims to bridge the gap between observation and control and help guide toward a predictive framework of fine-tuned selectivity for CO 2 R.

Fonseca Guzman, Maria V. [University of California↗
DOE OSTI2026

Microcanonical Kinetics of Water-Mediated Proton Transfer in 4ABAH + ·(H 2 O) n = 4–6 Clusters (ABA = Aminobenzoic Acid): A Model System for Size-Dependent Relaxation to Ergodic Behavior

Here, we leverage the unique properties of the 4ABAH + · (H 2 O) n clusters (ABA = 4-aminobenzoic acid, n = 4−6) to quantitatively address how a finite, isolated system evolves into an ergodic condition starting from localized arrangements in configuration space. This system adopts two distinct structural isomers in which water molecules cluster around the cationic centers of its two protomers with widely separated positive charge centers. These isomers arise from excess proton attachment to either the acid (O) or amino (N) group on opposite sides of the benzene ring. Both forms are captured and kinetically trapped using cryogenic ion methods and then selectively vibrationally excited through their mutually exclusive IR bands involving NH and OH stretching fundamentals. Because the IR excitation lies below the water binding energy, the system can evolve to explore slow, rare events that lead to the interconversion between the two isomers. The rates of these intracluster reactions are determined by using a pump−probe scheme involving ∼5 ns IR pump and UV probe lasers. The rates occur on the microsecond time scale, leading to steady state populations of the isomers, thus revealing the cluster size-dependent fractionation between the two species at microcanonical equilibrium. The steady state distributions are correlated with the expected trend in the cluster size-dependent reaction energetics, which are in turn consistent with changes in the relative densities of states of the two species. These results thus provide an unusually clear example in which complex, protic-solvent-mediated chemical transformations are captured within a finite system at a precisely determined internal energy.

Rana, Abhijit [Yale Univ., New Haven, CT (United S↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗