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 OSTI2014

Molecular dynamics simulation studies of the wild type and E92Q/N155H mutant of Elvitegravir-resistance HIV-1 integrase

Although Elvitegravir (EVG) is a newly developed antiretrovirals drug to treat the acquired immunodeficiency syndrome (AIDS), drug resistance has already been found in clinic, such as E92Q/N155H and Q148H/G140S. Several structural investigations have already been reported to reveal the molecular mechanism of the drug resistance. As full length crystal structure for HIV-1 integrase is still unsolved, we use in this paper the crystal structure of the full length prototype foamy virus (PFV) in complex with virus DNA and inhibitor Elvitegravir as a template to construct the wild type and E92Q/N155H mutant system of HIV-1 integrase. Molecular dynamic simulations was used to revel the binding mode and the drug resistance of the EVG ligand in E92Q/N155H. Several important interactions were discovered between the mutated residues and the residues in the active site of the E92Q/N155H double mutant pattern, and cross correlation and clustering methods were used for detailed analysis. The results from the MD simulation studies will be used to guide the experimental efforts of developing novel inhibitors against drug-resistant HIV integrase mutants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Alumina(0001)/Water Interface: Structural Properties and Infrared Spectra from First-Principles Molecular Dynamics Simulations

We investigated the atomic structure and infrared spectra of the alumina(0001)/water interface, using first-principles molecular dynamics simulations based on density functional theory within the generalized gradient approximation. The computed structural properties of the interface are in good agreement with the results of synchrotron X-ray experiments. In this work detailed analyses of the computed infrared spectra revealed two types of water molecules at the hydrophilic oxide/water interface: molecules participating in strong “ice-like” hydrogen bonding with the oxide surface and molecules involved in weaker “liquid-like” hydrogen bonding. Our results provide a molecular interpretation of the “ice-like” and “liquid-like” bands observed in sum-frequency vibrational spectroscopy experiments and underscore the significance of strong hydrogen-bonding interactions in determining the orientation of interfacial water molecules.

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

Spectral shapes of Ar-broadened HCl lines in the fundamental band by classical molecular dynamics simulations and comparison with experiments

Spectral shapes of isolated lines of HCl perturbed by Ar are investigated for the first time using classical molecular dynamics simulations (CMDS). Using reliable intermolecular potentials taken from the literature, these CMDS provide the time evolution of the auto-correlation function of the dipole moment, whose Fourier-Laplace transform leads to the absorption spectrum. In order to test these calculations, room temperature spectra of various lines in the fundamental band of HCl diluted in Ar are measured, in a large pressure range, with a difference-frequency laser spectrometer. Comparisons between measured and calculated spectra show that the CMDS are able to predict the large Dicke narrowing effect on the shape of HCl lines and to satisfactorily reproduce the shapes of HCl spectra at different pressures and for various rotational quantum numbers.

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

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

The Potential of Computational Fluid Dynamics Simulation on Serial Monitoring of Hemodynamic Change in Type B Aortic Dissection

PurposeWe aimed to assess the potential of computational fluid dynamics simulation (CFD) in detecting changes in pressure and flow velocity in response to morphological changes in type B aortic dissection.Materials and MethodsPressure and velocity in four morphological models of type B aortic dissection before and after closure of the entry tear were calculated with CFD and analyzed for changes among the different scenarios. The control model (Model 1) was patient specific and built from the DICOM data of CTA, which bore one entry tear and three re-entry tears. Models 2–4 were modifications of Model 1, with two re-entry tears less in Model 2, one re-entry tear more in Model 3, and a larger entry tear in Model 4.ResultsThe pressure and velocity pertaining to each of the morphological models were unique. Changes in pressure and velocity findings were accountable by the changes in morphological features of the different models. There was no blood flow in the false lumen across the entry tear after its closure, the blood flow direction across the re-entry tears was reversed after closure of the entry tear.ConclusionCFD simulation is probably useful to detect hemodynamic changes in the true and false lumens of type B aortic dissection in response to morphological changes, it may potentially be developed into a non-invasive and patient-specific tool for serial monitoring of hemodynamic changes of type B aortic dissection before and after treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Modelling the local atomic structure of molybdenum in nuclear waste glasses with ab initio molecular dynamics simulations

