Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Comparison of Chain Conformation of Poly(vinyl alcohol) in Solutions and Melts from Quantum Chemistry Based Molecular Dynamics Simulations

Confirmations of 2,4-dihydroxypentane (DHP), a model molecule for poly(vinyl alcohol), have been studied by quantum chemistry (QC) calculations and molecular dynamics (MD) simulations. QC calculations at the 6-311G MP2 level show the meso tt conformer to be lowest in energy followed by the racemic tg, due to intramolecular hydrogen bond between the hydroxy groups. The Dreiding force field has been modified to reproduce the QC conformer energies for DHP. MD simulations using this force field have been carried out for DHP molecules in the gas phase, melt, and CHCl3 and water solutions. Extensive intramolecular hydrogen bonding is observed for the gas phase and CHCl3 solution, but not for the melt or aqueous solution, Such a condensed phase effect due to intermolecular interactions results in a drastic change in chain conformations, in agreement with experiments.

Jaffe, Richard↗
DOE OSTI2020

Fast neutrino cooling of nuclear pasta in neutron stars: Molecular dynamics simulations

The direct Urca process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crusts of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct Urca reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct Urca reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified Urca process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Molecular dynamics simulation of vacancy and void effects on strain-induced martensitic transformations in Fe-50 at.% Ni model concentrated solid solution alloy

Strain-induced martensitic transformations can improve the strength and ductility of some face centered cubic (FCC) metals and alloys. Irradiation-induced defects such as vacancies, dislocation loops, and voids activate martensitic transformations over a wider range of mechanical loading conditions than in pristine material. Furthermore, the mechanisms underlying irradiation-enabled martensite transformations remain unclear. In this work, we use molecular dynamics simulations to study the effect of vacancies and voids on strain-induced martensitic transformations in a model concentrated solid solution alloy Fe-50 at.% Ni. It has been found that single vacancies have no resolvable effect on the transformation because they reduce the stacking fault energy by a relatively insignificant margin and do so only if the vacancy is located on the stacking fault plane. Voids, however, activate the martensite transformation through shear strain accumulation around the void due to dislocation pileup. The larger the void, the more pronounced this effect.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

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

Sticking Coefficients of Fusion Reactor Impurities from Molecular Dynamics Simulations for the Design of Cryopumps

A cryopump can be utilized as an impurity removal component of a direct internal recirculation (DIR) system for the fusion fuel cycle. The DIR facilitates a low fuel inventory by continuously pumping unburnt fuel while removing impurities from the fusion exhaust stream. A cryopump can target multiple impurity species by maintaining a temperature lower than the gas triple-point temperature that promotes desublimation. The desublimation/condensation of gases in cryopumps can be characterized by the sticking coefficient, which is defined as the probability for a gas particle to stick to a (cryo-)surface upon collision. The sticking coefficient is one of the important design/operation parameters for cryopumps, and it depends on a variety of surface and gas properties. Here, in this study, molecular dynamics simulations were utilized to estimate the sticking coefficients of typical fusion gas impurity species N 2 , CO 2 , and CH 4 over a Cu surface for a range of gas temperatures and surface coverages. The molecular dynamics study showed that the sticking coefficients for gases decrease with an increase in gas temperature. The presence of a single full monolayer of condensate on the metallic surface showed an adverse effect on the sticking of gases; however the sticking improved with two full monolayers of condensate on the surface. The sticking of gases over the mixed condensate on a surface was more favorable than the condensate of the same species for N 2 and CH 4 , with an exception for CO 2 , which showed a decrease in sticking over the mixed condensate.

cryopump↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

All atom molecular dynamics (MD) simulations of planar Na+-counterion-neutralized polyacrylic acid (PAA) brushes are performed for varying degrees of ionization (and thereby varying charge density) and varying grafting density. Variation in the PE charge density (or degree of ionization) and grafting density leads to massive changes of the properties of the PE molecules (quantified by the changes in the height and the mobility of the PE brushes) as well as the local arrangement and distribution of the brush-supported counterions and water molecules within the brushes. The effect on the counterions is manifested by the corresponding variation of the counterion mobility, counterion concentration, extent of counterion binding to the charged site of the PE brushes, water-in-salt-like structure formation, and counterion-water–oxygen radial distribution function within the PE brushes. On the other hand, the effect on water molecules is manifested by the corresponding variation of water–oxygen–water–oxygen RDF, local water density, water–water and water–PE functional group hydrogen bond networks, static dielectric constant of water molecules, orientational tetrahedral order parameter, and water mobility. Enforcing such varying degree of ionization of weak polyelectrolytes is possible by changing the pH of the surrounding medium. Thus, our results provide insights into the changes in microstructure (at the atomistic level) of weak polyionic brushes at varying pH. Finally, we anticipate that this knowledge will prove to be vital for the efficient design of several nano-scale systems employing PE brushes such as nanomechanical gates, current rectifiers, etc.

42 ENGINEERING↗
DOE OSTI2022

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2020

Statistically averaged molecular dynamics simulations of hydrogen diffusion in magnesium and magnesium hydrides

