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

Harness the power of atomistic modeling and deep learning in biofuel separation

Biofuels offer a remarkable, sustainable energy source for a future of clean energy. The development of efficient biofuel separation plays a crucial role in achieving cost-effective utilization of biofuel. In this chapter, we provide an overview of the recent advancements in atomistic-level modeling and deep learning in the rational design of novel, efficient biofuel separation. The fundamental principles of quantum and statistical mechanics are covered in appropriate detail to highlight their underlying differences in theory. The methodologies of several molecular representations and deep learning algorithms applicable to biofuel separation are briefly demonstrated as well. The applications, successes, and risks of employing density functional theory, ab initio molecular dynamics, classical molecular dynamics, and deep learning are provided to showcase their recent accomplishments in biofuel separation as well as potential improvements in both methodology and application. Lastly, a vision for the future growth of these methods is illustrated.

deep learning, artificial intelligence, biofuels, ↗
DOE OSTI2024

Molecular Understanding of Nitrogen Oxide Fixation of Water-Lean Carbon Capture Solvents by Atomistic Modeling

Nitrogen oxides, present in flue gas, can cause negative impacts on amine carbon capture solvents by the formation of heat-stable salts and suspected carcinogens. Thus, to maximize the performance of water-lean solvents, a better understanding of this process in these systems is necessary. Here, a computational study for the fixation of the CO 2 capture solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (EEMPA) to nitramine/nitrosamine was conducted. The first step involves the dissociation of the NH bond of EEMPA, in which the homolytic mechanism is energetically more favorable than the heterolytic mechanism. The second step involves radical recombination to form N–N bonds. While NO 2 directly reacts with EEMPA, NO has almost no effect. However, in the presence of O 2 , fixation of EEMPA by NO is enhanced via the formation of N 2 O 4 species. Finally, low reaction energies indicate that the formation of nitramine/nitrosamine may be a reversible process, suggesting that EEMPA could be recovered under thermal stripping conditions.

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

Atomistic modeling of tensile deformation and fracture of carbon fibers: Nanoscale stress redistribution, effect of local structural characteristics and nanovoids

The mechanisms of tensile deformation and fracture of carbon fibers are investigated in large-scale molecular dynamics simulations. The simulations performed for computational samples with different structural characteristics enable analysis of the structural sensitivity of the mechanical properties. It is found that the carbon fiber nanostructure can be represented by a fine (nanoscale) mixture of stiff and soft regions. Despite the nanoscale structural heterogeneity of the computational fibers, their fracture strain and tensile strength exceed the experimental values by almost an order of magnitude. Here, this observation is explained by the ability of the carbon fiber nanostructure, largely consisting of curved and folded layers of turbostratic carbon, to exhibit continuous nanoscale load redistribution during the deformation. The transient, fluid nature of a percolating cluster of stress-holding structural elements prevents stress localization at the early stage of the deformation and enables the fiber nanostructure to adapt to the increasing load. Overall, the results of the simulations suggest that the design of advanced carbon fiber manufacturing and post-processing methods focused on the elimination of critical mesoscopic structural defects (internal voids, foreign inclusions, and surface flaws) can yield a substantial increase in the strength and fracture strain of carbon fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An atomistic model for the thermal resistance of a liquid–solid interface

The thermal resistance associated with the interface between a solid and a liquid is analysed from an atomistic point of view. Partial evaluation of the associated Green–Kubo integral elucidates the various factors governing heat transport across the interface and leads to a quantitative model for the thermal resistance in terms of atomistic-level system parameters. The model is validated using molecular dynamics simulations.

Mechanics↗
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scalable training of trustworthy and energy-efficient predictive graph foundation models for atomistic materials modeling: a case study with HydraGNN

We present our work on developing and training scalable, trustworthy, and energy-efficient predictive graph foundation models (GFMs) using HydraGNN, a multi-headed graph convolutional neural network architecture. HydraGNN expands the boundaries of graph neural network (GNN) computations in both training scale and data diversity. It abstracts over message passing algorithms, allowing both reproduction of and comparison across algorithmic innovations that define nearest-neighbor convolution in GNNs. This work discusses a series of optimizations that have allowed scaling up the GFMs training to tens of thousands of GPUs on datasets consisting of hundreds of millions of graphs. Our GFMs use multitask learning (MTL) to simultaneously learn graph-level and node-level properties of atomistic structures, such as energy and atomic forces. Using over 154 million atomistic structures for training, we illustrate the performance of our approach along with the lessons learned on two state-of-the-art US Department of Energy (US-DOE) supercomputers, namely the Perlmutter petascale system at the National Energy Research Scientific Computing Center and the Frontier exascale system at Oak Ridge Leadership Computing Facility. The HydraGNN architecture enables the GFM to achieve near-linear strong scaling performance using more than 2000 GPUs on Perlmutter and 16,000 GPUs on Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Platinum nanoparticle compression: Combining in situ TEM and atomistic modeling

The mechanical behavior of nanoparticles governs their performance and stability in many applications. However, the small sizes of technologically relevant nanoparticles, with diameters in the range of 10 nm or less, significantly complicate experimental examination. These small nanoparticles are difficult to manipulate onto commercial test platforms and deform at loads that are below the typical noise floor of the testing instruments. Here, we synthesized small platinum nanoparticles directly onto a mechanical tester and used a modified nanomanipulator to enhance load resolution to the nanonewton scale. We demonstrated the in situ compression of an 11.5-nm platinum nanoparticle with simultaneous high-resolution measurements of load and particle morphology. Molecular dynamics simulations were performed on similarly sized particles to achieve complementary measurements of load and morphology, along with atomic resolution of dislocations. The experimental and simulation results revealed comparable values for the critical resolved shear stress for failure, 1.28 and 1.15 GPa, respectively. Altogether, this investigation demonstrated the promise of, and some initial results from, the combination of atomistic simulations and in situ experiments with an unprecedented combination of high spatial resolution and high load resolution to understand the behavior of metal nanoparticles under compression.

