Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation framework”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

ToPolyAgent: AI agents for coarse-grained bead-spring topological polymer simulations

We introduce ToPolyAgent, a multi-agent AI framework for performing coarse-grained molecular dynamics (MD) simulations of topological polymers through natural language instructions. By integrating large language models (LLMs) with domain-specific computational tools, ToPolyAgent supports both interactive and autonomous simulation workflows across diverse polymer architectures, including linear, ring, brush, and star polymers, as well as dendrimers. The system consists of four LLM-powered agents: a Config Agent for generating initial polymer–solvent configurations, a Simulation Agent for executing LAMMPS-based MD simulations and conformational analyses, a Report Agent for compiling markdown reports, and a Workflow Agent for streamlined autonomous operations. Interactive mode incorporates user feedback loops for iterative refinements, while autonomous mode enables end-to-end task execution from detailed prompts. We demonstrate ToPolyAgent's versatility through case studies involving diverse polymer architectures under varying solvent conditions, thermostats, and simulation lengths. Furthermore, we highlight its potential as a research assistant by directing it to investigate the effect of interaction parameters on the linear polymer conformation, and the influence of grafting density on the persistence length of the brush polymer. By coupling natural language interfaces with rigorous simulation tools, ToPolyAgent lowers barriers to complex computational workflows and advances AI-driven materials discovery in polymer science. It lays the foundation for autonomous and extensible multi-agent scientific research ecosystems.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2022

DIGS: deep inference of galaxy spectra with neural posterior estimation

