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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2026

Momentum-Resolved Direct Observation of Chiral Phonons in Elemental Tellurium

Chiral phonons carry finite phonon angular momentum, yet their momentum-resolved behavior in acoustic branches remains largely unexplored. In this work, we report the direct detection of chiral acoustic phonons in elemental Te. Through a combination of high-energy-resolution inelastic x-ray scattering experiments and atomistic modeling, we confirm that longitudinal acoustic modes develop pronounced circular atomic motion by mixing with transverse modes. Our results establish a direct link between unexpected scattering intensity in nominally forbidden geometries and circular phonons, offering a framework for characterizing momentum-resolved phonon chirality and highlighting the emergence of intrinsic phonon angular momentum in acoustic branches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Effect of surface steps on chemical ordering in the subsurface of Cu(Au) solid solutions

Atomic steps are typically assumed to only influence surface phenomena of solids. In contrast, in this paper we show, using in-situ atomic-resolution electron microscopy observations and atomistic modeling, the pronounced effect from surface stepsin inducing compositional and structural evolution in the subsurface of a Cu(Au) solid solution. We find that Au surface segregation results in a stacked sequence of Cu-Au ordered phases in the subsurface. The presence of a monatomic step at the surface induces the formation of an anti-phase boundary that extends from the surface step to deeper layers by maintaining the same composition profile associated with each terrace. The bunching of surface steps induces chemical disordering of the Cu-Au ordered phases in the subsurface region of the bunched steps. These results demonstrate the instant propagation of surface dynamics of atomic steps into the subsurface region and can find broader applicability in utilizing surface defects to tune the composition and structure in the subsurface of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Phonon-Assisted Intertube Electronic Transport in an Armchair Carbon Nanotube Film

The electrical conductivity of a macroscopic assembly of nanomaterials is determined through a complex interplay of electronic transport within and between constituent nano-objects. Phonons play dual roles in this situation: their increased populations tend to reduce the conductivity via electron scattering, while they can boost the conductivity by assisting electrons to propagate through the potential-energy landscape. Here, we identified a phonon-assisted coherent electron transport process between neighboring nanotubes in temperature-dependent conductivity measurements on a macroscopic film of armchair single-wall carbon nanotubes. Through atomistic modeling of electronic states and calculations of both electronic and phonon-assisted junction conductances, we conclude that phonon-assisted conductance is the dominant mechanism for observed high-temperature transport in armchair carbon nanotubes. The unambiguous manifestation of coherent intertube dynamics proves a single-chirality armchair nanotube film to be a unique macroscopic solid-state ensemble of nano-objects promising for the development of room-temperature coherent electronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Shear strain alters the structure and migration mechanism of self-interstitial atoms in copper

Here, we use atomistic modeling to show that externally applied shear strain causes the lowest energy self-interstitial atom (SIA) structure in copper (Cu) to change from a <100>-type dumbbell to a <110>-type dumbbell. Concurrently, SIA migration switches from the 3-D random walk characteristic of <100>-type dumbbells to a 1-D mechanism analogous to that of crowdion SIAs. Furthermore, the relative energies of these two dumbbell structures as a function of strain are well predicted using elastic dipole tensors computed at zero strain, indicating that examination of these tensors may be used to assess the likelihood of strain-induced SIA structure transitions in other materials. Changes in lowest energy SIA structures and associated migration mechanisms stand to impact predictions of SIA behavior in irradiated solids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structure and migration of heavily irradiated grain boundaries and dislocations in Ni in the athermal limit

The microstructural evolution at and near preexisting grain boundaries (GBs) and dislocations in materials under high radiation doses is still poorly understood. In this work, we use the creation relaxation algorithm (CRA) developed for atomistic modeling of high-dose irradiation in bulk materials to probe the athermal limit of saturation of GB and dislocation core regions under irradiation in fcc Ni. We find that, upon continuously subjecting a single dislocation or GB to Frenkel pair creation in the athermal limit, a local steady-state disordered defect structure is reached with excess properties that fluctuate around constant values. Case studies are given for a straight screw dislocation which elongates into a helix under irradiation and several types of low- and high-angle GBs, which exhibit coupled responses such as absorption of extrinsic dislocations, roughening, and migration. A positive correlation is found between the initial GB energy and the local steady-state GB energy under irradiation across a wide variety of GB types. Metastable GB structures with similar density in the defect core region but different initial configurations are found to converge to the same limiting structure under CRA. The mechanical responses of pristine and irradiated dislocations and GB structures are compared under an applied shear stress. Irradiated screw and edge dislocations are found to exhibit a hardening response, migrating at larger flow stresses than their pristine counterparts. Mobile GBs are found to exhibit softening or hardening responses depending on GB character. Although some GBs recover their initial pristine structures upon migration outside of the radiation zone, many GBs sustain different flow stresses corresponding to altered mobile core structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Thermochemical Equilibria with Interacting Defects: Sr 1 − x Ce x Mn O 3 − δ Alloys for Water Splitting

