Search NASA⌕ Search

SEARCH · Search NASA

Results for “materials and molecular modeling”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bayesian Conavigation: Dynamic Designing of the Material Digital Twins via Active Learning

Scientific advancement is universally based on the dynamic interplay between theoretical insights, modeling, and experimental discoveries. However, this feedback loop is often slow, including delayed community interactions and the gradual integration of experimental data into theoretical frameworks. This challenge is particularly exacerbated in domains dealing with high-dimensional object spaces, such as molecules and complex microstructures. Hence, the integration of theory within automated and autonomous experimental setups, or theory in the loop-automated experiment, is emerging as a crucial objective for accelerating scientific research. The critical aspect is to use not only theory but also on-the-fly theory updates during the experiment. Furthermore, we introduce a method for integrating theory into the loop through Bayesian conavigation of theoretical model space and experimentation. Our approach leverages the concurrent development of surrogate models for both simulation and experimental domains at the rates determined by latencies and costs of experiments and computation, alongside the adjustment of control parameters within theoretical models to minimize epistemic uncertainty over the experimental object spaces. This methodology facilitates the creation of digital twins of material structures, encompassing both the surrogate model of behavior that includes the correlative part and the theoretical model itself. While being demonstrated here within the context of functional responses in ferroelectric materials, our approach holds promise for broader applications, such as the exploration of optical properties in nanoclusters, microstructure-dependent properties in complex materials, and properties of molecular systems.

Microscopy↗
DOE OSTI2020

Scale and rate in CdS pressure-induced phase transition

Here, we describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.

42 ENGINEERING↗
NASA NTRS2007

Measurement of Damage Profiles from Solar Wind Implantation

NASA's Genesis Mission launched from Cape Canaveral in August of 2001 with the goal of collecting solar wind in ultra-pure materials. The samples were returned to Earth more than three years later for subsequent analysis. Although the solar wind is comprised primarily of protons, it also contains ionized species representing the entire periodic table. The Genesis mission took advantage of the natural momentum of these ionized species to implant themselves in specialized collectors including single crystal Si and SiC. The collectors trapped the solar wind species of interest and sustained significant damage to the surface crystal structure as a result of the ion bombardment. In this work, spectroscopic ellipsometry has been used to evaluate the extent of this damage in Si and SiC samples. These results and models are compared for artificially implanted samples and pristine non-flight material. In addition, the flown samples had accumulated a thin film of molecular contamination as a result of outgassing in flight, and we demonstrate that this layer can be differentiated from the material damage. In addition to collecting bulk solar wind samples (continuous exposure), the Genesis mission actually returned silicon exposed to four different solar wind regimes: bulk, high speed, low speed, and coronal mass ejections. Each of these solar wind regimes varies in energy, but may vary in composition as well. While determining the composition is a primary goal of the mission, we are also interested in the variation in depth and extent of the damage layer as a function of solar wind regime. Here, we examine flight Si from the bulk solar wind regime and compare the results to both pristine and artificially implanted Si. Finally, there were four samples which were mounted in an electrostatic "concentrator" designed to reject a large fraction (>85%) of incoming protons while enhancing the concentration of ions mass 4-28 amu by a factor of at least 20. Two of these samples were single crystal 6H silicon carbide. (The others were polycrystalline CVD diamond and amorphous carbon that were not examined in the work.) The ion damaged SiC samples from the concentrator were studied in comparison to the flight Si from the bulk array to understand differences in the extent of the damage.

McNamara, K. M.↗
DOE OSTI2021

Porous Covalent Organic Polymers for efficient Fluorocarbon-based Adsorption Cooling

Adsorption-based cooling is an energy-efficient renewable-energy technology that can be driven using low-grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear-shaped isotherms are revealed for the materials, namely COP-2 and COP-3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect-containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect-containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect-containing porous COPs are promising for adsorption-based cooling applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

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↗
DOE OSTI2026

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhancement of Photovoltaic Current through Dark States in Donor–Acceptor Pairs of Tungsten–Based Transition Metal Di–Chalcogenides

As several photovoltaic materials experimentally approach the Shockley– Queisser limit, there has been a growing interest in unconventional materials and approaches with the potential to cross this efficiency barrier. One such candidate is dark state protection induced by the dipole–dipole interaction between molecular excited states. This phenomenon has been shown to significantly reduce carrier recombination rate and enhance photon-to-current conversion, in elementary models consisting of few interacting chromophore centers. Atomically thin 2D transition metal di-chalcogenides (TMDCs) have shown great potential for use as ultrathin photovoltaic materials in solar cells due to their favorable photon absorption and electronic transport properties. TMDC alloys exhibit tunable direct bandgaps and significant dipole moments. In this work, the dark state protection mechanism has been introduced to a TMDC based photovoltaic system with pure tungsten diselenide (WSe 2 ) as the acceptor material and the TMDC alloy tungsten sulfo-selenide (WSeS) as the donor material. Furthermore, our numerical model demonstrates the first application of the dark state protection mechanism to a photovoltaic material with a photon current enhancement of up to 35% and an ideal photon-to-current efficiency exceeding the Shockley–Queisser limit.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Landau-Lifshitz-Bloch equation for ferrimagnets with higher-order interaction

