SEARCH · Engineering Papers

Results for “native state model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects

Acetylcholinesterase (EC 3.1.1.7), a key acetylcholine-hydrolyzing enzyme in cholinergic neurotransmission, is present in a variety of states in situ, including monomers, C-terminally disulfide-linked homodimers, homotetramers, and up to three tetramers covalently attached to structural subunits. Could oligomerization that ensures high local concentrations of catalytic sites necessary for efficient neurotransmission be affected by environmental factors? Using small-angle X-ray scattering (SAXS) and cryo-EM, we demonstrate that homodimerization of recombinant monomeric human acetylcholinesterase (hAChE) in solution occurs through a C-terminal four-helix bundle at micromolar concentrations. We show that diethylphosphorylation of the active serine in the catalytic gorge or isopropylmethylphosphonylation by the R P enantiomer of sarin promotes a 10-fold increase in homodimer dissociation. We also demonstrate the dissociation of organophosphate (OP)-conjugated dimers is reversed by structurally diverse oximes 2PAM, HI6, or RS194B, as demonstrated by SAXS of diethylphosphoryl-hAChE. However, binding of oximes to the native ligand-free hAChE, binding of high-affinity reversible ligands, or formation of an S P -sarin-hAChE conjugate had no effect on homodimerization. Dissociation monitored by time-resolved SAXS occurs in milliseconds, consistent with rates of hAChE covalent inhibition. OP-induced dissociation was not observed in the SAXS profiles of the double-mutant Y337A/F338A, where the active center gorge volume is larger than in wildtype hAChE. These observations suggest a key role of the tightly packed acyl pocket in allosterically triggered OP-induced dimer dissociation, with the potential for local reduction of acetylcholine-hydrolytic power in situ. Computational models predict allosteric correlated motions extending from the acyl pocket toward the four-helix bundle dimerization interface 25 Å away.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

DOE OSTI2025

A Post-translational Histidine–Histidine Cross-Link Enhances Enzymatic Oxygen Reduction Activity with Greater pH Adaptability

Cross-linked protein residues exist as enzyme cofactors to enable or enhance catalytic activities. Despite their importance in nature, the chemical identity of the cross-links is limited to certain amino acid combinations, whose function and the formation mechanism remain insufficiently understood due to the difficulty in isolating native enzymes without the cross-links. Herein, we report the formation and characterization of both His-Tyr and His-His cross-links under oxidative enzymatic turnover conditions in L29H/F33Y/F43H Mb, a structural and functional model of heme-copper oxidase (HCO). The connectivity of the cross-link was characterized as Ne ε2 (His29)-Cd δ2 (His43) by mass spectrometry (LC-MS/MS) and nuclear magnetic resonance (NMR). Interestingly, formation of the cross-link significantly enhances the oxygen reduction activity of the enzyme at neutral or basic pH with higher product specificity. X-ray crystallography has identified a novel Tyr-His cross-link through a Tyr-O-His linkage. Our mechanistic studies indicate the involvement of high-valent heme-iron and the neighboring tyrosine in an oxidative self-processing pathway to generate the cross-link. In conclusion, this work serves as a new example while providing insights into the enzyme cross-link formation, allowing the design of artificial biocatalysts containing these novel cross-links with higher activity and pH adaptability.

