Search NASA⌕ Search

SEARCH · Search NASA

Results for “semi-empirical 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 37 records · Page 2

Integration of the fundamental knowledge on solvent-packing interactions into the multiscale framework for column scale design and optimization

The interfacial area, also known as the effective mass transfer area, is a key factor for determining the mass transfer for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Most of the widely used empirical and semi-empirical models for interfacial area were derived indirectly through absorption mass transfer with simplifications based on fast chemical kinetics. This report presents the comprehensive unique multiscale approach to develop a surrogate model for effective mass transfer area in structured packed columns that accounts local hydrodynamics as well as variation in physical properties, and changes in solid surface characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-Domain Vortex Induced Vibration Modeling of Reference Dynamic Power Cable for the Gulf of Maine

Vortex-induced vibration (VIV) is a phenomenon known to decrease the fatigue life of dynamic power cables used in floating offshore wind systems through increased bending loading cycles. A variety of modeling solutions have been proposed to study this, but none enable fully-coupled simulations leveraging open-source tools. To address this, a time-domain VIV model has been successfully implemented into the open-source mooring dynamics model MoorDyn. The semi-empirical VIV model describes a lift force acting on segments of a flexible cylinder, making it well suited for implementation into MoorDyn. The new capability successfully predicted the frequencies and magnitudes of strain in both steady and oscillating flows when compared to the original validation results and matched peak spectral response when compared to experimental results, verifying successful implementation. MoorDyn with VIV was then leveraged to simulate a 15 MW floating turbine with a dynamic power cable in the Gulf of Maine using OpenFAST for six return periods of combined wind, wave, and current conditions. Increased curvature and tension fluctuations in higher flow speeds were observed when simulating VIV. Maximum tensions and curvatures also increased, with the 500 year conditions violating the curvature factor of safety of 2.0. These results highlight the importance of considering VIV when designing dynamic power cables for the Gulf of Maine, and indicate that cost effective mitigation strategies should be explored for the area. They also demonstrate the utility of this new modeling capability, which provides the first open-source tool for fully coupled time-domain floating offshore wind simulations with dynamic power cable VIV.

17 WIND ENERGY↗

Bayesian parameter estimation and evaluation of the K -ω shear stress transport model for plane impinging jets

Numerical simulations with semi-empirical turbulence models are commonly used to model impinging jets, often used for cooling solid surfaces. In this work, the constants in the k-ω shear stress transport model in ANSYS FLUENT are calibrated to experimental velocity and heat transfer data for a plane turbulent impinging air jet to determine if Kennedy-O'Hagan calibration (Kennedy and O'Hagan 2001 J. R. Stat. Soc. B 63 425–64) can improve predictions of near-surface velocities and surface Nusselt numbers for similar flows. Impinging jets have been proposed to cool the target plates of the divertor in future magnetic fusion energy reactors, where simulations are used to estimate divertor performance. The flat-plate divertor (Wang et al 2009 Fusion Sci. Technol .56 1023–7) uses a plane jet of helium issuing from a B = 0.5 mm slot to cool a surface with radius of curvature of 44 B at a distance 4 B from the slot. Predictions from the calibrated numerical model are compared with independent experimental data at different flow conditions, as well as surface temperature data for a flat plate divertor test section. The contribution of this work is evaluation of the accuracy of a calibrated turbulence model for modest extrapolations in flow geometry and flow conditions for a plane impinging jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust Combined Heat and Hybrid Power (CHHP) for High Electrical Efficiency Cogeneration

