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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Accurate noncovalent interactions in atomistic systems via quantum Drude oscillators

Accurately modeling polarization and van der Waals (vdW) interactions in atomistic systems typically requires high-level quantum-mechanical methods that are computationally expensive, hence limited in applicability. To address this challenge, efficient yet physically grounded models are needed—ones that not only enable accurate predictions but also provide insight into how noncovalent interactions scale in complex molecular and material systems. This review highlights the quantum Drude oscillator (QDO) model, a physically motivated and computationally efficient framework that captures the essential features of electronic response, including polarization and dispersion forces, across a wide range of chemical and material systems. We discuss how the QDO model quantitatively reproduces the polarization response of many-electron atoms and how key components of noncovalent interactions—exchange-repulsion, polarization, and dispersion—emerge naturally in QDO dimers. Furthermore, the model provides predictive scaling laws that elucidate trends in polarizability and dispersion across the periodic table and in molecular assemblies. By uniting interpretability, accuracy, and efficiency, the QDO model offers a versatile approach for modeling noncovalent interactions in systems ranging from isolated molecules to complex condensed phases and nanostructured materials.

Khabibrakhmanov, Almaz [Univ. of Luxembourg, Luxem↗
DOE OSTI2020

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling Exciton Trap Dynamics and Nonradiative Loss Pathways in Quantum Dots via Atomistic Simulations

Surface defects in colloidal quantum dots are a major source of nonradiative losses, yet the microscopic mechanisms underlying exciton trapping and recombination remain elusive. Here, we develop a model Hamiltonian based on atomistic electronic calculations to investigate exciton dynamics in CdSe/CdS core/shell QDs containing a single hole trap introduced by an unpassivated sulfur atom. By systematically varying the defect depth and reorganization energy, we uncover how defect-induced excitonic states mediate energy relaxation pathways. Our simulations reveal that a single localized defect can induce a rich spectrum of excitonic states, leading to multiple dynamical regimes, from slow, energetically off-resonant trapping to fast, cascaded relaxation through in-gap defect states. Crucially, we quantify how defect-induced polaron shifts and exciton-phonon couplings govern the balance between efficient radiative emission and rapid nonradiative decay. These insights clarify the microscopic origin of defect-assisted loss channels and suggest pathways for tailoring QD optoelectronic properties via surface and defect engineering.

Defects↗
DOE OSTI2023

Lightweight and effective tensor sensitivity for atomistic neural networks

Atomistic machine learning focuses on the creation of models that obey fundamental symmetries of atomistic configurations, such as permutation, translation, and rotation invariances. In many of these schemes, translation and rotation invariance are achieved by building on scalar invariants, e.g., distances between atom pairs. There is growing interest in molecular representations that work internally with higher rank rotational tensors, e.g., vector displacements between atoms, and tensor products thereof. Here, we present a framework for extending the Hierarchically Interacting Particle Neural Network (HIP-NN) with Tensor Sensitivity information (HIP-NN-TS) from each local atomic environment. Crucially, the method employs a weight tying strategy that allows direct incorporation of many-body information while adding very few model parameters. We show that HIP-NN-TS is more accurate than HIP-NN, with negligible increase in parameter count, for several datasets and network sizes. As the dataset becomes more complex, tensor sensitivities provide greater improvements to model accuracy. In particular, HIP-NN-TS achieves a record mean absolute error of 0.927 $\frac{\text{kcal}}{\text{mol}}$ for conformational energy variation on the challenging COMP6 benchmark, which includes a broad set of organic molecules. We also compare the computational performance of HIP-NN-TS to HIP-NN and other models in the literature.

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

Atomistic Study of the Plastic Deformation of Transition Metals and High Entropy Alloys (Final Report)

The overarching goal of the research performed in the studies supported by the DOE Basic Energy Sciences has been investigation of mechanical properties of materials employed in technologies related to energy production and conversion using atomic level computer modeling. The materials studied are body-centered-cubic (BCC) transition metals, including α-iron, their alloys and most recently also compositionally complex, so called high entropy alloys (HEAs), which are believed to become future materials in high temperature applications. The goal of the research has been to achieve atomic level understanding of the deformation and fracture of these materials starting from the electronic structure, coarse grain the material energetics into potentials describing interaction between atoms and proceed via atomistic modeling and mesoscopic models to investigation of the glide of dislocations and thus advance the theory of plastic deformation of studied materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reweighting configurations generated by transferable, machine learned models for protein sidechain backmapping

Multiscale modeling requires the linking of models at different levels of detail, with the goal of gaining accelerations from lower fidelity models while recovering fine details from higher resolution models. Communication across resolutions is particularly important in modeling soft matter, where tight couplings exist between molecular-level details and mesoscale structures. While multiscale modeling of biomolecules has become a critical component in exploring their structure and self-assembly, backmapping from coarse-grained to fine-grained, or atomistic, representations presents a challenge, despite recent advances through machine learning. A major hurdle, especially for strategies utilizing machine learning, is that backmappings can only approximately recover the atomistic ensemble of interest. We demonstrate conditions for which backmapped configurations may be reweighted to exactly recover the desired atomistic ensemble. By training separate decoding models for each sidechain type, we develop an algorithm based on normalizing flows and geometric algebra attention to autoregressively propose backmapped configurations for any protein sequence. Critical for reweighting with modern protein force fields, our trained models include all hydrogen atoms in the backmapping and make probabilities associated with atomistic configurations directly accessible. We also demonstrate, however, that reweighting is extremely challenging despite state-of-the-art performance on recently developed metrics and generation of configurations with low energies in atomistic protein force fields. Through detailed analysis of configurational weights, we show that machine-learned backmappings must not only generate configurations with reasonable energies, but also correctly assign relative probabilities under the generative model. These are broadly important considerations in generative modeling of atomistic molecular configurations.

Monroe, Jacob I. [Univ. of Arkansas, Fayetteville,↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗
DOE OSTI2021

Consortium for Production of Affordable Carbon Fibers in the United States (Final Report)

Various feedstocks from coal, petroleum and bio-based acrylonitrile were investigated for their ability to produce carbon fiber within DOE strength (250 ksi), modulus (25 Msi), elongation (>1%) and cost (≤$5/lb) goals for wider utilization in light-weight vehicles for the automotive industry. The overarching approach was to chemically identify and map the various feedstocks, understand the molecular transformation to produce carbon fiber precursors (mesophase for pitch materials of polyacrylonitrile polymer) and how these impact carbon fiber quality. Predictive models were built to relate the feedstocks and precursors to carbon fiber properties using atomistic and micro-scale modeling within ±15%. Using tow-level carbon fiber epoxy composites, carbon fibers tow-level predictions of elastic components of the carbon fibers were also successfully modeled within ±15%, using finite element analysis to confirm single filament carbon fiber properties. These models were validated with tow-level epoxy composite experiments from carbon fibers produced in this study.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Structural modeling of high-entropy oxides battery anodes using x-ray absorption spectroscopy

High-entropy oxides (HEOs) are single phase solid solutions where five or more metals share the same sublattice, giving rise to unexpected features in various fields of applications. Recently, HEOs have emerged as an alternative conversion electrode anode material for next-generation Li-ion batteries, where the combination of several different elements in a single solid solution can synergistically act to overcome some of its main drawbacks, improving performance. Due to their chemical complexity, x-ray absorption spectroscopy (XAS) emerges as an appropriate technique to study the electronic (x-ray absorption near edge structure, XANES) and local structure (extended x-ray absorption fine structure, EXAFS) of these compounds as a function of cycling. This work aims to highlight the capabilities of XAS as an element-specific probe to understand a material’s structure at the atomistic level through EXAFS modeling of (MgFeCoNiCuZn)O high-entropy system and how to extract valuable information about the bond distance, number of near neighbors, and local disorder, which are crucial to a full understanding of the electrochemical reaction mechanisms of such battery electrodes.

25 ENERGY STORAGE↗
NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗
DOE OSTI2021

Modeling Disorder in Pyrochlores and Other Anion-Deficient Fluorite Structural Derivative Oxides

Their very flexible chemistry gives oxide materials a richness in functionality and wide technological application. A specific group of oxides that have a structure related to fluorite but with less oxygen, termed anion-deficient fluorite structural derivatives and with pyrochlores being the most notable example, has been shown to exhibit a diversity of useful properties. For example, the possibility to undergo a transition from an ordered to disordered state allows these oxides to have high radiation tolerance. Atomistic-scale calculations in the form of molecular dynamics (MD) and density functional theory (DFT) have been extensively used to understand what drives this order/disorder transition. Here we give a brief overview of how atomistic-scale calculations are utilized in modeling disorder in pyrochlores and other anion-deficient fluorite structural derivatives. We discuss the modeling process from simple point defects to completely disordered structures, the dynamics during the disordering process, and the use of mathematical models to generate ordered solid-solution configurations. We also attempt to identify the challenges in modeling short range order and discuss future directions to more comprehensive models of the disordered structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