Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic 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
DOE OSTI2025

Combined molecular and spin dynamics simulation of BCC iron with vacancy defects

Utilizing an atomistic computational model, which handles both translational and spin degrees of freedom, combined molecular and spin dynamics simulations have been performed to investigate the effect of vacancy defects on spin wave excitations in ferromagnetic iron. Fourier transforms of space- and time-displaced correlation functions yield the dynamic structure factor, providing characteristic frequencies and lifetimes of the spin wave modes. A comparison of the system with a 5% vacancy concentration with pure lattice data shows a decrease in frequency and a decrease in lifetime for all transverse spin wave excitations observed. In addition, the clearly defined transverse spin wave excitations are distorted with the introduction of vacancy defects, and we observe reduced excitation lifetimes due to increased magnon–magnon scattering. We observe further evidence of increased magnon–magnon scattering, as the peaks in the longitudinal spin wave spectrum become less distinct. Finally, similar impacts are observed in the vibrational subsystem, with a decrease in characteristic phonon frequency and flattening of lattice excitation signals due to vacancy defects.

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

Insights into the distinct lithiation/sodiation of porous cobalt oxide by in operando synchrotron x-ray techniques and ab initio molecular dynamics simulations

Sodium-ion batteries (SIBs) have been considered as one of the promising power source candidates for the stationary storage industries owing to the much lower cost of sodium than lithium. It is well-known that the electrode materials largely determine the energy density of the battery systems. However, recent discoveries on the electrode materials showed that most of them present distinct lithium and sodium storage performance, which is not yet well understood. In this work, we performed a comparative understanding on the structural changes of porous cobalt oxide during its electrochemical lithiation and sodiation process by in operando synchrotron small angel X-ray scattering, X-ray diffraction, and X-ray absorption spectroscopy. It was found that compared to the lithiation process, the porous cobalt oxide undergoes less pore structure changes, oxidation state, and local structure changes as well as crystal structure evolution during its sodiation process, which is attributed to the intrinsic low sodiation activity of cobalt oxide as evidenced by ab initio molecular dynamics simulations. Moreover, it was indicated that the sodiation activity of metal sulfides is higher than that of metal oxides, indicating a better candidate for SIBs. Furthermore, such understanding is crucial for future design and improvement of high-performance electrode materials for SIBs.

25 ENERGY STORAGE↗
DOE OSTI2016

Molecular Dynamics Simulations Reveal the Conformational Flexibility of Lipid II and Its Loose Association with the Defensin Plectasin in the Staphylococcus aureus Membrane

Lipid II is critical for peptidoglycan synthesis, which is the main component of the bacterial cell wall. Lipid II is a relatively conserved and important part of the cell wall biosynthesis pathway and is targeted by antibiotics such as the lantibiotics, which achieve their function by disrupting the biosynthesis of the cell wall. Given the urgent need for development of novel antibiotics to counter the growing threat of bacterial infection resistance, it is imperative that a thorough molecular-level characterization of the molecules targeted by antibiotics be achieved. To this end, we present here a molecular dynamics simulation study of the conformational dynamics of Lipid II within a detailed model of the Staphylococcus aureus cell membrane. We show that Lipid II is able to adopt a range of conformations, even within the packed lipidic environment of the membrane. Our simulations also reveal dimerization of Lipid II mediated by cations. In the presence of the defensin peptide plectasin, the conformational lability of Lipid II allows it to form loose complexes with the protein, via a number of different binding modes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Parallel replica dynamics simulations of reactions in shock compressed liquid benzene

The study of the long-term evolution of slow chemical reactions is challenging because quantum-based reactive molecular dynamics simulation times are typically limited to hundreds of picoseconds. In this work, the extended Lagrangian Born-Oppenheimer molecular dynamics formalism is used in conjunction with parallel replica dynamics to obtain an accurate tool to describe the long-term chemical dynamics of shock-compressed benzene. Langevin dynamics has been employed at different temperatures to calculate the first reaction times in liquid benzene at pressures and temperatures consistent with its unreacted Hugoniot. Our coupled engine runs for times on the order of nanoseconds (one to two orders of magnitude longer than traditional techniques) and is capable of detecting reactions that are characterized by rates significantly slower than we could study before. At lower pressures and temperatures, we mainly observe Diels-Alder metastable reactions, whereas at higher pressures and temperatures we observe stable polymerization reactions.

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

Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations

