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 271 records · Page 15

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↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

Data-driven modeling of dislocation mobility from atomistics using physics-informed machine learning

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

36 MATERIALS SCIENCE↗

Grain boundary segregation beyond the dilute limit: Separating the two contributions of site spectrality and solute interactions

Solute segregation at grain boundaries (GBs) is known to have a profound impact on material properties, and as such is becoming routinely used as an element in alloy design. Beyond the dilute limit, the extent of solute GB segregation is known to be concentration dependent. Using hybrid Monte Carlo/Molecular Statics simulations of Mg segregation in Al, in this work we decouple the two contributions to this composition dependence: (i) spectrality of atomic environments at the boundary and (ii) solute-solute interactions. Although only contribution (ii) is typically considered in the literature, we argue that both contributions are equally important to understand concentration dependence and correctly quantify GB solute segregation in a binary alloy. Finally, a thermodynamic segregation isotherm is outlined that accounts for both the spectrality of grain boundary sites and solute-solute interactions. Unlike classical isotherms like those of McLean or Fowler-Guggenheim, which can be successfully fitted to GB segregation data only over a limited range of composition and temperature, our proposed model is shown to be accurate across the composition-temperature space.

36 MATERIALS SCIENCE↗

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗

Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation

The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.

Materials science↗

Sampling lattices in semi-grand canonical ensemble with autoregressive machine learning

Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.

36 MATERIALS SCIENCE↗

Machine learning potential assisted exploration of complex defect potential energy surfaces

Abstract Atomic-scale defects generated in materials under both equilibrium and irradiation conditions can significantly impact their physical and mechanical properties. Unraveling the energetically most favorable ground-state configurations of these defects is an important step towards the fundamental understanding of their influence on the performance of materials ranging from photovoltaics to advanced nuclear fuels. Here, using fluorite-structured thorium dioxide (ThO 2 ) as an exemplar, we demonstrate how density functional theory and machine learning interatomic potential can be synergistically combined into a powerful tool that enables exhaustive exploration of the large configuration spaces of small point defect clusters. Our study leads to several unexpected discoveries, including defect polymorphism and ground-state structures that defy our physical intuitions. Possible physical origins of these unexpected findings are elucidated using a local cluster expansion model developed in this work.

36 MATERIALS SCIENCE↗

Nanocomposite structure of two-line ferrihydrite powder from total scattering

Ferrihydrite is one of the most important iron-containing minerals on Earth. Yet determination of its atomic-scale structure has been frustrated by its intrinsically poor crystallinity. The key difficulty is that physically-different models can appear consistent with the same experimental data. Using X-ray total scattering and a nancomposite reverse Monte Carlo approach, we evaluate the two principal contending models—one a multi-phase system without tetrahedral iron(III), and the other a single phase with tetrahedral iron(III). Our methodology is unique in considering explicitly the complex nanocomposite structure the material adopts: namely, crystalline domains embedded in a poorly-ordered matrix. The multi-phase model requires unphysical structural rearrangements to fit the data, whereas the single-phase model accounts for the data straightforwardly. Hence the latter provides the more accurate description of the short- and intermediate-range order of ferrihydrite. We discuss how this approach might allow experiment-driven (in)validation of complex models for important nanostructured phases beyond ferrihydrite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

Modeling Si/SiGe quantum dot variability induced by interface disorder reconstructed from multiperspective microscopy

SiGe heteroepitaxial growth yields pristine host material for quantum dot qubits, but residual interface disorder can lead to qubit-to-qubit variability that might pose an obstacle to reliable SiGe-based quantum computing. By convolving data from scanning tunneling microscopy and high-angle annular dark field scanning transmission electron microscopy, we reconstruct 3D interfacial atomic structure and employ an atomistic multi-valley effective mass theory to quantify qubit spectral variability. The results indicate (1) appreciable valley splitting (VS) variability of ~50% owing to alloy disorder and (2) roughness-induced double-dot detuning bias energy variability of order 1–10 meV depending on well thickness. For measured intermixing, atomic steps have negligible influence on VS, and uncorrelated roughness causes spatially fluctuating energy biases in double-dot detunings potentially incorrectly attributed to charge disorder. Our approach yields atomic structure spanning orders of magnitude larger areas than post-growth microscopy or tomography alone, enabling more holistic predictions of disorder-induced qubit variability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scale-invariant machine-learning model accelerates the discovery of quaternary chalcogenides with ultralow lattice thermal conductivity

We design an advanced machine-learning (ML) model based on crystal graph convolutional neural network that is insensitive to volumes (i.e., scale) of the input crystal structures to discover novel quaternary chalcogenides, AMM'Q 3 (A/M/M' = alkali, alkaline earth, post-transition metals, lanthanides, and Q = chalcogens). These compounds are shown to possess ultralow lattice thermal conductivity (κ l ), a desired requirement for thermal-barrier coatings and thermoelectrics. Upon screening the thermodynamic stability of ~1 million compounds using the ML model iteratively and performing density-functional theory (DFT) calculations for a small fraction of compounds, we discover 99 compounds that are validated to be stable in DFT. Taking several DFT-stable compounds, we calculate their κ l using Peierls–Boltzmann transport equation, which reveals ultralow κ l (<2 Wm -1 K -1 at room temperature) due to their soft elasticity and strong phonon anharmonicity. Our work demonstrates the high efficiency of scale-invariant ML model in predicting novel compounds and presents experimental-research opportunities with these new compounds.

36 MATERIALS SCIENCE↗