Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel computational chemistry”

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.

376 records · Page 21

Speeding up equation of motion coupled cluster theory with the chain of spheres approximation

In the present paper, the chain of spheres exchange (COSX) approximation is applied to the highest scaling terms in the equation of motion (EOM) coupled cluster equations with single and double excitations, in particular, the terms involving integrals with four virtual labels. It is found that even the acceleration of this single term yields significant computational gains without compromising the desired accuracy of the method. For an excitation energy calculation on a cluster of five water molecules using 585 basis functions, the four virtual term is 9.4 times faster using COSX with a loose grid than using the canonical implementation, which yields a 2.6 fold acceleration for the whole of the EOM calculation. For electron attachment calculations, the four virtual term is 15 times and the total EOM calculation is 10 times faster than the canonical calculation for the same system. The accuracy of the new method was tested using Thiel’s test set for excited states using the same settings and the maximum absolute deviation over the whole test set was found to be 12.945 cm{sup −1} (59 μHartree) for excitation energies and 6.799 cm{sup −1} (31 μHartree) for electron attachments. Using MP2 amplitudes for the ground state in combination with the parallel evaluation of the full EOM equations in the manner discussed in this paper enabled us to perform calculations for large systems. Electron affinity values for the two lowest states of a Zn protoporphyrine model compound (224 correlated electrons and 1120 basis functions) were obtained in 3 days 19 h using 4 cores of a Xeon E5-2670 processor allocating 10 GB memory per core. Calculating the lowest two excitation energies for trans-retinal (114 correlated electrons and 539 basis functions) took 1 day 21 h using eight cores of the same processor and identical memory allocation per core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reversible Olefin Addition to Extended Lattices of a Nickel–Selenium Framework

We performed density functional theory computations to study the structural and electronic properties as the basis of ethylene addition activity for [Ni(XC)4]n (X = Se, S)-extended lattices. We demonstrated that the mechanism of ethylene cycloaddition to a periodic [Ni(SeC)4]n two-dimensional (2D) network is analogous to that previously described for [Ni(SC)4]n 2D sheets and similar to the metal bis(dithiolene) molecular complexes [M(S2C2R2)2] (M = Ni, Pd, Pt, Co, Cu). These nanosheet materials avoid decomposition upon olefin addition, which is one of the main limitations of the molecular metal bis(dithiolene) complexes, as we find the decomposition processes to be thermodynamically unfavorable. Our calculations also suggest that the preferred conformation of the [Ni(SeC)4]n bilayer lattice is parallel displaced, with the Se atoms positioned above the Ni atoms, which is different from the eclipsed conformation found for [Ni(SC)4]n. We also managed to optimize an adduct of [Ni(SC)4]n in the bilayer form, which exceed the ethylene coverage of molecular complexes. We calculate that the preferred three-dimensional geometry of the stacked sheets is eclipsed because of strong van der Waals interactions. Such an arrangement of the sheets indicates that these materials should be highly porous, pointing to the high capacity for olefin bindings. Indeed, a few moderately stable ethylene adducts have been located. Owing to their unique structures and chemical reactivity, these newly predicted materials can be potentially developed as electrocatalysts for olefin purification. © 2018 American Chemical Society.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers

