Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics simulation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Elucidating Lithium Transport Mechanisms in Disordered LiF from Machine-Learning Molecular Dynamics Simulations

Lithium fluoride (LiF) is a ubiquitous component of solid- and cathode–electrolyte interphases, yet its functional role remains unclear under the structural and chemical heterogeneity typical of cycling batteries. Here, we systematically quantify how structural disorder, off-stoichiometry, and strain govern Li-ion transport in LiF. Using a machine-learning potential to enable extensive molecular-dynamics sampling, we compare crystalline and amorphous LiF, Li 0.95 F, and LiF 0.95 , and evaluate the impact of small homogeneous deformations. Defect-free crystalline LiF is effectively ion-blocking at 300−500 K, whereas amorphization generates free-volume–assisted percolation pathways that facilitates Li-ion diffusion. At elevated temperatures, thermodynamically driven crystallization disrupts these pathways, leading to non-Arrhenius behavior. In crystalline phases, Li deficiency activates vacancy-mediated diffusion, while in amorphous LiF, transport is governed primarily by network connectivity. Strain is found to have only a marginal effect on Li mobility in both crystalline and amorphous structures.

Batteries↗
DOE OSTI2024

Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations

Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.

42 ENGINEERING↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2021

All-Atom Molecular Dynamics Simulations of the Temperature Response of Densely Grafted Polyelectrolyte Brushes

Water acts a “good” solvent for most polyelectrolyte (PE) systems due to its high dielectric strength. Theoretical and experimental studies have highlighted some remarkably unconventional properties of strongly water-swollen PE brushes such as a reduction in brush height with increase in temperature. However, an understanding of the atomistic-scale response of the grafted PE chains as well as the local arrangement and organization of the brush-supported counterions and water molecules to changes in temperature is missing. Here, we conduct an all-atom molecular dynamics (MD) study to probe the influence of temperature on the “microstructure” of densely grafted, fully ionized polyacrylic acid (PAA) brushes. The atomistic insights obtained from our study shed light on prior experimental observations and elucidate the key role played by changes in hydrogen bonding network in dictating the hydrophilicity of the brushes. In addition to the temperature-mediated reduction in height of the PE brush layer, we study the influence of temperature on several properties of the brush-supported water molecules such as the number distribution and kinetics of water-water and water-PE hydrogen bonds (HBs), mass density, orientational tetrahedral order parameter, and the activation energy associated with the self-diffusion of water. The effect of temperature on the properties of the brush-trapped counterions is also quantified via changes in their solvation structure, counterion-water radial distribution function (RDF), and the translational mobility. Our findings unravel several hitherto unknown phenomena such as the enhanced propensity for attaining PE brush-induced “water-in-salt”-like scenarios at elevated temperatures, and an increase in the activation energy for the self-diffusion of water with increase in the degree of PE brush-induced nanoconfinement (from ‘bulk water’ to ‘sparsely grafted PE brushes’ to ‘densely grafted PE brushes’). Further, we anticipate that this work will be helpful in improving our predictive capabilities of the temperature response of water-swollen PE brushes, which would be instrumental for better design of several PE brush systems used in a multitude of applications.

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

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2000

Molecular Dynamics Simulations of the Thermal Conductivity of Single-Wall Carbon Nanotubes

Carbon nanotubes (CNT) have very attractive electronic, mechanical. and thermal properties. Recently, measurements of thermal conductivity in single wall CNT mats showed estimated thermal conductivity magnitudes ranging from 17.5 to 58 W/cm-K at room temperature. which are better than bulk graphite. The cylinderical symmetry of CNT leads to large thermal conductivity along the tube axis, additionally, unlike graphite. CNTs can be made into ropes that can be used as heat conducting pipes for nanoscale applications. The thermal conductivity of several single wall carbon nanotubes has been calculated over temperature range from l00 K to 600 K using non-equilibrium molecular dynamics using Tersoff-Brenner potential for C-C interactions. Thermal conductivity of single wall CNTs shows a peaking behavior as a function of temperature. Dependence of the peak position on the chirality and radius of the tube will be discussed and explained in this presentation.

