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 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↗
DOE OSTI2026

HydraGNN_GFM_FineTuning4Materials v1.0

This repository enables fine-tuning of the HydraGNN Predictive GFM 2026 — an open-source ensemble of pre-trained graph foundation models for atomistic materials modeling, developed at Oak Ridge National Laboratory. The GFM 2026 is freely available and downloadable via Globus from the OLCF Data Constellation (DOI: 10.13139/OLCF/2562660). Starting from these pre-trained weights, this repository provides a complete transfer learning pipeline for adapting the GFM ensemble to domain-specific molecular and materials property prediction tasks. It includes: 1) Utilities for ensemble fine-tuning with task-specific output heads 2) Example pipelines for eight widely-used materials and molecular datasets 3) Tools for model adaptation and head configuration 4) Data preprocessing utilities for each supported dataset 5) Benchmarking and evaluation scripts

Ungerboeck, Linda↗
DOE OSTI2020

Functional Photoresists for Energy Applications. Final Report

Monolithic ultralow-density porous bulk materials have recently attracted much interest due to many emerging applications in the areas of catalysis, energy storage and conversion, and thermal insulation. They are also important components of high energy density (HED) and inertial confinement fusion (ICF) targets. However, despite tremendous progress that has been made in the synthesis of porous materials, deterministic and independent control over microscopic architecture, density and composition remain key issues, and their integration in high precision devices requires cost and time-intensive mechanical machining that not only reduces reproducibility by generating debris but also limits the complexity of the 3D shapes that can be realized. In this project, we overcame these limitations by developing a universal templating capability that provides deterministic and independent control over density, composition, architecture, and macroscopic sample shape. This was achieved by developing the technology to 1) 3D print ultrahigh resolution, ultra-high precision polymeric micro-lattice templates, 2) coat these templates with the desired materials, and 3) removing the template (Fig. 1a). Atomic layer deposition (ALD) provides the atomic scale coating thickness accuracy required for precisely controlling density. While this templating approach had been demonstrated in prior work, limitations in suitable photoresists, 3D print technologies, print design, and template removal techniques did not allow the fabrication of millimeter-sized high-precision parts with sub-micron resolution. To enable this technology, we developed 1) two-photon polymerization (TPP) print designs that enable the fabrication of millimeter-sized, mechanically robust polymeric templates with sub-micron resolution and 2) a continuum level TPP printing simulation capability for additional print design guidance; 3) atomistic models to study photoresist polymerization kinetics and network topography, 4) refractive index matched polymeric and preceramic TPP photoresists, and 5) functional TPP photoresists including porous voxel structures and self-immolative polymer photoresist chemistries; and 6) damage free template removal techniques that enable the fabrication of defect-free high-precision low-density foam components. We also developed a templating approach for pure carbon microlattices with a unique tube-in-tube ligament morphology. As a test platform, we pursued the fabrication of foam liners that promise to further increase the neutron yield in indirect drive ICF experiments by improving implosion symmetry control and coupling between the laser and the deuterium-tritium fuel. This application requires fabrication and integration of a ultra-high precision, millimeter-sized, thin-walled (200-400 micrometer thick), low-density (10-30 mg/cc), high atomic number (high Z) cylindrical foam tube into the gold hohlraum of an indirect drive ICF target (Fig. 1b). While our hohlraum liner test case will mainly find application in HED and ICF experiments, the underlying science will also directly apply to previously developed nanoparticle and additive manufacturing technologies and will advance those techniques as well.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theoretical Characterization of Correlated Redox Activity and Nanostructure in Hybrid Supercapacitors

Electrochemical energy storage devices such as supercapacitors, batteries, and fuel cells fundamentally rely on complex and correlated redox processes at electrode/electrolyte interfaces. Furthermore, the nanostructured morphology of high-surface area electrodes can significantly modulate this physics, and further complicates mechanistic understanding. Focusing on supercapacitor devices composed of carbon/polymer based electrodes and organic electrolytes, we will employ first-principles-based computer simulations to understand the coupling of these essential energy storage pathways occurring at the electrode/electrolyte interface. Of critical importance will be the development of theoretical methodologies to allow for the first-principles description of redox processes within electrochemical systems at fixed external voltage. Based on our hypothesis that redox activity and ion structuring is dramatically altered under nanoconfinement, we will then conduct a systematic study utilizing atomistic models of nanoporous carbon electrodes to understand the coupling of redox-activity, energy-storage, and nanostructuring for a variety of electrode morphologies.

25 ENERGY STORAGE↗
DOE OSTI2023

Nanosecond Laser Drying in Wafer Foundries for Significant Energy Efficiency Improvement

Surface drying is a critical process in microelectronics wafer fabrication. In order to reduce and eliminate watermarks, many surface drying processes have been studied and used. Most of the wafer foundries currently use drying processes based on Isopropyl Alcohol (IPA), which causes serious safety (fire), health, environmental, and energy efficiency issues. The long-term goal of this project is to develop a fast, effective, chemical-free, and extremely energy efficient wafer drying technology that employs nanosecond laser-induced sub-surface evaporation and explosion. The objectives of this project are to 1) experimentally develop and verify the laser surface drying method, explore and identify how and to what extent various physical parameters influence the drying efficiency and effectiveness; and 2) understand and optimize the underlying thermodynamics in the proposed laser surface drying, including micro/nanoscale sub-surface superheating, evaporation, explosion, and photo/thermochemical bond breaking via atomistic modeling and experimental characterization. Toward these objectives, five tasks have been accomplished successfully, including 1) develop in-lab laser surface drying setup and study the effect of varied laser parameters, 2) use multiscale hybrid-modeling to study and understand the physics of water behavior during drying toward process optimization, 3) investigate the effect of drying conditions on the drying effectiveness, 4) develop a large-scale laser drying scanning system with automatic scanning, and 5) conduct large-scale simulation to study the effect of dopant level and develop strategies for system development toward industry deployment.

42 ENGINEERING↗
DOE OSTI2025

High temperature nuclear data measurements of SiC, ZrC, and MgO [Slides]

Performed temperature dependent measurements of SiC and ZrC at ARCS instrument at SNS. Performed temperature dependent measurements of SiC, ZrC, and MgO at VISION instrument at SNS. Performed initial atomistic modeling of these materials using various techniques, including machine learned potentials. Future work includes temperature dependent transmission measurements of these materials, as well as improving the machine learned potentials with more training data and different machine learned frameworks.

ARCS↗