Georgia Tech (Prime Recipient), the University of Texas at El Paso (UTEP) and the National Energy Technology Laboratory (NETL) investigated a hybrid fuel cell/ gas turbine system concept as a combined heat and hybrid power (CHHP) system for both robust and high power-to-process heat ratio cogeneration. The novelty of the proposed system entailed the distinct, elevated electrical efficiencies it maintains while simultaneously supporting a broad span of heating needs (e.g., supply temperatures) demanded across variable heat loads. The scope included: 1) leveraging a pioneering national lab facility configured for dynamic system operability development of hybrid fuel cell/gas turbine cycles; 2) enabling technology development to adjust and modulate the quality and quantity of thermal supply to bottoming heat loads via novel extreme temperature gas bypass valves. Hybrid fuel cell/gas turbine systems have primarily been reduced-to-practice in a constrained (e.g., initial proof-of-concept) manner and have still demonstrated considerable electrical efficiencies. However, these pre-pilot systems have focused upon electrical efficiencies and electrical power generation as the exclusive energy demand. Such hybrid systems had not been extensively researched or developed for flexible and variable operation consisting of both power and heat demands; however, these variable combined power and heat demands are characteristic of many types of manufacturers such as animal/poultry processing, bakeries and milk/flour/pastry manufacturing, textile mills, and electrochemical processing. Commercially, developing the system into a working combined heat and power system benefits these types of manufacturers by allowing them to meet their power and heat demands at a lower cost, higher efficiency, and/or through onsite generation. Therefore, the technical scope of this project was largely to study and facilitate these hybrid systems as combined heat and hybrid power (CHHP) systems that include dynamic operability for variable heat and power loads and/or grid dynamics for various types of manufacturers. Simulation results were used to predict the performance of the CHHP system and conceptually develop it to achieve desired dynamic operability. Experimentally, the primary goal was to design, manufacture, and experiment upon a high-temperature bypass valve. Experimental data included air mass flow rates through the valve orifice when the valve was changed to variable extent between fully closed and fully open. The experimental data was then used to create a semi-empirical computational model of the bypass valve. Concluded simulation goals for the research included developing computational heat exchanger models for the hybrid system inclusive of the bottoming heat exchanger and the recuperative heat exchanger, and then combining the computational recuperator model with the computational valve model. Afterwards, the computational models were then integrated to predict the dynamic operation of hybrid fuel cell/gas turbine cycles throughout a design space and reporting such. The scope stated in the preceding paragraph was packaged into five specific goals: 1) enabling the simulation of dynamic combined heat and power through the creation of computational, modular heat exchanger models; 2) simulation and exploration of the CHHP system’s performance by integrating the heat exchanger models with the national lab’s pre-existing hybrid system (computational) simulation, but without the recuperator bypass valve concept in order to initially determine how the (baseline) system behaves and can be controlled in order to meet variable heat and power demands; 3) development and initial deployment of the high-temperature recuperator bypass valve technology in order to confirm and characterize the approach; 4) usage of the experimental data for the valve to create a semi-empirical computational model for the bypass valve which could then be combined with the heat exchanger computational models; 5) repeat of the second task of simulating and exploring the system’s performance, but this time including the bypass valve to resolve its efficacy. Tasks were successfully completed, and the general notion of flexibly operating, high electrical efficiency CHHP was further corroborated. Supportive details are provided in the report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Daily operational impacts on battery degradation in heavy-duty electric drayage trucks

Battery aging is a critical factor influencing the performance, longevity, and cost of ownership of battery electric trucks (BETs). This paper presents a comprehensive evaluation of battery aging for two Li-ion battery chemistries, Nickel-Manganese-Cobalt (NMC) and Lithium-Iron-Phosphate (LFP), accounting for both cycling and calendar aging. In contrast to traditional methods that rely on simplified linear degradation models based on manufacturer-provided data, this study employs semi-empirical aging models calibrated to experimentally collected data. The models are integrated into a detailed vehicle simulation environment, enabling a comprehensive assessment of battery degradation under realistic operating conditions. A case study focusing on heavy-duty electric drayage truck operations in the Port of Savannah, GA, is presented to illustrate the impact on battery pack lifespan of: seasonal variations, daily operational activities, charging strategies, and battery storage conditions. The results illuminate the significance of the battery pack’s state of charge during stationary periods, such as overnight storage or weekend parking, on battery degradation and its potential implications for long-term vehicle viability. Additionally, the study explores how different operational and environmental factors affect battery degradation, offering critical insights into best battery charging and storage practices. Our results demonstrate that LFP outperforms NMC in terms of years of useful life; however, by utilizing charging strategies that minimize the amount of time the battery spends resting at high levels of state-of-charge, the lifespan of the battery pack that uses NMC can nonetheless be increased by more than a factor of two.

25 ENERGY STORAGE↗

Modeling the effects of active wake mixing on wake behavior through large-scale coherent structures

The use of active wake mixing (AWM) to mitigate downstream turbine wakes has created new opportunities for reducing power losses in wind farms. However, many current analytical or semi-empirical wake models do not capture the flow instabilities that are excited through the blade pitch actuation. In this work, we develop a framework, which accounts for the impacts of the large-scale coherent structures and turbulence on the mean flow, for modeling AWM. The framework uses a triple-decomposition approach for the unsteady flow field and models the mean flow and fine-scale turbulence with a parabolized Reynolds-averaged Navier–Stokes (RANS) system. The wave components are modeled using a simplified spatial linear stability formulation that captures the growth and evolution of the coherent structures. Comparisons with high-fidelity large eddy simulations (LESs) of the turbine wakes showed that this framework was able to capture the additional wake mixing and faster wake recovery in the far-wake regions for both the pulse and helix AWM strategies with minimal computational expense. In the near-wake region, some differences are observed in both the RANS velocity profiles and initial growth of the large-scale structures, which may be due to some simplifying assumptions used in the model.

