Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive molecular dynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Precipitation of calcium–alumino–silicate–hydrate gels: The role of the internal stress

Concrete gains its strength from the precipitation of a calcium–alumino–silicate–hydrate (C–A–S–H) colloidal gel, which acts as its binding phase. However, despite concrete’s ubiquity in the building environment, the atomic-scale mechanism of C–A–S–H precipitation is still unclear. Here, we use reactive molecular dynamics simulations to model the early-age precipitation of a C–A–S–H gel. We find that, upon gelation, silicate and aluminate precursors condensate and polymerize to form an aluminosilicate gel network. Notably, we demonstrate that the gelation reaction is driven by the existence of a mismatch of atomic-level internal stress between Si and Al polytopes, which are initially experiencing some local tension and compression, respectively. The polymerization of Si and Al polytopes enables the release of these competitive stresses.

Zhao, Cheng↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
NASA NTRS2010

Modeling Initial Stage of Ablation Material Pyrolysis: Graphitic Precursor Formation and Interfacial Effects

Reactive molecular dynamics simulations are used to study initial stage of pyrolysis of ablation materials and their composites with carbon nanotubes and carbon fibers. The products formed during pyrolysis are characterized and water is found as the primary product in all cases. The water formation mechanisms are analyzed and the value of the activation energy for water formation is estimated. A detailed study on graphitic precursor formation reveals the presence of two temperature zones. In the lower temperature zone (less than 2000 K) polymerization occurs resulting in formation of large, stable graphitic precursors, and in the high temperature zone (greater than 2000 K) polymer scission results in formation of short polymer chains/molecules. Simulations performed in the high temperature zone on the phenolic resin composites (with carbon nanotubes and carbon fibers) shows that the presence of interfaces had no substantial effect on the chain scission rate or the activation energy value for water formation.

Desai, Tapan G.↗
NASA NTRS2019

Toward Ultralight High Strength Structural Materials via Collapsed Carbon Nanotube Bonding

The growing commercial availability of carbon nanotube (CNT) macro-assemblies such as sheet and yarn is making their use in structural composite components increasingly feasible. However, the mechanical properties of these materials continue to trail those of state-of-the-art carbon fiber composites due to relatively weak inter-tube load transfer. Forming covalent links between adjacent CNTs promises to mitigate this problem, but it has proven difficult in practice to introduce them chemically within densified and aligned CNT materials due to their low permeability. To avoid this limitation, this work explores the combination of pulsed electrical current, temperature, and pressure to introduce inter-CNT bonds. Reactive molecular dynamics simulations identify the most probable locations, configurations, and conditions for inter-nanotube bonds to form. This process is shown to introduce covalent linkages within the CNT material that manifest as improved macroscale mechanical properties. The magnitude of this effect increases with increasing levels of prealignment of the CNT material, promising a new synthesis pathway to ultralight structural materials with specific strengths and stiffnesses exceeding 1 and 100 GPa/(g/cu.cm), respectively.

Jensen, Benjamin D.↗
DOE OSTI2020

Lithium Metal Electrode — Understanding Its Unique Characteristics and Functions

Lithium metal electrode (LME) is attractive for the next generation of commercial high-energy rechargeable lithium battery (RLB) design and development. As a light element with a very reducing electrochemical potential, Li metal certainly has the unique characteristics that also present challenges for its application in RLBs. To overcome technical barriers with LME such as the dendrite formation that causes uneven current distributions, excessive SEI formation that consumes Li inventory, and formation of inactive Li metal in granular forms that induces safety hazards, more fundamental understanding of the unique characteristics of LME and its impact on functions in performance in RLBs is critical. Here, we present some recent work based on temporally-correlated reactive molecular dynamics simulations to understand the Li deposition and the associated nucleation and growth process on the LME or graphitic electrode surfaces. Such a study provides some intriguing perspectives on how Li reacts with electrolytes on LME and conventional graphitic electrodes in the RLB operation. Such an insight could help improving our strategy to deal with LME problems and provide functional improvements for LME’s application in RLBs. This presentation is dedicated to Professor Robert Huggins for his life-long teaching of the fundamentals on solid state ionics and its implications in understanding the electrochemical behavior of Li in RLB applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Self-Healing Ceramics and Mechanical Metamaterials

