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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Effect of Strain Field on Threshold Displacement Energy of Tungsten Studied by Molecular Dynamics Simulation

The influence of strain field on defect formation energy and threshold displacement energy (Ed) in body-centered cubic (BCC) tungsten (W) has been studied with molecular dynamics simulations. Two different W potentials (Fikar and Juslin) were compared and the results indicate that the connection distance and selected function linking the short-range and long-range portions of the potentials affects the threshold displacement energy and its direction-specific values. The minimum Ed direction calculated with the Fikar-potential is <100> and with the Juslin-potential is <111>. Nevertheless, the most stable self-interstitial configuration is found to be a <111>-crowdion for both potentials. This stable configuration does not change with applied strain. Varying the strain from compression to tension increases the vacancy formation energy but decreases the self-interstitial formation energy. The formation energy of a self-interstitial changes more significantly than a vacancy such that Ed decreases with applied hydrostatic strain from compression to tension.

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

Modeling of Grid-Forming and Grid-Following Inverters for Dynamic Simulation of Large-Scale Distribution Systems

Historically, distribution system planning studies mainly focused on steady state and quasi-steady state analysis, with limited attention paid to dynamic analysis. This paper develops three-phase, electromechanical models for both grid-forming and grid-following inverters, and integrates them into an open source, three-phase distribution network solver, thereby enabling dynamic simulation of large-scale, three-phase unbalanced distribution systems with high penetration of inverter-based DERs. The proposed inverter models are validated against electromagnetic simulation and field test data from the CERTS/AEP microgrid testbed, and simulated in an islanded 5252 node distribution system in the GridLAB-D simulation environment. Simulation verifies the effectiveness of the proposed inverter models for large-scale distribution system analysis. Finally, study results show that compared to traditional grid-following inverters, the high penetration of grid-forming inverters can improve the voltage and frequency stability of islanded distribution systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
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 OSTI2017

Different dynamics and pathway of disulfide bonds reduction of two human defensins, a molecular dynamics simulation study

Human defensins are a class of antimicrobial peptides that are crucial components of the innate immune system. Both human α defensin type 5 (HD5) and human β defensin type 3 (hBD-3) have 6 cysteine residues which form 3 pairs of disulfide bonds in oxidizing condition. Disulfide bond linking is important to the protein structure stabilization, and the disulfide bond linking and breaking order have been shown to influence protein function. In this project, microsecond long molecular dynamics simulations were performed to study the structure and dynamics of HD5 and hBD-3 wildtype and analogs which have all 3 disulfide bonds released in reducing condition. The structure of hBD-3 was found to be more dynamic and flexible than HD5, based on RMSD, RMSF, and radius of gyration calculations. The disulfide bridge breaking order of HD5 and hBD-3 in reducing condition was predicted by two kinds of methods, which gave consistent results. It was found that the disulfide bonds breaking pathways for HD5 and hBD-3 are very different. The breaking of disulfide bonds can influence the dimer interface by making the dimer structure less stable for both kinds of defensin. In order to understand the difference in dynamics and disulfide bond breaking pathway, hydrophilic and hydrophobic accessible surface areas (ASA), buried surface area between cysteine pairs, entropy of cysteine pairs, and internal energy were calculated. Comparing to the wildtype, hBD-3 analog is more hydrophobic, while HD5 is more hydrophilic. For hBD-3, the disulfide breaking is mainly entropy driven, while other factors such as the solvation effects may take the major role in controlling HD5 disulfide breaking pathway. Proteins 2017; 85:665–681. © 2016 Wiley Periodicals, Inc.

Biochemistry & Molecular Biology↗
DOE OSTI2016

Molecular Dynamics Simulation Study of Solvent and State of Charge Effects on Solid-Phase Structure and Counterion Binding in a Nitroxide Radical Containing Polymer Energy Storage Material

Here we performed molecular dynamics simulations to understand the effects of solvent swelling and state of charge (SOC) on the redox active, organic radical cathode material poly(2,2,6,6-tetramethylpiperidinyloxy methacrylate) (PTMA). We show that the polar solvent acetonitrile primarily solvates the nitroxide radical without disrupting the packing of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) pendant groups of PTMA. We also simulated bulk PTMA in different SOC, 25%, 50%, 75%, and 100%, by converting the appropriate number of TEMPO groups to the cation charge state and adding BF 4 - counterions to the simulation. At each SOC the packing of PTMA, the solvent, and the counterions were examined. The binding of the anion to the nitroxide cation site was examined using the potential of mean force and found to be on the order of tens of meV, with a binding energy that decreased with increasing SOC. Additionally, we found that the cation state is stabilized by the presence of a nearby anion by more than 1 eV, and the implications of this stabilization on charge transport are discussed. Finally, we describe the implications of our results for how the SOC of an organic electrode affects electron and anion charge transport during the charging and discharging processes.

36 MATERIALS SCIENCE↗
DOE OSTI2014

All-atom molecular dynamics simulation of a photosystem I/detergent complex

All-atom molecular dynamics (MD) simulation was used to investigate the solution structure and dynamics of the photosynthetic pigment protein complex photosystem I (PSI) from Thermosynechococcus elongatus embedded in a toroidal belt of n -dodecyl-β-d-maltoside (DDM) detergent. Evaluation of root-mean-square deviations (RMSDs) relative to the known crystal structure show that the protein complex surrounded by DDM molecules is stable during the 200 ns simulation time, and root-mean-square fluctuation (RMSF) analysis indicates that regions of high local mobility correspond to solvent-exposed regions such as turns in the transmembrane α-helices and flexible loops on the stromal and lumenal faces. Comparing the protein detergent complex to a pure detergent micelle, the detergent surrounding the PSI trimer is found to be less densely packed but with more ordered detergent tails, contrary to what is seen in most lipid bilayer models. We also investigated any functional implications for the observed conformational dynamics and protein detergent interactions, discovering interesting structural changes in the psaL subunits associated with maintaining the trimeric structure of the protein. Moreover, we find that the docking of soluble electron mediators such as cytochrome c 6 and ferredoxin to PSI is not significantly impacted by the solubilization of PSI in detergent.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Ab initio molecular dynamics simulation of supercritical aqueous ionic solutions: Spectral diffusion of water in the vicinity of Br¯ and I¯ ions

