Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Further, experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Electronic stopping in molecular dynamics simulations of cascades in 3C–SiC

Here, we investigate the effect of the electronic stopping power on defect production due to ion irradiation of cubic silicon carbide using molecular dynamics simulations. We simulate 20 keV and 30 keV Si and C ions, with and without the electronic energy loss. The results show that the electronic stopping effects are more profound in the case of C irradiation, where the ratio of the electronic energy loss S e to the nuclear energy loss S n is much larger compared to the ratio for Si ions. These findings indicate that this ratio plays a role in the effect of the electronic stopping on ion irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Dynamics Simulations of Complexation of Am(III) with a Preorganized Dicationic Ligand in an Ionic Liquid

Preorganized ligands with imidazolium arms have been found to be highly selective in extracting Am(III) into ionic liquids (ILs), but the detailed structure and mechanism of the complexation process in the ionic solvation environment are unclear. Here, we carry out molecular dynamics simulation of the complexation of Am(III) with a preorganized 1,10-phenanthroline-2,9-dicarboxamide complexant (L) functionalized with alkyl chains and imidazolium cations in the butylmethylimidazolium bistriflimide ([BMIM][NTf 2 ]) IL. Both Am:L (1:1) and Am:L 2 (1:2) complexes are examined. In the absence of the ligand, Am(III) is found to be coordinated by six NTf 2 anions via nine O donors in the first solvation shell. In the Am:L complex, Am(III) is coordinated to the ligand via two O donors and four NTf 2 anions via seven O donors in the first coordination shell. In the Am:L 2 complex, Am(III) is coordinated to the two ligands via four O donors and four NTf 2 anions via five O donors. The imidazolium arms of the ligands play an important role in the secondary solvation environment by attracting NTf 2 anions closer to the metal center. As a result, we find that the binding free energy for the second L 2+ ligand is twice that for the first L 2+ ligand, making the Am:L 2 complex significantly more stable than the Am:L complex. This study highlights the multiple factors and tunability in using preorganized ligands with charged functional groups in an ionic solvation environment, which could hold the key to achieving desired selectivity in ion extraction efficiency.

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

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗
DOE OSTI2021

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Pressure-induced phase transformations in Nb 2 O 5 from ab initio molecular dynamics simulations

Here, the equations of state for the H-, L- and T-Nb 2 O 5 polymorphs were investigated with density functional theory. The shock Hugoniot curves were predicted from ab initio molecular dynamics simulations along 300–2500 K isotherms. Enthalpies derived from simulations at 300 K show that the H → L transition occurs at 6.2 GPa, in agreement with the measured value of ~5.9 GPa, and is more favorable than the H → T transition. Pressure-induced amorphization (PIA) of L-Nb 2 O 5 is in the range ~18-25 GPa, consistent with experiments reporting PIA above 21.4 GPa. The simulated PIA of T-Nb 2 O 5 between ~16 and 26 GPa reproduces PIA observed above 19.2 GPa.

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

Interpretable artificial intelligence and exascale molecular dynamics simulations to reveal kinetics: Applications to Alzheimer's disease

The rapid increase in computing power, especially with the integration of graphics processing units, has dramatically increased the capabilities of molecular dynamics simulations. To date, these capabilities extend from running very long simulations (tens to hundreds of microseconds) to thousands of short simulations. However, the expansive data generated in these simulations must be made interpretable not only by the investigator who performs them but also by others as well. Here, we demonstrate how integrating learning techniques, such as artificial intelligence, machine learning, and neural networks, into analysis pipelines can reveal the kinetics of Alzheimer's disease (AD) protein aggregation. Finally, we review select AD targets, describe current simulation methods, and introduce learning concepts and their application in AD, highlighting limitations and potential solutions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

Ion Pairing, Clustering and Transport in a LiFSI-TMP Electrolyte as Functions of Salt Concentration using Molecular Dynamics Simulations

Battery capacity is highly related to ion-pairing mechanisms in electrolytes, since a cluster formation can lead to dead Li formation, reducing the number of charge carriers and leading to capacity fading. We use molecular dynamics simulations to model an electrolyte comprising trimethyl phosphate (TMP) solvent and a lithium bis(fluorosulfonyl)imide (LiFSI) salt, exploring effects of salt concentration on solvation and ion-transport. We simulate the LiFSI-TMP electrolyte for salt concentrations of 0.7, 1.43 and 3.82 molar. A statistical analysis was performed to study ion-pairing, clustering, diffusivity, conductivity, and coordination of Li-ions, providing insights into relations between molecular structures and transport properties. Molecular structure of ionic components changes as concentration increases, from a predominant solvent separated ion pair (SSIP) and contact ion pair (CIP) to aggregate salt (AGG) and ionic cluster formation. Given the formation of the ionic cluster, the diffusion mechanism followed by Li-ions changes from a hopping/exchange to a vehicular mechanism as concentration increases; this is reflected in a decrease of ionic conductivities. Ionicity was also calculated to reveal how the ionic motion changes from an uncorrelated to a correlated one as the salt concentration increases. Furthermore, we also compared our results with experimental calculations performed for similar electrolyte systems