Osman, M.↗
DOE OSTI2020

Horizons of modern molecular dynamics simulation in digitalized solid freeform fabrication with advanced materials

Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid free form fabrication. Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror and thereafter an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is becoming ubiquitous. This article is intended to introduce the most recent advances in the computation benefit specific to the area of solid free form fabrication as these systems are traversing through the journey of digitalisation and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations

Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. Here, we present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2021

Molecular dynamics simulations of ultrafast radiation induced melting at metal–semiconductor interfaces

Metal–semiconductor contacts in silicon carbide (SiC) diodes endure damages at the interface when exposed to harsh radiation environments. Due to the rapid rise in temperature and ultrafast cooling that follows the radiation impact, the structural properties of the materials can be altered through melting, recrystallization, and amorphization. A detailed understanding of the material failure modes at the interface is lacking, specifically at the nanoscale. We use molecular simulations to investigate the ultrafast melting at tungsten (W)–SiC interfaces following radiation damage and apply deep learning techniques to track the transient evolution of the local molecular structures. We show that W near the radiation track undergoes melting and, eventually, most of it recrystallizes with a noticeable degree of undercooling, while SiC is rendered permanently amorphous. The observation of local undercooling in W films is important as it can affect the device performance even before the bulk melting temperature of the material is reached. We also show that at high temperatures, the interface undergoes a fracture-like failure. The results presented here are significant in understating the different failure modes of SiC diode materials.

Ashwin Ravichandran↗
DOE OSTI2020

A comparative study of the effectiveness of empirical potentials for molecular dynamics simulations of borosilicate glasses

Despite their practical importance, atomistic modeling of B 2 O 3 containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO 2 ]/[B 2 O 3 ]) =2 and varying R (ratio of [Na 2 O]/[B 2 O 3 ] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment was established on the basis of the short- and medium-range structure features, such as boron N 4 values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

Carbon fiber synthesis from pitch: Insights from ReaxFF based molecular dynamics simulations

Despite the potential to lowering fabrication costs of pitch-based carbon fiber (CF) compared to polyacrylonitrile-based CF, large scale utilization is currently limited by the complexity in manufacturing. This is mainly due to the lack of consistency between the pitch structure/chemistry and CF properties under different processing conditions. Here, we employed full atomistic simulations based on reactive force fields to investigate the conversion of pitch precursors to CFs. Additionally, it was found that at the carbonization stage, the consumption rates of precursors, gas formations, and solid yields (carbonized CFs) are all dependent on the composition and arrangement of molecular moieties in the precursor. Specifically, precursors with sp 3 carbons in the cata-condensed core area have limited decomposition pathways (i.e., limited reaction sites), which leads to the formation of carbonized CFs with similar compositions. Furthermore, carbonized CFs possess distinct structural characteristics, and finally lead to the formation of CFs with different moduli. By analyzing CF crystallinity and elemental composition, general rationales are proposed for the selection of optimal precursors and processing conditions, and thus shed lights on parameter screening during industrial CF synthesis. For instance, inclusion of sp 3 carbons in cata- condensed rings can be preferred over the presence of sp 3 carbons on side chains.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Mechanism of Solid Electrolyte Interphase Formation on Lithium Anodes Studied by Molecular Dynamics Simulations

We investigated the initial formation of the SEI on lithium-metal anodes in presence of several electrolyte compositions using ReaxFF, in absence of a bias potential. We observed the formation of well-defined regions with different density and oxidation states of Li. Dissolution of Li metal atoms results in an uneven distribution of Li in the liquid phase. The electrolyte decomposition mechanism was studied for various electrolyte mixtures. Overall, the reaction mechanisms observed with ReaxFF are in good agreement with previous reports using DFT. However, we found some discrepancies in specific predictions by the ReaxFF that should be addressed in future work.

25 ENERGY STORAGE↗