Solar thermochemical hydrogen is one of the few potential routes towards direct fuel production from renewable energy sources, but the thermodynamic boundary conditions for efficient and economic energy conversion are challenging. Success or failure of a given oxide working material depends on the subtle balance between enthalpy and entropy contributions in the redox processes. Developing a mechanistic understanding of the behavior of materials on the basis of atomistic models and first-principles calculations is an important part of advancing the technology. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox-active defects start to interact with each other, thereby impeding the formation of additional defects. This problem is of more general importance to applications that rely on high levels of off-stoichiometry or doping, including, for example, batteries, thermoelectrics, and ceramic fuel cells. To account for such repulsive defect interactions, we introduce a statistical mechanics approach, defining an expression for the free energy of defect interaction based on limited sampling of defect configurations in density functional theory supercell calculations. The parameterization of this energy contribution as a function of defect concentration and temperature allows on-the-fly simulation of thermochemical equilibria. The approach consistently incorporates finite temperature effects by including the leading contributions to the temperature-dependent free energy for the case at hand, i.e., the ideal gas and configurational enthalpies and entropies. We demonstrate the capability and utility of the approach by simulating the water splitting redox processes for Sr 1 − x Ce x Mn O 3 − δ alloys. Published by the American Physical Society 2024

08 HYDROGEN↗
DOE OSTI2024

Toward a quantitative description of solvation structure: a framework for differential solution scattering measurements

Appreciating that the role of the solute–solvent and other outer-sphere interactions is essential for understanding chemistry and chemical dynamics in solution, experimental approaches are needed to address the structural consequences of these interactions, complementing condensed-matter simulations and coarse-grained theories. High-energy X-ray scattering (HEXS) combined with pair distribution function analysis presents the opportunity to probe these structures directly and to develop quantitative, atomistic models of molecular systems in situ in the solution phase. However, at concentrations relevant to solution-phase chemistry, the total scattering signal is dominated by the bulk solvent, prompting researchers to adopt a differential approach to eliminate this unwanted background. Though similar approaches are well established in quantitative structural studies of macromolecules in solution by small- and wide-angle X-ray scattering (SAXS/WAXS), analogous studies in the HEXS regime—where sub-ångström spatial resolution is achieved—remain underdeveloped, in part due to the lack of a rigorous theoretical description of the experiment. To address this, herein we develop a framework for differential solution scattering experiments conducted at high energies, which includes concepts of the solvent-excluded volume introduced to describe SAXS/WAXS data, as well as concepts from the time-resolved X-ray scattering community. Our theory is supported by numerical simulations and experiment and paves the way for establishing quantitative methods to determine the atomic structures of small molecules in solution with resolution approaching that of crystallography.

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

Effect of relative humidity on the interlayer spacing of phosphate intercalated Mg, Al layered double hydroxide (hydrotalcite‐like) crystals

Abstract The availability of water and the type of interlayer anions present affect the arrangement and interlayer spacing of layered double hydroxide (LDH). A systematic study of the effect of relative humidity (RH) on the basal spacing of phosphate intercalated, magnesium aluminum LDHs (MgAl–PO 4 LDH) was conducted, confirming that MgAl‐PO 4 LDH exhibit distinct interlayer spacing corresponding to separate low‐ and high‐hydration states produced when exposed to low and high humidity conditions. The transition from the low‐ to high‐hydration form begins when RH exceeds 33% and is accompanied by the transition from one to two interlayers of water. An improved understanding of the effect of RH on MgAl–PO 4 LDHs will enable more accurate atomistic modeling and experimental characterization of MgAl–PO 4 LDHs and may lead to improvements in the tunability of MgAl‐PO 4 LDHs for commercial applications such as slow‐release fertilizer.

Reed, Titus↗
DOE OSTI2022

Hole in one: Pathways to deterministic single-acceptor incorporation in Si(100)-2 x 1

We report that stochastic incorporation kinetics can be a limiting factor in the scalability of semiconductor fabrication technologies using atomic-precision techniques. While these technologies have recently been extended from donors to acceptors, the extent to which kinetics will impact single-acceptor incorporation has yet to be assessed. To identify the precursor molecule and dosing conditions that are promising for deterministic incorporation, we develop and apply an atomistic model for the single-acceptor incorporation rates of several recently demonstrated molecules: diborane (B 2 H 6 ), boron trichloride (BCl 3 ), and aluminum trichloride in both monomer (AlCl 3 ) and dimer forms (Al 2 Cl 6 ). While all three precursors can realize single-acceptor incorporation, we predict that diborane is unlikely to realize deterministic incorporation, boron trichloride can realize deterministic incorporation with modest heating (50 °C), and aluminum trichloride can realize deterministic incorporation at room temperature. We conclude that both boron and aluminum trichloride are promising precursors for atomic-precision single-acceptor applications, with the potential to enable the reliable production of large arrays of single-atom quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Atomistic mechanisms of water vapor–induced surface passivation

