Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

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

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the dynamical correlation function S(q,t) increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, here, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44\% - 71\% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Molecular Dynamics Simulations of Aldol Condensation Catalyzed by Alkylamine-Functionalized Crystalline Silica Surfaces

Molecular dynamics simulations are performed to investigate the cooperatively catalyzed aldol condensation between acetone and 4-nitrobenzaldehyde on alkylamine (or alkylenamine)-grafted silica surfaces, focusing on the mechanism of the catalytic activation of the acetone and 4-nitrobenzaldehyde by the acidic surface silanols followed by the nucleophilic attack of the basic amine functional group toward the activated reactant. From the analysis of the correlations between the catalytically active acid-base sites and reactants, it is concluded that the catalytic cooperativity of the acid-base pair can be affected by two factors: 1) the competition between the silanol and the amine (or enamine) to form a hydrogen bond with a reactant and 2) the flexibility of the alkylamine (or alkylenamine) backbone. Increasing the flexibility of the alkylamine facilitates the nucleophilic attack of the amine on the reactants. From the molecular dynamics simulations, it is found that C 3 propylamine and C 4 butylamine linkers exhibit the highest probability of reaction, which is consistent with the experimental observation that the activity of the aldol reaction on mesoporous silica depends on the length of alkylamine grafted on the silica surface. Furthermore, this simulation work serves as a pioneering study demonstrating how the molecular simulation approach can be successfully employed to investigate the cooperative catalytic activity of such bifunctional acid-base catalysts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Nonadiabatic molecular dynamics simulation: An approach based on quantum measurement picture

Mixed-quantum-classical molecular dynamics simulation implies an effective quantum measurement on the electronic states by the classical motion of atoms. Based on this insight, we propose a quantum trajectory mean-field approach for nonadiabatic molecular dynamics simulations. The new protocol provides a natural interface between the separate quantum and classical treatments, without invoking artificial surface hopping algorithm. Moreover, it also bridges two widely adopted nonadiabatic dynamics methods, the Ehrenfest mean-field theory and the trajectory surface-hopping method. Excellent agreement with the exact results is illustrated with representative model systems, including the challenging ones for traditional methods.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of water on structure and dynamics of [BMIM][PF{sub 6}] ionic liquid: An all-atom molecular dynamics simulation investigation

Composition dependent structural and dynamical properties of aqueous hydrophobic 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF{sub 6}]) ionic liquid (IL) have been investigated by using all-atom molecular dynamics simulation. We observe that addition of water does not increase significant number of dissociated ions in the solution over the pure state. As a consequence, self-diffusion coefficient of the cation and anion is comparable to each other at all water concentration similar to that is observed for the pure state. Voronoi polyhedra analysis exhibits strong dependence on the local environment of IL concentration. Void and neck distributions in Voronoi tessellation are approximately Gaussian for pure IL but upon subsequent addition of water, we observe deviation from the Gaussian behaviour with an asymmetric broadening with long tail of exponential decay at large void radius, particularly at higher water concentrations. The increase in void space and neck size at higher water concentration facilitates ionic motion, thus, decreasing dynamical heterogeneity and IL reorientation time and increases self-diffusion coefficient significantly.

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

Adaptive variational quantum dynamics simulations with compressed circuits and fewer measurements