This report is the final scientific one for the award DE- FE0000988 entitled “Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers.” The work has been divided into six tasks. In task, “Development of a Three-Phase Non-Isothermal CO 2 Flow Module,” we developed a fluid property module for brine-CO 2 mixtures designed to handle all possible phase combinations of aqueous phase, sub-critical liquid and gaseous CO 2 , supercritical CO 2 , and solid salt. The thermodynamic and thermophysical properties of brine-CO 2 mixtures (density, viscosity, and specific enthalpy of fluid phases; partitioning of mass components among the different phases) use the same correlations as an earlier fluid property module that does not distinguish between gaseous and liquid CO 2 -rich phases. We verified the fluid property module using two leakage scenarios, one that involves CO 2 migration up a blind fault and subsequent accumulation in a secondary “parasitic” reservoir at shallower depth, and another investigating leakage of CO 2 from a deep storage reservoir along a vertical fault zone. In task, “Development of a Rock Mechanical Module,” we developed a massively parallel reservoir simulator for modeling THM processes in porous media brine aquifers. We derived, from the fundamental equations describing deformation of porous elastic media, a momentum conservation equation relating mean stress, pressure, and temperature, and incorporated it alongside the mass and energy conservation equations from the TOUGH2 formulation, the starting point for the simulator. In addition, rock properties, namely permeability and porosity, are functions of effective stress and other variables that are obtained from the literature. We verified the simulator formulation and numerical implementation using analytical solutions and example problems from the literature. For the former, we matched a one-dimensional consolidation problem and a two-dimensional simulation of the Mandel-Cryer effect. For the latter, we obtained a good match of temperature and gas saturation profiles, and surface uplift, after injection of hot fluid into a model of a caldera structure. In task, “Incorporation of Geochemical Reactions of Selected Important Species,” we developed a novel mathematical model of THMC processes in porous and fractured saline aquifers, simulating geo-chemical reactions associated with CO 2 sequestration in saline aquifers. Two computational frameworks, sequentially coupled and fully coupled, were used to simulate the reactions and transport. We verified capabilities of the THMC model to treat complex THMC processes during CO 2 sequestration by analytical solutions and we constructed reactive transport models to analyze the THMC process quantitatively. Three of these are 1D reactive transport under chemical equilibrium, a batch reaction model with equilibrium chemical reactions, and a THMC model with CO 2 dissolution. In task “Study of Instability in CO 2 Dissolution-Diffusion-Convection Processes,” We reviewed literature related to the study of density driven convective flows and on the instability of CO 2 dissolution-diffusion-convection processes. We ran simulations that model the density-driven flow instability that would occur during CO 2 sequestration. CO 2 diffused through the top of the system and dissolved in the aqueous phase there, increasing its density. Density fingers formed along the top boundary, and coalesced into a few prominent ones, causing convective flow that forced the fluid to the system bottom. These simulations were in two and three dimensions. We ran additional simulations of convective mixing with density contrast caused by variable dissolved CO 2 concentration in saline water, modeled after laboratory experiments in which supercritical CO2 was circulated in the headspace above a brine saturated packed sand in a pressure vessel. As CO 2 dissolved into the upper part of the saturated sand, liquid phase density increases causing instability and setting off convective mixing. We obtained good agreement with the laboratory experiments, which were characterized by finger development and associated mixing of dissolved CO 2 into the system. We then varied a wide range of parameters and conceptual models in order to analyze the possibility of convective mixing under different conditions, such as various boundary conditions, and chemical reaction conditions. The CO 2 fingers from different simulations showed great differences as time progressed, caused by permeability heterogeneity. The early time diffusive phenomenon was captured by fine grid resolution, and the permeability heterogeneity affected the pattern of the CO 2 fingers. In addition, the fingers from three-dimensional simulations tended to be larger and flatter than the two-dimensional ones. In task “Implementation of Efficient Parallel Computing Technologies,” we made enhancements and modifications to our code in order to substantially increase the grid size that could be run. We installed and ran it on various platforms, including a multi-core PC and a cluster, and verified the numerical implementation and parallel code using an example problem from the literature. This problem, with a grid size of sixty million, utilized the cluster’s entire memory and all of its processors. In task “Implementation of General Fracture Conceptual Models,” we used the MINC approach, a generalization of the double-porosity concept, to model flow through porous and fractured media. In this approach, flow within the matrix is described by subdividing the matrix into nested volumes, with flow occurring between adjacent nested matrix volumes as well as between the fractures and the outer matrix volume. We generalized Hooke’s law to a thermo-multi- poroelastic medium, and derived from the fundamental equations describing deformation of porous and fractured elastic media a momentum conservation equation for thermo-multi- poroelastic media. This equation is a generalization to multi-poroelastic media of the one derived in Task 3.0 for single porosity media. We describe two simulations to provide model verification and application examples. The first, one-dimensional consolidation of a double-porosity medium, is compared to an analytical solution. The second is a match of published results from the literature, a simulation of CO 2 injection into hypothetical aquifer-caprock systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determining polarizable force fields with electrostatic potentials from quantum mechanical linear response theory

