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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Conformational dynamics of a crystalline protein from microsecond-scale molecular dynamics simulations and diffuse X-ray scattering

Significance A major challenge of protein crystallography is to accurately describe protein structure variations. Whereas Bragg peaks only yield a picture in which the variations are superimposed, diffuse X-ray scattering (diffraction between the peaks) reports on the way different atoms move together and can be used to increase the accuracy of models. We have performed a molecular-dynamics (MD) simulation of crystalline staphylococcal nuclease, yielding a rich picture of functionally important motions. This simulation is 100-fold longer than previous studies and yields reproducible calculations of diffuse intensity, overcoming a major challenge in the field. Experimental data support the MD models and indicate that diffuse scattering can be used both to validate and to improve MD simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Molecular Dynamics Simulation and Experimental Studies of Gold Nanoparticle Templated HDL-like Nanoparticles for Cholesterol Metabolism Therapeutics

High-density lipoprotein (HDL) plays an important role in the transport and metabolism of cholesterol. Mimics of HDL are being explored as potentially powerful therapeutic agents for removing excess cholesterol from arterial plaques. Gold nanoparticles (AuNPs) functionalized with apolipoprotein A-I and with the lipids 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine-N-[3-(2-pyridyldithio)propionate] have been demonstrated to be robust acceptors of cellular cholesterol. However, detailed structural information about this functionalized HDL AuNP is still lacking. Here in this study, we have used X-ray photoelectron spectroscopy and lecithin/cholesterol acyltransferase activation experiments together with coarse-grained and all-atom molecular dynamics simulations to model the structure and cholesterol uptake properties of the HDL AuNP construct. By simulating different apolipoprotein-loaded AuNPs, we find that lipids are oriented differently in regions with and without apoA-I. We also show that in this functionalized HDL AuNP, the distribution of cholesteryl ester maintains a reverse concentration gradient that is similar to the gradient found in native HDL.

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

Effects of polarizability on the structural and thermodynamics properties of [C{sub n}mim][Gly] ionic liquids (n = 1–4) using EEM/MM molecular dynamic simulations

An extended electronegativity equalization method/molecular mechanics (EEM/MM) model for ionic liquids is used to investigate the structures and properties of 1-alkyl-3-methylimidazolium glycine ionic liquids [C{sub n}mim][Gly] (n = 1–4) with alkyl substituents of different lengths. The EEM/MM model describes the electrostatic interactions of atoms and their changes in different ambient environments. This property is the most outstanding characteristic of the model. EEM parameters (i.e., valence electronegativities and valence hardness parameters) are calibrated using linear regression and least-squares methods, which can accurately predict the gas-phase properties of [C{sub n}mim]{sup +}, [Gly]{sup −}, and [C{sub n}mim][Gly] ion pairs. We utilize the EEM/MM force field to systematically investigate the effects of polarizability on the accuracy of [C{sub n}mim][Gly] properties predicted through the molecular dynamic simulations. EEM/MM explicitly describes the atom-based polarizability of [C{sub n}mim][Gly]; thus, the densities, enthalpies of vaporization, self-diffusion coefficients, and conductivities of the [C{sub n}mim][Gly] are consistent with the experimental values. The calculated radial distribution functions provide a mechanistic understanding of the effects of polarizability on ionic aggregations in amino acid ionic liquids. The effects of alkyl chain length on the diffusion coefficient and conductivity are also discussed.

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

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

EXPLORING THE STRUCTURE OF GLASS-FORMING LIQUIDS USING HIGH ENERGY X-RAY DIFFRACTION, CONTAINERLESS METHODOLOGY AND STATE-OF-THE-ART MOLECULAR DYNAMICS SIMULATION

High energy X-ray diffraction can be combined with containerless techniques to provide information on the atomic arrangements in glass-forming liquids in stable and metastable regimes. The high incident energies provide bulk diffraction data to high values of scattering vector which enables significantly more robust analysis of the local and medium-range order that influences important physical properties such as viscosity and crystal nucleation. These combined techniques have been applied to a range of oxide liquids. In this contribution we illustrate addition of further dimensions to phase space by controlling the partial pressure of oxygen that permits the study liquids containing iron. The advantages of rapid data acquisition are also demonstrated in a study of tellurite glass-forming systems where a transition from ergodic to non-ergodic regimes in the deeply supercooled liquid is shown. Finally we demonstrate how descriptions of the liquid structure can be developed by combining HEXRD with molecular dynamics simulations.

containerless techniques↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗
NASA NTRS1991

System dynamic simulation of precision segmented reflector