We investigate helium flux effects on helium transport and surface evolution in plasma-facing tungsten using molecular dynamics. The simulations span two orders of magnitude, from ITER-relevant levels to those more typical of simulations published to date. Simulation times of up to 2.5 µs (corresponding to actual fluences of 1.4 x 10 21 m -2 ) are achieved, revealing concerted bubble-bursting events that are responsible for significant and very sudden changes in surface morphology. The depth distribution of helium depends very strongly on helium flux: helium self-trapping becomes more probable near the surface at high flux, and a layer of near-surface bubbles forms. Helium retention prior to the onset of bubble bursting is also substantially lower at low flux than it is at high flux. Surface features at low fluence are correlated with the positions of bubbles, but at high fluence, bubbles tend to coalesce, venting to the surface at one or more locations and leaving large interconnected cavities below the surface. Ruptured bubbles may serve as pathways deeper into the material, allowing helium to bypass the layer of near-surface bubbles and fill deeper, potentially much larger, bubbles that can produce more substantial surface features. Deeper bubbles also emit prismatic dislocation loops that can fill in cavities closer to the surface. Our results suggest that nearly all molecular dynamics simulations published to date are hampered by finite-size effects, and that helium flux is a very important parameter in determining the behavior of helium in plasma-facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Molecular dynamics simulations of the structure and single-particle dynamics of mixtures of divalent salts and ionic liquids

We report a molecular dynamics study of the structure and single-particle dynamics of mixtures of a protic (ethylammonium nitrate) and an aprotic (1-butyl-3-methylimidazolium hexaflurophosphate [BMIM][PF{sub 6}]) room-temperature ionic liquids doped with magnesium and calcium salts with a common anion at 298.15 K and 1 atm. The solvation of these divalent cations in dense ionic environments is analyzed by means of apparent molar volumes of the mixtures, radial distribution functions, and coordination numbers. For the protic mixtures, the effect of salt concentration on the network of hydrogen bonds is also considered. Moreover, single-particle dynamics of the salt cations is studied by means of their velocity autocorrelation functions and vibrational densities of states, explicitly analyzing the influence of salt concentration, and cation charge and mass on these magnitudes. The effect of the valency of the salt cation on these properties is considered comparing the results with those for the corresponding mixtures with lithium salts. We found that the main structural and dynamic features of the local solvation of divalent cations in ionic liquids are similar to those of monovalent salts, with cations being localized in the polar nanoregions of the bulk mixture coordinated in monodentate and bidentate coordination modes by the [NO{sub 3}]{sup −} and [PF{sub 6}]{sup −} anions. However, stronger electrostatic correlations of these polar nanoregions than in mixtures with salts with monovalent cations are found. The vibrational modes of the ionic liquid (IL) are seen to be scarcely affected by the addition of the salt, and the effect of mass and charge on the vibrational densities of states of the dissolved cations is reported. Cation mass is seen to exert a deeper influence than charge on the low-frequency vibrational spectra, giving a red shift of the vibrational modes and a virtual suppression of the higher energy vibrational modes for the heavier Ca{sup 2+} cations. No qualitative difference with monovalent cations was found in what solvation is concerned, which suggests that no enhanced reduction of the mobility of these cations and their complexes in ILs respective to those of monovalent cations is to be expected.

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

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2021

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2019

Transfer of an erbium ion across the water/dodecane interface: Structure and thermodynamics via molecular dynamics simulations

The thermodynamics and structural changes involved in the transfer of the Er 3+ ion across the water/dodecane interface are investigated by molecular dynamics simulations. We show that the Er 3+ ion is transferred as a highly conserved 8-water coordinated species and that the transfer involves significant perturbation of the interfacial water structure. Furthermore, several structural properties are used to quantify this process. Implications for the ion extraction process are discussed.

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

Distance Restraints from Crosslinking Mass Spectrometry: Mining a Molecular Dynamics Simulation Database to Evaluate Lysine-Lysine Distances

Integrative structural biology models the structures of protein complexes that are intractable by classical structural methods (because of extreme size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such method is chemical cross-linking mass spectrometry (XL-MS), in which cross-linked peptides, derived from a covalently cross-linked protein complex and identified by liquid chromatography-mass spectrometry, pinpoint protein residues close in three-dimensional space. The commonly used lysine-reactive N-hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS3) have a linker arm that is 11.4 Å long when fully extended. However, XL-MS studies on proteins of known structure frequently report cross-links that exceed this distance. Typically, a tolerance of ~3 Å is added to the theoretical maximum to account for this observation, with little justification for the value chosen. We used the Dynameomics database, a repository of high-quality molecular dynamics simulations of 807 proteins representative of all protein folds, to investigate the change in lysine-lysine distances resulting from native-state dynamics on the time-scale of tens of nanoseconds. We conclude that observed cross-links are consistent with a protein structure if the distance between cross-linked lysine Nζ atoms is less than the cross-linker length plus 11.3 Å. For DSS or BS3, this corresponds to a Cα to Cα distance of 30.4 Å. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL-MS results to structures, and indicates the appropriate values of an XLMS derived distance constraint to use in structural modeling.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2023

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein–solvent interface. To address this limitation, a set of MD–MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Direct comparison of elastic incoherent neutron scattering experiments with molecular dynamics simulations of DMPC phase transitions

