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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

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

Spin-lattice dynamics simulation of external field effect on magnetic order of ferromagnetic iron

Modeling of field-induced magnetization in ferromagnetic materials has been an active topic in the last dozen years, yet a dynamic treatment of distance-dependent exchange integral has been lacking. In view of that, we employ spin-lattice dynamics (SLD) simulations to study the external field effect on magnetic order of ferromagnetic iron. Our results show that an external field can increase the inflection point of the temperature. Also the model provides a better description of the effect of spin correlation in response to an external field than the mean-field theory. An external field has a more prominent effect on the long range magnetic order than on the short range counterpart. Furthermore, an external field allows the magnon dispersion curves and the uniform precession modes to exhibit magnetic order variation from their temperature dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Overcoming the Time Limitation in Molecular Dynamics Simulation of Crystal Nucleation: A Persistent-Embryo Approach

The crystal nucleation from liquid in most cases is too rare to be accessed within the limited time scales of the conventional molecular dynamics (MD) simulation. Here, we developed a “persistent embryo” method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. In conclusion, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO 2 -16.0B 2 O 3 -12.6Na 2 O-3.8Al 2 O 3 -5.7CaO-1.7ZrO 2 . Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additional short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO 4 ]-[SiO 4 ] and [BO 4 ]-[SiO 4 ] connections are the most abundant while [ZrO 6 ] octohedra having higher probability to connect to another [SiO 4 ], [BO 4 ] and [BO 3 ]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Coarse-grained molecular dynamics simulation of the interface behaviour and self-assembly of CTAB cationic surfactants

Here, we study the behaviour at interfaces and the micelle self-assembly of a cationic surfactant (CTAB) by Molecular Dynamics (MD) simulations of coarse-grained models. We consider both the standard (with explicit water) Martini force field and the implicit solvent version of the Martini force field (Dry Martini). First, we study the behaviour of CTAB at a water/vacuum interface, at a water/organic solvent interface and in a pre-assembled CTAB micelle using both standard and Dry Martini and all-atomic simulations. Our results indicate that there are significant quantitative differences between the predictions of the two models. Interestingly, implicit solvent simulations with Dry Martini show good quantitative agreement with all-atomic MD simulations, better than explicit solvent Martini MD simulations. The computational efficiency of the Martini and Dry Martini models allowed us to study the self-assembly of CTAB in a large system with many micelles. We observe the self-assembly of CTAB into micelles and also the exchange of CTAB molecules between micelles by events such as micelle fusion and fission which are difficult to observe in all-atomic MD simulations due to the time and length scales involved. Under the studied conditions, both Martini models predict a rather different self-assembly behaviour. The standard Martini model predicts a final equilibrium state with spherical micelles with an average size of ≈70 CTAB molecules. In contrast, the Dry Martini model predicts the formation of large tubular micelles with ≈330 CTAB molecules. Compared with experiments, standard Martini and Dry Martini underestimate and overestimate, respectively, the micelle size.

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

A classical molecular dynamics simulation method for the formation of “dry” gels from boro-aluminosilicate glass structures

In contact with water, glass transforms into amorphous and porous structures called gels. A simulation method based on classical molecular dynamics is proposed here to mimic “dry” gels forming from initial oxide glass structures. Six glass compositions were investigated. Two behaviours were evidenced depending on the initial glass composition, and in particular on the quantity of elements removed. If a large quantity of soluble elements (B, Na) was removed, it induced an increase in the average pore size within the gels, and the time needed to stabilise the gel structure increased because more local atomic rearrangements occurred. The gel network displayed a higher proportion of Si-Q4 at the expense of Si-Q3 and a lower average ring size compared to the glass network, irrespective of the glass composition. Surface effects were also highlighted in the dry gels, such as the presence of 3-coordinated Al and a decrease in the average angle Si-O-Si and Al-O-Al. Our findings will be compared to both wet gels and experimental data in further studies, to help find the best procedure to simulate such structures.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Does fluoride disrupt hydrogen bond network in cationic lipid bilayer? Time-dependent fluorescence shift of Laurdan and molecular dynamics simulations

Time-dependent fluorescence shift (TDFS) of Laurdan embedded in phospholipid bilayers reports on hydration and mobility of the phospholipid acylgroups. Exchange of H{sub 2}O with D{sub 2}O prolongs the lifetime of lipid-water and lipid-water-lipid interactions, which is reflected in a significantly slower TDFS kinetics. Combining TDFS measurements in H{sub 2}O and D{sub 2}O hydrated bilayers with atomistic molecular dynamics (MD) simulations provides a unique tool for characterization of the hydrogen bonding at the acylgroup level of lipid bilayers. In this work, we use this approach to study the influence of fluoride anions on the properties of cationic bilayers composed of trimethylammonium-propane (DOTAP). The results obtained for DOTAP are confronted with those for neutral phosphatidylcholine (DOPC) bilayers. Both in DOTAP and DOPC H{sub 2}O/D{sub 2}O exchange prolongs hydrogen-bonding lifetime and does not disturb bilayer structure. These results are confirmed by MD simulations. TDFS experiments show, however, that for DOTAP this effect is cancelled in the presence of fluoride ions. We interpret these results as evidence that strongly hydrated fluoride is able to steal water molecules that bridge lipid carbonyls. Consequently, when attracted to DOTAP bilayer, fluoride disrupts the local hydrogen-bonding network, and the differences in TDFS kinetics between H{sub 2}O and D{sub 2}O hydrated bilayers are no longer observed. A distinct behavior of fluoride is also evidenced by MD simulations, which show different lipid-ion binding for Cl{sup −} and F{sup −}.

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

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Performance and Feature Improvements in Parareal-based Power System Dynamic Simulation

