Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Defect-driven anomalous transport in fast-ion conducting solid electrolytes

We report solid-state ionic conduction is a key enabler of electrochemical energy storage and conversion. The mechanistic connections between material processing, defect chemistry, transport dynamics, and practical performance are of considerable importance, but remain incomplete. Here, inspired by studies of fluids and biophysical systems, we re-examine anomalous diffusion in the iconic two-dimensional fast-ion conductors, the β- and β"-aluminas. Using large-scale simulations, we reproduce the frequency dependence of alternating-current ionic conductivity data. We show how the distribution of charge-compensating defects, modulated by processing, drives static and dynamic disorder, which lead to persistent sub-diffusive ion transport at macroscopic timescales. We deconvolute the effects of repulsions between mobile ions, the attraction between the mobile ions and charge-compensating defects, and geometric crowding on ionic conductivity. Our quantitative framework based on these model solid electrolytes connects their atomistic defect chemistry to macroscopic performance with minimal assumptions and enables mechanism-driven 'atoms-to-device' optimization of fast-ion conductors.

25 ENERGY STORAGE↗
DOE OSTI2024

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core–Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core–shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 10 1 –10 3 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core–cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ~100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

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

Deep-learning atomistic semi-empirical pseudopotential model for nanomaterials

The semi-empirical pseudopotential method (SEPM) has been widely applied to provide computational insights into the electronic structure, photophysics, and charge carrier dynamics of nanoscale materials. We present “DeepPseudopot”, a machine-learned atomistic pseudopotential model that extends the SEPM framework by combining a flexible neural network representation of the local pseudopotential with parameterized non-local and spin-orbit coupling terms. Trained on bulk quasiparticle band structures and deformation potentials from GW calculations, the model captures many-body and relativistic effects with very high accuracy across diverse semiconducting materials, as illustrated for silicon and group III-V semiconductors. DeepPseudopot’s accuracy, efficiency, and transferability make it well-suited for data-driven in silico design and discovery of novel optoelectronic nanomaterials.

Lin, Kailai [University of California, Berkeley, C↗
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2021

Porous Covalent Organic Polymers for Efficient Fluorocarbon‐Based Adsorption Cooling

Abstract 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.

Zheng, Jian↗
DOE OSTI2026

The structure, composition, and performance impact of a YSZ-GDC interdiffusion layer in solid oxide electrolysis cells

This study provides a combined experimental and computational investigation into the structure and impact of the cation interdiffusion layer that appears at the gadolinium doped ceria (GDC)/yttria stabilized zirconia (YSZ) interface in solid oxide electrolysis cells (SOECs). Scanning transmission electron microscopy (STEM) illustrates that a ∼0.4 μm interdiffusion layer (IDL) with an intermixed cation distribution and fine grain size forms upon sintering. STEM identifies that the interdiffusion layer exists in the cubic fluorite structure despite changes in cation composition. The interdiffusion layer microstructure formed during sintering does not change during SOEC testing at either 1.3V or heightened voltage pulse testing. Modeling predicts that ionic conductivity may decrease in the interdiffusion layer due to Coulombic trapping between mobile oxygen vacancies and excess Gd 3+ acceptor dopants. Yet, the density and continuous nature of the layer should benefit cell stability by substantially reducing the formation of SrZrO 3 , which is corroborated by STEM and Synchrotron X-ray diffraction (XRD). We conclude that the interdiffusion layer acts as a beneficial barrier to Sr diffusion, when operating in a regime where electrolyte void formation is not observed.

organic↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Model-free estimation of completeness, uncertainties, and outliers in atomistic machine learning using information theory

Abstract An accurate description of information is relevant for a range of problems in atomistic machine learning (ML), such as crafting training sets, performing uncertainty quantification (UQ), or extracting physical insights from large datasets. However, atomistic ML often relies on unsupervised learning or model predictions to analyze information contents from simulation or training data. Here, we introduce a theoretical framework that provides a rigorous, model-free tool to quantify information contents in atomistic simulations. We demonstrate that the information entropy of a distribution of atom-centered environments explains known heuristics in ML potential developments, from training set sizes to dataset optimality. Using this tool, we propose a model-free UQ method that reliably predicts epistemic uncertainty and detects out-of-distribution samples, including rare events in systems such as nucleation. This method provides a general tool for data-driven atomistic modeling and combines efforts in ML, simulations, and physical explainability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

Model ensembles are effective tools for estimating prediction uncertainty in deep learning atomistic force fields. However, their widespread adoption is hindered by high computational costs and overconfident error estimates. In this work, we address these challenges by leveraging distributions of per-sample errors obtained during training and employing a distance-based similarity search in the model latent space. Our method, which we call LTAU (Loss Trajectory Analysis for Uncertainty), efficiently estimates the full probability distribution function of errors for any test point using the logged training errors, achieving speeds that are 2–3 orders of magnitudes faster than typical ensemble methods and allowing it to be used for tasks where training or evaluating multiple models would be infeasible. We apply LTAU towards estimating parametric uncertainty in atomistic force fields (LTAU-FF), demonstrating that it produces well-calibrated confidence intervals and predicts errors that correlate strongly with the true errors for data near the training domain. Furthermore, we show that the errors predicted by LTAU-FF can be used in practical applications for detecting out-of-domain data, tuning model performance, and predicting failure during simulations. We believe that LTAU will be a valuable tool for uncertainty quantification in atomistic force fields and is a promising method that should be further explored in other domains of machine learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular dynamics studies of knotted polymers

Molecular dynamics calculations have been used to explore the influence of knots on the strength of a polymer strand. In particular, the mechanism of breaking 31, 41, 51, and 52 prime knots has been studied using two very different models to represent the polymer: (1) the generic coarse-grained (CG) bead model of polymer physics and (2) a state-of-the-art machine learned atomistic neural network (NN) potential for polyethylene derived from electronic structure calculations. While there is a broad overall agreement between the results on the influence of the pulling rate on chain rupture based on the CG and atomistic NN models, for the simple 31 and 41 knots, significant differences are found for the more complex 51 and 52 knots. Notably, in the latter case, the NN model more frequently predicts that these knots can break not only at the crossings at the entrance/exit but also at one of the central crossing points. The relative smoothness of the CG potential energy surface also leads to stabilization of tighter knots compared to the more realistic NN model.

DelloStritto, Mark (ORCID:0000000206785860)↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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 OSTI2023

Multi-scale modeling of the evolution of structure and properties in materials for nuclear energy applications

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. These modeling efforts make extensive of MOOSE (Multiphysics Object-Oriented Simulation Environment), a general-purpose open source finite element framework developed at INL. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Multi-scale modeling of the evolution of structure and properties in materials for nuclear energy applications

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. These modeling efforts make extensive of MOOSE (Multiphysics Object-Oriented Simulation Environment), a general-purpose open source finite element framework developed at INL. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