Neutron scattering techniques have been employed to investigate 1,2-dimyristoyl-sn-glycero-3- phosphocholine (DMPC) membranes in the form of multilamellar vesicles (MLVs) and deposited, stacked multilamellar-bilayers (MLBs), covering transitions from the gel to the liquid phase. Neutron diffraction was used to characterise the samples in terms of transition temperatures, whereas elastic incoherent neutron scattering (EINS) demonstrates that the dynamics on the sub-macromolecular length-scale and pico- to nano-second time-scale are correlated with the structural transitions through a discontinuity in the observed elastic intensities and the derived mean square displacements. Molecular dynamics simulations have been performed in parallel focussing on the length-, time- and temperature-scales of the neutron experiments. They correctly reproduce the structural features of the main gel-liquid phase transition. Particular emphasis is placed on the dynamical amplitudes derived from experiment and simulations. Two methods are used to analyse the experimental data and mean square displacements. They agree within a factor of 2 irrespective of the probed time-scale, i.e. the instrument utilized. Mean square displacements computed from simulations show a comparable level of agreement with the experimental values, albeit, the best match with the two methods varies for the two instruments. Consequently, experiments and simulations together give a consistent picture of the structural and dynamical aspects of the main lipid transition and provide a basis for future, theoretical modelling of dynamics and phase behaviour in membranes. The need for more detailed analytical models is pointed out by the remaining variation of the dynamical amplitudes derived in two different ways from experiments on the one hand and simulations on the other.

Aoun, Bachir↗
DOE OSTI2017

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes were used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Molecular dynamics simulation of seeded crystal growth in glass

The mechanism of non-congruent growth of a crystal from glass has been sought using molecular dynamics simulations. Specifically, as a model of this process, the growth of a lithium niobate (LiNbO 3 ) crystal seed sandwiched between two lithium niobosilicate (LNS) glass slabs has been simulated as a function of time and temperature. The growth of pre-existing crystal is strongly affected by the orientation of crystal seed, temperature, and the SiO 2 concentration in the surrounding LNS glass matrix. The orientation of LiNbO 3 seed surface that has inherently larger interplanar distance results in a relatively slower crystal growth. The addition of SiO 2 to LNS system significantly decreases the crystal growth, which primarily occurs in the region devoid of Si. The suppressive effect of SiO 2 on growth rate can be traced to the existence of defect complex comprising of Si substituted at the Nb site and a nearby Nb vacancy.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The FADE mass-stat: A technique for inserting or deleting particles in molecular dynamics simulations

The emergence of new applications of molecular dynamics (MD) simulation calls for the development of mass-statting procedures that insert or delete particles on-the-fly. In this paper we present a new mass-stat which we term FADE, because it gradually “fades-in” (inserts) or “fades-out” (deletes) molecules over a short relaxation period within a MD simulation. FADE applies a time-weighted relaxation to the intermolecular pair forces between the inserting/deleting molecule and any neighbouring molecules. The weighting function we propose in this paper is a piece-wise polynomial that can be described entirely by two parameters: the relaxation time scale and the order of the polynomial. FADE inherently conserves overall system momentum independent of the form of the weighting function. We demonstrate various simulations of insertions of atomic argon, polyatomic TIP4P water, polymer strands, and C{sub 60} Buckminsterfullerene molecules. We propose FADE parameters and a maximum density variation per insertion-instance that restricts spurious potential energy changes entering the system within desired tolerances. We also demonstrate in this paper that FADE compares very well to an existing insertion algorithm called USHER, in terms of accuracy, insertion rate (in dense fluids), and computational efficiency. The USHER algorithm is applicable to monatomic and water molecules only, but we demonstrate that FADE can be generally applied to various forms and sizes of molecules, such as polymeric molecules of long aspect ratio, and spherical carbon fullerenes with hollow interiors.

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

Molecular dynamics simulation of a binary mixture near the lower critical point

2,6-lutidine molecules mix with water at high and low temperatures but in a wide intermediate temperature range a 2,6-lutidine/water mixture exhibits a miscibility gap. We constructed and validated an atomistic model for 2,6-lutidine and performed molecular dynamics simulations of 2,6-lutidine/water mixture at different temperatures. We determined the part of demixing curve with the lower critical point. The lower critical point extracted from our data is located close to the experimental one. The estimates for critical exponents obtained from our simulations are in a good agreement with the values corresponding to the 3D Ising universality class.

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

Evaporative cooling of microscopic water droplets in vacuo : Molecular dynamics simulations and kinetic gas theory

In the present study, we investigate the process of evaporative cooling of nanometer-sized droplets in vacuum using molecular dynamics simulations with the TIP4P/2005 water model. The results are compared to the temperature evolution calculated from the Knudsen theory of evaporation which is derived from kinetic gas theory. The calculated and simulation results are found to be in very good agreement for an evaporation coefficient equal to unity. Lastly, our results are of interest to experiments utilizing droplet dispensers as well as to cloud micro-physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