In recent years, a novel Parareal-based approach has been developed for fast transient simulations of large power system interconnections. Parareal belongs to the class of Parallel-in-time algorithms for solution of systems of differential-algebraic equations in parallel over an interval of time. The selection of a reasonably fast and accurate coarse solution is crucial to improve the performance of Parareal algorithm. Semi-analytical solution methods are one promising approach to achieve this goal. They have been investigated, and some preliminary results are presented here. In addition, Parareal-based simulator has been expanded to enable co-simulation with OpenDSS, a widely used open-source distribution system simulator. Preserving the parallel nature of the Parareal approach and taking advantage of the parallel capabilities of the latest versions of OpenDSS, each distribution system can be solved in their entirety on different processors in parallel within the main Parareal simulator. This paper also presents the structure of the transmission and distribution co-simulation and some results with different dynamic models of inverter-based resources in the distribution systems.

Park, Byungkwon↗
DOE OSTI2019

An Overview of Recent Standard and Accelerated Molecular Dynamics Simulations of Helium Behavior in Tungsten

One of the most critical challenges for the successful adoption of nuclear fusion power corresponds to plasma-facing materials. Due to its favorable properties in this context (low sputtering yield, high thermal conductivity, high melting point, among others), tungsten is a leading candidate material. Nevertheless, tungsten is affected by the plasma and fusion byproducts. Irradiation by helium nuclei, in particular, strongly modifies the surface structure by a synergy of processes, whose origin is the nucleation and growth of helium bubbles. In this review, we present recent advances in the understanding of helium effects in tungsten from a simulational approach based on accelerated molecular dynamics, which emphasizes the use of realistic parameters, as are expected in experimental and operational fusion power conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.

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

Quantitative interpretation of molecular dynamics simulations for X-ray photoelectron spectroscopy of aqueous solutions

Over the past decade, energy-dependent ambient pressure X-ray photoelectron spectroscopy (XPS) has emerged as a powerful analytical probe of the ion spatial distributions at the vapor (vacuum)-aqueous electrolyte interface. These experiments are often paired with complementary molecular dynamics (MD) simulations in an attempt to provide a complete description of the liquid interface. There is, however, no systematic protocol that permits a straightforward comparison of the two sets of results. XPS is an integrated technique that averages signals from multiple layers in a solution even at the lowest photoelectron kinetic energies routinely employed, whereas MD simulations provide a microscopic layer-by-layer description of the solution composition near the interface. Here, we use the National Institute of Standards and Technology database for the Simulation of Electron Spectra for Surface Analysis (SESSA) to quantitatively interpret atom-density profiles from MD simulations for XPS signal intensities using sodium and potassium iodide solutions as examples. We show that electron inelastic mean free paths calculated from a semi-empirical formula depend strongly on solution composition, varying by up to 30% between pure water and concentrated NaI. The XPS signal thus arises from different information depths in different solutions for a fixed photoelectron kinetic energy. XPS signal intensities are calculated using SESSA as a function of photoelectron kinetic energy (probe depth) and compared with a widely employed ad hoc method. SESSA simulations illustrate the importance of accounting for elastic-scattering events at low photoelectron kinetic energies (<300 eV) where the ad hoc method systematically underestimates the preferential enhancement of anions over cations. Finally, some technical aspects of applying SESSA to liquid interfaces are discussed.

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

Correlation between time and angular alignment in molecular dynamics simulations of heavy ion collisions

Neutron-proton equilibration is a process which has been used to study the density dependence of the symmetry energy term in the nuclear equation-of-state. Here, this study utilizes constrained molecular dynamics (CoMD) simulations of 70 Zn + 70 Zn with collision energies of 35 and 45 MeV/nucleon. An algorithm is used which searches through CoMD events and identifies the PLF* after it separates from the target and determines its lifetime, Δt. It also determines the fragments that the PLF* breaks apart into and determines their angular alignment. This technique gives an opportunity to explore how the average alignment of dynamically produced fragments, $\langle$α $\rangle$ dyn , evolves with PLF* lifetime. An approximately linear relationship was determined with d$\langle$α $\rangle$ dyn /dΔt = 0.98 ± 0.08 rad / zs and 1.06 ± 0.09 rad / zs for the 35 and 45 MeV / nucleon, respectively, indicating a correlation with magnitude consistent with classically determined values which were used for prior experimental studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

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

Structural and mechanical properties of cardiolipin lipid bilayers determined using neutron spin echo, small angle neutron and X-ray scattering, and molecular dynamics simulations

We report that the detailed structural and mechanical properties of a tetraoleoyl cardiolipin (TOCL) bilayer were determined using neutron spin echo (NSE) spectroscopy, small angle neutron and X-ray scattering (SANS and SAXS, respectively), and molecular dynamics (MD) simulations. We used MD simulations to develop a scattering density profile (SDP) model, which was then utilized to jointly refine SANS and SAXS data. In addition to commonly reported lipid bilayer structural parameters, component distributions were obtained, including the volume probability, electron density and neutron scattering length density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