Liu, Yiwei [Univ. of Texas, Austin, TX (United Sta↗

DOE OSTI2025

A simple and highly efficient protocol for 13 C-labeling of plant cell wall for structural and quantitative analyses via solid-state nuclear magnetic resonance

Plant cell walls are made of a complex network of interacting polymers that play a critical role in plant development and responses to environmental changes. Thus, improving plant biomass and fitness requires the elucidation of the structural organization of plant cell walls in their native environment. The 13 C-based multi-dimensional solid-state nuclear magnetic resonance (ssNMR) has been instrumental in revealing the structural information of plant cell walls through 2D and 3D correlation spectral analyses. However, the requirement of enriching plants with 13 C limits the applicability of this method. To our knowledge, there is only a very limited set of methods currently available that achieve high levels of 13 C-labeling of plant materials using 13 CO 2 , and most of them require large amounts of 13 CO 2 in larger growth chambers. In this study, a simplified protocol for 13C-labeling of plant materials is introduced that allows ca 60% labeling of the cell walls, as quantified by comparison with commercially labeled samples. This level of 13 C-enrichment is sufficient for all conventional 2D and 3D correlation ssNMR experiments for detailed analysis of plant cell wall structure. The protocol is based on a convenient and easy setup to supply both 13 C-labeled glucose and 13 CO 2 using a vacuum-desiccator. The protocol does not require large amounts of 13 CO 2 . This study shows that our 13 C-labeling of plant materials can make the accessibility to ssNMR technique easy and affordable. The derived high-resolution 2D and 3D correlation spectra are used to extract structural information of plant cell walls. This helps to better understand the influence of polysaccharide-polysaccharide interaction on plant performance and allows for a more precise parametrization of plant cell wall models.

09 BIOMASS FUELS↗

DOE OSTI2024

4-Clique network minor embedding for quantum annealers

Quantum annealing is a quantum algorithm for computing solutions to combinatorial optimization problems. This study proposes a method for minor embedding optimization problems onto sparse quantum annealing hardware graphs called 4-clique network minor embedding. This method is in contrast to the standard minor embedding technique of using a path of linearly connected qubits in order to represent a logical variable state. The 4-clique minor embedding is possible on Pegasus graph connectivity, which is the native hardware graph for some of the current D-Wave quantum annealers. The Pegasus hardware graph contains many cliques of size 4, making it possible to form a graph composed entirely of paths of connected 4-cliques on which a problem can be minor-embedded. The 4-clique chains come at the cost of additional qubit usage on the hardware graph, but they allow for stronger coupling within each chain, thereby increasing chain integrity, reducing chain breaks, and allow for greater usage of the available energy scale for programming logical problem coefficients on current quantum annealers. The 4-clique minor embedding technique is compared with the standard linear path minor embedding with experiments on two D-Wave quantum annealing processors with Pegasus hardware graphs. We show proof-of-concept experiments where the 4-clique minor embeddings can use weak chain strengths while successfully carrying out the computation of minimizing random all-to-all spin glass problem instances. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Larval connectivity for European green crab management in the Salish Sea and surrounding waters

The presence of invasive species is a growing concern in coastal marine ecosystems because of their adverse effects on biodiversity. The European green crab Carcinus maenas (EGC) is a small crab inhabiting inshore areas. Although it is native to the Northeast Atlantic Ocean and Baltic Sea, its distribution has expanded to North America, where it is an invasive species. Its main food sources are small invertebrate species that support valuable fisheries in the USA. The first presence of EGC in northern Washington was observed around 20 yr ago, along the Pacific Coast of the USA. Recently, EGC has been detected throughout the Salish Sea (Washington, USA, and British Columbia, Canada) wherein spread dynamics are unknown. The overall distribution of EGC is mainly driven by larval dispersal and, in the Salish Sea and surrounding waters, the assessment of EGC population dynamics is essential to understand its migration patterns and prevent its future expansion. To investigate the dispersal patterns of EGC larvae, a larval dispersal model was developed which couples a regional model of hydrodynamic circulation with an individual-based model of ichthyoplankton dynamics. Simulations were performed over 9 yr (2013-2022) to analyze average larval transport trends in the Salish Sea and surrounding waters, interannual variability of EGC larval connectivity, and the influence of larval behavior on connectivity patterns. Lastly, areas were identified to inform invasive species management moving forward. The prediction of likely sources and settlement locations of EGC larvae from the model will help improve the management of the population in the Salish Sea and surrounding waters.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Combining MicroED and native mass spectrometry for structural discovery of enzyme–small molecule complexes

With the goal of accelerating the discovery of small molecule–protein complexes, we leverage fast, low-dose, event-based electron counting microcrystal electron diffraction (MicroED) data collection and native mass spectrometry. This approach, which we term electron diffraction with native mass spectrometry (ED-MS), allows assignment of protein target structures bound to ligands with data obtained from crystal slurries soaked with mixtures of known inhibitors and crude biosynthetic reactions. This extends to libraries of printed ligands dispensed directly onto TEM grids for later soaking with microcrystal slurries, and complexes with noncovalent ligands. ED-MS resolves structures of the natural product, epoxide-based cysteine protease inhibitor E-64, and its biosynthetic analogs bound to the model cysteine protease, papain. It further identifies papain binding to its preferred natural products, by showing that two analogs of E-64 outcompete others in binding to papain crystals, and by detecting papain bound to E-64 and an analog from crude biosynthetic reactions, without purification. ED-MS also resolves binding of the CTX-M-14 β-lactamase, a target of active drug development, to the non-β-lactam inhibitor, avibactam, alone or in a cocktail of unrelated compounds. These results illustrate the utility of ED-MS for natural product ligand discovery and for structure-based screening of small molecule binders to macromolecular targets, promising utility for drug discovery.

MicroED↗

DOE OSTI2026

Insights into Native Single-Atom Electrocatalyst Site Structures

Single-atom electrocatalysts consisting of metal atoms embedded in a carbon matrix are promising next-generation catalysts for green hydrogen production and utilization, CO2 reduction, low-temperature CO oxidation, ammonia production, plastic decomposition, and electrochemical energy storage. The origins of activity and stability for the single-atom sites are still debatable, however, because of constrained insights into their local structure resulting from idealized models and experiments derived from a large number of individual sites. Insights into structural variations around single atomic sites are therefore critical for the continued development of these next-generation catalysts. While electron microscopy commonly provides atomic-scale information about these materials, the beam sensitivity of individual sites makes structural determination by conventional low-voltage (60 keV) techniques challenging. Here, we introduce ultralow-voltage electron ptychography, performed at 30 keV, that enables determination of the lattice structure around individual metal sites in a well-defined single-atom electrocatalyst system while essentially eliminating knock-on structural modifications. Pairing these atomic-scale, site-specific measurements with computational methods will broaden our understanding of the activity and stability of these materials, which will accelerate the development of the next generation of catalysts.

Zachman, Michael [ORNL] (ORCID:0000000319101357)↗

DOE OSTI2026

Impacts of Pasture Conversion to Sugarcane on Water Fluxes and Water Use Efficiency in the Southeastern US

The expansion of sugarcane (cane), a high-yielding perennial crop, will likely reshape the bioenergy landscape in the Southeastern US. However, its ecohydrological implications, particularly following conversion from grazed pastures, a dominant land use in the region, remain highly uncertain. We investigated the impact of cane expansion on evapotranspiration (ET) and its partitioning, and the mechanisms influencing both ET components and water use efficiency (WUE) across multiple scales and growth cycles in subtropical Florida. We combined eddy covariance, biometric measurements, and process-based stomatal conductance (g s ) models. ET was 1.7% lower in cane than in improved pasture (IMP) but exceeded that in semi-native pasture (SN) by 21%. Transpiration (T) followed a similar pattern, consistent with lower g s in cane relative to IMP. Cane had more conservative water use and greater sensitivity of g s to vapor pressure deficit (VPD) compared to IMP pasture, suggesting cane may be more tolerant of increasing atmospheric water demand. In contrast, SN showed lower g s and weaker stomatal sensitivity to VPD compared to cane, resulting in lower T. In cane, stomatal regulation and T varied across growth cycles, with stomata becoming less water conservative as stands matured, highlighting the importance of incorporating stand age-dependent stomatal regulation into hydrological models. Evaporation (E) was higher in cane than pastures (19%–26%), partially offsetting WUE gains. Cane exhibited higher intrinsic WUE (GPP/g s ; Gross Primary Productivity), ecosystem WUE (GPP/ET), and harvest WUE (harvest/ET) than both pasture types. Large-scale pasture-to-cane conversion could produce widely contrasting hydrological outcomes. The net regional impact will depend on the proportion of each pasture type converted and on cane's high g s sensitivity to VPD, which triggers tight stomatal regulation and conservative water use, both of which will become increasingly consequential under intensifying atmospheric water demand.

bioenergy↗

DOE OSTI2025

Maximizing long-term biohydrogen production with Clostridium thermocellum for high solids conversion of lignocellulosic biomass

Biological hydrogen production from lignocellulosic biomass sustainably couples organic waste reduction with renewable energy generation. Efficient conversion is challenged by the structural complexity of lignocellulose and resulting recalcitrance to enzymatic degradation. Clostridium thermocellum natively breaks down biomass with highly effective hemi-/cellulases systems (i.e., cellulosomes) and generates hydrogen in anaerobic cultivation, creating a compelling platform for lignocellulosic biohydrogen production. Achieving commercially viable production rates requires balancing high biomass loading and throughput against uniform mixing conditions required for enzyme dispersion, pH and temperature control, and efficient hydrogen and metabolite removal in continuous operation. To address these barriers to process intensification, we implemented novel reactor and process designs for high-solids lignocellulosic biomass fermentations using the C. thermocellum KJC19-9 strain, genetically engineered for co-utilization of cellulose and hemicellulose sugars (i.e., xylose). Via computational fluid dynamics (CFD) modeling and experimental validation, we achieved a >50% improvement in biohydrogen production with an improved anchor-type impeller morphology, coupled to a threefold reduction in agitation rate. To further reduce rheological constraints and accumulation of toxic metabolites, we then transitioned the process to sequencing fed-batch operation. The resulting process generated 24.87 L H 2 L −1 from 160 g L −1 of deacetylated and mechanically refined (DMR)-pretreated corn stover biomass over 16 days while solubilizing >95% of influent cellulose and hemicellulose, setting a new performance benchmark for continuous production of biohydrogen from lignocellulose.

08 HYDROGEN↗

DOE OSTI2020

A Comparison of RESRAD and GoldSim Models for Assessing Radiological Dose for a RCRA Landfill - 20418

US Ecology Idaho (USEI) operates a Resource Conservation and Recovery Act (RCRA) Subtitle C, Hazardous Waste Treatment, Storage, and Disposal (TSD) Facility on a 640-acre property in Owyhee County, Idaho. USEI accepts a wide variety of RCRA-exempt low activity radioactive wastes, including naturally occurring radioactive material (NORM) and technologically enhanced NORM (TENORM) (USEI 2009). The performance of the landfill with regard to radiological operating permit requirements was previously assessed using a model constructed with the residual radioactivity (RESRAD) computer program, developed by Argonne National Laboratory for the U.S. Department of Energy (DOE). RESRAD supports the evaluation of several environmental transport pathways related to a radionuclide-contaminated soil source term, but it was not conceived as a model for the evaluation of landfill radiological performance. A more comprehensive and realistic model is desired in order to better support submissions to both the State of Idaho's Department of Environmental Quality as well as the U.S. Nuclear Regulatory Commission. A Performance Assessment (PA) computer model was developed using GoldSim software for the USEI RCRA Subtitle C landfill. The RESRAD computer model was specifically developed for calculating soil cleanup criteria and radiological dose and cancer risk from residual radioactive material in soil. To evaluate site-specific conditions, users may select from among a number of available environmental transport and exposure pathways and modify 'default' parameter values. In GoldSim, a user must construct the model 'from scratch.' However, this modeling is facilitated by a number of specialized elements available in GoldSim to support a radiological mass transport model. These include elements for defining radiological decay and ingrowth, container failure and radiological release, advective and diffusive transport, and other processes. Advantages afforded by the use of RESRAD for modeling radiological dose for a RCRA Subtitle C landfill include ease of use and reasonable flexibility in specifying site-specific conditions. GoldSim allows for considerably more flexibility and site-specificity than RESRAD, including evaluation of potentially relevant environmental transport processes not supported in RESRAD. The probabilistic modeling capabilities of GoldSim also far exceed those of RESRAD. The pros and cons of RESRAD and GoldSim for this modeling problem are explored with a focus on identifying approaches and critical factors in identifying the appropriate platform. The capabilities of RESRAD and GoldSim for mathematically modeling the disposal system will be contrasted and compared. For example, the existing RESRAD radiological safety assessment evaluated potentially complete exposure pathways related to infiltration to groundwater and upwards diffusion of radon. Additional transport pathways identified in the USEI Idaho facility Conceptual Site Model include deposition of radon decay products in cover material, root uptake of radionuclides in disposed wastes by native plants, and mixing of cover material by burrowing animals. Results of the two models and how the different transport and dose pathways affect results will be discussed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Selective recovery of native copper from basalt tailings using alkaline glycinate solution

Mafic and ultramafic rocks present an intriguing pathway for CO 2 capture through strategic enhancements to the natural silicate weathering cycle. Simultaneously, these rock types are often hosts to appreciable amounts of metals critical to U.S. energy independence, particularly in the context of alkaline mine tailings from historical metal mining. The research and scientific prerogative, then, is to identify promising mafic and ultramafic feedstocks and conduct exploratory studies to effectively recover these critical minerals while preserving the potential for mineral carbonation. The Keweenaw Peninsula of Michigan hosts the largest native Cu reservoir within basalt in the world and experienced a rich history of Cu production spanning from pre-historic times to the mid-1900s. Today, much of the Cu mining legacy remains as mine waste tailings. In this study, we examined the Cu extractability from Keweenaw Basalt tailings using a sodium glycinate solution, in comparison with acid and sodium hydroxide leaching. Experimental results showed an 85 % Cu extraction rate using sodium glycinate as the extraction solution with negligible release of other cations from the basalt. The kinetic and extraction mechanisms of Cu selective recovery using glycinate solution were discussed using time-resolved experimental data and kinetic geochemical modeling. Theoretical estimation of carbon mineralization potential of all the existing basalt waste tailings (∼500 million tons) can reach 85.5 MMT CO 2 . A total of 0.786 MMT Cu can be recovered with sodium glycinate, with a value of 7.7 billion USD. In conclusion, this novel application of alkaline glycinate for selective Cu recovery from basalt mine tailings demonstrates the viability of selective metal recovery using a non-hazardous chemical while preserving CO 2 capture potential and presents a potential pathway toward reducing energy-related emissions and providing an unconventional domestic source of critical minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗

NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

DOE OSTI2026

Contrasting Carbon–Water–Energy Dynamics in Perennial and Annual Bioenergy Agroecosystems Using Eddy Covariance and Interpretable Machine Learning

Understanding how agroecosystems respond to environmental variability is fundamental to predicting productivity and sustainability under a changing climate. We analyzed 55 site-years of high-frequency eddy covariance observations from five agroecosystems—two perennial grasses (miscanthus and switchgrass), two annual rotation systems (maize–soybean and sorghum–soybean), and a restored native prairie—to examine ecosystem-scale carbon, water, and energy fluxes. Using an interpretable machine-learning framework with regression tree ensembles, Shapley Additive Explanations, and Accumulated Local Effects, we quantified how environmental and temporal factors regulate gross primary productivity (GPP), evapotranspiration (ET), water-use efficiency, and the Bowen ratio. Perennials exhibited stronger physiological buffering and maintained fluxes across a broader range of temperature and moisture conditions, reflecting deeper rooting and persistent canopy cover. Annuals, in contrast, showed greater short-term variability and stronger coupling to atmospheric demand, with GPP and ET declining rapidly under low humidity or soil moisture. Differences in temperature sensitivity of Bowen ratio further revealed that perennials sustained proportionally greater sensible heat flux under cool conditions, whereas annuals exhibited constrained energy exchange when evaporative demand was low. Together, these results demonstrate that crop life cycle and canopy structure are fundamental determinants of ecosystem-scale carbon–water–energy coupling. By integrating long-term flux observations with interpretable machine learning, this study identifies the environmental drivers that shape agroecosystem function and highlights how conversion from annual to perennial feedstocks can enhance climatic resilience and alter land–atmosphere energy feedbacks. These findings provide a data-driven basis for improving crop and Earth-system models and for guiding bioenergy landscape design under future climate scenarios.

Accumulated Local Effects↗

DOE OSTI2026

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization↗

DOE OSTI2025

Digital Droplet PCR and Mesocosm-Based Methods to Evaluate Biocontainment Strategies in a Native Soil Ecosystem

Genetically modified industrial production microbes and their associated bioproducts have emerged as an integral component of a sustainable bioeconomy. However, the rapid development of these innovative technologies raises biosecurity concerns, namely, the risk of environmental escape. Thus, the realization of a bioeconomy hinges not only on the development and deployment of microbial production hosts, but also on the development of secure biosystems and biocontainment designs. Current laboratory-based biocontainment testing systems do not accurately reflect the complexities found in natural environments, necessitating an environmentally relevant analysis pipeline that allows for the detection of rare escapees within a complex soil microbiome and differentiation between closely related strains. To this end, we have developed an approach that utilizes soil mesocosms and integrated digital droplet PCR (ddPCR) system to evaluate the efficacy of novel biocontainment strategies. We demonstrate the utility of this approach by modeling contamination with industrial microbial chasses versus their biocontained counterparts. Here we demonstrate the broad utility of this system by highlighting findings from strains of Saccharomyces cerevisiae that are contained with an inducible toxin anti-toxin system, strains of Synechocystis sp. PCC 6803 contained via gene knockout or toxin anti-toxin system, and strains of Escherichia coli that are contained via genomic recoding. We also show that ddPCR can be used to detect gene copies from E. coli equal to those counted by traditional spot plating assays. The resultant data demonstrates that this system has broad utility across diverse microbial chassis and biocontainment strategies and enables researchers to track the fate of our contaminating microbe with high sensitivity in the soil. The findings presented here support the use of this mesocosm-based approach to assess the environmental impact of industrial microbes and to validate biocontainment strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