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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Te/CdTe and Al/CdTe Interfacial Energy Band Alignment by Atomistic Modeling

A synergistic approach that incorporates first-principles atomistic modeling with numerical device simulations is used to systematically evaluate the role of heterointerfaces within metal-chalcogenide-based photovoltaic technologies. Two interfaces involving either a tellurium back contact or aluminum back electrode combined with a cadmium telluride absorber layer within cadmium-telluride-based thin-film solar cells are investigated on an atomic scale to determine the mechanisms contributing to variations in device performance. Electronic structures and predicted charge transport behavior with respect to cadmium and tellurium termination of the absorber layer are studied along the polar oriented CdTe{111} facets. The computational methodology reveals a noticeable contrast between the Schottky barrier forming Al/CdTe interface versus the Type I Te/CdTe heterojunction. Greater band bending features are exhibited by the cadmium termination as opposed to the tellurium termination for each interface case. Subsequent device modeling suggests that 3.6% higher photovoltaic conversion efficiency is achievable for the cadmium termination relative to the tellurium termination of the Te/CdTe interface. Based strictly on an idealistic representation, both interface models show the importance of atomic-scale interfacial properties for cadmium telluride solar cell device performance with their bulk properties being validated in comparison to published experimental data. Here, the synergistic approach offers a suitable method to analyze solar cell interfaces through a predictive computational framework for the engineering and optimization of metal-chalcogenide-based thin-film photovoltaic technologies.

Green’s function↗
DOE OSTI2021

Atomistic modeling of meso -timescale processes with $\mathrm{SEAKMC}$: A perspective and recent developments

On-the-fly kinetic Monte Carlo (kMC) methods have recently garnered significant attentions after successful applications to various atomic-scale problems using a timescale outside the reach of classical molecular dynamics. These methods play a critical role in modeling atomistic meso-timescale processes, and it is therefore essential to further improve their capabilities. In this report we review one of the on-the-fly kMC methods, Self-Evolving Atomistic kinetic Monte Carlo (SEAKMC) and propose two schemes that considerably enhance the efficiency of saddle point searches (SPSs) during the simulations. The performance of these schemes is tested using the diffusion of point defects in bcc Fe. In addition, we discuss approaches to significantly mitigate limitations of these schemes, which further improves their efficiencies. Importantly, these schemes improve the SPS efficiency not only for SEAKMC but also for other on-the-fly kMC methods, broadening the applications of on-the-fly kMC simulations to complex meso-timescale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scaling Laws of Graph Neural Networks for Atomistic Materials Modeling

Atomistic materials modeling is a critical task with wide-ranging applications, from drug discovery to materials science, where accurate predictions of the target material property can lead to significant advancements in scientific discovery. Graph Neural Networks (GNNs) represent the state-of-the-art approach for modeling atomistic material data thanks to their capacity to capture complex relational structures. While machine learning performance has historically improved with larger models and datasets, GNNs for atomistic materials modeling remain relatively small compared to large language models (LLMs), which leverage billions of parameters and terabyte-scale datasets to achieve remarkable performance in their respective domains. To address this gap, we explore the scaling limits of GNNs for atomistic materials modeling by developing a foundational model with billions of parameters, trained on extensive datasets in terabytescale. Our approach incorporates techniques from LLM libraries to efficiently manage large-scale data and models, enabling both effective training and deployment of these large-scale GNN models. This work addresses three fundamental questions in scaling GNNs: the potential for scaling GNN model architectures, the effect of dataset size on model accuracy, and the applicability of LLM-inspired techniques to GNN architectures. Specifically, the outcomes of this study include (1) insights into the scaling laws for GNNs, highlighting the relationship between model size, dataset volume, and accuracy, (2) a foundational GNN model optimized for atomistic materials modeling, and (3) a GNN codebase enhanced with advanced LLM-based training techniques. Our findings lay the groundwork for large-scale GNNs with billions of parameters and terabyte-scale datasets, establishing a scalable pathway for future advancements in atomistic materials modeling.

Li, Chaojian [ORNL] (ORCID:0000000340309777)↗
DOE OSTI2026

Multi-task Parallelism for Robust Pre-training of Graph Foundation Models on Multi-source, Multi-fidelity Atomistic Modeling Data