We developed a new method to calculate the atomic polarizabilities by fitting to the electrostatic potentials (ESPs) obtained from quantum mechanical (QM) calculations within the linear response theory. This parallels the conventional approach of fitting atomic charges based on electrostatic potentials from the electron density. Our ESP fitting is combined with the induced dipole model under the perturbation of uniform external electric fields of all orientations. QM calculations for the linear response to the external electric fields are used as input, fully consistent with the induced dipole model, which itself is a linear response model. The orientation of the uniform external electric fields is integrated in all directions. The integration of orientation and QM linear response calculations together makes the fitting results independent of the orientations and magnitudes of the uniform external electric fields applied. Another advantage of our method is that QM calculation is only needed once, in contrast to the conventional approach, where many QM calculations are needed for many different applied electric fields. The molecular polarizabilities obtained from our method show comparable accuracy with those from fitting directly to the experimental or theoretical molecular polarizabilities. Since ESP is directly fitted, atomic polarizabilities obtained from our method are expected to reproduce the electrostatic interactions better. Our method was used to calculate both transferable atomic polarizabilities for polarizable molecular mechanics’ force fields and nontransferable molecule-specific atomic polarizabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Theoretical Study on Reaction Mechanisms and Kinetics of Thiophene Hydrodesulfurization over MoS 2 Catalysts

In the present work, thiophene hydrodesulphurization (HDS) over the Mo-edge, the S-edge, and the Mo-S connection edge of MoS 2 catalyst with 50% sulfur coverage was studied using first-principles based microkinetic modeling. Two parallel HDS routes, i.e., direct desulfurization (DDS) and hydrogenation (HYD) were taken into account. It has been found that the major reaction route of thiophene HDS on the Mo- and the Mo-S edges is temperature dependent. In the low temperature range of 500~600 K, the HYD route is dominant, leading to the C 4 H 8 formation. As the temperature increases, the DDS route becomes competitive with the HYD route. At the temperature above 650 K, the DDS route will be the dominant HDS reaction route on the Mo- and the Mo-S edges. While the DDS route that leads to the formation of C4H6 is the major thiophene HDS reaction route on the S-edge in the entire temperature range of 500~750 K. The microkinetic modeling results show the overall HDS activity on the S-edge is lower than it on the Mo- and the Mo-S edges. The Mo-S edge also provides a preferential reaction pathway, which facilitates 2-hydrothiophene migration from the Mo-edge to the S-edge, followed by remaining elementary steps with lower activation barriers in the DDS route. The work was financially supported by the National Natural Science Foundation of China (No. 21476012, 21571012, 21407007 and 91534201). The computing time was granted by a scientific theme user proposal in the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL), which is a U.S. Department of Energy national scientific user facility located at Pacific Northwest National Laboratory (PNNL) in Richland, Washington. Part of computing time was provided the National Energy Research Scientific Computing Center (NERSC). D. Mei is supported by the US Department of Energy, Office of Sciences, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences.

Thiophene, hydrodesulphurization, MoS2, DFT, micro↗

BF{sub 3} PIII modeling: Implantation, amorphisation and diffusion