A joint effort was undertaken on a Precision Segmented Reflector (PSR) Project. The missions in which the PSR is to be used will use large (up to 20 m in diameter) telescopes. The essential requirement for the telescopes is that the reflector surface of the primary mirror must be made extremely precise to allow no more than a few microns of errors and, additionally, this high surface precision must be maintained when the telescope is subjected to on-orbital mechanical and thermal disturbances. Based on the mass, size, and stability considerations, reflector surface formed by segmented, probably actively or passively controlled, composite panels are regarded as most suitable for future space based astronomical telescope applications. In addition to the design and fabrication of composite panels with a surface error of less than 3 microns RMS, PSR also develops related reflector structures, materials, control, and sensing technologies. As part of the planning effort for PSR Technology Demonstration, a system model which couples the reflector, consisting of panels, support truss and actuators, and the optical bench was assembled for dynamic simulations. Random vibration analyses using seismic data obtained from actual measurements at the test site designated for PSR Technology Demonstration are described.

Shih, Choon-Foo↗
DOE OSTI2024

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗
DOE OSTI2019

Molecular Dynamics Simulation of Hydration and Swelling of Mixed-Layer Clays in the Presence of Carbon Dioxide

Swelling of clay minerals is a prime anomaly that is relevant to many problems in geoscience, oil and gas reservoirs, geophysics, and engineering. The phenomenon is even more complex when carbon dioxide is also present. Since sequestration of CO 2 in sedimentary rock has been under active consideration as a way of alleviating the increasing concerns over climate change for which CO 2 is a main culprit, studying swelling of clays minerals in the presence of CO 2 has taken on added urgency. The swelling changes the permeability and porosity of porous formations, which in turn may lead to a reactivation of dormant fractures, triggering seismic activities, and opening up new pathways for CO 2 to return to the porous formations’ surface. Although the large majority of sedimentary rocks contain various types of mixed-layer clays (MLCs) with intermixed stacking sequences of several types of distinct layers within a single crystal, the vast majority of the past experimental and computational studies of the swelling was focused on pure clays. We present the results of extensive molecular dynamics simulation of hydration and swelling of illite-montmorillonite (I-MMT) MLCs, the most common type of mixed clays, in the presence of both water and CO 2 and various combinations of interlayer cations, Na + and K + . To understand the differences with pure clays, we also report on the same phenomenon in the MMT only. At low CO 2 concentrations in Na-MMT, which has its layers’ charge concentrated in its octahedral sheet, weak ion–surface interactions result in fully hydrated ions and, therefore, more extensive swelling than in I-MMT. Without CO 2 , adsorption of the cations at the illite surface increases the hydrophobicity of its surface. Thus, in the asymmetric interlayer of I-MMT, illite with stronger surface–ion interactions causes accumulation of cations near its surface, which limits their hydration and, therefore, controls swelling of the MLCs. Further inhibition of swelling of Na-MLC can be achieved by increasing the concentration of K + in the interlayer. At higher CO 2 concentrations, however, intercalation of water and CO 2 results in a completely different behavior, since in the 2W and 3W hydration states, the hydrophobicity of the MMT surface is stronger than that of illite. Thus, the density distributions of the intercalated molecules are considerably different from that with water only, which together with disruption of the water network by CO 2 reduces the difference between the extent of swelling of the MLC and pure MMT.

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

Dynamic Simulation with on the fly Doppler Broadening treatment in RMC

In conventional transient analysis, thermal hydraulics calculations are generally coupled to deterministic neutronics solvers. Due to the limited accuracy of deterministic methods and the rapid development of computer power, Monte Carlo neutronics calculations have received increasing attention nowadays. One of challenges for Monte Carlo methods coupled to thermal hydraulics calculations is the temperature dependent cross sections. Taking the detailed temperature distribution feedback into account, each nuclide may be in a large range of temperatures, pre-generated libraries at every temperature point are no longer applicable owing to unaffordable memory footprint. Different solutions have been proposed to deal with this problem. TMS, accounting for Doppler broadening effect by the explicit treatment of thermal motion, was used to model fast transients in Serpent 2. Stochastic mixing, using two sets of libraries for each nuclide and mixing them to get the correct temperature, was used in transient analysis by Sjenitzer and Hoogenboom. Besides, there are other on the fly Doppler broadening methods being developed, including series expansion of cross sections, multipole representation, different integration algorithms based on Doppler broadening equation. In this paper, on the fly Doppler broadening treatment based on integration algorithms was applied to transient calculation in Reactor Monte Carlo code RMC . In the following section, the integration algorithms implemented and dynamic simulation method in RMC is briefly introduced. Finally, a pure physical transient fuel assembly model was chosen for numerical validation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Molecular dynamics simulations of Zika virus NS3 helicase: Insights into RNA binding site activity