We performed large scale molecular dynamics (MD) simulations and reactive & quantum molecular dynamics simulations to investigate: (1) faceting, grain growth, and crack healing in ceramics; (2) stress-release mechanism and healing in fibers; and (3) mechanical properties of two-dimensional metamaterials. We also designed highly accurate interatomic potentials and forces using artificial neural networks (ANN) and first-principles molecular dynamics simulation data for use in our MD and reactive MD simulations. We demonstrated the efficacy of modeling quantum materials atomistic simulation data using unsupervised deep learning approaches on a quantum computer. These achievements are elaborated in the following sections.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2023

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An examination of the performance of molecular dynamics force fields: Silicon and silicon dioxide reactive ion etching

Classical molecular dynamics (MD) simulations of plasma-surface interactions were performed of physical sputtering and reactive ion etching (RIE), with predictions based on several force fields. In this paper, we focus mainly on SiO 2 but include some results for Si substrates as well. Here, we compare predictions from these MD simulations to experimental studies of SiO 2 physical sputtering (by Ar + ions), RIE of Si, and RIE of SiO 2 (both using F atoms and Ar + ions). MD results using different published force fields are compared to reported yields from published vacuum beam experiments. The near-surface depth profiles predicted using different force fields are compared. One motivation for the present study is to document the nature and magnitude of differences in the predictions for selected systems and conditions of practical interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Lithiation of Sulfur-Graphene Compounds Using Reactive Force-Field Molecular Dynamics Simulations

We performed molecular dynamics simulations of lithium-sulfur-graphene compounds using reactive force fields, providing a time scale to observe atomistic features relevant to the microscopic behavior of the of the bulk of sulfur-based cathodes to be used beyond our present Li-ion batteries. The samples we used were set to realistic geometries through sophisticated protocols to simulate ultrafast reactions that occur within the picosecond range, thus allowing us to get some insights into the characteristics of the bulk material in working cathodes of Li–S batteries, which are mixed with carbon to increase the poor electronic conductivity of S. We report chemical speciation and geometrical data at atomistic levels. We observed that slowly lithiated cathodes were more stable and with higher density than those that were suddenly fully-lithiated. We did not observe molecular Li 2 S formation; however, we observed an amorphous solid arrangement with the same stoichiometry of Li and S, with S–Li–S angles of ~111° and smaller ones due to the interaction between polysulfides that did not reacted totally. In addition, graphene keeps its planar shape; however, S8 changes its shape from rings to chains. Lithiated structures are more stable with lower energies, and more close-packed structures than structures with Li already inserted.

25 ENERGY STORAGE↗
DOE OSTI2023

Isolating Chemical Reaction Mechanism as a Variable with Reactive Coarse-Grained Molecular Dynamics: Step-Growth versus Chain-Growth Polymerization

Here we present a general approach to isolate chemical reaction mechanism as an independently controllable variable across chemically distinct systems. Modern approaches to reduce the computational expense of molecular dynamics simulations often group multiple atoms into a single “coarse-grained” interaction site, which leads to a loss of chemical resolution. In this work we convert this shortcoming into a feature and use identical coarse-grained models to represent molecules that share nonreactive characteristics but react by different mechanisms. As a proof of concept, we use this approach to simulate and investigate distinct, yet similar, trifunctional isocyanurate resin formulations that polymerize by either chain- or step-growth. Because the underlying molecular mechanics of these models are identical, all emergent differences are a function of the reaction mechanism only. We find that the microscopic morphologies resemble related all-atom simulations and that simulated mechanical testing reasonably agrees with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