Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Low-temperature carbonization of polyacrylonitrile/graphene carbon fibers: A combined ReaxFF molecular dynamics and experimental study

Graphene inclusion in a polymer matrix is a promising route to significantly enhance the mechanical properties of low-grade carbon fibers (CFs). Using ReaxFF molecular dynamics simulation, the atomistic mechanism leading to this enhancement is investigated. We demonstrate that the graphene edges along with the nitrogen and oxygen functional groups play a catalytic role and act as seeds to expedite alignment of the all-carbon rings, which are starting sites for the growth of graphitic structures. To examine the role of this proposed mechanism that enhances the graphitic structure of PAN/graphene CFs, we discuss the experimental results wherein the PAN/graphene CFs carbonized at 1250 C demonstrate 91% (from 632 to 1207 MPa) increase in strength and 101% (from 88 to 177 GPa) enhancement in Young’s modulus compared to PAN-based CFs carbonized at 1500 C. In conclusion, these enhanced mechanical properties of low-grade carbon fibers achieved via graphene inclusion at decreased carbonization temperature provide a means to realize both energy savings and cost reduction.

42 ENGINEERING↗
DOE OSTI2023

Ab Initio Molecular Dynamics Spectra for Characterization of Hydrated Al 2 O 3 Supported MoO x

The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Effect of Aromaticity on Soot Formation: A Reactive Molecular Dynamics Study

This study investigates the formation of soot from n-dodecane and its blend with m-xylene under pyrolytic and low-oxygen conditions with reactive molecular dynamics (MD) simulations. ReaxFF force field is used to simulate the interactions of fuel and oxygen molecules at engine-relevant conditions (2500-3500 K and 70-75 bar). Fuel molecules with carbon atom density of 0.03 g/cm 3 are dispersed into a cube and constant temperature simulations are carried out with a Nosé-Hoover thermostat. Pyrolysis of n-dodecane and the chemical pathway towards the formation of aromatic hydrocarbons are investigated. The initial pyrolysis rate of n-dodecane is observed to be slightly higher with the addition of aromatic compounds and oxygen. While the first aromatics formed from the decomposed n-dodecane are phenyl radicals at 2500 K, decomposition of the aliphatic fuel compound is observed to be necessary to build a carbon cluster, and eventually soot. Morphological properties of the created soot particles are also analyzed. The maturity, size and aromatics content of the soot particle are slightly affected with the addition of m-xylene and oxygen at 3500 K within the 3 ns simulation time possibly due the formation of larger and more condensed aromatic rings, resulting in larger soot molecules. The outcomes of this study provide important insights into the formation mechanism of soot and its morphological characteristics depending on varying fuel concentrations at pyrolysis and low oxygen engine conditions. Results of this analysis will further contribute to the understanding of contrail formation on generated soot particles originated from the combustion of different chemical components with particle chemistry matching companion experiments.

soot formation↗
DOE OSTI2021

High-salinity brine desalination with amine-based temperature swing solvent extraction: A molecular dynamics study

Solvent extraction desalination (SED) is a promising technique for efficiently treating high salinity brines. However, only a few studies in the literature have investigated the influence of the molecular characteristics of the organic solvents on the emergent performance of the separation process. Thus, in this work we investigate the influence of the molecular structure of four different secondary amines on the desalination process of a concentrated aqueous solution of NaCl. Molecular dynamics simulations and free energy calculations are used to obtain the thermodynamic and structural properties of the pure solvent bulk phase, as well as a homogeneous solvent-water bulk mixture. Furthermore, solvent-water and solvent-brine interfacial systems are used to explicitly analyze the water partitioning behavior in the relevant temperature range for SED. The water solubilization in these secondary amines is found to be an exothermic process. Additionally, although the hydrogen-bonding interaction between the water and solvent molecules are beneficial, this interaction is not necessarily the dominant factor for ranking the water solubility performance of the different solvents. Finally, the excess volume of mixing of the bulk solvent-water mixture, as well as the electrostatic surface features of the solvents are also found to provide important clues for screening solvent behavior.

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

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗
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↗
NASA NTRS2021

Molecular Dynamics Investigation of the Structural and Mechanical Properties of Off-Stoichiometric Epoxy Resins

Molecular dynamic (MD) simulations were performed to compute the mechanical properties of off-stoichiometric epoxy resins as a function of hardener/epoxy mixture ratio (r). Properties were characterized in relation to their microscopic structures. Such resins have been used recently for adhesive-free bonding of large-scale composite structures using the co-curing-ply method. In this process, two partially precured composite panels with hardener-poor (HP) off-stoichiometric resins are coupled with ply(ies) of complementary hardener-rich (HR) formulations and then cured simultaneously. This bonding process has the potential to produce reliable and certifiable composite joints without the need for additional fasteners, which are often required for many conventional bonding methods because even small amounts of contamination can cause a weak bond. The reflow and mixing of the HP/HR resin in this bonding process result in a joint with no discernable interface that should not be susceptible to surface contamination. However, incomplete mixing of the two offset resins may result in chemical heterogeneity of the cured polymeric joint. Thus, different r values may be obtained across the joint. Classical MD simulations were performed to compute the Young’s modulus of polymers with different r values and correlate their properties to network structures. High stiffness was associated with molecular packing due to chemical crosslinking, leading to a single network structure. Moreover, the networks became denser as the ratio approached the stoichiometric value r = 1. Thus, the r = 1 systems were single clusters, with high stiffness, high molecular weight, and a high degree of crosslinking. Structural properties such as radius of gyration and mean square displacement were determined to investigate the variation in the stiffness with respect to r. This MD simulation study was validated with experimental measurements.

