Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical Simulations”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Molecular Dynamic Simulations of Interaction of an AFM Probe with the Surface of an SCN Sample

Molecular dynamic (MD) simulations is conducted in order to estimate forces of probe-substrate interaction in the Atomic Force Microscope (AFM). First a review of available molecular dynamic techniques is given. Implementation of MD simulation is based on an object-oriented code developed at the University of Delft. Modeling of the sample material - succinonitrile (SCN) - is based on the Lennard-Jones potentials. For the polystyrene probe an atomic interaction potential is used. Due to object-oriented structure of the code modification of an atomic interaction potential is straight forward. Calculation of melting temperature is used for validation of the code and of the interaction potentials. Various fitting parameters of the probe-substrate interaction potentials are considered, as potentials fitted to certain properties and temperature ranges may not be reliable for the others. This research provides theoretical foundation for an interpretation of actual measurements of an interaction forces using AFM.

Bune, Adris↗
DOE OSTI2015

Utilizing fast multipole expansions for efficient and accurate quantum-classical molecular dynamics simulations

Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10{sup 3}-10{sup 5} molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.

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

Review of the fundamental theories behind small angle X-ray scattering, molecular dynamics simulations, and relevant integrated application

In this paper, the fundamental concepts and equations necessary for performing small angle X-ray scattering (SAXS) experiments, molecular dynamics (MD) simulations, and MDSAXS analyses were reviewed. Furthermore, several key biological and non-biological applications for SAXS, MD, and MD-SAXS are presented in this review; however, this article does not cover all possible applications. SAXS is an experimental technique used for the analysis of a wide variety of biological and non-biological structures. SAXS utilizes spherical averaging to produce one- or two-dimensional intensity profiles, from which structural data may be extracted. MD simulation is a computer simulation technique that is used to model complex biological and non-biological systems at the atomic level. MD simulations apply classical Newtonian mechanics’ equations of motion to perform force calculations and to predict the theoretical physical properties of the system. This review presents several applications that highlight the ability of both SAXS and MD to study protein folding and function in addition to non-biological applications, such as the study of mechanical, electrical, and structural properties of non-biological nanoparticles. Lastly, the potential benefits of combining SAXS and MD simulations for the study of both biological and non-biological systems are demonstrated through the presentation of several examples that combine the two techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Molecular dynamics simulations of NMR relaxation and diffusion of bulk hydrocarbons and water

Molecular dynamics (MD) simulations are used to investigate 1 H nuclear magnetic resonance (NMR) relaxation and diffusion of bulk n-C 5 H 12 to n-C 17 H 36 hydrocarbons and bulk water. The MD simulations of the 1H NMR relaxation times T 1,2 in the fast motion regime where T 1 = T 2 agree with measured (de-oxygenated) T2 data at ambient conditions, without any adjustable parameters in the interpretation of the simulation data. Likewise, the translational diffusion DT coefficients calculated using simulation con gurations are well-correlated with measured diffusion data at ambient conditions. The agreement between the predicted and experimentally measured NMR relaxation times and diffusion coefficient also validate the force elds used in the simulation. The molecular simulations naturally separate intramolecular from intermolecular dipole-dipole interactions helping bring new insight into the two NMR relaxation mechanisms as a function of molecular chain-length (i.e. carbon number). Comparison of the MD simulation results of the two relaxation mechanisms with traditional hard-sphere models used in interpreting NMR data reveals important limitations in the latter. With increasing chain length, there is substantial deviation in the molecular size inferred on the basis of the radius of gyration from simulation and the tted hard-sphere radii required to rationalize the relaxation times. This deviation is characteristic of the local nature of the NMR measurement, one that is well-captured by molecular simulations.

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

Molecular dynamics simulations of grain boundary thermal resistance in UO2

By means of molecular dynamics (MD) simulations, we have calculated Kaptiza resistance of UO2 with or without radiation damage. For coincident site lattice boundaries of different configurations, the boundary thermal resistance of unirradiated UO2 can be well described by a parameter-reduced formula by using boundary energies as variables. We extended the study to defect-loaded UO2 by introducing damage cascades in close vicinity to the boundaries. Following cascade annealing and defect migrations towards grain boundaries, the boundary energy increases and so does Kaptiza resistance. The correlations between these two still follow the same formula extracted from the unirradiated UO2. The finding will benefit multi-scale modeling of UO2 thermal properties under extreme radiation conditions by combining effects from boundary configurations and damage levels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations

To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. Overall, the results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.

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

Dynamic Simulation of a Wave Rotor Topped Turboshaft Engine