In the race for highly doped ultra-shallow junctions (USJs) in complementary metal oxide semi-conductor (CMOS) technologies, plasma immersion ion implantation (PIII) is a promising alternative to traditional beamline implantation. Currently, no commercial technology computer aided design (TCAD) process simulator allows modeling the complete USJ fabrication process by PIII, including as-implanted dopant profiles, damage formation, dopant diffusion and activation. In this work, a full simulation of a p-type BF{sub 3} PIII USJ has been carried out. In order to investigate the various physical phenomena mentioned above, process conditions included a high energy/high dose case (10 kV, 5 Multiplication-Sign 10{sup 15} cm{sup -2}), specifically designed to increase damage formation, as well as more technology relevant implant conditions (0.5 kV) for comparison. All implanted samples were annealed at different temperatures and times. As implanted profiles for both boron and fluorine in BF{sub 3} implants were modeled and compared to Secondary Ion Mass Spectrometry (SIMS) measurements. Amorphous/crystalline (a/c) interface depths were measured by transmission electron microscopy (TEM) and successfully simulated. Diffused profiles simulations agreed with SIMS data at low thermal budgets. A boron peak behind the a/c interface was observed in all annealed SIMS profiles for the 10 kV case, indicating boron trapping from EOR defects in this region even after high thermal budgets. TEM measurements on the annealed samples showed an end of range (EOR) defects survival behind the a/c interface, including large dislocation loops (DLs) lying on (001) plane parallel to the surface. In the last part of this work, activation simulations were compared to Hall measurements and confirmed the need to develop a (001) large BICs model.

36 MATERIALS SCIENCE↗

Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method

A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation. Two methods for computing the conformational free energy (FE) along the string are compared, a restrained method, and a tessellation method introduced by E. Vanden-Eijnden and M. Venturoli [J. Chem. Phys. 130, 194103 (2009)]. It is found that obtaining meaningful free energy profiles using the present atom-based coordinates requires restricting sampling to a vicinity of the converged path, where the hyperplanar approximation to the isocommittor surface is sufficiently accurate. This sampling restriction can be easily achieved using restraints or constraints. The endpoint FE differences computed from the FE profiles are validated by comparison with previous calculations using a path-independent confinement method. The FE profiles are decomposed into the enthalpic and entropic contributions, and it is shown that the entropy difference contribution can be as large as 10 kcal/mol for intermediate regions along the path, compared to 15–20 kcal/mol for the enthalpy contribution. This result demonstrates that enthalpic barriers for transitions are offset by entropic contributions arising from the existence of different paths across a barrier. The possibility of using systematically coarse-grained representations of amino acids, in the spirit of multiple interaction site residue models, is proposed as a means to avoid ad hoc sampling restrictions to narrow transition tubes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: Difference between atomistic and coarse-grained simulations

Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT{sup OPLS}). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT{sup OPLS} PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT{sup CHARMM}) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16–22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT{sup OPLS} PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA){sub 3} peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}: Structural, spectroscopic and computational studies on a sorosilicate

Synthesis experiments in the system Li{sub 2}O–CaO–SiO{sub 2} resulted in the formation of single-crystals of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}. Structural investigations were based on single-crystal diffraction. At ambient conditions the compound has the following basic crystallographic data: hexagonal symmetry, space group P6{sub 1}22, a=5.0961(2) Å, c=41.264(2) Å, V=928.07(6) Å{sup 3}, Z=6. Structure solution was performed using direct methods. The final least-squares refinement calculations converged at a residual of R(|F|)=0.0260. From a structural point the lithium calcium silicate belongs to the group of pyrosilicates containing [Si{sub 2}O{sub 7}]-groups. Additional lithium and calcium cations are incorporated between the silicate dimers and are coordinated by four and six nearest oxygen neighbours, respectively. Each [LiO{sub 4}]-tetrahedron shares two common corners with directly neighboring tetrahedra forming zweier single-chains which are running parallel to 〈1 0 0〉 in z-levels defined by the presence of the 6{sub 1}{sup [0} {sup 0} {sup 1]}-screw axes. From the corner-sharing [LiO{sub 4}]- and [SiO{sub 4}]-moieties a three dimensional framework can be constructed. An interesting feature of this framework is the presence of an O{sup [3]}-type bridging oxygen linking three tetrahedra (one [LiO{sub 4}]- and two [SiO{sub 4}]-units). Structural similarities with other silicates are discussed in detail. The high-temperature behavior of the Si–O, Ca–O and Li–O bond distances in Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} was investigated by in{sub -}situ single-crystal X-ray diffraction in the range between 65 and 700 °C. From the evolution of the lattice parameters, the thermal expansion tensor α{sub ij} has been determined. The structural characterization has been supplemented by micro-Raman spectroscopy. Interpretation of the spectroscopic data including the allocation of the bands to certain vibrational species has been aided by DFT-calculations. - Graphical abstract: Framework of [SiO{sub 4}]- and [LiO{sub 4}]-tetrahedra in the crystal structure of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}. - Highlights: • We report the crystal structure of the sorosilicate Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} at ambient conditions. • The thermal expansion tensor as well as the evolution of the structure between 25 and 700 °C was studied. • A topological analysis of the tetrahedral framework based on [SiO{sub 4}]- and [LiO{sub 4}]-units is presented. • The crystal structure of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} is discussed with respect to related compounds. • Interpretation of the Raman spectra of the crystals has been aided by DFT-calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exceptionally Complex Electronic Structures of Lanthanide Oxides and Small Molecules