Graph foundation models using graph neural networks promise sustainable, efficient atomistic modeling. To tackle challenges of processing multi-source, multi-fidelity data during pre-training, recent studies employ multi-task learning, in which shared message passing layers initially process input atomistic structures regardless of source, then route them to multiple decoding heads that predict data-specific outputs. This approach stabilizes pre-training and enhances a model’s transferability to unexplored chemical regions. Preliminary results on approximately four million structures are encouraging, yet questions remain about generalizability to larger, more diverse datasets and scalability on supercomputers. We propose a multi-task parallelism method that distributes each head across computing resources with GPU acceleration. Implemented in the open-source HydraGNN architecture, our method was trained on over 24 million structures from five datasets and tested on the Perlmutter, Aurora, and Frontier supercomputers, demonstrating efficient scaling on all three highly heterogeneous super-computing architectures.

Lupo Pasini, Massimiliano [ORNL] (ORCID:0000000249↗
DOE OSTI2020

Atomistic modeling of energy band alignment in CdSeTe surfaces

Atomistic modeling based on Desnity Functional Theory (DFT)- 1/2 method coupled with surface Green’s function has been employed to investigate the energy band alignment results in cadmium selenium telluride (CdSeTe) surfaces. The structural and electronic properties of the bulk ternary alloy were established before cleaving the CdSeTe low index surface facets. The dependency of atomic-scale electronic features on different plane orientations was explored by studying the energy band alignment diagrams in the unreconstructed and reconstructed surfaces. While the low index CdSeTe surface geometry reconstructions replicate the CdTe surface geometries, the energy band alignment features of the unreconstructed and reconstructed CdSeTe surfaces differ from those observed in CdTe low index surfaces. The structurally relaxed reconstructed CdSeTe surface results in purely flat bands as opposed to an internal cusp feature observed in CdTe surfaces. The presence of an internal energy cusp in CdTe(111) surfaces is expected to play a key role in enhancing the hole charge transport towards the back of CdTe-only solar cell device and the absence of such feature in CdSeTe surface may explain one of the reason for the lower performance of CdSeTe-only solar cell seen experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Persistence of negative vacancy and self-interstitial formation energies in atomistic models of amorphous silicon

We use classical potential atomistic modeling to demonstrate that negative vacancy and selfinterstitial formation energies are inherently present in amorphous silicon (a-Si). Our conclusions are based on comprehensive calculations of vacancy and self-interstitial formation energies in nine different models, ranging in size from 1000 to 2744 atoms and formed by quenching liquid Si at rates between 10 10 and 10 12 K/s. We account for a range of possible reference energies corresponding to different choices of reservoirs for atoms inserted into or removed from our models. Moreover, we show that successive addition of either vacancies or self-interstitials does not deplete the number of negative formation energy sites. Instead, it further lowers defect formation energies. Furthermore, these findings lead us to conclude that a-Si is not a structure at or near a local enthalpy minimum nor can it be brought to a local enthalpy minimum through the introduction of constitutional defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mapping TpPa-1 covalent organic framework (COF) molecular interactions in mixed solvents via atomistic modeling and experimental study

Complex solvent environments continue to limit the widespread adoption of organic solvent nanofiltration (OSN) in many chemical industry applications. In this paper we employ a commercially available covalent organic framework (COF), TpPa-1, and force field models to molecularly map separation performance of TpPa-1 membrane in mixed solvents. To minimize time and length scale mismatch between atomistic modeling and experiments, solvent permeance was normalized with water in modeling and experimental results to enable direct comparison. Model outputs, such as organic solvent permeance and solute rejection rate, matched well with filtration results. Since the atomistic models assume that all mass transfer is via through-pore transport, the discrepancies between modeling and experimental results provide insights on the effect of linear polymer defects, adsorption and interstitial mass transfer on polycrystalline COF membrane performance. In sum, force field models can serve as digital twins of COF membranes to simulate separation processes while capturing the effects of COF structure, chemistry, and crystallinity on membrane performance in complex organic solvent environments. Finally, this approach will provide insight into future COF design and synthesis for persisting separation challenges.

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

Atomistic modeling of radiation damage in crystalline materials

This review discusses atomistic modeling techniques used to simulate radiation damage in crystalline materials. Radiation damage due to energetic particles results in the formation of defects. The subsequent evolution of these defects over multiple length and time scales requiring numerous simulations techniques to model the gamut of behaviors. This work focuses attention on current and new methodologies at the atomistic scale regarding the mechanisms of defect formation at the primary damage state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Large-scale atomistic model construction of subbituminous and bituminous coals for solvent extraction simulations with reactive molecular dynamics

Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.

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

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic modeling of energy band alignment in CdTe(1 0 0) and CdTe(1 1 1) surfaces

An atomic-scale perspective of energy band alignment in CdTe surfaces has not been spatially studied despite the major role surfaces play in forming interfaces within CdTe-based thin film photovoltaic devices. Atomistic modeling based on density functional theory coupled with surface Green’s function is used for calculating energy band alignment of CdTe surfaces. The CdTe(1 0 0) ((1×1) and c(2×2) reconstruction) and CdTe(1 1 1) ((1×1) and (2×2) reconstruction) facets without and with surface relaxation provide insightful band bending characteristics that influence charge carrier transport. Results show that unrelaxed (1×1) CdTe(1 0 0) and CdTe(1 1 1) surfaces bend the valence band downward with surface polarity dictating the surface potential magnitude. The reconstructed CdTe(1 0 0) c(2×2) and CdTe(1 1 1) (2×2) surfaces result in favorable surface electronic features in relation to their unreconstructed variants. In addition, the structurally relaxed CdTe(1 1 1) surfaces develop an internal energy cusp potential that may enhance hole charge transport toward the back of CdTe solar cell devices. In conclusion, energy band alignments calculated within the study lead to a detailed understanding of how CdTe surfaces may affect CdTe-based thin film photovoltaic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic modeling of metal–nonmetal interphase boundary diffusion

Atomistic computer simulations are applied to investigate the atomic structure, thermal stability, and diffusion processes in Al–Si interphase boundaries as a prototype of metal–ceramic interfaces in composite materials. Some of the most stable orientation relationships between the phases found in this work were previously observed in epitaxy experiments. Here, a non-equilibrium interface can transform to a more stable state by a mechanism that we call interface-induced recrystallization. Diffusion of both Al and Si atoms in stable Al–Si interfaces is surprisingly slow compared with diffusion of both elements in Al grain boundaries but can be accelerated in the presence of interface disconnections. A qualitative explanation of the sluggish interphase boundary diffusion is proposed. Atomic mechanisms of interphase boundary diffusion are similar to those in metallic grain boundaries and are dominated by correlated atomic rearrangements in the form of strings and rings of collectively moving atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DDStore: Distributed Data Store for Scalable Training of Graph Neural Networks on Large Atomistic Modeling Datasets

Graph neural networks (GNNs) are a class of Deep Learning models used in designing atomistic materials for effective screening of large chemical spaces. To ensure robust prediction, GNN models must be trained on large volumes of atomistic data on leadership class supercomputers. Even with the advent of modern architectures that consist of multiple storage layers that include node-local NVMe devices in addition to device memory for caching large datasets, extreme-scale model training faces I/O challenges at scale.We present DDStore, an in-memory distributed data store designed for GNN training on large-scale graph data. DDStore provides a hierarchical, distributed, data caching technique that combines data chunking, replication, low-latency random access, and high throughput communication. DDStore achieves near-linear scaling for training a GNN model using up to 1000 GPUs on the Summit and Perlmutter supercomputers, and reaches up to a 6.15x reduction in GNN training time compared to state-of-the-art methodologies.

Choi, Jong Youl↗
DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2023

Atomistic modeling of a superconductor–transition metal dichalcogenide–superconductor Josephson junction

Using an atomistic tight-binding model, we investigate the characteristics of a Josephson junction formed by monolayers of MoS 2 sandwiched between Pb superconducting electrodes. We derive and apply Green's function–based formulation to compute the Josephson current as well as the local density of states in the junction. Our analysis of diagonal and off-diagonal components of the local density of states reveals the presence of triplet superconducting correlations in the MoS 2 monolayers and spin-polarized subgap (Andreev bound) states. Our formulation can be extended to other systems where atomistic details and large scales are needed to obtain accurate modeling of Josephson junction physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Emerging Atomistic Modeling Methods for Heterogeneous Electrocatalysis

Heterogeneous electrocatalysis lies at the center of various technologies that could help enable a sustainable future. However, its complexity makes it challenging to accurately and efficiently model at an atomic level. Herein, we review emerging atomistic methods to simulate the electrocatalytic interface with special attention devoted to the components/effects that have been challenging to model, such as solvation, electrolyte ions, electrode potential, reaction kinetics, and pH. Additionally, we review relevant computational spectroscopy methods. Then, we showcase several examples of applying these methods to understand and design catalysts relevant to green hydrogen. We also offer experimental views on how to bridge the gap between theory and experiments. Finally, we provide some perspectives on opportunities to advance the field.

36 MATERIALS SCIENCE↗