The nature of chemical bonding of molybdenum in high level nuclear waste glasses has been elucidated by ab initio molecular dynamics simulations. Two compositions, (SiO2)57.5 – (B2O3)10 – (Na2O)15 – (CaO)15 – (MoO3)2.5 and (SiO2)57.3 – (B2O3)20 – (Na2O)6.8 – (Li2O)13.4 – (MoO3)2.5 , were considered in order to investigate the effect of ionic and covalent components on the glass structure and the formation of the crystallisation precursors (Na2MoO4 and CaMoO4). The coordination environments of Mo cations and the corresponding bond lengths calculated from our model are in excellent agreement with experimental observations. The analysis of the first coordination shell reveals two different types of molybdenum host matrix bonds in the lithium sodium borosilicate glass. Based on the structural data and the bond valence model, we demonstrate that the Mo cation can be found in a redox state and the molybdate tetrahedron can be connected with the borosilicate network in a way that inhibits the formation of crystalline molybdates. These results significantly extend our understanding of bonding in Mo-containing nuclear waste glasses and demonstrate that tailoring the glass composition to specific heavy metal constituents can facilitate incorporation of heavy metals at high concentrations. K.K. was supported through the Impact Studentship scheme at UCL co-funded by the IHI Corporation and UCL. P.V.S. thanks the Royal Society, which supported preliminary work on this project, and the Laboratory Directed Research and Development program at PNNL, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy. Via our membership of the UK's HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202), this work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

Laser Tracker Utilization Methodology in Measuring Truth Trajectories for INS Testing on 6 Degree of Freedom Table at the Marshall Space Flight Center's Contact Dynamics Simulation Laboratory with Lessons Learned

When performing Inertial Navigation System (INS) testing at the Marshall Space Flight Center's (MSFC) Contact Dynamics Simulation Laboratory (CDSL) early in 2017, a Leica Geosystems AT901 Laser Tracker system (LLT) measured the twist & sway trajectories as generated by the 6 Degree Of Freedom (6DOF) Table in the CDSL. These LLT measured trajectories were used in the INS software model validation effort. Several challenges were identified and overcome during the preparation for the INS testing, as well as numerous lessons learned. These challenges included determining the position and attitude of the LLT with respect to an INS-shared coordinate frame using surveyed monument locations in the CDSL and the accompanying mathematical transformation, accurately measuring the spatial relationship between the INS and a 6DOF tracking probe due to lack of INS visibility from the LLT location, obtaining the data from the LLT during a test, determining how to process the results for comparison with INS data in time and frequency domains, and using a sensitivity analysis of the results to verify the quality of the results. While many of these challenges were identified and overcome before or during testing, a significant lesson on test set-up was not learned until later in the data analysis process. It was found that a combination of trajectory-dependent gimbal locking and environmental noise introduced non-negligible noise in the angular measurements of the LLT that spanned the evaluated frequency spectrum. The lessons learned in this experiment may be useful for others performing INS testing in similar testing facilities.

Leggett, Jared O.↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

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

Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge–Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call “monomer PIMD”, correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Dynamic Simulation Probabalistic Risk Assessment Model for an Enceladus Sample Return Mission

Enceladus, a moon of Saturn, has geyser-like jets that spray plumes of material into orbit. These jets could enable a free-flying spacecraft to collect samples and return them to Earth for study to determine if they contain the building blocks of life. The Office of Planetary Protection at NASA requires containment of any unsterilized samples and prohibits destructive impact of the spacecraft upon return to Earth, with a sample release probability of less than 1 in 1,000,000 as a recommended goal. This paper describes a probabilistic risk assessment model that uses dynamic simulation techniques to capture the physics-based, time- and state-dependent interactions between the sample return system and the environment, which drive the risk of sample release. The dynamic approach uses a Monte Carlo-style simulation to integrate the many phases and sources of risk for a sample return mission. The model is used to assess the achievability of the planetary protection reliability goal. This is accomplished by performing sensitivity studies assessing the impact of modeling assumptions to identify where uncertainties drive the risk. These results, in turn, are used to examine the feasibility of meeting key design and performance parameters that are needed to achieve the reliability goal for a given architecture with existing technologies.

Dynamic Simulation↗