Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tensor factorizations”

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.

393 records · Page 22

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗

Simulating spin dynamics of supersolid states in a quantum Ising magnet

Motivated by a recent experimental study on the quantum Ising magnet K 2 Co(SeO 3 ) 2 that presented spectroscopic evidence of zero-field supersolidity (Chen et al., arXiv:2402.15869), we simulate the excitation spectrum of the corresponding microscopic XXZ model for the compound, using the recently developed excitation ansatz for infinite projected entangled-pair states. Here, we map out the ground state phase diagram and compute the dynamical spin structure factors across a range of magnetic field strengths, focusing especially on the two supersolid phases found near zero and saturation fields. Our simulated excitation spectra for the zero-field supersolid “Y” phase are in excellent agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TU-CD-BRB-05: Radiation Damage Signature of White Matter Fiber Bundles Using Diffusion Tensor Imaging (DTI)

Purpose: To develop an automated and scalable approach and identify temporal, spatial and dosimetric patterns of radiation damage of white matter (WM) fibers following partial brain irradiation. Methods: An automated and scalable approach was developed to extract DTI features of 22 major WM fibers from 33 patients with low-grade/benign tumors treated by radiation therapy (RT). DTI scans of the patients were performed pre-RT, 3- and 6-week during RT, and 1, 6 and 18 months after RT. The automated tractography analysis was applied to 198 datasets as: (1) intra-subject registration of longitudinal DTI, (2) spatial normalization of individual-patient DTI to the Johns Hopkins WM Atlas, (3) automatic fiber tracking regulated by the WM Atlas, and (4) segmentation of WM into 22 major tract profiles. Longitudinal percentage changes in fractional anisotropy (FA), and mean, axial and radial diffusivity (MD/AD/RD) of each tract from pre-RT were quantified and correlated to 95%, 90% and 80% percentiles of doses and mean doses received by the tract. Heatmaps were used to identify clusters of significant correlation and reveal temporal, spatial and dosimetric signatures of WM damage. A multivariate linear regression was further carried out to determine influence of clinical factors. Results: Of 22 tracts, AD/MD changes in 12 tracts had significant correlation with doses, especially at 6 and 18 months post-RT, indicating progressive radiation damage after RT. Most interestingly, the DTI-index changes in the elongated tracts were associated with received maximum doses, suggesting a serial-structure behavior; while short association fibers were affected by mean doses, indicating a parallel-structure response. Conclusion: Using an automated DTI-tractography analysis of whole brain WM fibers, we reveal complex radiation damage patterns of WM fibers. Damage in WM fibers that play an important role in the neural network could be associated with late neurocognitive function declines after brain irradiation. NIH NS064973.

60 APPLIED LIFE SCIENCES↗

Synchronization effects on rest frame energy and momentum densities in the proton

We obtain two-dimensional relativistic densities and currents of energy and momentum in a proton at rest. These densities are obtained at surfaces of fixed light front time, which physically corresponds to using an alternative synchronization convention. Mathematically, this is done using tilted light front coordinates, which consist of light front time and ordinary spatial coordinates. In this coordinate system, all 16 components of the energy-momentum tensor (EMT) obtain clear physical interpretations, and the nine Galilean components reproduce results from standard light front coordinates. We find angular modulations in several densities that are absent in the corresponding instant form results, which are explained as optical effects arising from using fixed light front time when motion is present within the target. Additionally, transversely polarized spin-half targets exhibit an energy dipole moment—which evaluates to −1/4 for all targets if the Belinfante EMT is used, but which is target dependent and vanishes for pointlike fermions if the asymmetric EMT is instead used.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Requirements on the gain calibration for LiteBIRD polarisation data with blind component separation

The detection of primordial B modes of the cosmic microwave background (CMB) could provide information about the early stages of the Universe's evolution. The faintness of this signal requires exquisite calibration accuracy and control of instrumental systematic effects which otherwise could bias the measurements. In this work, we study the impact of an imperfect relative polarisation gain calibration on the recovered value of the tensor-to-scalar ratio r for the LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We derive requirements on the relative calibration accuracy of the overall polarisation gain (Δg ν ) for each LiteBIRD frequency channel. Our results show that minimum variance techniques, as NILC, are less sensitive to systematic gain calibration uncertainties compared to a parametric approach, if the latter is not equipped with a proper modelling of these instrumental effects. In this study, the most stringent requirements are found in the channels where the CMB signal is relatively brighter, with the tightest constraints at 166 GHz (Δg ν ≈ 0.16%). This differs from the outcome of an analogous analysis performed with a parametric method, where the tightest requirements are obtained for the foreground-dominated channels. Gain calibration uncertainties, corresponding to the derived requirements, are then simultaneously propagated into all frequency channels. By doing so, we find that the overall impact on estimated r is lower than the total gain systematic budget for LiteBIRD approximately by a factor 5, due to the correlations of the impacts of gain calibration uncertainties in different frequency channels. In order to decouple the systematic effect from the specific choice of the model, we derive the requirements assuming constant spectral parameters for the foreground emission. To assess the robustness of the obtained results against more realistic scenarios, we repeat the analysis assuming sky models of intermediate and high complexity. In these further cases, we adopt an optimised NILC pipeline, called the Multi-Clustering NILC (MC-NILC). We find that the impact of gain calibration uncertainties on r is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor 1.8.