The adaptive variational quantum dynamics simulation (AVQDS) method performs real-time evolution of quantum states using automatically generated parametrized quantum circuits that often contain substantially fewer gates than Trotter circuits. Here we report an improved version of the method, which we call AVQDS(T), by porting the tiling efficient trial circuits with rotations implemented simultaneously technique. The algorithm adaptively adds layers of disjoint unitary gates to the ansatz circuit so as to keep the McLachlan distance, a measure of the accuracy of the variational dynamics, below a fixed threshold. Here we perform benchmark noiseless AVQDS(T) simulations of quench dynamics in local spin models and compare with an alternative adaptive variational approach on quantum resource requirement. Quantum dynamical simulations implementing realistic noise channels are also reported. Finally, we propose a way to substantially alleviate the measurement overhead of AVQDS(T) while maintaining high accuracy by synergistically integrating quantum circuit calculations on quantum processing units with classical calculations using, e.g., tensor networks to evaluate the quantum geometric tensor. We showcase that this approach enables AVQDS(T) to deliver more accurate results than simulations using a fixed ansatz of comparable final depth for a significant time duration with fewer quantum resources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Trends in Ln(III) Sorption to Quartz Assessed by Molecular Dynamics Simulations and Laser Induced Flourescence Studies

Molecular dynamics simulations were performed to examine trends in trivalent lanthanide [Ln(III)] sorption to quartz surface SiOH0 and SiO- sites across the 4f period. Complementary laser induced fluorescence studies examined Eu(III) sorption to quartz at varying ionic strength such that the surface sorbed species could be extrapolated at zero ionic strength, the conditions under which the simulations are performed. This allowed for direct comparison of the data, enabling a molecular understanding of the surface sorbed species and the role of the ion surface charge density upon the interfacial reactivity. Thus, this combined theoretical and experimental approach aids in the prediction of the fate of trivalent radioactive contaminants at temporary and permanent nuclear waste storage sites. Potential of mean force molecular dynamics, as well as simulations of pre-sorbed Ln(III) species agrees with the spectroscopic study of Eu(III) sorption, indicating that strongly bound inner-sphere complexes are formed upon sorption to an SiO- site. The coordination shell of the ion contains 6-7 waters of hydration and it is predicted that surface OH groups dissociate from the quartz and bind within the inner coordination shell of Eu(III). Molecular simulations predict less-strongly bound inner2 sphere species in early lanthanides and more strongly bound species in late lanthanides, following trends in the ionic radius of the 4f ions. The participation of surface dissociated OHgroups within the inner coordination shell of the Ln(III) ion is, however, consistent across the series studied. Sorption to a fully protonated quartz surface is not predicted to be favorable by any Ln(III), except perhaps Lu.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2012

The electronic and optical properties of warm dense nitrous oxide using quantum molecular dynamics simulations

First-principles molecular-dynamics simulations based on density-functional theory have been used to study the electronic and optical properties of fluid nitrous oxide under extreme conditions. Systematic descriptions of pair-correlation function, atomic structure, and the charge density distribution are used to investigate the dissociation of fluid nitrous oxide. The electrical and optical properties are derived from the Kubo-Greenwood formula. It is found that the nonmetal-metal transition for fluid nitrous oxide can be directly associated to the dissociation and has significant influence on the optical properties of the fluid.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Clustering effects in ionic polymers: Molecular dynamics simulations

Ionic clusters control the structure, dynamics, and transport in soft matter. Incorporating a small fraction of ionizable groups in polymers substantially reduces the mobility of the macromolecules in melts. Furthermore, these ionic groups often associate into random clusters in melts, where the distribution and morphology of the clusters impact the transport in these materials. Here, using molecular dynamic simulations we demonstrate a clear correlation between cluster size and morphology with the polymer mobility in melts of sulfonated polystyrene. We show that in low dielectric media ladderlike clusters that are lower in energy compared with spherical assemblies are formed. Reducing the electrostatic interactions by enhancing the dielectric constant leads to morphological transformation from ladderlike clusters to globular assemblies. Finally, decrease in electrostatic interaction significantly enhances the mobility of the polymer.

36 MATERIALS SCIENCE↗
DOE OSTI2013

On the structural affinity of macromolecules with different biological properties: Molecular dynamics simulations of a series of TEM-1 mutants

