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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
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↗
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↗
DOE OSTI2025

Quantum dynamics simulation of the advection-diffusion equation

The advection-diffusion equation is simulated via several quantum algorithms. Three formulations are considered: (1) Trotterization, (2) variational quantum time evolution (VarQTE), and (3) adaptive variational quantum dynamics simulation (AVQDS). These schemes were originally developed for the Hamiltonian simulation of many-body quantum systems. The finite-difference discretized operator of the transport equation is formulated as a Hamiltonian and solved without the need for ancillary qubits. Computations are conducted on a quantum simulator (IBM Qiskit Aer) and a superconducting quantum hardware (IBM Fez). The former emulates the latter without the noise. The actual hardware implementation experiences significant noise. The results of the quantum simulator are compared with data from direct numerical simulation (DNS) with infidelities of the order 10 −5 . In the quantum simulator, Trotterization is observed to have the lowest infidelity and is suitable for fault-tolerant computation. The AVQDS algorithm requires the lowest gate count and circuit depth. The VarQTE algorithm is the next best in terms of gate counts, but the number of its optimization variables is directly proportional to the number of qubits. Due to current hardware limitations, Trotterization cannot be implemented, as it has an overwhelmingly large number of operations. Meanwhile, AVQDS and VarQTE can be executed at the hardware level. These algorithms present a new paradigm for computational transport phenomena on quantum computers.

Alipanah, Hirad [Univ. of Pittsburgh, PA (United S↗
DOE OSTI2018

Topological Effects in Isolated Poly[ n ]catenanes: Molecular Dynamics Simulations and Rouse Mode Analysis

Poly[n]catenanes are mechanically interlocked polymers consisting of interlocking ring molecules. We report that over the years, researchers have speculated that the permanent topological interactions within the poly[n]catenane backbone could lead to unique dynamical behaviors. To investigate these unusual polymers, molecular dynamics simulations of isolated poly[n]catenanes have been conducted, along with a Rouse mode analysis. Owing to the mechanical bonds within the molecule, the dynamics of poly[n]catenanes at short length scales are significantly slowed and the distribution of relaxation times is broadened; these same behaviors have been observed in melts of linear polymers and are associated with entanglement. Lastly, despite these entanglement-like effects, at large length scales poly[n]catenanes do not relax much slower than isolated linear polymers and are less strongly impacted by increased segmental stiffness.

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

Uncertainty Quantification in Dynamic Simulations of Large-scale Power System Models using the High-Order Probabilistic Collocation Method on Sparse Grids

This paper proposes a probabilistic collocation method (PCM) to quantify the uncertainties with dynamic simulations in power systems. The appraoch was tested on a single-machine-infinite-bus system and the over 15,000 -bus Western Electricity Coordinating Council (WECC) system. Comparing to classic Monte-Carlo (MC) method, the proposed PCM applies the Smolyak algorithm to reduce the number of simulations that have to be performed. Therefore, the computational cost can be greatly reduced using PCM. The algorithm and procedures are described in the paper. Comparison was made with MC method on the single machine as well as the WECC system. The simulation results shows that using PCM only a small number of sparse grid points need to be sampled even when dealing with systems with a relatively large number of uncertain parameters. PCM is, therefore, computationally more efficient than MC method.

Uncertainty quantification, parameter estimation, ↗
DOE OSTI2015

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

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

Spin dynamics simulation of electron spin relaxation in Ni{sup 2+}(aq)

The ability to quantitatively predict and analyze the rate of electron spin relaxation of open-shell systems is important for electron paramagnetic resonance and paramagnetic nuclear magnetic resonance spectroscopies. We present a combined molecular dynamics (MD), quantum chemistry (QC), and spin dynamics simulation method for calculating such spin relaxation rates. The method is based on the sampling of a MD trajectory by QC calculations, to produce instantaneous parameters of the spin Hamiltonian used, in turn, to numerically solve the Liouville-von Neumann equation for the time evolution of the spin density matrix. We demonstrate the approach by simulating the relaxation of electron spin in an aqueous solution of Ni{sup 2+} ion. The spin-lattice (T{sub 1}) and spin-spin (T{sub 2}) relaxation rates are extracted directly from the simulations of the time dependence of the longitudinal and transverse magnetization, respectively. Good agreement with the available, indirectly obtained experimental data is obtained by our method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