79 ASTRONOMY AND ASTROPHYSICS↗

Several versions of the compensated compactness principle

The convergence of the product of a solenoidal vector w{sub {epsilon}} and a gradient {nabla}u{sub {epsilon}} in L{sup 1}({Omega}) (where {Omega} is a region in R{sup d}) is investigated in the case when the factors converge weakly in the spaces L{sup {gamma}({Omega})}{sup d} and L{sup {alpha}({Omega})}{sup d}, respectively, with 1/{gamma}+1/{alpha}>1, which means that the main assumption of the classical div-curl lemma fails. Nevertheless, the same convergence (in the sense of distributions in {Omega}) as in the framework of the div-curl lemma, survives under certain additional assumptions. The new versions of the compensated compactness principle proved in the paper can be used in homogenization and in the theory of G-convergence of monotone operators with non-standard coercivity and growth properties, for instance, some degenerate operators. Bibliography: 20 titles.

97 MATHEMATICS AND COMPUTING↗

Polarization angle requirements for CMB B-mode experiments. Application to the LiteBIRD satellite

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio r parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., ≈ 150 GHz) and of few tens of arcmin at the lowest (i.e., ≈ 40 GHz) and highest (i.e., ≈ 400 GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are ≈ 5 times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on r due to systematics below the limit established by the LiteBIRD collaboration.

79 ASTRONOMY AND ASTROPHYSICS↗

Improvements to Bison Validation Base for LWR Fuel

The Bison fuel performance code is being developed by Idaho National Laboratory (INL) to analyze a wide range of fuel forms and operating conditions with the flexibility to conduct simulations in 1D through 3D geometries. Bison plays an important role in the development of advanced fuel with improved accident tolerance for existing light water reactors (LWRs), the improved understanding of mechanisms in fuel designs in a wider range of operating conditions, and the facilitation of future development of fuel for advanced reactor designs. Additionally, industry interest in advanced fuels has also driven Bison development for advanced fuels in the last few years. All of these factors combine to create a potential increase of new Bison users. These new users will potentially be interested in the Bison LWR validation base for examples of Bison's capabilities and how to use Bison. Current users also rely on the Bison validation base as a demonstration of the best practices settings and material behavioral models. Bison, both through direct updates to the Bison source code and to the nearly daily changes and improvements in the underlying MOOSE framework source code, undergoes nearly continuous changes. It is necessary that these source code changes are reflected with updated input files in the Bison validation base. Periodic updates of the validation case input files should also be accompanied by a review and revision of the validation case documentation on which our users rely. Therefore our efforts under this milestone were guided by two objectives: (1) Update the Bison LWR validation input files to use current Bison simulation settings, and (2) Enhance the ease-of-use for Bison through improved documentation These objectives build on the goals identified FY18 CASL milestone reports. During FY18 we addressed issues with the robustness of Bison and identified a need to migrate the assessment suite to the more reliable tensor mechanics module system and set a goal of standardizing the Bison input files. The FY18 CASL milestone report on updates to the Bison documentation system listed continued development of user-focused guides and examples along with revisions to the validation case documentation as future development goals. During FY19 we have separated these two objectives into four different focus areas for our work on this milestone. Our first two focus areas target the existing Bison validation suite: update the existing LWR validation cases and the accompanying documentation. The second two focus areas center on adding new capabilities: a new Bison LWR validation case and expansion of user support documentation. In the following report each of these four focus areas are discussed as a separate section. We conclude the report with a summary of potential future work to build on the efforts of this year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dimerization and spin decoupling in a two-leg Heisenberg ladder with frustrated trimer rungs

We study the antiferromagnetic spin-half Heisenberg ladder in the presence of an additional frustrating rung spin which is motivated and relevant also for the description of real two-dimensional materials such as the two-dimensional trimer magnet Ba 4 Ir 3 O 10 . We study the zero-temperature phase diagram, where we combine numerical and analytical methods into an overall consistent description. All numerical simulations are also accompanied by studies of the dynamical spin structure factor obtained via the density matrix renormalization group. Overall, we find in the regime of strong rung coupling a gapped dimerized phase related to competing symmetry sectors in Hilbert space that ultimately results in frustration-driven spin-Peierls transition. In the weak rung-coupling regime, the system is uniform, yet shows a gapped spinon continuum together with a sharp coherent low-energy branch which renders the system critical overall. In either case, the additional rung spin quickly get sidelined and nearly decouple in the regime when their bare coupling to the ladder is weakened

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

BICEP Array: 150 GHz Detector Module Development

The Background Imaging of Cosmic Extragalactic Polarization (BICEP)/Keck (BK) collaboration is currently leading the quest for the highest-sensitivity measurements of the polarized cosmic microwave background (CMB) anisotropies on a degree scale with a series of cryogenic telescopes, of which BICEP Array (BA) is the latest Stage-3 upgrade with a total of ~ 32,000 detectors. The instrument comprises 4 receivers spanning 30–270 GHz, with the low-frequency 30/40 GHz deployed to the South Pole Station in late 2019. The full complement of receivers is forecast to set the most stringent constraints on the tensor-to-scalar ratio r. Building on these advances, the overarching small-aperture telescope concept is already being used as the reference for further Stage-4 experiment design. This paper describes the development of the BICEP Array 150 GHz detector module and its fabrication requirements, with highlights on the high-density time division multiplexing (TDM) design of the cryogenic circuit boards. The low-impedance wiring required between the detectors and the frst stage of superconducting quantum interference device amplifers is crucial to maintaining a stable bias current on the detectors. Here, a novel multilayer FR4 Printed Circuit Board with superconducting traces, capable of reading out up to 648 detectors, is detailed along with its validation tests. An ultra-high-density TDM detector module concept we developed for a CMB-S4-like experiment that allows up to 1920 detectors to be read out is also presented. TDM has been chosen as the detector readout technology for the Cosmic Microwave Background Stage-4 (CMB-S4) experiment based on its proven low-noise performance, predictable costs, and overall maturity of the architecture. The heritage for TDM is rooted in mm- and sub-mm-wave experiments dating back 20 years and has since evolved to support a multiplexing factor of 64x in Stage-3 experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Partons as unique ground states of quantum Hall parent Hamiltonians: The case of Fibonacci anyons

We present microscopic, multiple Landau level, (frustration-free and positive semi-definite) parent Hamiltonians whose ground states, realizing different quantum Hall fluids, are parton-like and whose excitations display either Abelian or non-Abelian braiding statistics. We prove ground state energy monotonicity theorems for systems with different particle numbers in multiple Landau levels, demonstrate S-duality in the case of toroidal geometry, and establish complete sets of zero modes of special Hamiltonians stabilizing parton-like states, specifically at filling factor \nu=2/3 ν = 2 / 3 . The emergent Entangled Pauli Principle (EPP), introduced in [Phys. Rev. B 98, 161118(R) (2018)] and which defines the “DNA” of the quantum Hall fluid, is behind the exact determination of the topological characteristics of the fluid, including charge and braiding statistics of excitations, and effective edge theory descriptions. When the closed-shell condition is satisfied, the densest (i.e., the highest density and lowest total angular momentum) zero-energy mode is a unique parton state. We conjecture that parton-like states generally span the subspace of many-body wave functions with the two-body M M -clustering property within any given number of Landau levels, that is, wave functions with M M th-order coincidence plane zeroes and both holomorphic and anti-holomorphic dependence on variables. General arguments are supplemented by rigorous considerations for the M=3 M = 3 case of fermions in four Landau levels. For this case, we establish that the zero mode counting can be done by enumerating certain patterns consistent with an underlying EPP. We apply the coherent state approach of [Phys. Rev. X 1, 021015 (2011)] to show that the elementary (localized) bulk excitations are Fibonacci anyons. This demonstrates that the DNA associated with fractional quantum Hall states encodes all universal properties. Specifically, for parton-like states, we establish a link with tensor network structures of finite bond dimension that emerge via root level entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic dipole interactions in crystals

The influence of magnetic dipole interactions (MDIs) on the magnetic properties of local-moment Heisenberg spin systems is investigated. A general formulation is presented for calculating the eigenvalues λ and eigenvectors μ ˆ of the MDI tensor of the magnetic dipoles in a line (one dimension, 1D), within a circle (2D) or a sphere (3D) of radius r surrounding a given moment μ → i for given magnetic propagation vectors k for collinear and coplanar noncollinear magnetic structures on both Bravais and non-Bravais spin lattices. Results are calculated for collinear ordering on 1D chains, 2D square and simple-hexagonal (triangular) Bravais lattices, 2D honeycomb and kagomé non-Bravais lattices, and 3D cubic Bravais lattices. The λ and μ ˆ values are compared with previously reported results. Calculations for collinear ordering on 3D simple tetragonal, body-centered tetragonal, and stacked triangular and honeycomb lattices are presented for c/a ratios from 0.5 to 3 in both graphical and tabular form to facilitate comparison of experimentally determined easy axes of ordering on these Bravais lattices with the predictions for MDIs. Comparisons with the easy axes measured for several illustrative collinear antiferromagnets (AFMs) are given. The calculations are extended to the cycloidal noncollinear 120 ° AFM ordering on the triangular lattice where λ is found to be the same as for collinear AFM ordering with the same k. The angular orientation of the ordered moments in the noncollinear coplanar AFM structure of GdB 4 with a distorted stacked 3D Shastry-Sutherland spin-lattice geometry is calculated and found to be in disagreement with experimental observations, indicating the presence of another source of anisotropy. Similar calculations for the undistorted 2D and stacked 3D Shastry-Sutherland lattices are reported. The thermodynamics of dipolar magnets are calculated using the Weiss molecular field theory for quantum spins, including the magnetic transition temperature T m and the ordered moment, magnetic heat capacity, and anisotropic magnetic susceptibility χ versus temperature T . The anisotropic Weiss temperature θ p in the Curie-Weiss law for T>T m is calculated. A quantitative study of the competition between FM and AFM ordering on cubic Bravais lattices versus the demagnetization factor in the absence of FM domain effects is presented. The contributions of Heisenberg exchange interactions and of the MDIs to T m and to θ p are found to be additive, which simplifies analysis of experimental data. Some properties in the magnetically-ordered state versus T are presented, including the ordered moment and magnetic heat capacity and, for AFMs, the dipolar anisotropy of the free energy and the perpendicular critical field. The anisotropic χ for dipolar AFMs is calculated both above and below the Néel temperature T N and the results are illustrated for a simple tetragonal lattice with c/a>1, c/a=1 (cubic), and c/a<1 , where a change in sign of the χ anisotropy is found at c/a=1 . Finally, following the early work of Keffer [Phys. Rev. 87, 608 (1952)], the dipolar anisotropy of χ above T N =69 K of the prototype collinear Heisenberg-exchange-coupled tetragonal compound MnF 2 is calculated and found to be in excellent agreement with experimental single-crystal literature data above 130 K, where the smoothly increasing deviation of the experimental data from the theory on cooling from 130 K to T N is deduced to arise from dynamic short-range collinear c -axis AFM ordering in this temperature range driven by the exchange interactions.

36 MATERIALS SCIENCE↗

Correlating Protein Dynamics and Catalytic Activity of a Model Hydrogenase Using Paramagnetic and Biological Nuclear Magnetic Resonance Spectroscopy

Rational catalyst design remains a significant challenge, with electronic structure, steric, and electrostatic effects known to contribute to activity. Recently, dynamics has been recognized as another factor that impacts catalysis, though identifying and predicting these effects has remained out of reach. Nickel-substituted rubredoxin (NiRd), a protein-based mimic of a hydrogenase enzyme, serves as a model catalytic system in which dynamics can be systematically investigated with respect to activity. While over 30 secondary-sphere mutants of NiRd have been shown to be catalytically active, no significant correlation was observed between the rates and catalytic overpotential or electronic structure, prompting questions about the protein-derived factors that modulate activity. Here, in this work, NMR spectroscopy was used to investigate the roles of substrate accessibility, protein dynamics, and protein stability in controlling catalysis. Significant paramagnetic effects from the nickel center (S = 1) isolate the methylene proton resonances of the metal-coordinating cysteine residues. The sensitivity of resonance positions and linewidths to local environment offers an opportunity to study dynamical molecular changes around the metal center with high resolution. Machine learning algorithms were employed to identify correlations between the catalytic activity and the paramagnetic NMR spectra. These analyses revealed spectroscopic features of specific cysteine protons that report on catalytic overpotential and increased turnover rates, which are further supported by the results obtained using high-field NMR techniques. Collectively, these studies indicate the potential for multifrequency NMR techniques to resolve key contributors to catalytic activity and highlight the importance of local and outer-sphere dynamics.

Protein Engineering↗

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" 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 "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