42 ENGINEERING↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Advances in atomistic modeling and understanding of drying shrinkage in cementitious materials

Despite more than a century of research, drying shrinkage still causes unwanted stresses and cracks in our concrete infrastructures. This is partly due to the lack of understanding of the underlying mechanisms governing the drying shrinkage in cementitious materials. Inspired by the tremendous impact of atomistic simulations on materials science, we expect such innovative simulation techniques to hold a key to unlocking the enigmatic nature of drying shrinkage at the nanoscale. In this treatise, we first introduce basic physics concepts and atomistic simulation methods in detail. Afterward, we discuss the progress brought by atomistic simulations in understanding the drying shrinkage in cementitious materials and other nanoporous and nanolayered systems. This review also highlights important remaining fundamental questions and practical issues regarding drying shrinkage and delineates how atomistic simulations can help resolve them.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of surface and grain boundary dislocation nucleation in FCC metals

Dislocation nucleation plays a critical role in the plastic deformation of crystalline materials. However, it is challenging to predict the active mode and associated rate of dislocation nucleation under typical experimental loading conditions through molecular dynamics simulation due to timescale limitations. In this report we use the free-end nudged elastic band method to determine the activation energies and activation volumes of dislocation nucleation in four typical face-centered cubic metals of Au, Al, Cu and Ni. We focus on the representative processes of surface and grain boundary dislocation nucleation. The atomistically determined activation volumes of these dislocation nucleation processes are larger than 10b3 (with b being the Burgers vector length) under typical experimental loading conditions. These results are compared with experimentally measured activation volumes in ultrafine-grained and nanocrystalline metals, thereby providing mechanistic insight into their rate-controlling deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel

Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding (La,Sr)(Co,Fe)O$_{3-δ}$ Phase Instability within SOECs Using a Combined Experimental and Atomistic Modeling Approach

Understanding the onset of degradation in the air electrode within solid oxide electrolysis cells (SOECs), and the subsequent impact on cell performance, is a critical step in mitigating the performance losses and stability issues of SOECs. In an effort to identify early onset degradation phenomena, SOECs were characterized as fabricated and after testing potentiostatically at 1.3 V for 1000 h at 750 °C. SOEC air electrodes composed of a 1:1 composite of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (6428-LSCF) and Gd 0.1 Ce 0.9 O 1.95 (GDC) were studied using synchrotron X-ray diffraction (XRD), scanning transmission electron microscopy coupled with energy dispersive X-ray spectroscopy (STEM-EDS), and X-ray absorption near-edge spectroscopy (XANES) to evaluate the changes in the air electrode structurally and chemically. These techniques show the migration of Sr species from the air electrode through pores in the GDC barrier layer, progressing to the electrolyte boundary, where it accumulates and reacts with (Zr 0.84 Y 0.16 )O 2–δ (YSZ) to form SrZrO 3 . Microscopy results are paired with atomistic simulations to better understand the relationship between the thermodynamic instability of 6428-LSCF and cell fabrication/testing conditions. First-principles calculations reveal that LSCF-6428 is not stable during cell manufacturing and testing conditions, which supports the experimental identification of secondary phases in both as-fabricated and tested cells. Together, these results demonstrate that the challenging environments encountered by SOECs during cell manufacturing and operation lead to instabilities of the target 6428-LSCF anode material and underscore the need for more durable, high-performing SOEC components.

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

Atomistic modeling of interface strengthening in Al-Si eutectic alloys

Al-Si cast alloys are usually composed of α-Al and Al-Si eutectic. Si flakes and Al matrix generally hold cube-on-cube orientation relationship with the primary interface (111) Al ∥(111) Si . Extensive experimental studies demonstrated that Si flakes cannot significantly improve mechanical properties of Al-Si cast alloys. We hypothesize that the weak strengthening effect associated with Si flakes might be attributed to thermomechanical properties of Al-Si interfaces besides their morphologies. To characterize Al-Si interfaces with a large lattice mismatch (> 30%), we proposed the quasi-coincident site lattice (Q-CSL) as reference lattice, and demonstrated that the Q-CSL Al-Si coherent interface has three characteristic coherent structures, one stable and low energy structure and two metastable and high energy structures. The translation vectors for the same type of coherent Q-CSL structures are consistent with three displacement shift complete (DSC) vectors. The two metastable structures can be obtained by shifting the low energy structure with three partial DSC vectors. Semi-coherent interface is composed of the low energy Q-CSL patches and three sets of interface misfit dislocations with Burgers vectors same as the DSC vectors. Atomistic simulations revealed that Al-Si interface exhibits low shear resistance. Ideal shear strength of the Q-CSL coherent interface is 110 MPa and semi-coherent interface is 20 MPa. The low shear resistance is attributed to the glide of interface misfit dislocations. Al-Si interface also exhibits low formation and migration energies of point defects. Owing to low shear strength and low formation and migration energies of point defects, interface sliding or shear readily happen under mechanical loading or during dislocation-interface interactions. Lattice dislocations can cross slip onto or climb along Al-Si interfaces. These reactions decrease the number of accumulated dislocation loops around Si flakes and promote nucleation and emission of lattice dislocations from Al-Si interfaces to matrix, consequently reduce the repulsive force on approaching dislocations and weaken Si flakes strengthening effect. In situ tension and compression tests in a scanning electron microscope reveal relatively weak strengthening effect due to Si flakes, consistent with the computed dislocation interaction with interfaces and shear behavior of interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