Lanthanide (Ln) oxide clusters and molecular systems provide a bottom-up look at the electronic structures of the bulk materials because of close parallels in the patterns of Ln 4fN subshell occupancy between the molecular and bulk Ln 2 O 3 size limits. At the same time, these clusters and molecules offer a challenge to the theory community to find appropriate and robust treatments for the 4f N patterns across the Ln series. Anion photoelectron (PE) spectroscopy provides a powerful experimental tool for studying these systems, mapping the energies of the ground and low-lying excited states of the neutral relative to the initial anion state, providing spectroscopic patterns that reflect the Ln 4f N occupancy. In this Account, we review our anion PE spectroscopic and computational studies on a range of small lanthanide molecules and cluster species. The PE spectra of LnO – (Ln = Ce, Pr, Sm, Eu) diatomic molecules show spectroscopic signatures associated with detachment of an electron from what can be described as a diffuse Ln 6s-like orbital. While the spectra of all four diatomics share this common transition, the fine structure in the transition becomes more complex with increasing 4f occupancy. This effect reflects increased coupling between the electrons occupying the corelike 4f and diffuse 6s orbitals with increasing N. Understanding the PE spectra of these diatomics sets the stage for interpreting the spectra of polyatomic molecular and cluster species. In general, the results confirm that the partial 4f N subshell occupancy is largely preserved between molecular and bulk oxides and borides. However, they also suggest that surfaces and edges of bulk materials may support a low-energy, diffuse Ln 6s band, in contrast to bulk interiors, in which the 6s band is destabilized relative to the 5d band. We also identify cases in which the molecular Ln centers have 4f N+1 occupancy rather than bulklike 4f N , which results in weaker Ln–O bonding. Specifically, Sm centers in mixed Ce–Sm oxides or in Sm x O y – (y ≤ x) clusters have this higher 4f N+1 occupancy. The PE spectra of these particular species exhibit a striking increase in the relative intensities of excited-state transitions with decreasing photon energy (resulting in lower photoelectron kinetic energy). Furthermore, this is opposite of what is expected on the basis of the threshold laws that govern photodetachment. We relate this phenomenon to strong electron–neutral interactions unique to these complex electronic structures. The time scale of the interaction, which shakes up the electronic configuration of the neutral, increases with decreasing electron momentum. From a computational standpoint, we point out that special care must be taken when considering Ln cluster and molecular systems toward the center of the Ln series (e.g., Sm, Eu), where treatment of electrons explicitly or using an effective core potential can yield conflicting results on competing subshell occupancies. However, despite the complex electronic structures associated with partially filled 4f N subshells, we demonstrate that inexpensive and tractable calculations yield useful qualitative insight into the general electronic structural features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