Highlights: •We have performed molecular dynamics simulations of TEM-1 mutants. •Mutations effects on the mechanical properties are considered. •Mutants do not significantly alter the average enzymes structure. •Mutants produce sharp alterations in enzyme conformational repertoire. •Mutants also produce changes in the active site volume. -- Abstract: Molecular Dynamics simulations have been carried out in order to provide a molecular rationalization of the biological and thermodynamic differences observed for a class of TEM β-lactamases. In particular we have considered the TEM-1(wt), the single point mutants TEM-40 and TEM-19 representative of IRT and ESBL classes respectively, and TEM-1 mutant M182T, TEM-32 and TEM-20 which differ from the first three for the additional of M182T mutation. Results indicate that most of the thermodynamic, and probably biological behaviour of these systems arise from subtle effects which, starting from the alterations of the local interactions, produce drastic modifications of the conformational space spanned by the enzymes. The present study suggests that systems showing essentially the same secondary and tertiary structure may differentiate their chemical–biological activity essentially (and probably exclusively) on the basis of the thermal fluctuations occurring in their physiological environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Tensile testing of Fe and FeCr nanowires using molecular dynamics simulations

Using molecular dynamics, we have studied the behaviour of cylindrical [001]-oriented Fe and FeCr nanowires under uniaxial tensile strain with both an embedded atom method (EAM) and a Tersoff-like bond order potential. The mechanical properties were analysed and the deformation mechanism was studied and compared between the potentials. The effects of chromium content and size of the wire were studied. Both potentials show elongation by deformation twinning in the 〈111〉/(211) system resulting in a significantly stiffer and stronger [110]-axial nanowire. The pure iron nanowires are elastically softer than bulk iron and an addition of chromium has both a softening and weakening effect. The bond order potential shows a strong dependence on chromium concentration, while the dependence is considerably weaker for the EAM potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Structure and Dynamics of N, N-diethyl-N-methylammonium Triflate Ionic Liquid, Neat and with Water, from Molecular Dynamics Simulations

We investigated by means of molecular dynamics simulations the properties (structure, thermodynamics, ion transport, and dynamics) of the protic ionic liquid N,N-diethyl-N-methyl-ammonium triflate (dema:Tfl) and of selected aqueous mixtures of dema:Tfl. This ionic liquid, a good candidate for a water-free proton exchange membrane, is shown to exhibit high ion mobility and conductivity. For bulk melts in the temperature range of 303-453K, both liquid densities and enthalpies of vaporization are found to decrease roughly linearly with increasing temperature. The radial distribution functions reveal a significant long-range structural correlation. The ammonium cations [dema]+ are found to diffuse slightly faster than the triflate anions [Tfl]-, and both types of ions exhibit enhanced mobility at higher temperatures, leading to higher ionic conductivity of these ionic liquids. Analysis of the dynamics of ion pairing clearly points to the existence of long-lived contact ion pairs in this ionic liquid. We also examined the effects of water on the ionic properties of dema:Tfl-water mixtures. From the structural analysis it was found that water molecules tend to replace counter ions in the coordination shell and hydrogen bond to both ions, thus weakening their mutual association. As water concentration increases, water molecules start to connect with each other and eventually form a large network that percolates through the system. It is also found that water has a strong influence on the ion dynamics in the mixtures. As the concentration of water increases, both translational and rotational motion of [dema]+ and [Tfl]- are significantly enhanced. As a result, higher ionic conductivity is observed with increased hydration level. This work was supported by the US Department of Energy Basic Energy Sciences' Chemical Sciences, Geosciences & Biosciences Division. Pacific Northwest National Laboratory is operated by Battelle for the US Department of Energy.

08 HYDROGEN↗
DOE OSTI2016

Nanoindentation of the pristine and irradiated forms of a sodium borosilicate glass: Insights from molecular dynamics simulations

