Engineering Papers⌕ Search

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.

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
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↗
DOE OSTI2025

Akiachak Energy Efficiency Retrofit Project

The goal of the project is to reduce the overall energy use of the Akiachak Native Community (ANC) by implementing energy efficiency measures in five high-use Tribal buildings. This project will have the following outcomes: Projected annual energy savings of $17,369; projected annual reduction in fuel oil #1 of 1,200 gallons and electricity of 17,751 kWh; annual reduction in carbon dioxide emissions of approximately 60,340 pounds/year. ANC will install energy efficiency measures in the Laundry, Tribal Indian Reorganization Act (IRA) Office, Clinic, Daycare, and Police Station. ANC obtained energy audits on these buildings in 2018, and this project will implement high-payback recommendations such as replacing lighting with LEDs, installing setback thermostats and occupancy sensors, replacing furnaces with more efficient models, replacing the circulation pumps with variable speed ones, air tightening, and adding insulation. Buildings will see energy cost reductions from 15% to 40%. These retrofits will help build ANC’s long-term vision for sustainable energy usage and address the first goal of the Tribal IRA Council’s Energy Efficiency and Conservation Strategy, to “create and maintain functionally appropriate, sustainable, accessible, high quality tribal infrastructure and facilities.” ANC intends to replicate this project by using the resulting energy savings to address audit recommendations in other buildings as well as to demonstrate the value of energy efficiency to community members. Other outcomes will include an increase in community resiliency, reduced dependence on outside shipments of fuel oil, training for maintenance staff, and no-touch control of building appliances to reduce transmission of diseases such as COVID-19.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Coupled Roles of Surface Chemistry and Hydrogen-Assisted Cycling in Ruthenium Atomic Layer Deposition on Silicon Oxides

Ruthenium (Ru) is a promising interconnect material for advanced semiconductor technologies due to its favorable scaling characteristics, including a short electron mean free path and strong electromigration resistance. In semiconductor integration, silicon oxide-based dielectrics serve as dominant insulating materials and constitute ubiquitous interfaces for metallization; however, their formation-dependent surface chemistry and its impact on Ru growth remain insufficiently explored. Here, we investigate Ru ALD on native oxide SiO x (N-SiO x ) and thermally grown SiO 2 (T-SiO 2 ) as model substrates using bis(ethylcyclopentadienyl)ruthenium(II) [Ru(EtCp) 2 ] under two distinct reactant-sequence environments: AB-type (Ru(EtCp) 2 /O 2 ) and hydrogenassisted ABC-type (Ru(EtCp) 2 /O 2 /H 2 ). Under the AB-type process, both N-SiO x and T-SiO 2 exhibit pronounced nucleation delay. N-SiO x shows earlier nucleation and higher nucleation density than T-SiO 2 , plausibly attributed to differences in surface hydroxyl populations. Similar temperature-dependent phase evolution is observed on both substrates, with mixed Ru and RuO 2 phases at 250 °C and predominantly metallic Ru at 300 °C accompanied by increased morphological roughening. In contrast, incorporating an H 2 subpulse (ABC-type) mitigates nucleation delay, particularly on hydroxyl-deficient T-SiO 2 , thereby reducing the substratedependent disparity observed under AB cycling. Moreover, RuO 2 formation is suppressed even at 250 °C on both substrates, shifting growth toward more metallic Ru with reduced resistivity (∼20 μΩ·cm at ∼ 20 nm on N-SiO x ). These trends suggest that H 2 influences the surface reaction pathway, contributing to enhanced metallic stabilization and altered early stage growth kinetics. Overall, this work clarifies the coupled roles of substrate chemistry and reactant-sequence design in governing Ru nucleation and early stage film evolution, providing insight relevant to next-generation interconnect integration and future area-selective deposition strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Project Phase 1 Report: Reducing Data Center Peak Cooling Demand and Energy Costs With Cold Underground Thermal Energy Storage (Cold UTES)

Cold Underground Thermal Energy Storage (Cold UTES) is an ultra-long duration grid energy storage technology. With Cold-UTES, low-cost grid power is converted to cold thermal energy and stored in the native subsurface rock at the point of use. Cold UTES is one approach within the general category of engineered geothermal systems. Cold UTES for peak-hour cooling of data centers (DCs) was studied for deployment in Maricopa County, Arizona and Loudoun County, Virgina using thermal storage capacities from 4 GWh-th to over 1,000 GWh-th (>1 Terawatt-hour). The two sites have different power and transmission systems, available grid energy resources, daily and seasonal load profiles, weather conditions, and grid regulatory requirements. The study results indicate high value for both locations and because of this, likely indicates value across most of the US and the world. The basis of the study was a 1,000 MW-e hourly electric use of DC computing and auxiliary loads, which was modeled as 1 GW-th of thermal load to a dry cooled heat rejection system - i.e. a cooling system that does not consume water. The electric power required for cooling the DC varies as the air temperature changes. In cold weather only the dry-coolers are used, with an electrical load for cooling load as low as 10 MW-e. In hot summer hours, chillers and dry- coolers are required, which raises the electrical load for cooling load to as much as 300 MW-e. The continuous and peak cooling electrical loads result in a grid interconnection requirement of no less than 1,300 MW-e. From both a grid and thermal design modeling perspective the 1.3 GW-e could either be a single facility or result from the total load at multiple sites.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗
DOE OSTI2024

STITCHES: a Python package to amalgamate existing Earth system model output into new scenario realizations

Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Iowa Tribe of Kansas and Nebraska: Advancing Clean, Resilient, and Sovereign Energy (Summary Report of Communities LEAP Activities)

The Iowa Tribe of Kansas and Nebraska (ITKN) is a federally recognized Native American Tribe located along the Missouri River on the border of northeast Kansas and southeastern Nebraska. There are over 800 residents (Tribal citizens and non-Tribal) who live on the reservation, as well as more than 500 people who visit or work on the reservation on a daily basis. The ITKN faces many energy challenges, including rising service costs and dozens of power outages annually that impact resident well-being and business activities on Tribal lands. Power service issues are made more challenging by the remoteness of the reservation, which is 20 miles from the nearest town. Despite this, the ITKN has a long history of cultural and economic resilience: Local self-reliance, environmental stewardship, respecting the carrying capacity of the land, and strengthening the community are Tribal communities' traditional strengths. Long-term energy goals for the ITKN are centered around achieving energy sovereignty. Priority actions include: (1) Establishing a Tribal Utility Authority (TUA) to promote social welfare and community development.; (2) Deploying renewable community microgrids with ground-mount solar arrays and sustainable energy storage systems to advance energy sovereignty, resilience, and reliability. To advance these goals, the ITKN partnered with the U.S. Department of Energy's (DOE's) Communities LEAP (Local Energy Action Program) pilot. From August 2022 to March 2024, the ITKN community coalition collaborated with technical assistance providers at DOE's National Renewable Energy Laboratory (NREL) and Sandia National Laboratories to evaluate TUA planning needs and microgrid deployment scenarios. This final report details the Communities LEAP technical assistance process, models and analysis performed, and results.

24 POWER TRANSMISSION AND DISTRIBUTION↗