Magnesium has a different crystal structure from its dihydride with hydrogenation leading to a phase transition from the hexagonal closely packed Mg into a tetragonal α-MgH 2 rutile type structure. Such materials exhibit complex hydrogen uptake and release kinetics because hydrogen diffusivities significantly change when the crystal structure changes. To provide a foundational understanding of (de)hydrogenation kinetics that is applicable to all stages of the reaction, we performed statistically averaged molecular dynamics simulations to derive hydrogen diffusivities as a function of temperature and hydrogen content for both magnesium and magnesium hydride. Furthermore, our studies confirm that hydrogen diffusivities in magnesium hydride are much lower than in magnesium, in agreement with experimental data. Additionally, we observe that in either magnesium or magnesium hydride, higher hydrogen compositions result in reduced diffusivities. The latter was not revealed by prior experiments, which were conducted at fixed hydrogen composition. Finally, we discover a non-Arrhenius behavior in magnesium hydride. The physical origin of this behavior is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗
DOE OSTI2025

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗
DOE OSTI2022

On the equivalence of nonequilibrium and equilibrium measurements of slip in molecular dynamics simulations

Here, we show that nonequilibrium and equilibrium measurements of slip are consistent, provided the hydrodynamic wall location associated with the equilibrium measurement is properly taken into account. The latter is a strong function of the fluid state and wall-fluid interaction and cannot be neglected as it typically has been. Our results are based on an alternative approach for calculating the hydrodynamic wall location which alleviates most of the difficulties associated with its calculation via a Green-Kubo integral that appear to have contributed to its neglect. Extensive molecular dynamics simulations are used to validate our approach including a model for calculating the slip length that does not involve a Green-Kubo integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of Zwitterionic Additives on Solvation and Transport of Sodium and Potassium Cations in (Ethylene Oxide)10: A Molecular Dynamics Simulation Study

Sodium- (Na + ) and potassium- (K + ) ion batteries are cost-effective alternatives to lithium-ion (Li + ) batteries due to the abundant sodium and potassium resources. Solid polymer electrolytes (SPEs) are essential for safer and more efficient Na + and K + batteries because they often exhibit low ionic conductivity at room temperature. While zwitterionic (ZW) materials enhance Li+ battery conductivity, their potential for Na + and K + transport in batteries remains unexplored. In this study, we investigated the effect of three ZW molecules (ChoPO4, i.e., 2-methacryloyloxyethyl phosphorylcholine, ImSO3, i.e., sulfobetaine ethylimidazole, and ImCO2, i.e., carboxybetaine ethylimidazole) on the dissociation of Na + and K + coordination with ethylene oxide (EO) chains in EO-based electrolytes through molecular dynamics simulations. Our results showed that ChoPO4 possessed the highest cation–EO10 dissociation ability, while ImSO3 exhibited the lowest. Such dissociation ability correlated with the cation–ZW molecule coordination strength: ChoPO4 and ImSO3 showed the strongest and the weakest coordination with cations. However, the cation–ZW molecule coordination could slow the cationic diffusion. The competition of these effects resulted in accelerating or decelerating cationic diffusion. Our simulated results showed that ImCO2 enhanced Na + diffusion by 20%, while ChoPO4 and ImSO3 led to a 10% reduction. For K + , ChoPO4 reduced its diffusion by 40%, while ImCO2 and ImSO3 caused a similar decrease of 15%. These findings suggest that the ZW structure and the cationic size play an important role in the ionic dissociation effect of ZW materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Advanced analysis of laser-driven pulsed magnetic diffusion based on quantum molecular dynamics simulation

Magnetic diffusion plays an important role in inertial confinement fusion with strong magnetic fields. In this paper, we improve a previous analysis of the generation and diffusion of the magnetic field [Morita et al., Phys. Plasmas 25, 094505 (2018)]. For the generation process, we calculate the temporal evolution of the coil current using a self-consistent circuit model. The results show that the peak of the calculated magnetic field is delayed by 1.2 ns compared with that of the incident laser pulse. For the diffusion process, we evaluate the electrical conductivity of warm dense gold over a wide temperature range (300 K–100 eV) by combining the Kubo–Greenwood formula based on a quantum molecular dynamics simulation with the modified Spitzer model. Our simulation shows that the maximum magnetic field (530 T) that penetrates the cone is delayed by 2.5 ns compared with the laser peak. This result is consistent with experiments [Sakata et al., Nat. Commun. 9, 3937 (2018)] that showed that applying a strong magnetic field improved the heating efficiency of fusion fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Investigating structure and dynamics of unentangled poly(dimethyl- co -diphenyl)siloxane via molecular dynamics simulation

Polysiloxane is one of the most important polymeric materials in technological use. Polydimethylsiloxane displays glass-like mechanical properties at low temperatures. Incorporation of phenyl siloxane, via copolymerization for example, improves not only the low-temperature elasticity but also enhances its performance over a wide range of temperatures. Copolymerization with the phenyl component can significantly change the microscopic properties of polysiloxanes, such as chain dynamics and relaxation. However, despite much work in the literature, the influence of such changes is still not clearly understood. In this work, we systematically study the structure and dynamics of random poly(dimethyl- co -diphenyl)siloxane via atomistic molecular dynamics simulations. As the molar ratio Φ of the diphenyl component increases, we find that the size of the linear copolymer chain expands. At the same time, the chain-diffusivity slows down by over an order of magnitudes. In conclusion, the reduced diffusivity appears to be a result of a complex interplay between the structural and dynamic changes induced by phenyl substitution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