America is still suffering with the outbreak of Zika virus (ZIKV) infection. Congenital ZIKV syndrome has already caused a public health emergency of international concern. However, there are still no vaccines to prevent or drugs to treat the infection caused by ZIKV. The ZIKV NS3 helicase (NS3h) protein is a promising target for drug discovery due to its essential role in viral genome replication. NS3h unwinds the viral RNA to enable the replication of the viral genome by the NS5 protein. NS3h contains two important binding sites: the NTPase binding site and the RNA binding site. Here, we used molecular dynamics (MD) simulations to study the molecular behavior of ZIKV NS3h in the presence and absence of ssRNA and the potential implications for NS3h activity and inhibition. Although there is conformational variability and poor electron densities of the RNA binding loop in various apo flaviviruses NS3h crystallographic structures, the MD trajectories of NS3h-ssRNA demonstrated that the RNA binding loop becomes more stable when NS3h is occupied by RNA. Our results suggest that the presence of RNA generates important interactions with the RNA binding loop, and these interactions stabilize the loop sufficiently that it remains in a closed conformation. This closed conformation likely keeps the ssRNA bound to the protein for a sufficient duration to enable the unwinding/replication activities of NS3h to occur. In addition, conformational changes of this RNA binding loop can change the nature and location of the optimal ligand binding site, according to ligand binding site prediction results. These are important findings to help guide the design and discovery of new inhibitors of NS3h as promising compounds to treat the ZIKV infection. - Highlights: • Ligand binding site predictions could provide useful insights into design of compounds to block ZIKV helicase activities. • Insights into dynamic features of NS3h and its interactions with ssRNA are provided through molecular dynamics simulations. • Presence of RNA generates stabilizing interactions with RNA binding loop of NS3h, which holds it in a closed conformation. • Closed conformation of the RNA binding loop keeps the ssRNA bound to the protein to enable the activity of NS3h.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

Digital system for structural dynamics simulation

State-of-the-art digital hardware and software for the simulation of complex structural dynamic interactions, such as those which occur in rotating structures (engine systems). System were incorporated in a designed to use an array of processors in which the computation for each physical subelement or functional subsystem would be assigned to a single specific processor in the simulator. These node processors are microprogrammed bit-slice microcomputers which function autonomously and can communicate with each other and a central control minicomputer over parallel digital lines. Inter-processor nearest neighbor communications busses pass the constants which represent physical constraints and boundary conditions. The node processors are connected to the six nearest neighbor node processors to simulate the actual physical interface of real substructures. Computer generated finite element mesh and force models can be developed with the aid of the central control minicomputer. The control computer also oversees the animation of a graphics display system, disk-based mass storage along with the individual processing elements.

Krauter, A. I.↗
DOE OSTI2013

Solvent effects on the ultrafast nonradiative deactivation mechanisms of thymine in aqueous solution: Excited-state QM/MM molecular dynamics simulations

On-the-fly excited-state quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations of thymine in aqueous solution are performed to investigate the role of solvent water molecules on the nonradiative deactivation process. The complete active space second-order perturbation theory (CASPT2) method is employed for a thymine molecule as the QM part in order to provide a reliable description of the excited-state potential energies. It is found that, in addition to the previously reported deactivation pathway involving the twisting of the C-C double bond in the pyrimidine ring, another efficient deactivation pathway leading to conical intersections that accompanies the out-of-plane displacement of the carbonyl group is observed in aqueous solution. Decay through this pathway is not observed in the gas phase simulations, and our analysis indicates that the hydrogen bonds with solvent water molecules play a key role in stabilizing the potential energies of thymine in this additional decay pathway.

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

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures

Performance improvements in biomolecular simulations based on molecular dynamics (MD) codes are widely desired. Unfortunately, the factors, which allowed past performance improvements, particularly the microprocessor clock frequencies, are no longer increasing. Hence, novel software and hardware solutions are being explored for accelerating performance of widely used MD codes. In this paper, we describe our efforts on porting, optimizing and tuning of Large-scale Atomic/Molecular Massively Parallel Simulator, a popular MD framework, on heterogeneous architectures: multi-core processors with graphical processing unit (GPU) accelerators. Our implementation is based on accelerating the most computationally expensive non-bonded interaction terms on the GPUs and overlapping the computation on the CPU and GPUs. This functionality is built on top of message passing interface that allows multi-level parallelism to be extracted even at the workstation level with the multi-core CPUs and allows extension of the implementation on GPU-enabled clusters. We hypothesize that the optimal benefit of heterogeneous architectures for applications will come by utilizing all possible resources (for example, CPU-cores and GPU devices on GPU-enabled clusters). Benchmarks for a range of biomolecular system sizes are provided, and an analysis is performed on four generations of NVIDIA's GPU devices. On GPU-enabled Linux clusters, by overlapping and pipelining computation and communication, we observe up to 10-folds application acceleration in multi-core and multi-GPU environments illustrating significant performance improvements. Detailed analysis of the implementation is presented that allows identification of bottlenecks in algorithm, indicating that code optimization and improvements on GPUs could allow microsecond scale simulation throughput on workstations and inexpensive GPU clusters, putting widely desired biologically relevant simulation time-scales within reach of a large user community. In order to systematically optimize simulation throughput and to enable performance prediction, we have developed a parameterized performance model that will allow developers and users to explore the performance potential of future heterogeneous systems for biological simulations.

Agarwal, Pratul K.↗