Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomistic 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 55 records · Page 3

Emerging Atomistic Modeling Methods for Heterogeneous Electrocatalysis

Heterogeneous electrocatalysis lies at the center of various technologies that could help enable a sustainable future. However, its complexity makes it challenging to accurately and efficiently model at an atomic level. Herein, we review emerging atomistic methods to simulate the electrocatalytic interface with special attention devoted to the components/effects that have been challenging to model, such as solvation, electrolyte ions, electrode potential, reaction kinetics, and pH. Additionally, we review relevant computational spectroscopy methods. Then, we showcase several examples of applying these methods to understand and design catalysts relevant to green hydrogen. We also offer experimental views on how to bridge the gap between theory and experiments. Finally, we provide some perspectives on opportunities to advance the field.

36 MATERIALS SCIENCE↗

Configurational-Bias Monte Carlo Back-Mapping Algorithm for Efficient and Rapid Conversion of Coarse-Grained Water Structures into Atomistic Models

Coarse-grained molecular dynamics (MD) simulations represent a powerful approach to simulate longer time scale and larger length scale phenomena than those accessible to all-atom models. The gain in efficiency, however, comes at the cost of atomistic details. The reverse transformation, also known as back-mapping, of coarse grained beads into their atomistic constituents represents a major challenge. Most existing approaches are limited to specific molecules or specific force-fields and often rely on running a long time atomistic MD of the back-mapped configuration to arrive at an optimal solution. Such approaches are problematic when dealing with systems with high diffusion barriers. Here, we introduce a new extension of the configurational-bias-Monte-Carlo (CBMC) algorithm, which we term the crystalline-configurational-bias-Monte-Carlo (C-CBMC) algortihm, that allows rapid and efficient conversion of a coarse-grained model back into its atomistic representation. Although the method is generic, we use a coarse-grained water model as a representative example and demonstrate the back-mapping or reverse transformation for model systems ranging from the ice-liquid water interface to amorphous and crystalline ice configurations. A series of simulations using the TIP4P/Ice model are performed to compare the new CBMC method to several other standard Monte Carlo and Molecular Dynamics based back-mapping techniques. In all the cases, the C-CBMC algorithm is able to find optimal hydrogen bonded configuration many thousand evaluations/steps sooner than the other methods compared within this paper. For crystalline ice structures such as a hexagonal, cubic, and cubic-hexagonal stacking disorder structures, the C-CBMC was able to find structures that were between 0.05 and 0.1 eV/water molecule lower in energy than the ground state energies predicted by the other methods. Detailed analysis of the atomistic structures show a significantly better global hydrogen positioning when contrasted with the existing simpler back-mapping methods. The error in the RDFs of backmapped relative to reference configuration for the CCBMC, MD, and MC were found to be 6.9, 8.7, and 12.9 respectively for the hexagonal system. For the cubic system the relative error of the RDFs for the CCBMC, MD, and MC were found to be 18.2, 34.6, and 39.0 respectively. Finally, our results demonstrate the efficiency and efficacy of our new back-mapping approach, especially for crystalline systems where simple force-field based relaxations have a tendency to get trapped in local minima.

74 ATOMIC AND MOLECULAR PHYSICS↗

Phase transformations at interfaces: Observations from atomistic modeling

Here, we review the recent progress in theoretical understanding and atomistic computer simulations of phase transformations in materials interfaces, focusing on grain boundaries (GBs) in metallic systems. Recently developed simulation approaches enable the search and structural characterization of GB phases in single-component metals and binary alloys, calculation of thermodynamic properties of individual GB phases, and modeling of the effect of the GB phase transformations on GB kinetics. Atomistic simulations demonstrate that the GB transformations can be induced by varying the temperature, loading the GB with point defects, or varying the amount of solute segregation. The atomic-level understanding obtained from such simulations can provide input for further development of thermodynamics theories and continuous models of interface phase transformations while simultaneously serving as a testing ground for validation of theories and models. They can also help interpret and guide experimental work in this field.

36 MATERIALS SCIENCE↗

Harness the power of atomistic modeling and deep learning in biofuel separation

Biofuels offer a remarkable, sustainable energy source for a future of clean energy. The development of efficient biofuel separation plays a crucial role in achieving cost-effective utilization of biofuel. In this chapter, we provide an overview of the recent advancements in atomistic-level modeling and deep learning in the rational design of novel, efficient biofuel separation. The fundamental principles of quantum and statistical mechanics are covered in appropriate detail to highlight their underlying differences in theory. The methodologies of several molecular representations and deep learning algorithms applicable to biofuel separation are briefly demonstrated as well. The applications, successes, and risks of employing density functional theory, ab initio molecular dynamics, classical molecular dynamics, and deep learning are provided to showcase their recent accomplishments in biofuel separation as well as potential improvements in both methodology and application. Lastly, a vision for the future growth of these methods is illustrated.