17 WIND ENERGY↗

Modelling Tritium Production and Release at High-Energy Proton Accelerators

Tritium is a well-known byproduct of particle accelerator operations. To keep levels of tritium below regulatory limits, tritium production is actively monitored and managed at Fermilab. We plan to study tritium production in the targets, beamline components, and shielding elements of the Fermilab facilities such as NuMI, BNB, and MI-65. To facilitate the analysis, we construct a simple model and use three Monte Carlo radiation codes, FLUKA, MARS, and PHITS, to estimate the amount of tritium produced in these facilities. The analysis could also serve as an intercomparison between these code results related to tritium production. To assess the actual amounts of tritium that would be released from various materials, we employ a semi-empirical diffusion model. The results of this analysis are compared to experimental data whenever possible. This approach also helps to optimize proposed target materials with respect to the tritium production and release.

Georgobiani, Dali [Fermilab]↗

Next Generation Heat Transfer Fluids for Two-Phase Immersion Cooling of Data Centers

The purpose of this study is to evaluate the performance of next generation dielectric fluids in a Two-Phase Immersion Cooling (2PIC) system, which was designed for use in data centers. Hence, this report contains the performance evaluations of a new developmental dielectric fluid, Opteon™ 2P50, in a commercially available small-scale 2PIC system under typical and off-design range of operating conditions. Accordingly, ambient temperature and thermal loads were varied to simulate different ambient conditions. Additionally, this research report describes the development of a semi-empirical lumped model to predict the energy efficiency of the 2PIC system using Opteon™ 2P50 across a wide range of conditions. The model aims to offer a comprehensive understanding of the system’s efficiency and potential improvements. The outcomes of this study are expected to contribute to the adoption of sustainable 2PIC cooling technologies in data centers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An experimentally validated electro-thermal EV battery pack model incorporating cycle-life aging and cell-to-cell variations

Lithium-ion batteries are used in a wide variety of applications. To meet the power and energy demands of these applications battery packs are composed of hundreds to thousands of cells. The electrical and thermal interactions between cells introduce additional complexity in the pack dynamics. To capture these effects, a battery pack model composed of 192 cells based on a first-generation (2012) Nissan Leaf battery pack is developed in MATLAB/Simulink/Simscape. Here, with this model, we simulate the electrical dynamics (using a first-order equivalent-circuit model), the thermal dynamics (using a first-order lumped-parameter thermal model), and the aging dynamics (using a semi-empirical severity factor-based model) of every cell in the pack and we also create a pack thermal model that explicitly captures the heat exchange between the modules, and the cells contained within, during operation. The models are calibrated and validated, both at the cell and pack level, with experimental data. Two different case studies of this pack model are investigated. In the first case study, an initial, normally-distributed, cell-to-cell capacity variation is introduced and its effect on the pack voltage and module temperatures is studied. In the second case study, we deliberately insert cells with lower than nominal capacity into the pack and we investigate how this type of initial cell-to-cell capacity variation affects the pack’s ability to deliver energy over time. Finally, we also study how parallel-connected cells can reduce the effects of cell-to-cell variations at the expense of increased aging of the pack overall.

25 ENERGY STORAGE↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗

Learning thermodynamically constrained equations of state with uncertainty

Numerical simulations of high energy-density experiments require equation of state (EOS) models that relate a material’s thermodynamic state variables—specifically pressure, volume/density, energy, and temperature. EOS models are typically constructed using a semi-empirical parametric methodology, which assumes a physics-informed functional form with many tunable parameters calibrated using experimental/simulation data. Since there are inherent uncertainties in the calibration data (parametric uncertainty) and the assumed functional EOS form (model uncertainty), it is essential to perform uncertainty quantification (UQ) to improve confidence in EOS predictions. Model uncertainty is challenging for UQ studies since it requires exploring the space of all possible physically consistent functional forms. Thus, it is often neglected in favor of parametric uncertainty, which is easier to quantify without violating thermodynamic laws. This work presents a data-driven machine learning approach to constructing EOS models that naturally captures model uncertainty while satisfying the necessary thermodynamic consistency and stability constraints. We propose a novel framework based on physics-informed Gaussian process regression (GPR) that automatically captures total uncertainty in the EOS and can be jointly trained on both simulation and experimental data sources. A GPR model for the shock Hugoniot is derived, and its uncertainties are quantified using the proposed framework. We apply the proposed model to learn the EOS for the diamond solid state of carbon using both density functional theory data and experimental shock Hugoniot data to train the model and show that the prediction uncertainty is reduced by considering thermodynamic constraints.