25 ENERGY STORAGE↗
DOE OSTI2020

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Determination of the structure of lithium niobosilicate glasses by molecular dynamics simulation with a new Nb-O potential

A partial charge empirical potential for the Nb-O pair has been developed based on existing parameters of other pairs in a Teter potential to enable molecular dynamics simulations of both lithium niobate crystal and lithium niobosilicate glasses. The developed potential is capable of describing structural features of lithium niobosilicate glasses in a wide composition range, including niobium coordination number and bond length, density, bond angle distribution, polyhedral distribution and fraction of bridging oxygen. Furthermore, the results obtained using this new potential show good agreement with experimental data of density and structure from Raman spectroscopy of lithium niobosilicate glasses and X-ray absorption spectroscopy near Nb K-edge.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

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↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

The identity of metal ions surrounding DNA is key to its biological function and materials applications. In this work, we compare atomistic molecular dynamics simulations of double strand DNA (dsDNA) with four alkaline earth metal ions (Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ ) to elucidate the physical interactions that govern DNA–ion binding. Simulations accurately model the ion–phosphate distance of Mg 2+ and reproduce ion counting experiments for Ca 2+ , Sr 2+ , and Ba 2+ . Our analysis shows that alkaline earth metal ions prefer to bind at the phosphate backbone compared to the major groove and negligible binding occurs in the minor groove. Larger alkaline earth metal ions with variable first solvation shells (Ca 2+ , Sr 2+ , and Ba 2+ ) show both direct and indirect binding, where indirect binding increases with ion size. Mg 2+ does not fit this trend because the strength of its first solvation shell predicts indirect binding only. Ions bound to the phosphate backbone form fewer contacts per ion compared to the major groove. Within the major groove, metal ions preferentially bind to guanine–cystosine base pairs and form simultaneous contacts with the N7 and O6 atoms of guanine. Overall, we find that the interplay among ion size, DNA–ion interaction, and the size and flexibility of the first solvation shell are key to predicting how alkaline earth metal ions interact with DNA.

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

Ab-Initio molecular dynamics simulations of binary NaCl-ThCl 4 and ternary NaCl-ThCl 4 -UCl 3 molten salts

Molten salt reactor (MSR) with Thorium (Th) fuel cycle has attracted growing attention due to its merits such as safety, low radioactive waste production, and non-proliferation. The design and safe operation of MSR rely on a thorough understanding of the thermophysical properties of molten salts over a wide range of composition and temperature. However, experimental data of Th-based molten salts are limited due to the inherent challenges of dealing with corrosive molten salts and the radioactive nature of actinides. Here, in this work, thermophysical properties including density, heat capacity, thermal expansion coefficient, and mixing energy of binary NaCl-ThCl 4 and ternary NaCl-ThCl 4 -UCl 3 molten salts are explored using ab-initio molecular dynamic simulations (AIMD) with the dDsC dispersion correction. The calculated mixing energy of binary NaCl-ThCl 4 exhibits a minimum close to the eutectic composition. The heat capacity of the mixtures is linearly dependent on the mole fraction of each component. 7-fold and 6-fold coordinated Th complexes are dominant in the mixtures. 8-fold coordinated Th complex increases with ThCl 4 fraction due to the formation of network structures. The average coordination number of Th exhibits a minimum near the eutectic composition. The minimum mixing energy for the ternary mixtures is observed in systems with a composition close to [NaCl] 0.5 [ThCl 4 ] 0.25 [UCl 3 ] 0.25 . The density positively deviates from the ideal solution in mixtures near this composition. These results are important to fill the data and knowledge gap of ThCl 4 molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Structures of fluoride containing aluminosilicate low activity nuclear waste glasses: A molecular dynamics simulations study

Fluorine and other halides commonly exist in nuclear waste forms, and due to their volatile nature, halide retention poses an issue affecting waste loading during vitrification. The compositional effect on fluorine incorporation in aluminosilicate glasses is investigated through molecular dynamics simulations. Oxygen and fluorine coordination numbers around glass former and modifier cations, bond angle distributions, and medium range structure features such as Q n distributions, ring size distributions and neutron diffraction structure factors were calculated. It was found that fluorine has higher preference to bond to Ca 2+ than to Na + , both in the melt and the glass, and there is no Si-F bond formation in the glass but they do exist in the melt. Consequently, CaO for Na 2 O substitution can be an effective way to help fluorine retention without significantly changing the glass chemistry. Furthermore, these results thus provide insights on fluorine incorporation in the aluminosilicate nuclear waste glasses and the strategy on how to improve fluorine retention both in the glass and the melt.

36 MATERIALS SCIENCE↗