Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemistry 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 199 records · Page 11

A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests

Abstract Dry deposition is the second largest tropospheric ozone (O 3 ) sink and occurs through stomatal and nonstomatal pathways. Current O 3 uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O 3 fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O 3 fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O 3 deposition and its drivers. We show that GEOS‐Chem overpredicts observed O 3 fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O 3 deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O 3 flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O 3 deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O 3 fluxes; and (c) the CTM big‐leaf treatment overestimates O 3 ‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O 3 predictions.

Vermeuel, Michael P. [Department of Soil, Water, a↗

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Updated Production Inventory for Lithium-Ion Battery Anodes for the GREET® Model, and Review of Advanced Battery Chemistries

The Greenhouse gases, Regulated Emissions, and Energy use in Technologies (GREET®) model considers lithium-ion batteries with multiple anode materials. Synthetic graphite is the primary anode material used in the previous GREET versions, even as the model offered options to choose a lithium anode and/or a blended anode (blend of synthetic graphite and silicon). Yet, the inventory (material and energy flows) considered for these anodes is dated, and the anode options do not consider natural graphite, which is another important anode material for lithium-ion batteries. This report documents the material and energy flows for natural graphite anode production from raw material extraction to anode production – as incorporated in the updated GREET model. We also present a brief literature review on the current state of inventory for the other three anodes (synthetic graphite, silicon, and lithium), as well as updates made in the recent GREET model on material and energy flows associated with their respective production. Finally, this study provides a summary of advanced battery systems that may be alternatives to LIBs for use in future electric vehicles.

25 ENERGY STORAGE↗

Cloud Atlas: Unraveling the Vertical Cloud Structure with the Time-series Spectrophotometry of an Unusually Red Brown Dwarf

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. Additionally, in the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH 4 chemistry, our models also reproduce the redder $J-{K}_{{\rm{s}}}$ color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Droplet burning of n-butyl acetate/n-heptane mixtures in the standard atmosphere: experiments and numerical modeling

Blending petroleum fuels with oxygenates is a common approach to stem the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that n-butyl acetate (BA, C6H12O2 , boiling point of 399K) could be a viable oxygenate additive to diesel fuel. In this application it is important to determine the influence of the fractional amount of BA on burning of the mixture. This presentation considers this problem from an experimental and computational approach using n-heptane (C7H14 , boiling point of 372K) as an essentially single component surrogate for diesel fuel. The burning configuration considered was ostensibly spherical symmetry as promoted by burning droplets under conditions where forced and natural convection effects were minimized. In this configuration the droplet and flame are concentric spheres and soot aggregates are trapped in a shell structure between the droplet and flame. Measurements were made of the droplet, flame and soot shell diameters through video imaging of the burning process. The experimental results were used to validate a detailed numerical model of the spherically symmetric droplet burning process that incorporated a consistent combustion kinetic mechanism for BA/heptane mixtures, comprised of 402 species and 16,872 reactions inclusive of soot chemistry. A model for soot formation was also included in the simulation to predict formation of the soot shell. The simulations agreed well with measured droplet and flame diameters. It was found both experimentally and computationally that adding BA to heptane had a minimal effect on the mixture burning rate, while the flame was positioned closer to the droplet surface and the simulated soot volume fraction decreased. The implications of these results are discussed for using BA as a potential additive to diesel fuel.

09 BIOMASS FUELS↗

DRACO: An Overview [Slides]

DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗

Flame Chemistry Workshop: a perspective on challenges and strategic actions in combustion experiments and chemical kinetics modeling

Continued progress in the development of predictive models for combustion chemistry—including ignition and flame behavior, species evolution, and combustion system performance—relies on overcoming enduring and emerging challenges in experimental measurements, theoretical formulations, and chemical kinetics mechanism construction. As combustion science continues to coincide with advances in sustainable fuels development, plasma technologies, and automated modeling capabilities, the need for coordinated, community-driven strategies is essential. The Flame Chemistry Workshop (FCWS), held biennially before the International Symposium on Combustion, serves as a dedicated platform to identify, consolidate, and address these challenges in a structured and collaborative manner. This perspective arises from discussions at the 7th FCWS in Milan, Italy (2024), and presents a collective view of the critical barriers currently limiting progress. Across the five technical domains discussed during the 7th FCWS – sustainable fuels combustion, advanced diagnostics for combustion measurements, experiments and modeling in plasma combustion, artificial intelligence and automated methods for theory and mechanisms generation, and chemical kinetic models—a series of persistent and emerging scientific challenges were identified, highlighting the need for deeper integration between three areas: theory, experiments, and modeling. In conclusion, the present article concisely describes present challenges that were identified in each of the technical domains in an effort to streamline and coordinate solutions to accelerate progress in combustion science.

Chemical kinetics↗

Slag Management of Carbon Feedstock Used in Gasification FY17 TCF between Eastman Chemical and NETL (Final Report)

NETL, in collaboration with Eastman Chemical Corp (Kingsport, TN), targeted the joint evaluation and refinement of a slag management program developed by NETL (NETL IP, patent applied for) to control gasifier slag properties formed from ash in individual or combinations of different coal-based carbon feedstock materials and specific slag modifiers used in gasification. Slag forms from mineral or organic-metallic impurities in the carbon feedstock that melt at the high temperatures of gasification, coalescing to form a liquid (slag) that must flow from the gasifier. Slag additives, also known as slag modifiers, are often added to the carbon feedstock to modify the ash chemistry and achieve desired slag flow properties at target temperatures used in gasification. If slag does not flow from a gasifier, its buildup within the gasification chamber will cause either gasifier shutdown to physically remove the slag (consuming time and damaging the protective refractory lining), or the gasification temperature must be increased to lower slag viscosity so it can flow from the gasifier, leading to greater refractory wear/corrosion, along with increased fouling and downstream component wear. Because of these issues, gasification is typically targeted at the lowest temperature possible where slag will flow in a controlled manner (a temperature range called T100 for the lowest temperature, with T250 being the highest temperature where the slag has achieved the desired fluidity). That temperature range also needs to yield a desired syngas composition, highlighting the importance of knowing and being able to control carbon feedstock slag viscosity and the temperature of thermal breakdown of the carbon feedstock. Slag viscosity varies with changes in the gasification chamber temperature, oxygen partial pressure, and slag chemistry. The slag model work with Eastman Chemical targeted predicting and controlling molten coal ash slag behavior at specific gasification conditions. The final goal of the refined model was to develop/refine the NETL slag model so a gasifier operator would be able to predict slag behavior prior to purchasing or using a carbon feedstock and control slag properties once a carbon feedstock is being gasified, allowing a gasifier user to determine what mineral additives or blends of different carbon materials would be needed to achieve desired slag properties at targeted gasification conditions. For a given carbon feedstock and desired gasification condition, the actual gasification cost and slag behavior would be able to be predicted using the slag management program in advance of purchasing or using the carbon feedstock – allowing accurate control of slag ash behavior in real time. As part of the TCF project work with Eastman Chemical, the viscosity of several ash compositions of interest to them would be measured and used to populate the slag model, along with a large database of historical slag temperature, chemistry, and viscosity information they possess – information that would in-crease the predictive accuracy of the slag management model with slag ash chemistry of interest to them. As part of the model evaluation and refinement, the predictive model and its database will be installed on computers at Eastman Chemical and used to manage slag properties in a gasifier during real-time operation of their gasifier. In summary, success of the slag management program at Eastman Chemical would support its use by other gasifier operators to manage slag chemistry and properties over other desired gasification temperature ranges, allowing greater process control and efficiency of their carbon feedstock conversion during gasification.

01 COAL, LIGNITE, AND PEAT↗

Sensitivity of Regional WRF‐Chem Air Quality and Weather Simulations to Biomass‐Burning Emission Data Sets: A Case Study of the Impact of Canadian Wildfire on the US°

This study focuses on the period from June 26 to 29, 2023, when record‐breaking Canadian wildfires severely impacted air quality in the Midwest United States. Using the Weather Research and Forecasting Model with Chemistry (WRF‐Chem) and four biomass‐burning data sets (Fire Inventory from NCAR version 1, Fire Inventory from NCAR version 2.5, Quick Fire Emissions Data set [QFED], and Regional ABI‐VIIRS Emission), we analyzed aerosol transport from Canada to the US and assessed the model's accuracy in predicting PM 2.5 , O 3 , CO and aerosol weather feedback. Model simulations were compared with ground‐based and remote sensing observations as well as field measurements from the Community Research on Climate and Urban Science (CROCUS) project. Our findings show that the movement of a low‐pressure system from the Great Lakes to the Atlantic, combined with the high‐pressure system over the Atlantic, caused the transport of aerosols from Canadian wildfires to the US. Results show WRF‐Chem significantly underestimated key atmospheric components: aerosol optical depth (AOD) by over 50%, PM 2.5 by 65%–90% and peak O 3 concentrations by 50%–55% across four biomass burning data sets. Additionally, CO and NO 2 concentrations were underpredicted. The substantial underestimation of PM 2.5 led to an overestimation of temperature by up to 3.6 °C primarily due to excessive downward shortwave radiation, which resulted from the underestimation of direct aerosol effects and an increase in sensible heat flux. Among the biomass‐burning data sets, QFED produced the most accurate AOD and PM 2.5 predictions due to improved wildfire emission estimates, leading to a 1.0 to 1.5 °C reduction in temperature overestimation during the daytime. These findings underscore the need for improving wildfire emission estimates for trace gases and aerosols to enhance air quality and weather feedback predictions.

WRF-chem model↗

A Comprehensive Machine Learning Model for Metal–Ligand Binding Prediction: Applications in Chemistry and Biology

A machine-learning (ML) model that predicts metal–ligand binding constants was developed using the open-source Chemprop software. The model was trained on over 30,000 experimental log K 1 values, which include both protonation and metal–ligand stability constants, comprising over 3500 ligands and 10 2 metal ions from 73 total elements, thus generalizing beyond existing limited approaches, which focus only on specific metals or ligand families. The best-performing model included a combination of SMILES-based molecular representations along with descriptors for the metal ion and experimental conditions. It had an external test R 2 value of 0.942, and MAE value of 0.834. A “SMILES-only” simpler version also produced accurate predictions and preserved the binding trends, serving as a quick and easily accessible alternative for users without computational expertise. The SMILES-only model performed comparably to density functional theory (DFT) calculations but utilized a fraction of the computational resources. The model was successfully applied across diverse domains, including bioinorganic chemistry, heavy metal remediation, and sensor development and demonstrated its effectiveness as a rapid and reliable screening tool for both academic and industrial uses.

Ligands↗

Global nitrogen and sulfur deposition mapping using a measurement–model fusion approach

Abstract. Global reactive nitrogen (N) deposition has more than tripled since 1860 and is expected to remain high due to food production and fossil fuel consumption. Global sulfur emissions have been decreasing worldwide over the last 30 years, but many regions are still experiencing unhealthily high levels of deposition. We update the 2010 global deposition budget for reactive nitrogen and sulfur components with new regional wet deposition measurements from Asia, improving the ensemble results of 11 global chemistry transport models from the second phase of the United Nations Economic Commission for Europe's Task Force on Hemispheric Transport of Air Pollution (HTAP II). The observationally adjusted global N deposition budget is 114.5 Tg N, representing a minor increase of 1 % from the model-only derived values, and the adjusted global sulfur deposition budget is 88.9 Tg S, representing a 6.5 % increase from the modeled values, using an interpolation distance of 2.5∘. Regionally, deposition adjustments can be up to ∼ 73 % for nitrogen and 112 % for sulfur. Our study demonstrates that a global measurement–model fusion approach can improve N and S deposition model estimates at a regional scale, with sufficient availability of observations; however, in large parts of the world, alternative approaches need to be explored. The analysis presented here represents a step forward toward the World Meteorological Organization's goal of global fusion products for accurately mapping harmful air pollution deposition.

Rubin, Hannah J. (ORCID:0000000182226924)↗

Validating a microphysical prognostic stratospheric aerosol implementation in E3SMv2 using observations after the Mount Pinatubo eruption

Abstract. This paper describes the addition of a stratospheric prognostic aerosol (SPA) capability – developed with the goal of accurately simulating sulfate aerosol formation and evolution in the stratosphere – in the Department of Energy (DOE) Energy Exascale Earth System Model, version 2 (E3SMv2). The implementation includes changes to the four-mode Modal Aerosol Module microphysics in the stratosphere to allow for larger particle growth and more accurate stratospheric aerosol lifetime following the Pinatubo eruption. E3SMv2-SPA reasonably reproduces stratospheric aerosol lifetime, burden, aerosol optical depth, and top-of-atmosphere flux when compared to remote sensing observations. E3SMv2-SPA also has close agreement with the interactive chemistry–climate model CESM2-WACCM (Community Earth System Model version 2–Whole Atmosphere Community Climate Model) – which has a more complete chemical treatment – and the observationally constrained, prescribed volcanic aerosol treatment in E3SMv2. Global stratospheric aerosol size distributions identify the nucleation and growth of sulfate aerosol from volcanically injected SO2 from both major and minor volcanic eruptions from 1991 to 1993. The modeled aerosol effective radius is consistently lower than satellite and in situ measurements (max differences of ∼ 30 %). Comparisons with in situ size distribution samples indicate that this simulated underestimation in both E3SMv2-SPA and CESM2-WACCM is due to overly small accumulation and coarse-mode aerosols 6–18 months post-eruption, with E3SMv2-SPA simulating ∼ 50 % of the coarse-mode geometric mean diameters of observations 11 months post-eruption. Effective radii from the models and observations are used to calculate offline scattering and absorption efficiencies to explore the implications of smaller simulated aerosol size for the Pinatubo climate impacts. Scattering efficiencies at wavelengths of peak solar irradiance (∼ 0.5 µm) are 10 %–80 % higher for daily samples in models relative to observations through 1993, suggesting higher diffuse radiation at the surface and a larger cooling effect in the models due to the smaller simulated aerosol; absorption efficiencies at the peak wavelengths of outgoing terrestrial radiation (∼ 10 µm) are 15 %–40 % lower for daily samples in models relative to observations, suggesting an underestimation in stratospheric heating in the models due to the smaller simulated aerosol. These potential biases are based on aerosol size alone and do not take into account differences in the aerosol number. The overall agreement of E3SMv2-SPA with observations and its similar performance to the well-validated CESM2-WACCM makes E3SMv2-SPA a viable alternative to simulating climate impacts from stratospheric sulfate aerosols.

Brown, Hunter York (ORCID:0000000218147874)↗

The ballad of LLM agents: philosophical reasoning for chemistry

Large language models (LLMs) show remarkable potential for scientific reasoning but often produce unreliable or scientifically unactionable outputs when faced with multi-step logic, domain grounding, and interpretability challenges, especially in complex fields like chemistry and materials science. Here, we introduce a framework of philosophical reasoning agents, inspired by canonical thinkers such as Socrates, Descartes, Kant, and Hume, to guide LLM behavior via structured prompt engineering. These agents embody distinct reasoning paradigms (dialectical inquiry, deductive logic, rule-based judgment, and empirical validation) and are evaluated across multiple chemistry subdomains, physical, analytical, general, inorganic, and organic chemistry, using the ChemBench benchmark. Our agentic prompting approach yields substantial accuracy gains on open-ended numerical chemistry questions, with gains of +11.5 percentage points for GPT-4o with Hume, +4.5 percentage points for GPT-5 with Kant, and +21.8 percentage points for GPT-5.1 with Socrates at the strict 1% error threshold, relative to the corresponding base models. Beyond accuracy, we observe benchmark-level model–agent performance patterns, suggesting that different prompting styles interact differently with each base model. These findings demonstrate that embedding philosophy-of-science principles into multi-agent frameworks can improve and produce interpretable, adaptive, and domain-aligned scientific LLMs.

Harb, Hassan [Argonne National Laboratory (ANL), A↗

Basic Research Needs for Transformative Manufacturing (Brochure)

Manufacturing is central to the nation’s prosperity and security. Manufacturing currently represents about 12% of the gross domestic product, provides nearly 13 million jobs, and accounts for about 25% of energy use. The nation’s economy relies heavily on wide-ranging manufacturing sectors - all of which share common challenges including data issues, lack of physics and chemistry-based models across scales, and resource constraints in a global environment. Furthermore, there are many hurdles that must be overcome to move basic science innovations to market. Addressing broad-ranging challenges demands a basic-science strategy that underpins applied research activities. This strategy would accelerate innovation and transform manufacturing. A Basic Research Needs workshop for Transformative Manufacturing was held in March 2020. The focus of the workshop was to identify the basic science research priorities that could accelerate innovation to transform manufacturing in the future. This was the first workshop of its kind to examine how basic energy science can drive manufacturing forward and innovate new ways to manufacture goods. Five Priority Research Directions were identified that address these science challenges: (1) innovative synthetic approaches to enable scalable assembly of matter, (2) computational methods and theoretical models to transform how manufacturing processes are controlled, (3) new characterization tools that can handle the necessary complexity, scales, and processing speeds to meet manufacturing needs, (4) new science to address opportunities relevant to sustainable and energy-efficient manufacturing, and (5) foundational approaches to co-design of materials, process, and products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