Abstract With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of simulation-based inference (SBI) and amortized neural posterior estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is (a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; and (b) amortized calculations of posterior distributions of said galaxy parameters at the modest cost of a few galaxy fits with Markov chain Monte Carlo (MCMC) methods. We utilise the SED generator and inference framework Prospector to generate simulated spectra, and train a dataset of 2 × 10 6 spectra (corresponding to a five-parameter SED model) with NPE. We show that SBI—with its combination of fast and amortized posterior estimations—is capable of inferring accurate galaxy stellar masses and metallicities. Our uncertainty constraints are comparable to or moderately weaker than traditional inverse-modelling with Bayesian MCMC methods (e.g. 0.17 and 0.26 dex in stellar mass and metallicity for a given galaxy, respectively). We also find that our inference framework conducts rapid SED inference (0.9–1.2 × 10 5 galaxy spectra via SBI/NPE at the cost of 1 MCMC-based fit). With this work, we set the stage for further work that focuses of SED fitting of galaxy spectra with SBI, in the era of JWST galaxy survey programs and the wide-field Roman Space Telescope spectroscopic surveys.

spectroscopy↗
DOE OSTI2019

Insights into the solvent-assisted degradation of organophosphorus compounds by a Zr-based metal–organic framework

The degradation of a chemical warfare agent simulant using a catalytically active Zr-based metal-organic framework (MOF) as a function of different solvent systems was explored here. Complementary molecular modelling studies suggest that the differences in the degradation rates are related to the increasing size in the nucleophile, which hinders the rotation of the product molecule during degradation. Methanol was identified as an appropriate solvent for non-aqueous degradation applications and demonstrated to support the MOF-based destruction of both Sarin and Soman.

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

Simulation of Compound Flooding Using River‐Ocean Two‐Way Coupled E3SM Ensemble on Variable‐Resolution Meshes

Abstract Coastal zone compound flooding (CF) can be caused by the interactive fluvial and oceanic processes, particularly when coastal backwater propagates upstream and interacts with high river discharge. The modeling of CF is limited in existing Earth System Models (ESMs) due to coarse mesh resolutions and one‐way coupled river‐ocean components. In this study, we present a novel multi‐scale coupling framework within the Energy Exascale Earth System Model (E3SM), integrating global atmosphere and land with interactively coupled river and ocean models using different meshes with refined resolutions near the coastline. To evaluate this framework, we conducted ensemble simulations of a CF event (Hurricane Irene in 2011) in a Mid‐Atlantic estuary. The results demonstrate that the novel E3SM configuration can reasonably reproduce river discharge and sea surface height variations. The two‐way river‐ocean coupling improves the representation of coastal backwater effects at the terrestrial‐aquatic interface that are caused by the combined actions of tide and storm surge during the CF event, thus providing a valuable modeling tool for better understanding the river‐estuary‐ocean dynamics in extreme events under climate change. Notably, our results show that the most significant CF impacts occur when the highest storm surge generated by a tropical cyclone meets with a moderate river discharge. This study highlights the state‐of‐the‐art advancements developed within E3SM for simulating multi‐scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Data-Driven Probabilistic Anomaly Detection for Electricity Market under Cyber Attacks

Information and communication technologies have been widely used in smart grid for efficient operation. However, these technologies are vulnerable to malicious cyber attacks, which may lead to severe reliability and economic issues. Recently, a variety of data-driven anomaly detection approaches have been explored to detect potential cyber attacks in smart grids. In this paper, we researched on the electricity market data aiming to identify anomalies from the locational marginal prices (LMPs) and provide a new indicator for potential cyber attacks in power grids. Specifically, a novel data-driven probabilistic anomaly detection framework is proposed for electricity market, which consists of three major components: long short-term memory (LSTM) based deterministic electricity price forecasting, probabilistic electricity price forecasting and anomaly detection. This framework is tested on a model-based electricity market simulator under two types of cyber attacks, i.e., load redistribution attack (LRA) and price responsive attack (PRA). Numerical results on the simulated LMPs show that the proposed framework is capable of detecting data anomalies over these attacks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Challenges and opportunities in verification and validation of low temperature plasma simulations and experiments

This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Multiscale Investigation of Thermal Transport in UO{sub 2} Fuel using Modeling and Simulation

The thermal conductivity of the UO{sub 2} fuel in light water reactors (LWRs) defines how efficiently heat is transported from the fuel to the coolant, and thus directly influences the reactor efficiency. Unfortunately, the same fission events that generate heat also cause radiation damage and fission product formation that decreases the thermal conductivity. This reduction in thermal conductivity reduces the safety of the reactor in addition to its efficiency, because excess heat generated during an accident cannot be transported from the fuel, resulting in high temperatures that can result in melting. Because thermal conductivity plays such an important role in reactor efficiency and safety, many experiments have been conducted to measure the thermal conductivity of UO{sub 2} and how it changes during reactor operation. Experiments on unirradiated fuel have resulted in a large set of data indicating how the bulk thermal conductivity changes as a function of temperature, which was used to develop a semiempirical model. Additional experiments have been conducted to determine how the thermal conductivity of UO{sub 2} changes as a function of burnup, leading to burnup dependent models fitted to the data. While these empirical models perform well within the conditions at which the data was collected, they cannot be extrapolated outside of those conditions. In addition, while burnup is a measure of the amount of fissioning that has occurred, it cannot account for the impact of the operating conditions on the thermal conductivity. While experimental data provides an empirical view on the value of the thermal conductivity in the fuel during reactor operation, a more mechanistic understanding of how temperature and burnup impact the transport of heat through the fuel would be extremely valuable. Such an understanding would allow for a more accurate prediction of fuel performance, but it could also be used to direct modifications to the fuel that would raise the thermal conductivity or reduce the impact of burnup. In this work we summarize recent efforts to use modeling and simulation at the atomistic and mesoscales to obtain a mechanistic understanding of thermal transport in UO{sub 2}. The thermal conductivity of LWR fuel is one of the most critical material parameters impacting the efficiency and safety of the reactor. Thus, a significant effort is being made to develop a mechanistic understanding of the impact of temperature and burnup on it. In this work, we have provided a short summary of modeling and simulation efforts being carried out at the atomistic and mesoscale to understand this critical property. These research efforts are providing a fundamental understanding of the thermal transport in unirradiated UO{sub 2}, as well as the impact of radiation damage and fission products on thermal transport. These efforts have been focused on UO{sub 2}, but they also are laying the framework for conducted similar simulations to investigate new accident tolerant fuel concepts. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2022

Supporting Hazard Analysis for Wildfire Response Using fmdtools and MIKA

The System Wide Safety (SWS) Safety Demonstrator (SD) Series drives development of an increasingly capable In-Time Aviation Safety Management System (IASMS) focusing on humanitarian applications, starting with wildfire response (SD-1). The goals of this report are to (1) provide an early hazard analysis and mitigation evaluation of wildfire response to support these efforts and (2) provide a demonstration of capabilities of the Fault Model Design Tools (fmdtools) and Manager for Intelligent Knowledge Access (MIKA) tools. fmdtools provides a modeling, simulation, and resiliency analysis framework in which a wildfire response model, the System Modeling and Analysis of Resiliency in Scalable Traffic Management for Emergency Response Operations (SMARt-STEReO), is built. MIKA is an intelligent knowledge manager with several capabilities, including assisting in hazard analysis by extracting and analyzing hazards from historical incident reports. The following topics are covered in the report: Understanding Wildfire Hazard Dynamics. We provide a description and simulated examples of how hazards occur in the SMARt-STEReO model of wildfire response and their effect on its outcome. This provides a common mental model and focuses the analysis presented in the remainder of the report. Wildfire Hazard Identification. MIKA identifies wildfire hazards from three relevant datasets: the ICS-209-PLUS, SAFECOM, and SAFENET. Hazards are manually organized into a taxonomy and MIKA analyzes each hazard’s effects, likelihood, severity, and risk. Evaluating Mitigation Strategies. The SMARt-STEReO wildfire response model built in fmdtools evaluates a subset of identified hazards. Specifically, we simulate the effect of communications faults and equipment faults on operator safety, the effect of changing winds and flammability, and a scenario with multiple ignition points and heavy smoke. Tool Limitations and Usage Considerations. We provide a discussion of appropriate tool use cases as well as limitations and considerations for usage. The tool findings are used to synthesize recommendations for wildfire response operations, which can be captured as part of an IASMS. Key recommendations are as follows: Hazards are identified from a broad spectrum of sources including aircraft subsystems, operational sources, and ground crew operations. Highest risk operational environment hazards identified are Evacuations. The highest risk manned aerial operations hazard categorized is Jumper Operations Mishap. Ground crew hazards that are highest risk are Burns, Cargo Operations Overhead, Dehydration, Entrapment, Falling Objects, Heart Attacks, Heat Exhaustion, Inadequate Training or Certification, Vehicle Breakdown, and Vehicle Collision. Modelled containment failures arise from a mismatch between the difficulty of the firefighting scenario and the capacity (e.g., speed, effectiveness, awareness) of the response. In firefighting scenarios where containment is possible (e.g., because the fire does not spread too quickly), these mismatches can occur because of a change in environmental conditions (e.g., wind, flammability, etc) or because of planning, equipment, or communications faults. Improvements to communications increase the capacity of the firefighting response by reducing the time needed to respond to the fire. While surveillance does not increase this capacity by itself, it increases operator safety by increasing state awareness, enabling firefighters to evade approaching fires. Increasing both has a synergistic effect. In general, these performance and resilience increases generalize over fault scenarios as well as unforeseen changes to circumstances (i.e., wind, aridity, etc.). However, these improvements need to be designed so as not to make the system prone to persistent large-scale communications outages, which can reduce performance.

Hazard analysis↗
DOE OSTI2017

Impact of the strength and spatial distribution of adsorption sites on methane deliverable capacity in nanoporous materials

The methane deliverable capacity of adsorbent materials is a critical performance metric that will determine the viability of using adsorbed natural gas (ANG) technology in vehicular applications. ARPA-E recently set a target deliverable capacity of 315 cc(STP)/cc that a viable adsorbent material should achieve to yield a driving range competitive with incumbent fuels. However, recent computational screening of hundreds of thousands of materials suggests that the target is unattainable. In this work, we aim to determine whether the observed limits in deliverable capacity (similar to 200 cc(STP)/cc) are fundamental limits arising from thermodynamic or material design constraints. Our efforts focus on simulating methane adsorption isotherms in a large number of systems, resulting in a broad exploration of different combinations of spatial distributions and energetics of adsorption sites. All systems were classified into five adsorption scenarios with varying degrees of realism in the manner that adsorption sites are created and endowed with energetics. The scenarios range from methane adsorption on discrete idealized lattice sites to adsorption in metal-organic frameworks with coordinatively unsaturated sites (CUS) provided by metalated catechol groups. Our findings strongly suggest that the ARPA-E target is unattainable, although not due to thermodynamic constraints but due to material design constraints. On the other hand, we also find that the currently observed deliverable capacity limits may be moderately surpassed. For instance, incorporation of CUS in IRMOF-10 is predicted to yield a 217 cc(STP)/cc deliverable capacity. The modified material has a similar to 0.85 void fraction and a heat of adsorption of similar to 15 kJ/mol. This suggests that similar, moderate improvements over existing materials could be achieved as long as CUS incorporation still maintains a relatively large void fraction. Nonetheless, we conclude that more significant improvements in deliverable capacity will require changes in the currently proposed operation conditions.

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

Adsorbents and adsorption models for capture of Kr and Xe gas mixtures in fixed-bed columns

Off-gases produced during the reprocessing of used nuclear fuel (UNF) include 129-I 2 , 3-HHO, 14-CO 2 , 85-Kr, and 135-Xe, which are volatilized out into the off-gas. In order to meet regulatory requirements for reprocessing plant emissions, these gases must be captured and removed from the o?-gas stream prior to o?-gas emission. Of particular interest are the noble gases, Kr and Xe, which can be fairly difficult to remove from the off-gas due to their low chemical reactivity. Thus, this work is focused on utilizing engineered adsorbents, AgZ-PAN and HZ-PAN, to capture Kr and Xe from a mixed-gas stream at relatively low temperatures (191–295 K) and various flow rates (50–2000 mL/min). Isothermal data for Kr and Xe on each adsorbent are analyzed to produce the Langmuir parameters needed to model the mixture adsorption capacities at relevant temperatures using the Extended Langmuir model. Those parameters are then incorporated into a fixed-bed adsorption model developed in this work using the Mulitphysics Object-Oriented Simulation Environment (MOOSE). That model is used to simulate breakthrough times for Kr and Xe in packed columns of AgZ-PAN and HZ-PAN, ranging in length from 6 to 20 in., at relevant temperatures and flow rates. Breakthrough times varied from nearly instantaneous for Kr in AgZ-PAN to 30 h for Xe in HZ-PAN. After the developed model was validated by comparisons with experimental break-through data, the model framework was used to simulate the performance of multiple fixed-bed columns connected in series.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Applying Particle Swarm Optimization and Extended Kalman Filtering to Model Kaplan Generation Dynamics for Hydropower Systems

Variable renewable generation is increasing the need for hydropower plants to provide fast and flexible grid support, which places new demands on plant-level dynamic models used for monitoring, control, and operational decision-making. This need is especially important for hydroelectric systems, where turbine and generator dynamics are strongly coupled, nonlinear, and time-varying, making accurate real-time representation difficult. To address this problem, this paper develops a digital twin (DT) framework for a synchronous generator–Kaplan turbine system using an explicit separation of slow turbine dynamics and fast generator dynamics. The turbine subsystem is represented by a six-coefficient model, whose parameters are identified offline using particle swarm optimization, while the generator subsystem is updated online through an extended Kalman filter for real-time state and parameter estimation. These models are integrated within a closed-loop simulation that includes a proportional–integral–derivative–double-derivative governor and excitation system, allowing the DT to track plant behavior under realistic operating conditions. Unlike prior studies that treat turbine and generator modeling separately or rely mainly on simulated inputs, the proposed framework is validated using real operational data from a hydropower plant. Results show that the DT reproduces terminal voltage, active power, and reactive power with a normalized root mean square error of approximately 5%. This hybrid offline–online formulation constitutes the main contribution of the work, providing an adaptive and practically deployable DT for hydropower systems with direct relevance to control improvement, performance monitoring, and grid-support applications under high renewable penetration.

13 HYDRO ENERGY↗
DOE OSTI2026

Three-Dimensional Grid Visualization for Planning Activities: A Dubai Case Study

National Laboratory of the Rockies (NLR), in collaboration with the Dubai Electricity and Water Authority (DEWA) and Infra-X, has undertaken the Energy Visualization Analysis Project. The aim of this project is to enhance analytical and 3D visualization capabilities for distribution network planning and renewable energy integration. As modern grid continues to evolve with large-scale solar PV deployment and emerging distributed energy resources (DERs), the ability to effectively analyze, visualize, and communicate complex grid behaviors has become increasingly critical. The project focuses on developing empirical use cases based on real distribution feeder data and engineering workflows, ensuring the outcomes are directly aligned with operational environment. Through time-series power flow simulations and nodal hosting capacity analysis, the study quantifies the impacts of high PV penetration on voltage and thermal limits within representative 11 kV feeders. These analyses identify specific nodes and conditions where DER integration challenges arise. Furthermore, a Battery Energy Storage System (BESS) optimization algorithm was applied to determine the optimal size and placement of storage systems that can mitigate network constraints and enhance hosting capacity. The comparative results between base-case and BESS-augmented scenarios clearly demonstrate improvements in network stability and load management efficiency. In parallel, the NLR team developed an immersive 3D visualization framework, enabling interactive exploration of grid simulations using commodity head-mounted display (HMD) systems. This framework transforms conventional 2D simulation data into spatially intuitive visual environments - allowing engineers to analyze feeder conditions, PV hosting potential, and BESS effects in real time. This report represents the first foundational phase in establishing a visualization-driven analytical ecosystem. It provides a methodological foundation for data integration, visualization architecture, and simulation-based decision support, paving the way for large-scale adoption of immersive visualization across DEWA's Smart Grid Initiative, R&D activities, and future network resilience studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2017

Observations and Modeling of Thermal Structure in the Lower Atmosphere and the Upward Propagation of Tides into the Thermosphere

Thermal tides are the atmospheric response to diurnally varying thermal forcing resulting from radiative and convective heat transfer from the surface and from aerosol and gaseous heating within the atmosphere. Tides include sun-synchronous (migrating) waves driven in response to solar heating and additional non-migrating waves resulting from longitudinal variations in the distributions of topography, dust aerosol and water ice clouds. The systematic spatial mapping of temperature over 5 Mars years by the Mars Climate Sounder (MCS) has yielded a well-defined climatology of seasonally-varying temperature structures in the lower atmosphere, from 5 to ~80 km. Tide theory and Mars global circulation model (MGCM) simulations are a fruitful framework for relating temperature observations to thermal forcing by aerosol fields [1]. The analysis of density and temperature fields derived from MAVEN IUVS and NGIMS observations have revealed the presence of predominantly zonal wave 2 and 3 features at altitudes of 100-170 km that are almost certainly non-migrating tides propagating upward from the lower atmosphere [2,3]. In this presentation we will use the MCS climatology and MGCM simulations to relate the density variations seen by MAVEN with the seasonally varying tide activity in the lower atmosphere. Large amplitude perturbations in density are most sensitive to the tide components with the longest vertical wavelengths in temperature, which are well resolved in MCS observations.

Mars Climate Sounder (MCS)↗
DOE OSTI2018

Comparing Genome-Wide Association Study Results from Different Measurements of an Underlying Phenotype

Increasing popularity of high-throughput phenotyping technologies, such as image-based phenotyping, offer novel ways for quantifying plant growth and morphology. These new methods can be more or less accurate and precise than traditional, manual measurements. Many large-scale phenotyping efforts are conducted to enable genome-wide association studies (GWAS), but it is unclear exactly how alternative methods of phenotyping will affect GWAS results. In this study we simulate phenotypes that are controlled by the same set of causal loci but have differing heritability, similar to two different measurements of the same morphological character. We then perform GWAS with the simulated traits and create receiver operating characteristic (ROC) curves from the results. The areas under the ROC curves (AUCs) provide a metric that allows direct comparisons of GWAS results from different simulated traits. We use this framework to evaluate the effects of heritability and the number of causative loci on the AUCs of simulated traits; we also test the differences between AUCs of traits with differing heritability. We find that both increasing the number of causative loci and decreasing the heritability reduce a trait's AUC. We also find that when two traits are controlled by a greater number of causative loci, they are more likely to have significantly different AUCs as the difference between their heritabilities increases. When simulation results are applied to measures of tassel morphology, we find no significant difference between AUCs from GWAS using manual and image-based measurements of typical maize tassel characters. This finding indicates that both measurement methods have similar ability to identify genetic associations. These results provide a framework for deciding between competing phenotyping strategies when the ultimate goal is to generate and use phenotype-genotype associations from GWAS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Effectiveness of Geophysical Monitoring Methods in Detecting Brine and CO 2 Leakage in Underground Sources of Drinking Water

The US DOE National Risk Assessment Partnership (NRAP), funded through the Office of Fossil Energy and the National Energy Technology Laboratory (NETL), is developing methods to evaluate the effectiveness of monitoring techniques to detect brine and CO 2 leakage from legacy wells into underground sources of drinking water (USDW) overlying a CO 2 storage reservoir. As part of the NRAP Strategic Monitoring task group, we have generated Kimberlina Rev. 1.1 of the simulated aquifer impact data (140 simulations) from a model framework based on a hypothetical, compartmentalized, CO 2 storage reservoir near Kimberlina in the southern San Joaquin Basin, California. Brine and CO 2 leakage along damaged wellbores results in subsurface changes in pressure, CO 2 saturation and total dissolved solids (TDS). CO 2 buoyancy allows a significant fraction of leaked CO 2 to reach shallow permeable zones, which is more detectable by surface geophysical sensors. We generated synthetic geophysical data from the simulated aquifer impact data to assess electrical resistivity tomography (ERT), magnetotellurics (MT), gravity and seismic methods for remote and indirect detection of brine and CO 2 leakage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Accurate and scalable graph neural network force field and molecular dynamics with direct force architecture

Abstract Recently, machine learning (ML) has been used to address the computational cost that has been limiting ab initio molecular dynamics (AIMD). Here, we present GNNFF, a graph neural network framework to directly predict atomic forces from automatically extracted features of the local atomic environment that are translationally-invariant, but rotationally-covariant to the coordinate of the atoms. We demonstrate that GNNFF not only achieves high performance in terms of force prediction accuracy and computational speed on various materials systems, but also accurately predicts the forces of a large MD system after being trained on forces obtained from a smaller system. Finally, we use our framework to perform an MD simulation of Li 7 P 3 S 11 , a superionic conductor, and show that resulting Li diffusion coefficient is within 14% of that obtained directly from AIMD. The high performance exhibited by GNNFF can be easily generalized to study atomistic level dynamics of other material systems.

Chemistry↗
DOE OSTI2024

Initial energy-momentum to final flow: A general framework for heavy-ion collisions

The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Furthermore, we confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this framework, it is possible to separate the effects of early time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as 𝑝−𝑝, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.

Relativistic heavy-ion collisions↗