Uncertainty Quantification and Active Learning in Atomistic Computations.
Abstract not provided.
Najm, Habib↗
SEARCH · Search NASA
Results for “atomistic computations”
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.
Atomistic computer simulations of dissolution of inorganic glasses.
Abstract not provided.
Rimsza, Jessica↗
Deciphering the PMI Surface Chemistry of Lithium-Based PFCs and its Effects on High Performance Plasmas in NSTX
The primary aims or objectives of the proposed work here focus on the dynamic measurement (Illinois PI Allain) and atomistic-based multi-scale computational atomistic simulations (UTK PI: Wirth) of: 1) lithium coatings with re-deposited hydrogenated carbon surfaces and high-Z substrates, 2) lithium coatings on boronized ATJ graphite surfaces and high-Z substrates and 3) lithium coatings on ATJ graphite and high-Z substrates with variation in temperature. Three main primary tasks are proposed to address each of the primary objectives listed above: 1) dynamic in-situ irradiation with D+, He+ species of lithium coatings with variation in interface substrate morphology (e.g. smooth vs rough, fuzz vs nanostructured); 2) validation and coupling with multi-scale computational simulations connecting the irradiated surface to the plasma edge; 3) extrapolation to long-pulse conditions (e.g. function of flux vs fluence studies). In primary objective #3, temperature-based effects will be studied to evaluate three primary mechanisms: a) D uptake and recycling, b) erosion, and c) surface melting and evaporation as a function of substrate morphology and chemistry.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Computational Approaches for Li-O2 Battery Design
Threshold energy densities for general aviation electric aircraft are 400 Wh/kg with more ambitious air vehicles having significantly higher requirements. Li-O2 batteries, with the highest theoretical capacity, are one of the few “beyond Li-ion” chemistries that might satisfy the extraordinary specific capacity as well as specific power requirements of electric aircraft. However, side reactions at interfaces, in particular at the cathode, over the charge-discharge cycles result in very short cycle-life and dramatic reduction of capacity. Addressing these issues, in addition to others, are crucial for realizing practical, high performance Li-O2 batteries. In this talk, we discuss atomistic computational work to understand and mitigate some of the issues, including those at interfaces, that affect the Li-O2 electrochemistry. To start, we discuss the deposition mechanisms, both surface and solution based, of Li2O2 and their dependence on external potential. Next, we explore molten salt electrolytes as a stable alternative to organic electrolytes and approaches taken to develop new practical molten salt eutectic mixtures. Finally, we address the issue of reactive carbon-cathodes and possible cathode-candidates that were identified throughput high-throughput computations.
Li-Air battery↗
Heats of Segregation and Segregation Profiles of BCC Metals
The composition of metal alloy surfaces is often different from that of the bulk. Some alloys exhibit surface segregation, where one or more species reside preferentially at or near the surface. A detailed understanding of this behavior is necessary to correctly model such phenomena as adhesion or catalysis. Several phenomenological approaches to the problem have been put forward, falling into two broad categories: Thermodynamic approaches, where the equilibrium distribution of chemical species is computed. Atomistic approaches, where the tendency of a species to segregate is determined by computation of the energies of single atoms of that species in bulk and surface environments.
Good, Brian S.↗
Thermodynamics of grain-boundary phases
The past decade has seen a significant increase in research efforts aimed at understanding the thermodynamics of low-dimensional phases existing in many materials systems, ranging from two-dimensional materials to core regions of extended defects in crystalline solids. We review the current status of theoretical, computational, and experimental research on the “defect phases,” focusing on grain boundaries (GBs) in elemental and multicomponent polycrystalline materials. After reviewing the generalized concept of a phase of any dimensionality, we discuss recent progress in atomistic computer simulations of GB phase transformations and phase coexistences, including the observation of one-dimensional defects separating GB phases (defects in defects). Computational predictions compare well with experimental observations of multiple GB phases and segregation-induced phase transformations. An intriguing open question of GB thermodynamics is whether the GB free energy can be driven to a zero value by increasing solute segregation. We review recent efforts to understand this ultimate thermodynamic stabilization of GB phases and the possible polycrystalline microstructures that may arise. An outlook for future research in the field is discussed.
Materials science↗
Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials’ design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum–niobium–titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Atomic-level mechanisms of short-circuit diffusion in materials
This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.
36 MATERIALS SCIENCE↗
Atomistic modeling of metal–nonmetal interphase boundary diffusion
Atomistic computer simulations are applied to investigate the atomic structure, thermal stability, and diffusion processes in Al–Si interphase boundaries as a prototype of metal–ceramic interfaces in composite materials. Some of the most stable orientation relationships between the phases found in this work were previously observed in epitaxy experiments. Here, a non-equilibrium interface can transform to a more stable state by a mechanism that we call interface-induced recrystallization. Diffusion of both Al and Si atoms in stable Al–Si interfaces is surprisingly slow compared with diffusion of both elements in Al grain boundaries but can be accelerated in the presence of interface disconnections. A qualitative explanation of the sluggish interphase boundary diffusion is proposed. Atomic mechanisms of interphase boundary diffusion are similar to those in metallic grain boundaries and are dominated by correlated atomic rearrangements in the form of strings and rings of collectively moving atoms.
36 MATERIALS SCIENCE↗
Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces
As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Combined molecular and spin dynamics simulation of BCC iron with vacancy defects
Utilizing an atomistic computational model, which handles both translational and spin degrees of freedom, combined molecular and spin dynamics simulations have been performed to investigate the effect of vacancy defects on spin wave excitations in ferromagnetic iron. Fourier transforms of space- and time-displaced correlation functions yield the dynamic structure factor, providing characteristic frequencies and lifetimes of the spin wave modes. A comparison of the system with a 5% vacancy concentration with pure lattice data shows a decrease in frequency and a decrease in lifetime for all transverse spin wave excitations observed. In addition, the clearly defined transverse spin wave excitations are distorted with the introduction of vacancy defects, and we observe reduced excitation lifetimes due to increased magnon–magnon scattering. We observe further evidence of increased magnon–magnon scattering, as the peaks in the longitudinal spin wave spectrum become less distinct. Finally, similar impacts are observed in the vibrational subsystem, with a decrease in characteristic phonon frequency and flattening of lattice excitation signals due to vacancy defects.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020
Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
An idealized model for the crystal structure of intermetallic compounds isostructural with Mg 3 Cr 2 Al 18
Mg 3 Cr 2 Al 18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\overline{3}m$. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, $A, B, C, D,$ and $E$. With this, a general description for MCA structures can be written as $A^{8a}_{1}B^{16c}_{2}C^{16d}_{2}D^{48f}_{6}E^{96g}_{12}$, where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the $D$ and $E$ sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds.
36 MATERIALS SCIENCE↗
Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma
Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D
Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.
74 ATOMIC AND MOLECULAR PHYSICS↗
Boundary Condition for Modeling Semiconductor Nanostructures
A recently proposed boundary condition for atomistic computational modeling of semiconductor nanostructures (particularly, quantum dots) is an improved alternative to two prior such boundary conditions. As explained, this boundary condition helps to reduce the amount of computation while maintaining accuracy.
Lee, Seungwon↗
Uncertainty quantification for neural network potential foundation models
Abstract For neural network potentials (NNPs) to gain widespread use, researchers must be able to trust model outputs. However, the blackbox nature of neural networks and their inherent stochasticity are often deterrents, especially for foundation models trained over broad swaths of chemical space. Uncertainty information provided at the time of prediction can help reduce aversion to NNPs. In this work, we detail two uncertainty quantification (UQ) methods. Readout ensembling, by finetuning the readout layers of an ensemble of foundation models, provides information about model uncertainty, while quantile regression, by replacing point predictions with distributional predictions, provides information about uncertainty within the underlying training data. We demonstrate our approach with the MACE-MP-0 model, applying UQ to the foundation model and a series of finetuned models. The uncertainties produced by the readout ensemble and quantile methods are demonstrated to be distinct measures by which the quality of the NNP output can be judged.
36 MATERIALS SCIENCE↗