We have carried out classical molecular dynamics simulations in order to get insight into the atomistic mechanisms of the deformation during nanoindentation of the pristine and irradiated forms of a sodium borosilicate glass. In terms of the glass hardness, we have found that the primary factor affecting the decrease of hardness after irradiation is depolymerization rather than free volume, and we argue that this is a general trend applicable to other borosilicate glasses with similar compositions. We have analyzed the changes of the short- and medium-range structures under deformation and found that the creation of oxygen triclusters is an important mechanism in order to describe the deformation of highly polymerized borosilicate glasses and is essential in the understanding of the folding of large rings under stress. We have equally found that the less polymerized glasses present a higher amount of relative densification, while the analysis of bond-breaking during the nanoindentation has showed that shear flow is more likely to appear around sodium atoms. The results provided in this study can be proven to be useful in the interpretation of experimental results.

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

Ab initio molecular dynamics simulations of AlN responding to low energy particle radiation

Ab initio molecular dynamics simulations of low energy recoil events in wurtzite AlN have been performed to determine threshold displacement energies, defect production and evolution mechanisms, role of partial charge transfer during the process, and the influence of irradiation-induced defects on the properties of AlN. Here, the results show that the threshold displacement energies, E d , along the direction parallel to the basal planes are smaller than those perpendicular to the basal planes. The minimum E d values are determined to be 19 eV and 55 eV for N and Al atom, respectively, which occur along the [$\overline{11}20$] direction. In general, the threshold displacement energies for N are smaller than those for Al atom, indicating the N defects would be dominant under irradiation. The defect production mechanisms have been analyzed. It is found that charge transfer and redistribution for both the primary knock-on atom and the subsequent recoil atoms play a significant role in defect production and evolution. Similar to the trend in oxide materials, there is a nearly linear relationship between E d and the total amount of charge transfer at the potential energy peak in AlN, which provides guidance on the development of charge-transfer interatomic potentials for classic molecular dynamics simulations. Finally, the response behavior of AlN to low energy irradiation is qualitatively investigated. The existence of irradiation-induced defects significantly modifies the electronic structure, and thus affects the magnetic, electronic and optical properties of AlN. In conclusion, these findings further enrich the understanding of defects in the wide bandgap semiconductor of AlN.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

The VENUS/NWChem Software Package. Tight Coupling Between Chemical Dynamics Simulations and Electronic Structure Theory

The interface for VENUS and NWChem, and the resulting software package for direct dynamics simulations are described. The coupling of the two codes is considered to be a tight coupling. The two codes are compiled and linked together and act as one executable with data being passed between the two codes through routine calls. The advantages of this type of coupling are discussed. The interface has been designed to have as little interference as possible with the core codes of both VENUS and NWChem. VENUS is the code that propagates the direct dynamics trajectories and, therefore, is the program that drives the overall execution of VENUS/NWChem. VENUS has remained an essentially sequential code, which uses the highly parallel structure of NWChem. Subroutines of the interface which accomplish the data transmission and communication between the two computer programs are described. Recent examples of the use of VENUS/NWChem for direct dynamics simulations are summarized.

chemical dynamics simulation↗
DOE OSTI2024

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Molecular Dynamics Simulation of Binary Fluid in a Nanochannel

This paper presents the results from a molecular dynamics simulation of binary fluid (mixture of argon and krypton) in the nanochannel flow. The computational software LAMMPS is used for carrying out the molecular dynamics simulations. Binary fluids of argon and krypton with varying concentration of atom species were taken for two densities 0.65 and 0.45. The fluid flow takes place between two parallel plates and is bounded by horizontal walls in one direction and periodic boundary conditions are imposed in the other two directions. To drive the flow, a constant force is applied in one direction. Each fluid atom interacts with other fluid atoms and wall atoms through Week-Chandler-Anderson (WCA) potential. The velocity profile has been looked at for three nanochannel widths i.e for 12{sigma}, 14{sigma} and 16{sigma} and also for the different concentration of two species. The velocity profile of the binary fluid predicted by the simulations agrees with the quadratic shape of the analytical solution of a Poiseuille flow in continuum theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