Sharma, Himanshu (ORCID:000900050235718X)↗

A computational approach to complex junctions and interfaces

Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.

2D materials↗

A model for oblique impacts on material surfaces

Many practical situations of material damage, wear, and erosion involve collisions between small particles and surfaces at inclined angles. While there are many well-validated models of normal incidence impact situations, elastic-plastic models for oblique incidence impact events are lacking. Here the finite element method is used to predict the normal and tangential coefficient of restitution in oblique impacts for hard, elastic spheres impacting an elastic-perfectly plastic material surface. The proposed model covers various impact angles ranging from 0° to 45°, within a limiting impact velocity below which the effects of heating are negligible. The normal coefficient of restitution follows power-laws with respect to normalized values of the impact velocity. Interestingly, the tangential coefficient of restitution follows a linear relationship with impact velocity. Together, these results provide a semi-empirical set of equations predicting oblique impact rebounds (both velocity and trajectory) for a wide range of conditions and material properties, with which experimental results can be rapidly interpreted. As a result, Laser-Induced Particle Impact Test (LIPIT) data are also presented for aluminum particles impacting aluminum substrates, at impact angles of 25° and 40°; the results compare favorably with the model and validate the general use of such models for the analysis of experimental data.

Coefficient of restitution↗

Predictive numerical modeling of plasma-induced surface roughness and wettability evolution in LM-PAEK/CF tape

Plasma surface modification effectively enhances adhesion in thermoplastic composites, yet its impacts on high-performance polymers like low-melting polyaryletherketone (LM-PAEK) remain inadequately quantified. Here, this study integrates experimental analysis and numerical modeling to characterize surface roughness and wettability changes in LM-PAEK/carbon fiber composites treated with atmospheric plasma. Atomic Force Microscopy quantified surface topography (n = 10 per condition for contact angles), while static contact-angle assessments measured wettability. Roughness rapidly increased from ∼0.2 nm to 1.6 nm, and contact angle reduced from ∼90° to 24°, both stabilizing after 25–30 s of exposure. A semi-empirical, physics-informed framework was calibrated to these data, coupling surface chemistry via the Owens–Wendt decomposition with topography via the Wenzel roughness factor, and evaluated using out-of-sample (cross-validated) tests, while static contact angle assessments measured wettability. Numerical predictions matched experimental results closely (RMSE <5%, R 2 > 0.95). Incorporating material-specific parameters, the calibrated model supports plasma-treatment optimization and provides quantitative guidance for improving interfacial adhesion in thermoplastic composite manufacturing.

Atomic Force Microscopy↗

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗

Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization

The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.

Bayesian optimization↗

Scaling Ensembles of Data-Intensive Quantum Chemical Calculations for Millions of Molecules

Deep learning models are efficient computational tools that can accelerate the inverse design of molecules with desired functional properties by generating predictions at a fraction of the time required by traditional quantum chemical approaches. To ensure that a model maintains accuracy and transferability across broad regions of the chemical space explored during the inverse design, it must be trained on massively large volumes of simulation data. This requires running large-scale ensemble quantum chemical calculations on high-performance computing (HPC) systems for data collection. However, the efficient execution of such large ensemble calculations and the management of large volumes of output data require tools that can judiciously utilize computational resources and manage metadata overhead on the file system. Therefore, we present a high-performance, scalable, ensemble management framework for performing data-intensive quantum chemical electronic structure calculations for organic molecules. This framework provides abstractions to plug different ab initio, first principles, and first principles-based semi-empirical methods and executes them efficiently at large scale on HPC systems. It dynamically distributes tasks to resources and uses tiered storage for managing large collections of files. We employed this framework to process over ten million organic molecules and generate open-source datasets that provide UV-vis absorption spectra by running time-dependent density-functional tight-binding calculations. It is the largest database containing molecular optical spectra that were simulated with quantum chemical methods in a consistent manner.

Mehta, Kshitij↗