In this work, we present a micromagnetic formulation for modeling the magnetization dynamics and thermal equilibrium in ferrimagnetic materials at low and elevated temperatures. The formulation is based on a mean-field approximation (MFA). In this formulation, the ferrimagnet is described micromagnetically by two coupled sublattices with corresponding interactions, including inter- and intrasublattice micromagnetic exchange as well as four-spin interactions described as an intersublattice molecular field with a cubic dependence of the magnetization. The MFA is used to derive a Landau-Lifshitz-Bloch type equation for ferrimagnetic material, including cases with a ferromagnetic-antiferromagnetic phase transitions. For validation, the results obtained via the presented model are compared with recent experimental data for phase transitions in FeRh.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Polyalkenamers as Drop-In Additives for Ring-Opening Metathesis Polymerization: A Promising Upcycling Paradigm

Here we report a distinct strategy to upcycle waste polyalkenamers such as polybutadiene into new, performance-advantaged materials by using them as drop-in additives for ring-opening metathesis polymerization (ROMP). The polyalkenamers serve as competent chain-transfer agents in ROMPs of common classes of cyclic olefin monomers, facilitating good molecular weight control, allowing low Ru catalyst loadings, and enabling efficient incorporation of the polyalkenamer into the synthesized polymeric material. We successfully demonstrate ROMP using model polyalkenamers and translate these learnings to leverage commercial polybutadiene and acrylonitrile butadiene styrene (ABS) as chain transfer agents for ROMP copolymerizations. Critically, our strategy is shown to be highly efficient and operationally simple, quantitatively incorporating the polyalkenamer and inheriting aspects of its thermomechanical performance. Our results highlight a promising pathway for the upcycling of polyalkenamers and provide an alternative to existing deconstruction and functional upcycling strategies.

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

Epitaxial registry and transport properties in V4O7 thin films

We report the growth of V4O7 thin films deposited simultaneously on amorphous SiO2 and crystalline c-cut Al2O3 substrates. x-ray diffraction shows that films grown on SiO2 are polycrystalline with no preferred orientation, while films grown on Al2O3 exhibit a well-defined out-of-plane orientation and an in-plane registry consistent with epitaxial growth. Transmission electron microscopy confirms the crystallographic relationship between the V4O7 film and the Al2O3 substrate. Atomic force microscopy indicates substantially lower surface roughness for films grown on Al2O3 (~6 nm) compared to those grown on SiO2 (~22 nm), and Raman spectroscopy confirms stabilization of the same V4O7 phase on both substrates. Electrical transport measurements reveal a metal–insulator transition near 244 K for both substrates, with thermal hysteresis not exceeding ~1 K. Although the transition temperature remains essentially unchanged, films grown on Al2O3 exhibit higher conductivity over the entire temperature range, exceeding that of films grown on SiO2 by approximately two orders of magnitude at 100 K and by a factor of five at 300 K. These results indicate that the conductivity differences are consistent with variations in microstructural connectivity associated with crystallographic order.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Crystallographic registry in epitaxial V3O5 thin films

Crystallographic orientation plays a central role in determining transport properties in correlated oxides with low-hysteresis metal–insulator transitions. Here, V3O5 thin films were grown on Al2O3 (0001) and their crystallographic registry was evaluated using x-ray diffraction and cross-sectional transmission electron microscopy. θ–2θ measurements show a single out-of-plane orientation defined by the (202) reflection. Azimuthal φ-scans of the V3O5 (310) reflection exhibit twelve discrete maxima over 360°, indicating a finite set of in-plane orientations imposed by the substrate symmetry. High-resolution TEM and fast Fourier transform analysis confirm crystallographic coherence at the film–substrate interface. The electrical conductivity exhibits a metal–insulator transition at TMIT ≈ 423 K with no measurable thermal hysteresis and a total change of ~1.6 orders of magnitude between 300 and 480 K. These results establish crystallographic registry in epitaxial V3O5 thin films and provide a structurally well-defined system for future studies of orientation-dependent transport.\r\nT

36 MATERIALS SCIENCE↗
DOE OSTI2022

Discovery of Graphene-Water Membrane Structure: Toward High-Quality Graphene Process