Chang Woon Jang↗
DOE OSTI2022

Molecular dynamics study of silicon atomic layer etching by chlorine gas and argon ions

Classical molecular dynamics (MD) is used to simulate atomic layer etching processes of silicon by alternating exposure to chlorine gas and argon ions. In order to validate our model, a rigorous comparison is done with ion beam experiments found in the literature [Park et al., Jpn. J. Appl. Phys. 44, 389 (2005)]. Here it is shown that the etch per cycle (EPC) as a function of argon ion energy from simulations is in quantitative agreement with experimental results if the correct argon ion fluence is used in the simulations. The EPC as a function of ion irradiation time and amount of chlorine exposure also show good agreement with the experiment. The MD simulations also show the formation of an amorphous silicon region with chlorine atoms mixed uniformly throughout following ion bombardment. Finally, the etch products during the ion irradiation step are analyzed and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

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

Atomistic insights into the heat conductance across the interfaces between erythritol and different metals: A non-equilibrium molecular dynamics study

The heat conductance across the interface between a phase change material (PCM) and the metal shell is important for thermal design of latent heat storage systems. Non-equilibrium molecular dynamics simulations were carried out in this work to unravel, at the atomistic level, the interfacial heat conductance mechanisms between erythritol (the most promising medium-temperature polyol PCM) and two different metals (i.e., copper and aluminum). The interfacial contact ratio was introduced as a parameter to represent different contact conditions at the interfaces, for example, different surface roughness levels of the metals. It was found that with the increase of interfacial contact ratio, the temperature difference between the erythritol/metal interfaces decreases, whereas the interfacial heat transfer coefficient increases gradually. At the perfect contact condition, i.e., the interfacial contact ratio reaches 100%, the interfacial heat transfer coefficients across the erythritol/Cu and erythritol/Al interfaces were predicted to be 398 ± 26 MW/m 2 K and 544 ± 11 MW/m 2 K, respectively. The interfacial heat conductance of erythritol/Al is better than that of erythritol/Cu. The mismatch of atomic spectral density across the interfaces between erythritol and the two metals was analyzed and distinguished. These results could guide the selection of packaging materials and the thermal design of erythritol-based latent heat storage systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
DOE OSTI2020

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Graphic contrastive learning analyses of discontinuous molecular dynamics simulations: Study of protein folding upon adsorption

A comprehensive understanding of the interfacial behaviors of biomolecules holds great significance in the development of biomaterials and biosensing technologies. In this work, we used discontinuous molecular dynamics (DMD) simulations and graphic contrastive learning analysis to study the adsorption of ubiquitin protein on a graphene surface. Our high-throughput DMD simulations can explore the whole protein adsorption process including the protein structural evolution with sufficient accuracy. Contrastive learning was employed to train a protein contact map feature extractor aiming at generating contact map feature vectors. Subsequently, these features were grouped using the k-means clustering algorithm to identify the protein structural transition stages throughout the adsorption process. The machine learning analysis can illustrate the dynamics of protein structural changes, including the pathway and the rate-limiting step. Our study indicated that the protein–graphene surface hydrophobic interactions and the π–π stacking were crucial to the seven-stage adsorption process. Upon adsorption, the secondary structure and tertiary structure of ubiquitin disintegrated. The unfolding stages obtained by contrastive learning-based algorithm were not only consistent with the detailed analyses of protein structures but also provided more hidden information about the transition states and pathway of protein adsorption process and structural dynamics. Our combination of efficient DMD simulations and machine learning analysis could be a valuable approach to studying the interfacial behaviors of biomolecules.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Defect formation and bending properties in graphite under He atom implantation investigated by molecular dynamics method

The formation and evolution of radiation defects and their effects on bending property of a single graphite under low energy helium atom implantation have been studied at atomic scale through molecular dynamics simulations. Various simple interstitial and vacancy defects are observed during the implantation process, which is a function of incident energy and cross section of helium-carbon interaction. Large defects are observed after continuous helium atom implantations. The effects of incident energy on the density, size and distribution of defects are analyzed in detail. Furthermore, the kinetic evolution process of radiation defects at different annealing temperatures have also been simulated, including the migration, aggregation, transformation, coalescence, and self-healing. Here, the large defect clusters, crossing several carbon layers, have been confirmed to be the primary reasons that affect the bending property of graphite. All these results provide a new understanding for further applications of graphite in nuclear reactors.

Atomic implantation↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

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

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗
NASA NTRS1994

Dielectric Properties of Poly(ethylene oxide) from Molecular Dynamics Simulations

The order, conformations and dynamics of poly(oxyethylene) (POE) melts have been investigated through molecular dynamics simulations. The potential energy functions were determined from detailed ab initio electronic structure calculations of the conformational energies of the model molecules 1,2-dimethoxyethane (DME) and diethylether. The x-ray structure factor for POE from simulation will be compared to experiment. In terms of conformation, simulations reveal that chains are extended in the melt relative to isolated chains due to the presence of strong intermolecular O...H interactions, which occur at the expense of intramolecular O...H interactions. Conformational dynamics about the C-C bond were found to be significantly faster than in polymethylene, while conformational dynamics about the C-O bond even faster than the C-C dynamics. The faster local dynamics in POE relative to polymethylene is consistent with C-13 NMR spin-lattice relaxation experiments. Conformational transitions showed significant second-neighbor correlation, as was found for polymethylene. This correlation of transitions with C-C neighbors was found to be reduced relative to C-O neighbors. Dielectric relaxation from simulation will also be compared with experiment.

Smith, Grant D.↗