deep learning, artificial intelligence, biofuels, ↗

Atomistic modelling of scattering data in the Collaborative Computational Project for Small Angle Scattering (CCP-SAS)

The capabilities of current computer simulations provide a unique opportunity to model small-angle scattering (SAS) data at the atomistic level, and to include other structural constraints ranging from molecular and atomistic energetics to crystallography, electron microscopy and NMR. This extends the capabilities of solution scattering and provides deeper insights into the physics and chemistry of the systems studied. Realizing this potential, however, requires integrating the experimental data with a new generation of modelling software. To achieve this, the CCP-SAS collaboration (http://www.ccpsas.org/) is developing open-source, high-throughput and user-friendly software for the atomistic and coarse-grained molecular modelling of scattering data. Robust state-of-the-art molecular simulation engines and molecular dynamics and Monte Carlo force fields provide constraints to the solution structure inferred from the small-angle scattering data, which incorporates the known physical chemistry of the system. The implementation of this software suite involves a tiered approach in which GenApp provides the deployment infrastructure for running applications on both standard and high-performance computing hardware, and SASSIE provides a workflow framework into which modules can be plugged to prepare structures, carry out simulations, calculate theoretical scattering data and compare results with experimental data. GenApp produces the accessible web-based front end termed SASSIE-web , and GenApp and SASSIE also make community SAS codes available. Applications are illustrated by case studies: (i) inter-domain flexibility in two- to six-domain proteins as exemplified by HIV-1 Gag, MASP and ubiquitin; (ii) the hinge conformation in human IgG2 and IgA1 antibodies; (iii) the complex formed between a hexameric protein Hfq and mRNA; and (iv) synthetic `bottlebrush' polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Understanding of Nitrogen Oxide Fixation of Water-Lean Carbon Capture Solvents by Atomistic Modeling

Nitrogen oxides, present in flue gas, can cause negative impacts on amine carbon capture solvents by the formation of heat-stable salts and suspected carcinogens. Thus, to maximize the performance of water-lean solvents, a better understanding of this process in these systems is necessary. Here, a computational study for the fixation of the CO 2 capture solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (EEMPA) to nitramine/nitrosamine was conducted. The first step involves the dissociation of the NH bond of EEMPA, in which the homolytic mechanism is energetically more favorable than the heterolytic mechanism. The second step involves radical recombination to form N–N bonds. While NO 2 directly reacts with EEMPA, NO has almost no effect. However, in the presence of O 2 , fixation of EEMPA by NO is enhanced via the formation of N 2 O 4 species. Finally, low reaction energies indicate that the formation of nitramine/nitrosamine may be a reversible process, suggesting that EEMPA could be recovered under thermal stripping conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic modeling of tensile deformation and fracture of carbon fibers: Nanoscale stress redistribution, effect of local structural characteristics and nanovoids

The mechanisms of tensile deformation and fracture of carbon fibers are investigated in large-scale molecular dynamics simulations. The simulations performed for computational samples with different structural characteristics enable analysis of the structural sensitivity of the mechanical properties. It is found that the carbon fiber nanostructure can be represented by a fine (nanoscale) mixture of stiff and soft regions. Despite the nanoscale structural heterogeneity of the computational fibers, their fracture strain and tensile strength exceed the experimental values by almost an order of magnitude. Here, this observation is explained by the ability of the carbon fiber nanostructure, largely consisting of curved and folded layers of turbostratic carbon, to exhibit continuous nanoscale load redistribution during the deformation. The transient, fluid nature of a percolating cluster of stress-holding structural elements prevents stress localization at the early stage of the deformation and enables the fiber nanostructure to adapt to the increasing load. Overall, the results of the simulations suggest that the design of advanced carbon fiber manufacturing and post-processing methods focused on the elimination of critical mesoscopic structural defects (internal voids, foreign inclusions, and surface flaws) can yield a substantial increase in the strength and fracture strain of carbon fibers.

36 MATERIALS SCIENCE↗

An atomistic model for the thermal resistance of a liquid–solid interface

The thermal resistance associated with the interface between a solid and a liquid is analysed from an atomistic point of view. Partial evaluation of the associated Green–Kubo integral elucidates the various factors governing heat transport across the interface and leads to a quantitative model for the thermal resistance in terms of atomistic-level system parameters. The model is validated using molecular dynamics simulations.

Mechanics↗

Grain-Boundary Resistance in Copper Interconnects: From an Atomistic Model to a Neural Network

Orientation effects on the specific resistance of copper grain boundaries are studied systematically with two different atomistic tight-binding methods. A methodology is developed to model the specific resistance of grain boundaries in the ballistic limit using the embedded atom model, tight- binding methods, and nonequilibrium Green’s functions. The methodology is validated against first-principles calculations for thin films with a single coincident grain boundary, with 6.4% deviation in the specific resistance. A statistical ensemble of 600 large, random structures with grains is studied. For structures with three grains, it is found that the distribution of specific resistances is close to normal. Finally, a compact model for grain-boundary-specific resistance is constructed based on a neural network.

74 ATOMIC AND MOLECULAR PHYSICS↗

Atomistic model of xenon gas bubble re-solution rate due to thermal spike in uranium oxide

Atomistic simulations are performed to study the response of Xe gas bubbles in UO 2 to ionizing fission products through the thermal spike approximation. A portion of the total electronic stopping power (S e ) is taken as the thermal spike energy through a ratio variable ζ. The thermal spike energy causes extreme melting within the fission track cylindrical region. Molecular dynamics is employed to quantify the probability of a Xe gas atom to be re-solved (re-dissolved) back into the UO 2 matrix. Subsequently, a re-solution model is developed and parametrized as a function of bubble radius (R), off-centered distance (r), and thermal spike energy (ζS e ). The off-centered distance measures the shift of the thermal spike axis from the bubble center. To evaluate the re-solution model, independent fission product yield of U-235 fission due to thermal neutrons (0.0253 eV), taken from the JEFF-3.3 database, is used. The kinetic energy of the fission products is taken from the EXFOR database. Subsequently, the decay of Se over distance for each fission product is simulated. Finally, the evaluated re-solution rate (re-solution probability per second) is presented as a function of bubble radius for a range of ζ.

74 ATOMIC AND MOLECULAR PHYSICS↗

Reshaping and sintering of 3D fcc metal nanoclusters: Stochastic atomistic modeling with realistic surface diffusion kinetics

Far-from-equilibrium evolution of metallic nanocluster shapes is highly sensitive to the atomistic-level details of surface diffusion for diverse local surface configurations. A stochastic model was developed incorporating realistic values for the multiple diffusion barriers (contrasting previous unrealistic generic prescriptions) based upon insights from homoepitaxial film growth. Kinetic Monte Carlo simulation then elucidates the conversion of Ag nanocubes to Wulff polyhedra mediated by nucleation of new {100} facets, the pinch-off of sufficiently elongated Ag nanorods, and key aspects of sintering for orientationally aligned Ag and Au nanoclusters. In conclusion, the time scale for sintering of Au nanoclusters observed in high-resolution transmission electron microscopy studies was also recovered.

36 MATERIALS SCIENCE↗

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↗

Scalable training of trustworthy and energy-efficient predictive graph foundation models for atomistic materials modeling: a case study with HydraGNN

We present our work on developing and training scalable, trustworthy, and energy-efficient predictive graph foundation models (GFMs) using HydraGNN, a multi-headed graph convolutional neural network architecture. HydraGNN expands the boundaries of graph neural network (GNN) computations in both training scale and data diversity. It abstracts over message passing algorithms, allowing both reproduction of and comparison across algorithmic innovations that define nearest-neighbor convolution in GNNs. This work discusses a series of optimizations that have allowed scaling up the GFMs training to tens of thousands of GPUs on datasets consisting of hundreds of millions of graphs. Our GFMs use multitask learning (MTL) to simultaneously learn graph-level and node-level properties of atomistic structures, such as energy and atomic forces. Using over 154 million atomistic structures for training, we illustrate the performance of our approach along with the lessons learned on two state-of-the-art US Department of Energy (US-DOE) supercomputers, namely the Perlmutter petascale system at the National Energy Research Scientific Computing Center and the Frontier exascale system at Oak Ridge Leadership Computing Facility. The HydraGNN architecture enables the GFM to achieve near-linear strong scaling performance using more than 2000 GPUs on Perlmutter and 16,000 GPUs on Frontier.

97 MATHEMATICS AND COMPUTING↗

Atomistic modeling of nanoparticle generation in short pulse laser ablation of thin metal films in water

Laser ablation in liquids is actively used for generation of clean colloidal nanoparticles with unique shapes and functionalities. The fundamental mechanisms of the laser ablation in liquids and the key processes that control the nanoparticle structure, composition, and size distribution, however, are not yet fully understood. In this paper, we report the results of first atomistic simulations of laser ablation of metal targets in liquid environment. A model combining a coarse-grained representation of the liquid environment (parameterized for water), a fully atomistic description of laser interactions with metal targets, and acoustic impedance matching boundary conditions is developed and applied for simulation of laser ablation of a thin silver film deposited on a silica substrate. The simulations, performed at two laser fluences in the regime of phase explosion, predict a rapid deceleration of the ejected ablation plume and the formation of a dense superheated molten layer at the water-plume interface. The water in contact with the hot metal layer is brought to the supercritical state and transforms into an expanding low density metal-water mixing region that serves as a precursor for the formation of a cavitation bubble. Two distinct mechanisms of the nanoparticle formation are predicted in the simulations: (1) the nucleation and growth of small (mostly ≤10nm) nanoparticles in the metal-water mixing region and (2) the formation of larger (tens of nm) nanoparticles through the breakup of the superheated molten metal layer triggered by the emergence of complex morphological features attributed to the Rayleigh-Taylor instability of the interface between at the superheated metal layer and the supercritical water. The first mechanism is facilitated by the rapid cooling of the growing nanoparticles in the supercritical water environment, resulting in solidification of the nanoparticles located in the upper part of the mixing region on the timescale of nanoseconds. The computational prediction of the two mechanisms of nanoparticle formation yielding nanoparticles with different characteristic sizes is consistent with experimental observations of two distinct nanoparticle populations appearing at different stages of the ablation process.

Chemistry↗

Platinum nanoparticle compression: Combining in situ TEM and atomistic modeling

The mechanical behavior of nanoparticles governs their performance and stability in many applications. However, the small sizes of technologically relevant nanoparticles, with diameters in the range of 10 nm or less, significantly complicate experimental examination. These small nanoparticles are difficult to manipulate onto commercial test platforms and deform at loads that are below the typical noise floor of the testing instruments. Here, we synthesized small platinum nanoparticles directly onto a mechanical tester and used a modified nanomanipulator to enhance load resolution to the nanonewton scale. We demonstrated the in situ compression of an 11.5-nm platinum nanoparticle with simultaneous high-resolution measurements of load and particle morphology. Molecular dynamics simulations were performed on similarly sized particles to achieve complementary measurements of load and morphology, along with atomic resolution of dislocations. The experimental and simulation results revealed comparable values for the critical resolved shear stress for failure, 1.28 and 1.15 GPa, respectively. Altogether, this investigation demonstrated the promise of, and some initial results from, the combination of atomistic simulations and in situ experiments with an unprecedented combination of high spatial resolution and high load resolution to understand the behavior of metal nanoparticles under compression.

42 ENGINEERING↗

Atomistic Modeling of Semiconductors: Si, C, and 3C-SiC

An ongoing task of the Computational Materials Group (CMG) at the NASA Glenn Research Center is to enhance the role of atomistic simulations based on quantum-approximate methods in the study of new materials and their properties. One of the main goals of the activity continues to be breaching limitations that arise from the natural balance between accuracy, range of application, and computational simplicity. Whether that balance can be maintained while breaking new ground depends on the methods available with a minimum of constraints and limitations for the study of the energetics of arbitrary systems. The main tool used in CMG research, the Bozzolo- Ferrante-Smith (BFS) method for alloys, has no inherent constraint in its formulation, a feature that has allowed for successful research on various topics. In this article, we report on the latest development of the CMG program, namely, the extension and application of the BFS method to compound semiconductors, a departure from our previous research based primarily on metallic alloys.

Bozzolo, Guillermo↗

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗

Atomistic modeling of out-of-pile xenon diffusion by vacancy clusters in UO 2

We report that diffusion of Xe atoms in UO 2 fuel is important for nuclear fuel performance, and is enabled by interaction with U and O vacancies. Previous work using atomistic calculations based on density functional theory (DFT) and empirical potentials (EP) focused on the role of small vacancy clusters (Xe U2Oy , y = 0, 1, 2) for Xe transport, but the model predictions failed to reproduce the experimental observations by consistently underestimating Xe diffusivity. In this work, where we focused on out-of-pile conditions (i.e. in the absence of irradiation), we have explored two of the uncertainties associated with the model: the DFT methodology and the types of clusters considered. Additionally, we found that using the energy barriers obtained with GGA + U DFT allows to reconcile simulation results and experimental observations in the intrinsic regime (i.e. out-of-pile conditions). The Xe UO cluster is the preferred configuration at high temperature, while the mobile Xe U2O cluster takes over below 1720 K. The latter cluster is also the main contributor to Xe diffusion and good agreement is obtained with current models empirically fitted to experiments for diffusion at high temperature. A simple expression that captures most of Xe diffusivity is proposed, relying on Xe UO and Xe U2O clusters only. We also determined the out-of-pile concentration and diffusivity of extended clusters Xe U z O y , 3 ≤ x ≤ 9, and determined that these do not contribute significantly to the overall diffusivity under intrinsic conditions. However, we speculate that the Xe U 4 O 3 and Xe U 8 O 9 clusters, which have relatively low migration barrier, may play a role under irradiation conditions, because the radiation-induced enhanced vacancy concentration would rapidly increase the fraction of large Xe-vacancy clusters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