The microscopic mechanisms underpinning the spontaneous surface passivation of metals from ubiquitous water have remained largely elusive. Here, using in situ environmental electron microscopy to atomically monitor the reaction dynamics between aluminum surfaces and water vapor, we provide direct experimental evidence that the surface passivation results in a bilayer oxide film consisting of a crystalline-like Al(OH) 3 top layer and an inner layer of amorphous Al 2 O 3 . The Al(OH) 3 layer maintains a constant thickness of ~5.0 Å, while the inner Al 2 O 3 layer grows at the Al 2 O 3 /Al interface to a limiting thickness. On the basis of experimental data and atomistic modeling, we show the tunability of the dissociation pathways of H 2 O molecules with the Al, Al 2 O 3 , and Al(OH) 3 surface terminations. The fundamental insights may have practical significance for the design of materials and reactions for two seemingly disparate but fundamentally related disciplines of surface passivation and catalytic H 2 production from water.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identification of short-range ordering motifs in semiconductors

Chemical short-range ordering is expected to be a key factor for tuning the electronic structure of semiconductors. However, experimental evidence of short-range ordering is still lacking due to the challenge of characterizing atomic-scale ordering motifs. Here, we determined the presence of short-range order in a ternary GeSiSn semiconductor system using advanced energy-filtered four-dimensional scanning transmission electron microscopy and large-scale atomistic models generated by a machine learning neuroevolution potential of first-principles accuracy. This approach revealed preferred ordering of different atomic species with the dominant occurrence of Si–Ge–Sn triplets. Our findings not only confirmed the presence of short-range order but also directly revealed the actual atomic structure, demonstrating the potential for informed atomic order–based band engineering as a third degree of freedom beyond composition and strain tuning.

Vogl, Lilian M. [University of California Berkeley↗
DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
DOE OSTI2025

Impact of the LiPF6 Concentration on the Interfacial Charge Transfer and Fast-charging Capabilities of Lithium-Ion Batteries

Fast-charging lithium-ion batteries (LIB) demand optimized electrolyte formulations to balance ionic conductivity, viscosity, and interfacial charge transfer kinetics. This study examines how LiPF 6 concentration shapes solvation structure, desolvation energy, charge transfer activation energy, and solid electrolyte interphase (SEI) properties, which are critical for fast-charging performance. Using Raman spectroscopy, electrochemical cycling, X-ray photoelectron spectroscopy, and atomistic modeling, we analyze how varying LiPF 6 concentrations impact interfacial and bulk transport properties. Our findings show that increasing LiPF 6 concentration alters lithium solvation structures, reduces desolvation energy, and accelerates charge transfer at the electrode interface. Higher concentrations lower the activation energy for charge transfer and suppress excessive SEI growth, improving interfacial kinetics. However, concentrations above a certain threshold increase viscosity and reduce ionic conductivity, limiting transport efficiency. These results offer insights into electrolyte solvation and interfacial charge transfer mechanisms, providing guidelines for designing next-generation fast-charging LIB electrolytes with enhanced efficiency, stability, and longevity.

Son, Seoung-Bum [Argonne National Laboratory (ANL)↗
DOE OSTI2026

HydraGNN v5.0

HydraGNN v5.0 expands the code base into a more portable, scalable, and flexible framework for scientific graph learning, with particular strength in atomistic machine-learning interatomic potentials and large-scale distributed training. The release adds Fully Sharded Data Parallel (FSDP) support alongside existing DDP and DeepSpeed paths, including FSDP-aware checkpointing and optimizer integration, and introduces a configurable multi-precision training workflow supporting FP32, BF16, and FP64 across GPUs and Intel XPUs. For atomistic modeling, HydraGNN v5.0 strengthens its MLIP capabilities through dynamic graph construction at every forward pass, energy-conserving force prediction via automatic differentiation, and per-atom energy loss formulations, while extending EGNN models to properly handle periodic boundary conditions. The release also broadens model expressiveness through graph-level attribute conditioning, adds new multi-task and model-parallel extensions such as MACE support and encoder/decoder branch optimization, and expands application coverage with integrated examples for datasets including OC25, Nabla2-DFT, QCML, Open Polymers 2026, and OPF. In parallel, HydraGNN v5.0 improves production readiness through performance optimizations for large-scale runs, stratified sampling and linear-regression preprocessing utilities, and tested installation scripts for DOE supercomputers including Frontier, Aurora, Perlmutter, and Andes. Overall, the release advances HydraGNN as a robust software platform for scalable graph neural networks across materials science, chemistry, and scientific machine learning workflows

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