The dynamic behavior of a wave rotor topped turboshaft engine is examined using a numerical simulation. The simulation utilizes an explicit, one-dimensional, multi-passage, CFD based wave rotor code in combination with an implicit, one-dimensional, component level dynamic engine simulation code. Transient responses to rapid fuel flow rate changes and compressor inlet pressure changes are simulated and compared with those of a similarly sized, untopped, turboshaft engine. Results indicate that the wave rotor topped engine responds in a stable, and rapid manner. Furthermore, during certain transient operations, the wave rotor actually tends to enhance engine stability. In particular, there is no tendency toward surge in the compressor of the wave rotor topped engine during rapid acceleration. In fact, the compressor actually moves slightly away from the surge line during this transient. This behavior is precisely the opposite to that of an untopped engine. The simulation is described. Issues associated with integrating CFD and component level codes are discussed. Results from several transient simulations are presented and discussed.

Greendyke, R. B.↗
NASA NTRS1989

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗
DOE OSTI2014

Accelerating dissipative particle dynamics simulations for soft matter systems

Dissipative particle dynamics (DPD) is a coarse-grained particle-based simulation method that offers microscopic-scale insights into soft matter systems. We present an efficient implementation of a DPD model for graphical processing units (GPUs). As implemented in the LAMMPS molecular dynamics package, it can run effectively on current-generation supercomputers which often have hybrid nodes containing multi-core CPUs and (one or more) GPUs. Using efficient communication of information between the CPUs and GPUs, DPD interactions can be computed on the GPU while other portions of a full simulation model (boundary conditions, constraints, bonded interactions, diagnostic calculations, etc.) can be performed on the CPU. Our GPU-enhanced runs show a speedup of up to 9.5× versus many-core CPU simulations, and can run scalably across thousands of compute nodes. We briefly discuss how the new GPU implementation was validated against the CPU version for thermodynamics, diffusion, and hydrodynamic behavior. We also highlight large-scale models which the faster DPD implementation has enabled, for studies of monolayer self-assembly and thin-film instabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Scattering Density Profile Model of POPG Bilayers as Determined by Molecular Dynamics Simulations and Small-Angle Neutron and X-ray Scattering Experiments

