Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular dynamic simulations”

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 OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
DOE OSTI2022

Accelerators for Classical Molecular Dynamics Simulations of Biomolecules

Atomistic Molecular Dynamics (MD) simulations provide researchers the ability to model biomolecular structures such as proteins and their interactions with drug-like small molecules with greater spatiotemporal resolution than is otherwise possible using experimental methods. MD simulations are notoriously expensive computational endeavors that have traditionally required massive investment in specialized hardware to access biologically relevant spatiotemporal scales. Our goal is to summarize the fundamental algorithms that are employed in the literature to then highlight the challenges that have affected accelerator implementations in practice. We consider three broad categories of accelerators: Graphics Processing Units (GPUs), Field-Programmable Gate Arrays (FPGAs), and Application Specific Integrated Circuits (ASICs). These categories are comparatively studied to facilitate discussion of their relative trade-offs and to gain context for the current state of the art. We conclude by providing insights into the potential of emerging hardware platforms and algorithms for MD.

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

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Unraveling the Heterogeneous but Ordered Microstructure of the Nonionic Deep Eutectic Solvent Formed by Lauric Acid and N -Methylacetamide

The nonionic deep eutectic solvent, formed by lauric acid (LA) and N-methylacetamide (NMA), has been shown to have a heterogeneous molecular structure in which the LA and NMA form nonpolar and polar domains, respectively. Previous vibrational spectroscopy experiments demonstrated that the ability of the LA domains to solvate compounds was limited to long carbon chains, whereas other nonpolar molecules, such as W(CO) 6 , were found to be solvated by both LA and NMA. These experiments were not fully compatible with the previously proposed micelle-like structure of the nonpolar domains of the LA-NMA DES. In this work, the modeling of the DES molecular structure is pursued using classical molecular dynamics simulations. The new classical model reproduces both the SAXS structural factors and the previously experimentally derived interaction map for these LA-NMA DESs. In addition, the simulation also shows that LA-NMA DESs form highly organized LA aggregates that are difficult to disorganize. Further evidence of the correct description provided by the newly derived model is obtained using a moderately polar probe: chloroform-d. Computations using the classical model have a good agreement with the solvation behavior of the probe derived from experiments, in which the location of the probe is found to be mostly within the polar domain of the DES. The computational model also demonstrates that the probe solvation is a consequence of the tightly packed LA structure, which causes nonpolar molecules to be located at the interphase of the DES nonpolar domains.

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

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
DOE OSTI2024

Data and Code for Atomic Scale Etching of Diamond: Insights from Molecular Dynamics Simulations

This work investigates the effects of argon ions, hydrogen atoms, and hydrogen ions on the diamond (100) surface using classical molecular dynamics simulations. The purpose of this investigation was to asses plasma processing techniques for applications in quantum device manufacturing. The simulations suggest that combining argon ion smoothing with selective, near threshold energy H removal of amorphous C could be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices. Results were found to differ significantly with interatomic potential, and an analysis of these differences was also carried out. Included in this repository are LAMMPS source files, input scripts, and plotting scripts required to reproduce the data. Also included are the output data required to make all the plots included in the associated publication.

Brenner↗
DOE OSTI2022

Superhydrophilicity of $α$-alumina surfaces results from tight binding of interfacial waters to specific aluminols

Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. Here, contact angles at the alumina (0001) and ($11\bar{2}0$) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

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

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
NASA NTRS2020

Ultrafast Melting at Metal-Semiconductor Interfaces: A Molecular Dynamics Study

Quantifying the thermodynamics of interfacial melting is important for several applications such as laser machining and designing radiation hardened electronic devices. Due to the short time and length scales that are involved in the process, molecular simulations provide an ideal path for understating the interfacial melting phenomena. In this work we investigate tungsten (W) – silicon carbide (SiC) interfaces which are relevant in modern high-power electronic devices. Classical molecular dynamics simulations along with the two-temperature model (TTM) were used to study the radiation induced melting of the interface while a neural network-based machine learning framework was used to track the transient re-crystallization process. The damage endured by the semiconductor at the interface following the radiation event is discussed in detail while the re-crystallization of the metal is described in terms of the classical nucleation theory.

Ashwin Ravichandran↗
DOE OSTI2022

A molecular dynamics study of laser-excited gold

The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melting process, leading to large discrepancies between the predicted threshold energy density for complete melting, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. [Science 360, 1451–1455 (2018)]. We recreate the experimental electron diffraction data, employing both a constant and temperature-dependent electron–ion equilibration rate. In all cases, we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sulfonated polybenzimidazole membrane with graphene oxide additive for 2,3-butanediol/water separation: A molecular simulation

Membrane separation for 2,3-butanediol (2,3-BDO) recovery from fermentation broth is highly valued for sustainable and renewable processes, but it requires efficient membrane materials. Here, this work evaluates the sulfonated polybenzimidazole (sPBI) and its graphene oxide (GO) doped composite membrane for separating 2,3-BDO and water via atomistic simulations. Density functional theory calculations are applied to identify various forms of sPBI structures and quantify their binding interactions with 2,3-BDO and water. Classical molecular dynamic simulations are used to evaluate the structural changes, diffusivity, and selectivity of 2,3-BDO and water in different sPBI models, GO surfaces, and GO-doped sPBI composite models. Our results suggest that sPBI slightly increases the crystallinity of the membrane structures, enhances the adsorption strength for both 2,3-BDO and water, and improves the water/2,3-BDO selectivity by 2–3 times. The GO surfaces display a maximum selectivity at a surface coverage of 0.1–0.15 for both hydroxyl and epoxy surface groups. The addition of GO flakes to sPBI creates new interaction sites for 2,3-BDO and water at the interface of sPBI and GO, and the water/2,3-BDO selectivity of GO-doped sPBI models is further increased up to 3 times. This work illustrates how the integrated addition of sPBI and GO flakes offers a promising approach to selective separation of 2,3-BDO and water, providing theoretical guidance for polybenzimidazole-based membranes in the potential application of 2,3-BDO recovery.

2,3-butanediol↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tracking structural solvent reorganization and recombination dynamics following e – photoabstraction from aqueous I – with femtosecond x-ray spectroscopy and scattering

Here we present a sub-picosecond resolved investigation of the structural solvent reorganization and geminate recombination dynamics following 400 nm two-photon excitation and photodetachment of a valence p electron from the aqueous atomic solute, I–(aq). The measurements utilized time-resolved X-ray Absorption Near Edge Structure (TR-XANES) spectroscopy and X-ray Solution Scattering (TR-XSS) at the Linac Coherent Light Source x-ray free electron laser in a laser pump/x-ray probe experiment. The XANES measurements around the L 1 -edge of the generated nascent iodine atoms (I 0 ) yield an average electron ejection distance from the iodine parent of 7.4 ± 1.5 Å with an excitation yield of about 1/3 of the 0.1M NaI aqueous solution. The kinetic traces of the XANES measurement are in agreement with a purely diffusion-driven geminate iodine–electron recombination model without the need for a long-lived (I 0 :e–) contact pair. Nonequilibrium classical molecular dynamics simulations indicate a delayed response of the caging H 2 O solvent shell and this is supported by the structural analysis of the XSS data: We identify a two-step process exhibiting a 0.1 ps delayed solvent shell reorganization time within the tight H-bond network and a 0.3 ps time constant for the mean iodine–oxygen distance changes. The results indicate that most of the reorganization can be explained classically by a transition from a hydrophilic cavity with a well-ordered first solvation shell (hydrogens pointing toward I–) to an expanded cavity around I 0 with a more random orientation of the H 2 O molecules in a broadened first solvation shell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