It is widely accepted that solid-state membranes are indispensable media for the graphene process, particularly transfer procedures. But these membranes inevitably bring contaminations and residues to the transferred graphene and consequently compromise the material quality. This study reports a newly observed free-standing graphene-water membrane structure, which replaces the conventional solid-state supporting media with liquid film to sustain the graphene integrity and continuity. Experimental observation, theoretical model, and molecular dynamics simulations consistently indicate that the high surface tension of pure water and its large contact angle with graphene are essential factors for forming such a membrane structure. More interestingly, water surface tension ensures the flatness of graphene layers and renders high transfer quality on many types of target substrates. This report enriches the understanding of the interactions on reduced dimensional material while rendering an alternative approach for scalable layered material processing with ensured quality for advanced manufacturing.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Observations of ultraviolet-bright stars in globular clusters with the Hopkins Ultraviolet Telescope

Two UV-bright stars in globular clusters, UV5 in NGC 1851 and vZ 1128 in M3, were observed with the Hopkins Ultraviolet Telescope (HUT) during the Astro-1 space shuttle mission in 1990 December. The stars' spectra show weak absorption features and no significant emission features other than well known geocoronal lines. Detailed fitting of Kurucz (1991) stellar atmosphere models using a chi(sup 2) minimization technique indicates T(sub eff) = 16 000 K, log g = 2.5, and abundance (-1.0) (logarithm of abundance of elements heavier than helium relative to solar) for UV5, and T(sub eff) = 35 000 K, log g = 4.0, and abundance (-3.5) for vZ 1128. The Kurucz model which best fits vZ 1128 overpredicts the flux in the region below approximately 1000 A, an effect seen in previous models of O-type stars. Our results are robust (to within approximately 1000 K) with respect to uncertainties in interstellar reddening and atomic and molecular hydrogen column densities. We do not see significant molecular hydrogen absorption, which might have indicated material in a circumstellar shell, in either star's spectrum. We estimate the stellar luminosities to be log L/solar luminosity = 3.33 +/- 0.15 for UV5 and log L/solar luminosity = 3.21 +/- 0.12 for vZ 1128. These atmospheric parameters place both stars on the Schoenberner (post-AGB) tracks, though the stellar masses derived from the best-fitting Kurucz models are somewhat less than those predicted by the Schoenberner models. Examination of individual absorption line strengths reveals no significant abundance anomalies in either star.

Van Dyke Dixon, W.↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Next-Generation Materials Design: Quantum Mechanics and Data-Driven Modeling

The future of materials design is rapidly advancing through the combination of quantum mechanics and data-driven modeling. These approaches integrate quantum principles with advanced data analysis, enabling precise insights into material behavior. This talk will highlight recent progress in using these methods for computational design, particularly in high-entropy alloy catalysts, emphasizing the role of hierarchical machine-learning architectures for accurate predictions. Additionally, I will discuss our work on developing machine learning interatomic potentials (MLPs) for single-element metals, metal oxides, and alloys under extreme conditions, focusing on melting behavior and phase properties at high temperatures and pressures. We have also refined our MLP models to capture dynamic surface interactions, such as CO2 and CO adsorption on MgO, using both static and molecular dynamics simulations. These models maintain high accuracy while significantly reducing computational costs compared to first-principles calculations. By enabling efficient and accurate simulations, this work supports broader community adoption, optimizes datasets for materials discovery, and extends the accessible time, size, and environmental conditions beyond the limits of experiments and traditional simulations.

machine learning↗
DOE OSTI2020

Predicting 3D moisture sorption behavior of materials from 1D investigations

Abstract Moisture in materials can be a source of future outgassing and exacerbate unwanted changes in physical and chemical properties. Here, we investigate the effect of sample size and shape on the moisture transport phenomena through a combined experimental and modeling approach. Several different materials varying in size and shape were investigated over a wide range of relative humidities (0–90%) and temperatures ( $$30{-}70^{\,\circ } \hbox {C}$$ 30 - 70 ∘ C ) using gravimetric type dynamic vapor sorption (DVS). A dynamic triple-mode sorption model, developed previously, was employed to describe the experimental results with good success; the model includes absorption, adsorption, pooling (clustering) of species, and molecular diffusion. Here we show that the full triple-mode sorption model is robust enough to predict the dynamic uptake and outgassing of 3-dimensional (3D) samples using parameters derived from quasi-1D samples. This successful demonstration on three different materials (filled polydimethylsiloxane (PDMS), unfilled PDMS, and ceramic inorganic composite) illustrates that the model is robust at describing the scale-independent physics and chemistry of moisture sorption and diffusion materials. This work demonstrates that while sorption mechanisms manifest in testing of all sample sizes, some of these mechanisms were so subtle that they were overlooked in our initial modeling and assessment, illustrating the importance of multi-scale experiments in the development of robust predictive capabilities. Our study also outlines the challenges and viable solutions for global optimization of a multi-parameter model. The ability to quantify moisture sorption and diffusion, independent of scale, using 1D lab-scale experiments enables prediction of long-term bulk materials behavior in real applications.

36 MATERIALS SCIENCE↗