Supercritical aqueous ionic solutions containing bromide (Br¯) and iodide (I¯) ions have been studied using Car-Parrinello molecular dynamics simulation at two different concentrations and at two different densities. The time dependent trajectory has been used for calculating vibrational spectral diffusion, structural properties and, for other dynamical properties such as hydrogen bond dynamics, orientational relaxation and residence dynamics. The vibrational spectral diffusion of fluctuating OD stretch mode inside the solvation shell of the ions and in the bulk has been studied using frequency-frequency time correlation function. Hydrogen bond dynamics, orientational relaxation and residence dynamics of Sc-water have been investigated using population correlation function approach. The water molecules inside the solvation shell of Br¯ and I¯ ions are found to show complex molecular dynamics in comparison to the water-water dynamics. The concentration dependent spectral diffusion has been observed and it occurs in multiple time scales in the solvation shell of Br¯ and I¯ ions. Here, at supercritical condition, the ion-water hydrogen bond shows faster dynamics with respect to the same at ambient condition and slower dynamics with respect to the water-water hydrogen bond dynamics at supercritical condition.

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↗
NASA NTRS1996

Electromechanical Dynamics Simulations of Superconducting LSM Rocket Launcher System in Attractive-Mode

Further feasibility study on a superconducting linear synchronous motor (LSM) rocket launcher system is presented on the basis of dynamic simulations of electric power, efficiency and power factor as well as the ascending motions of the launcher and rocket. The advantages of attractive-mode operation are found from comparison with repulsive-mode operation. It is made clear that the LSM rocket launcher system, of which the long-stator is divided optimally into 60 sections according to launcher speeds, can obtain high efficiency and power factor.

Yoshida, Kinjiro↗
DOE OSTI2024

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

The Dynamic Conformational Cycle of the Group I Chaperonin C-Termini Revealed via Molecular Dynamics Simulation

Chaperonins are large ring shaped oligomers that facilitate protein folding by encapsulation within a central cavity. All chaperonins possess flexible C-termini which protrude from the equatorial domain of each subunit into the central cavity. Biochemical evidence suggests that the termini play an important role in the allosteric regulation of the ATPase cycle, in substrate folding and in complex assembly and stability. Despite the tremendous wealth of structural data available for numerous orthologous chaperonins, little structural information is available regarding the residues within the C-terminus. Herein, molecular dynamics simulations are presented which localize the termini throughout the nucleotide cycle of the group I chaperonin, GroE, from Escherichia coli. The simulation results predict that the termini undergo a heretofore unappreciated conformational cycle which is coupled to the nucleotide state of the enzyme. As such, these results have profound implications for the mechanism by which GroE utilizes nucleotide and folds client proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Molecular dynamics simulation of zirconium tungstate amorphization and the amorphous-crystalline interface

Classical molecular dynamics (MD) simulations were performed to provide a conceptual understanding of the amorphous-crystalline interface for a candidate negative thermal expansion (NTE) material, ZrW 2 O 8 . Simulations of pressure-induced amorphization at 300 K indicate that an amorphous phase forms at pressures of 10 GPa and greater, and this phase persists when the pressure is subsequently decreased to 1 bar. However, the crystalline phase is recovered when the slightly distorted 5 GPa phase is relaxed to 1 bar. Simulations were also performed on a two-phase model consisting of the high-pressure amorphous phase in direct contact with the crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the crystalline-amorphous interface. In conclusion, the differences in local atomic structure at the interface are quantified from the simulation trajectories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Molecular Dynamics Simulation of the Migration of Tilt Grain Boundaries in Ni and Ni{sub 3}Al

The migration of the 〈100〉 and 〈111〉 tilt boundaries in Ni and the Ni{sub 3}Al intermetallic compound is studied by molecular dynamics simulation. The low-angle 〈100〉 boundaries in Ni and Ni{sub 3}Al are shown to migrate much more slowly (by approximately two times at a temperature of 1700 K) that the 〈111〉 boundaries, which is related to different migration mechanisms of the 〈100〉 and 〈111〉 boundaries. The migration of the 〈100〉 boundaries is found to occur by the splitting of paired grain-boundary dislocations followed by a change of partner dislocations, which proceeds via the slip of split dislocations. The 〈111〉 tilt boundaries move according to a combined action of the following two mechanisms: the mechanism described above and the mechanism consisting in joint slip of paired grain-boundary dislocations, which have common slip planes in contrast to the grain-boundary dislocations in 〈100〉 boundaries. Since the joint slip of paired dislocations has a relatively low activation energy, the mobility of the 〈111〉 tilt boundaries is higher than that of the 〈100〉 boundaries. Our molecular dynamics studies demonstrate that the migration rate of analogous boundaries in Ni{sub 3}Al is significantly lower than in Ni (by approximately three times at a temperature of 1700 K), which is partly caused by the additional energy consumed for the formation of a disordered region behind a migrating boundary in Ni{sub 3}Al. Because of a relatively low boundary mobility in Ni{sub 3}Al, the contribution of diffusion atomic displacements during boundary migration is higher as compared to Ni.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 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 OSTI2016

Overcoming potential energy distortions in constrained internal coordinate molecular dynamics simulations

The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