We combine molecular dynamics (MD) simulations and experiment, both small-angle neutron (SANS) and small-angle X-ray scattering (SAXS), to determine the precise structure of bilayers composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), a lipid commonly encountered in bacterial membranes. Experiment and simulation are used to develop a one-dimensional scattering density profile (SDP) model suitable for the analysis of experimental data. The joint refinement of such data (i.e., SANS and SAXS) results in the area per lipid that is then used in the fixed-area simulations. In the final step, the direct comparison of simulated-to-experimental data gives rise to the detailed structure of POPG bilayers. From these studies we conclude that POPG s molecular area is 66.0 +/- 1.3 ^2, its overall bilayer thickness is 36.7 +/- 0.7 , and its hydrocarbon region thickness is 27.9 ( 0.6 , assuming a simulated value of 1203 ^3 for the total lipid volume.

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

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
DOE OSTI2017

Rapid sampling of stochastic displacements in Brownian dynamics simulations

Here, we present a new method for sampling stochastic displacements in Brownian Dynamics (BD) simulations of colloidal scale particles. The method relies on a new formulation for Ewald summation of the Rotne-Prager-Yamakawa (RPY) tensor, which guarantees that the real-space and wave-space contributions to the tensor are independently symmetric and positive-definite for all possible particle configurations. Brownian displacements are drawn from a superposition of two independent samples: a wave-space (far-field or long-ranged) contribution, computed using techniques from fluctuating hydrodynamics and non-uniform fast Fourier transforms; and a real-space (near-field or short-ranged) correction, computed using a Krylov subspace method. The combined computational complexity of drawing these two independent samples scales linearly with the number of particles. The proposed method circumvents the super-linear scaling exhibited by all known iterative sampling methods applied directly to the RPY tensor that results from the power law growth of the condition number of tensor with the number of particles. For geometrically dense microstructures (fractal dimension equal three), the performance is independent of volume fraction, while for tenuous microstructures (fractal dimension less than three), such as gels and polymer solutions, the performance improves with decreasing volume fraction. This is in stark contrast with other related linear-scaling methods such as the force coupling method and the fluctuating immersed boundary method, for which performance degrades with decreasing volume fraction. Calculations for hard sphere dispersions and colloidal gels are illustrated and used to explore the role of microstructure on performance of the algorithm. In practice, the logarithmic part of the predicted scaling is not observed and the algorithm scales linearly for up to 4×10 6 particles, obtaining speed ups of over an order of magnitude over existing iterative methods, and making the cost of computing Brownian displacements comparable to the cost of computing deterministic displacements in BD simulations. Finally, a high-performance implementation employing non-uniform fast Fourier transforms implemented on graphics processing units and integrated with the software package HOOMD-blue is used for benchmarking.

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

Molecular Dynamics Simulation of Tri-n-Butyl-Phophate Liquid: A Force Field Comparative Study

Molecular dynamics (MD) simulations were conducted to compare the performance of four force fields in predicting thermophysical properties of tri-n-butyl-phosphate (TBP) in the liquid phase. The intramolecular force parameters used were from the Assisted Model Building with Energy Refinement (AMBER) force field model. The van der Waals parameters were based on either the AMBER or the Optimized Potential for Liquid Simulation (OPLS) force fields. The atomic partial charges were either assigned by performing quantum chemistry calculations or utilized previously published data, and were scaled to approximate the average experimental value of the electric dipole moment. Canonical ensemble computations based on the aforementioned parameters were performed near the atmospheric pressure and temperature to obtain the electric dipole moment, mass density, and self-diffusion coefficient. In addition, the microscopic structure of the liquid was characterized via pair correlation functions between selected atoms. It has been demonstrated that the electric dipole moment can be approximated within 1% of the average experimental value by virtue of scaled atomic partial charges. The liquid mass density can be predicted within 0.5-1% of its experimentally determined value when using the corresponding charge scaling. However, in all cases the predicted self- diffusion coefficient is significantly smaller than a commonly quoted experimental measurement; this result is qualified by the fact that the uncertainty of the experimental value was not available.

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

Thermophysical properties of liquid Ni around the melting temperature from molecular dynamics simulation

Thermophysical properties of liquid nickel (Ni) around the melting temperature are investigated by means of classical molecular dynamics (MD) simulation, using three different embedded atom method potentials to model the interactions between the Ni atoms. Melting temperature, enthalpy, static structure factor, self-diffusion coefficient, shear viscosity, and thermal diffusivity are compared to recent experimental results. Using ab initio MD simulation, we also determine the static structure factor and the mean-squared displacement at the experimental melting point. For most of the properties, excellent agreement is found between experiment and simulation, provided the comparison relative to the corresponding melting temperature. We discuss the validity of the Hansen-Verlet criterion for the static structure factor as well as the Stokes-Einstein relation between self-diffusion coefficient and shear viscosity. The thermal diffusivity is extracted from the autocorrelation function of a wavenumber-dependent temperature fluctuation variable.

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

Is the Calcite-Water Interface Understood? Direct Comparisons of Molecular Dynamics Simulations with Specular X-ray Reflectivity Data

New insights into the structure of the calcite-water interface are obtained through direct model-independent comparison of multiple classical molecular dynamics (MD) simulations with high-resolution specular X-ray reflectivity (XR) data. This set of comparisons, with four different state-of-the-art force fields (including two non-polarizable, one polarizable, and one reactive force field), reveal new insights into the absolute accuracy of the simulated structures and the uniqueness of the XR-derived structural results. These four simulations, while qualitatively similar, have visibly distinct interfacial structure, and are distinguished through a quantitative comparison of the XR signals calculated from these simulations with experimental XR data. The results demonstrate that the simulated calcite-water interface structures, as a whole, are not consistent with the XR data (i.e., within their precision and accuracy). This disagreement is largely due to the simulation of the calcite lattice. The simulated interfacial water profiles show substantially different levels of agreement with the XR data. Of these, the rigid-ion model (RIM) simulations show the best consistency with the experimental XR data. Further model-dependent comparisons of the structural parameters that describe the interfacial structure (derived from both the MD simulations and the XR data) provide further insight into the sources of differences between these two approaches. Using the new insights from the RIM simulations, new structures of the calcite-water interface consistent with both the experimental data and the simulation are identified and compared to recent results.

Calcite-water interface↗
DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Hot spot and temperature analysis of shocked hydrocarbon polymer foams using molecular dynamics simulation

Hydrocarbon polymers, foams and nanocomposites are increasingly being subjected to extreme environments. Molecular scale modeling of these materials offers insight into failure mechanisms and complex response. Prior classical molecular dynamics (MD) simulations of the principal shock Hugoniot for two hydrocarbon polymers, polyethylene (PE) and poly (4-methyl-1-pentene) (PMP) have shown good agreement with density functional theory (DFT) calculations and experiments conducted at Sandia National Laboratories. We extended these results to include low-density polymer foams using nonequilibrium MD techniques and found good quantitative agreement with experiment. Here, we have measured the local temperature during void collapse to investigate the formation of hot spots and their relationship to polymer dissociation in foams.

36 MATERIALS SCIENCE↗