Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico experiments”

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.

1,516 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

PowerModelsGAT-AI: Physics-Informed Graph Attention for Multi-System Power Flow With Continual Learning

Solving the alternating current power flow equations in real time is essential for secure grid operation, yet classical Newton–Raphson solvers can be slow under stressed conditions. Existing graph neural networks for power flow are typically trained on a single system and often degrade on different systems. We present PowerModelsGAT-AI, a physics-informed graph attention network that predicts bus voltages and generator injections. The model uses bus-type-aware masking to handle different bus types and balances multiple loss terms, including a power-mismatch penalty, using learned weights. We evaluate the model on 14 benchmark systems (4 to 6,470 buses) and train a unified model on 13 of these under contingency conditions with up to two branch outages, achieving an average normalized mean absolute error of 0.89% for voltage magnitudes and R 2 >0.99 for voltage angles. We also show continual learning: when adapting a base model to a new 1,354-bus system, standard fine-tuning causes severe forgetting with error increases exceeding 1000% on base systems, while our experience replay and elastic weight consolidation strategy keeps error increases below 2% and in some cases improves base-system performance. Interpretability analysis shows that learned attention weights correlate with physical branch parameters (susceptance: r=0.38 ; thermal limits: r=0.22 ), and feature importance analysis supports that the model captures established power flow relationships.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗
DOE OSTI2025

Probing Operando Electrochemical Strain Generation in α-NaFeO 2 Composite Cathodes during Cycling of Na-Ion Batteries

The transition metal oxide (TMO) cathodes in Na-ion batteries suffer from low-capacity retention. Chemo-mechanical instabilities lead to the deterioration of the electrochemical performance of TMO cathodes in Li-ion batteries. However, there is not much known about the chemo-mechanical instabilities in the TMO cathodes for Na-ion batteries. Understanding the governing forces behind the interplay between the electrochemical performance and mechanical stability in TMO cathodes is critical for the development of Na-ion batteries. Here, we synchronize the digital image correlation (DIC) technique with electrochemical analysis to capture the real-time deformation behavior of the α-NaFeO 2 cathodes during cycling. When the charge cutoff voltage is 3.6 V, the cathode experiences reversible deformations (except for the first cycle). There is negative strain (shrinkage) generation during Na extraction and positive strain (expansion) generation during the subsequent Na insertion. A detailed analysis of the potential-dependent strain rate evolution points out complicated phase transformations and nonequilibrium conditions in the α-NaFeO 2 cathodes during cycling. When the charge cutoff voltage was increased to 4.2 V, there was a rapid capacity loss and large plastic deformations in the α-NaFeO 2 cathodes. We provide an in-depth discussion about the possible mechanisms behind the chemo-mechanical instabilities in the α-NaFeO 2 . In conclusion, the correlation is critical to develop material-based strategies to mitigate instability mechanisms in TMO cathodes for Na-ion batteries.

Wable, Minal [University of Maryland Baltimore Cou↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Direct Observation of Vortex Liquid Droplets in the Iron Pnictide Superconductor CaKFe 4 As 4 at 0.5T c

Type-II superconductors under magnetic fields remain in a quantum-coherent, non-dissipative state as long as vortices are pinned. Dissipation emerges when vortices depin, a process often driven by thermal fluctuations and commonly associated with a melting transition from a vortex solid to a vortex liquid. Macroscopic experiments almost always observe this transition close to the superconducting critical temperature 𝑇 𝑐 . However, how the vortex solid responds to thermal fluctuations at the scale of individual vortices, far below the melting transition, remains largely unexplored. Here, we use scanning tunneling microscopy (STM) to directly visualize vortices in the iron-based superconductor CaKFe 4 ⁢As 4 (𝑇 𝑐 ≈35 K ). We observe the formation of vortex liquid droplets—spatially localized regions where vortices exhibit strong thermal fluctuations—at temperatures as low as 0.5 𝑇 𝑐 . These results demonstrate that the onset of dissipation at the local scale occurs at temperatures significantly below 𝑇𝑐 in type-II superconductors, revealing a previously unrecognized regime of vortex dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗
DOE OSTI2021

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Roadmap and Benchmarking: Privacy in Federated Load Forecasting

Data-driven techniques for energy demand forecasting continue to emerge with promising impacts on distribution grid planning. However, the development of robust and generalizable machine learning models requires that representative high quality training data are available. Distributed energy resources have begun to embed intelligence, gathering large amounts of data on customer demand, behavior, and household devices that are connected to the grid. Though utilities aggregate meter-level demand data for load shaping, demand response, outage management, reliability planning, and billing applications, there lies an inherent privacy concern in sharing consumption data that may identify individual consumer behavioral patterns. Hence, while sharing the data is crucial, the private sensitive customer data must be safeguarded from being exposed or manipulated. In this study, we propose a roadmap for implementing a based privacy preserving framework to support the advancement of data-driven analytics in data-sensitive distributed energy resources environments. The roadmap incorporates federated learning–a distributed training framework, differential privacy–a statistical framework that provides guarantees to safeguard the leakage of sensitive data, secure multiparty computation and homomorphic encryption– techniques for encrypting model gradients and applying secure aggregation on the server. Moreover, we perform baseline experiments on the federated short-term load forecasting (STLF) task using open-source residential load profile datasets, offering insights into the challenges of integrating differential privacy into federated learning.

Abebe, Waqwoya [Oak Ridge National Laboratory (ORN↗
DOE OSTI2026

The U.S. Fusion Materials Community Roadmap: Near-term research priorities for the development of plasma-facing and structural materials for fusion power plants

In response to the needs of a rapidly growing private fusion industry, the U.S. Fusion Materials Coordinating Committee (FMCC) and the broader U.S. fusion materials research community undertook an extensive effort to create a comprehensive roadmap for fusion materials development. The result of this effort was the U.S. Fusion Materials Community Roadmap (US-FMCR), which describes the steps needed to advance the technical maturity of leading candidates for plasma-facing materials and structural materials for fusion power plants from laboratory-scale experiments to a point of sufficient technological readiness for industrial adoption and implementation. However, researchers face significant resource constraints as well as very aggressive pilot plant development timelines. Thus, the research strategies detailed in the US-FMCR require further assessment to downselect the specific tasks that must be prioritized within the next two to three years, in order to make the most efficient use of funding, human resources, and experimental facilities. This paper presents an overview of the US-FMCR and its development process. We also present the subset of research objectives that the FMCC identified as the most urgent research priorities for the U.S. fusion materials research community. The state-of-the-art of materials research is also highlighted for each class of materials considered in the US-FMCR. The recommendations presented here integrate an extensive evaluation of the current status of fusion materials research with a broad cross-section of opinion from the wider U.S. fusion community.

Ferry, Sara [Massachusetts Institute of Technology↗
DOE OSTI2023

Oxidation Chemistry of Bicarbonate and Peroxybicarbonate: Implications for Carbonate Management in Energy Storage

Carbonate formation presents a major challenge to energy storage applications based on low-temperature CO 2 electrolysis and recyclable metal–air batteries. While direct electrochemical oxidation of (bi)carbonate represents a straightforward route for carbonate management, knowledge of the feasibility and mechanisms of direct oxidation is presently lacking. Herein, we report the isolation and characterization of the bis(triphenylphosphine)iminium salts of bicarbonate and peroxybicarbonate, thus enabling the examination of their oxidation chemistry. Infrared spectroelectrochemistry combined with time-resolved infrared spectroscopy reveals that the photoinduced oxidation of HCO 3 – by an Ir(III) photoreagent results in the generation of the short-lived bicarbonate radical in less than 50 ns. The highly acidic bicarbonate radical undergoes proton transfer with HCO 3 – to furnish the carbonate radical anion and H 2 CO 3 , leading to the eventual release of CO 2 and H 2 O, thus accounting for the appearance of H 2 O and CO 2 in both electrochemical and photochemical oxidation experiments. Here, the back reaction of the carbonate radical subsequently oxidizes the Ir(II) photoreagent, leading to carbonate. In the absence of this back reaction, dimerization of the carbonate radical provides entry into peroxybicarbonate, which we show undergoes facile oxidation to O 2 and CO 2 . Together, the results reported identify tangible pathways for the design of catalysts for the management of carbonate in energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Machine learning pipeline for denoising low signal-to-noise ratio and out-of-distribution transmission electron microscopy datasets

High-resolution transmission electron microscopy (HRTEM) is crucial for observing material’s structural and morphological evolution at Angstrom scales, but the electron beam can alter these processes. Devices such as CMOS-based direct-electron detectors operating in electron-counting mode can be utilized to substantially reduce the electron dosage. However, the resulting images often lead to a low signal-to-noise ratio, which requires frame integration that sacrifices temporal resolution. Several machine learning (ML) models have been recently developed to successfully denoise HRTEM images. Yet, these models are often computationally expensive, and their inference speeds on GPUs are outpaced by the imaging speed of advanced detectors, precluding in situ analysis. Furthermore, the performance of these denoising models on datasets with imaging conditions that deviate from the training datasets has not been evaluated. To mitigate these gaps, we propose a new self-supervised ML denoising pipeline specifically designed for time-series HRTEM images. This pipeline integrates a blind-spot convolution neural network with pre-processing and post-processing steps, including drift correction and low-pass filtering. Results demonstrate that our model outperforms various other ML and non-ML denoising methods in noise reduction and contrast enhancement, leading to improved visual clarity of atomic features. Additionally, the model is drastically faster than U-Net-based ML models and demonstrates excellent out-of-distribution generalization. The model’s computational inference speed is in the order of milliseconds per image, rendering it suitable for application in in-situ HRTEM experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SILICON CARBIDE THERMOMETRY USING RAMAN SPECTROSCOPY ON IRRADIATED TRISO PARTICLES

Silicon carbide (SiC) passive thermometry has emerged as a promising post-irradiation examination (PIE) technique for estimating irradiation temperature near the end of irradiation. While dilatometry techniques have been traditionally used to analyze prismatic samples after irradiation, Raman spectroscopy has recently been shown to provide comparable results by analyzing Raman-active phonon modes. In this work, Raman spectroscopy has been applied to the SiC layer of cross-sectioned irradiation tristructural isotropic (TRISO) particles from the AGR-5/6/7 experiment to evaluate the feasibility of particle-scale passive thermometry. Two-dimensional Raman mapping was used to measure the position of the SiC longitudinal optical (LO) phonon, which was then converted to an apparent irradiation temperature using a previously established empirical correlation. Particles from AGR-5/6/7 Compacts 2-2-1 and 5-1-3 were selected as their calculated time-averaged, volume-averaged (TAVA) temperatures (828°C and 706°C, respectively) fall within the range of the sensitivity of the experimental approach. The use of SiC thermometry was anticipated to confirm or highlight potential deviations from calculated end-of-life TAVA temperature across compacts. Four particles from Compact 2-2-1 and two particles from Compact 5-1-3 were selected based on 110mAg inventory (either some measurable activity or below the minimum detection limit) which is commonly used as an indicator of in-pile temperature variation of particles within the same compact. Across every particle selected it was determined that the average LO peak position was located around 968 cm-1 to 969 cm-1. Using the previously determined empirical correlation, this corresponds to an irradiation temperature around 950°C to 966°C, which does not align with the reported TAVA temperature values. This discrepancy in apparent irradiation temperatures likely reflects a combination of uncertainties in the calculated particle temperatures and differences in irradiation history between the present specimens and those used to establish the empirical Raman calibration such as neutron flux (damage rate) and SiC microstructure (as fabricated and irradiated).

Vawdrey, Josh [ORNL]↗
DOE OSTI2026

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗
DOE OSTI2026

Effect of spacer grids on high-burnup fuel fragmentation, relocation, and dispersal

Increasing the fuel burnup limit in light-water reactors to improve fuel cycle economics requires a strong technical foundation. Experimental observations from the Halden and Studsvik programs have revealed severe fuel fragmentation during loss-of-coolant accident (LOCA) conditions, highlighting the need for additional technical evaluation. Consequently, further LOCA test data are needed to complement existing findings and improve the understanding of fuel fragmentation, relocation, and dispersal (FFRD) behavior. Oak Ridge National Laboratory’s Severe Accident Test Station has played a significant role in advancing the understanding of high-burnup fuel fragmentation, relocation, and dispersal phenomena. One remaining gap in the available experimental database is the effect of fuel assembly structural features on cladding deformation behavior during a LOCA, and more specifically, their impact on the fuel’s ability to fragment, relocate, and disperse. Recent analyses using the BISON fuel performance code suggest that cladding deformation near grid spacers will remain below the 3% threshold that has been reported in the NRC Research Information Letter, indicating that the cladding could remain mechanically constrained during the LOCA event. This paper builds upon the BISON analyses to design and conduct a series of out-of-cell tests aimed at further evaluating cladding deformation in and around grid spacers. In addition, these tests were used to assess local cladding temperature conditions and compare them against analytical predictions in order to better replicate expected in-reactor behavior. Finally, an in-cell high-burnup LOCA test was designed and performed to evaluate the effects of a grid spacer, or cladding restraint, on fuel fragmentation, relocation, and dispersal susceptibility. The high-burnup test results differed from those of historical LOCA experiments, with a recorded rupture temperature of 861°C. Two ballooned regions and corresponding rupture openings were observed, with rupture widths of approximately 0.64 mm for both ruptures and rupture lengths of 4.8 mm and 5.6 mm, respectively.

Capps, Nathan [ORNL]↗
DOE OSTI2026

Scan‐Path‐ and Initial‐State‐Dependent Superdomain Switching in (111)‐Oriented PZT

Polarization switching in ferroelectric materials arises from the collective evolution of complex domain hierarchies, yet deterministic control over these processes remains challenging. Here, we investigate scan-path- and initial-state-dependent switching in epitaxial (111)-oriented PbZr 0.2 Ti 0.8 O 3 thin films using automated AFM-based writing combined with quantitative 3D piezoresponse force microscopy. We show that the scan trajectory acts as an experimentally accessible control parameter for superdomain formation. Box-in-box raster scans reproducibly stabilize ordered stripe superdomains with a reduced subset of symmetry-allowed variants, whereas spiral trajectories generate frustrated mixed-variant states with a broader distribution of final microstructures. Automated pulsing experiments further show that the local superdomain configuration at the nucleation site strongly influences the final written morphology. Phase-field modeling qualitatively reproduces the contrast between representative initial-state geometries and supports the role of compatibility constraints among competing ferroelastic pathways. These findings establish scan-path and initial-state engineering as practical handles to program ferroic order in hierarchical ferroelectric domain structures.

Vasudevan, Rama K. [Oak Ridge National Laboratory ↗
DOE OSTI2026

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2026

Distributed Mafic Rock Resources for Carbon Mineralization in Arizona

Ex-situ carbon mineralization is a process by which CO2 is reacted with alkaline silicate minerals and rocks to produce stable carbonate materials, which can be used for other industrial processes. Arizona, U.S.A., hosts abundant surficial mafic rocks in three young volcanic fields, Geronimo-San Bernardino, San Francisco, and Springerville, and other distributed locations throughout the state. We created a Mafic Rock Resource Inventory (MRRI) that categorizes geochemical, physical, and textural characteristics of a diverse suite of surficial mafic rock samples and provide a benchmark reaction dataset parameterizing the temperature, pressure, and pH conditions best suited ex-situ mineralization in different rock types. MRRI data is publicly available online via a map-viewer. We establish two reaction condition sets, varied in temperature, pressure, and pH, where crystal-rich and glassy rocks reach maximum reaction extent and different carbonate phases are formed. Systematic ex-situ mineralization experiments on 21 diverse rock types show trends in geochemical, mineralogical, and reactivity behavior and establish maximum effective capture capacity. From this, scoria cones in three Arizona volcanic fields have a ~62 Gt effective CO₂ storage capacity with one of the fields having a ~42 Gt storage capacity in lava flows. Reactivity results have applications to alkaline mafic rock resources exposed globally, including producing additional effective storage capacity estimates and scaled commercialization of mafic rock ex-situ mineralization, should reaction extents be improved through advances in mineralization techniques. MRRI data were used to create a Direct Air Capture to Mineralization (DACM) systems model, technoeconomic analysis, and life-cycle assessment. These documents are presented as three appendices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