Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Understanding misfit strain releasing mechanisms via molecular dynamics simulations of CdTe growth on {112}zinc-blende CdS

Molecular dynamics simulations have been used to analyse microstructures of CdTe films grown on {112} surfaces of zinc-blende CdS. Interestingly, CdTe films grow in <331> orientations as opposed to <112> epitaxial orientations. At the CdTe-{331}/CdS-{112} interface, however, there exists an axis that is parallel to the <110> orientation of both CdS and CdTe. It is the direction orthogonal to this <110> that becomes different, being <116> for CdTe and <111> for CdS, respectively. Missing CdTe-{110} planes are found along the <110> axis, suggesting that the misfit strain is released by the conventional misfit dislocation mechanism along this axis. In the orthogonal axis, the misfit strain is found to be more effectively released by the new grain orientation mechanism. Our finding is supported by literature experimental observations of the change of growth direction when Cd 0.96 Zn 0.04 Te films are deposited on GaAs. Lastly the analyses of energetics clearly demonstrate the cause for the formation of the new orientation, and the insights gained from our studies can help understand the grain structures experimentally observed in lattice mismatched systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Domain-specific compilers for dynamic simulations of quantum materials on quantum computers

Abstract Simulation of the dynamics of quantum materials is emerging as a promising scientific application for noisy intermediate-scale quantum (NISQ) computers. Due to their high gate-error rates and short decoherence times, however, NISQ computers can only produce high-fidelity results for those quantum circuits smaller than some given circuit size. Dynamic simulations, therefore, pose a challenge as current algorithms produce circuits that grow in size with each subsequent time-step of the simulation. This underscores the crucial role of quantum circuit compilers to produce executable quantum circuits of minimal size, thereby maximizing the range of physical phenomena that can be studied within the NISQ fidelity budget. Here, we present two domain-specific (DS) quantum circuit compilers for the Rigetti and IBM quantum computers, specifically designed to compile circuits simulating dynamics under a special class of time-dependent Hamiltonians. The compilers outperform state-of-the-art general-purpose compilers in terms of circuit size reduction by around 25%–30% as well as wall-clock compilation time by around 40% (dependent on system size and simulation time-step). Drawing on heuristic techniques commonly used in artificial intelligence, both compilers scale well with simulation time-step and system size. Code for both compilers is open-source and packaged into a full-stack quantum simulation software with tutorials included for ease of use for future researchers wishing to perform dynamic simulations of quantum materials on quantum computers. As our DS compilers provide significant improvements in both compilation time and simulation fidelity, they provide a building block for accelerating progress toward physical quantum supremacy.

Physics↗
DOE OSTI2015

Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Aggregation of cyclic polypeptoids bearing zwitterionic end-groups with attractive dipole–dipole and solvophobic interactions: a study by small-angle neutron scattering and molecular dynamics simulation

Aggregation behavior of cyclic polypeptoids bearing zwitterionic end-groups in methanol has been studied using a combination of experimental and simulation techniques. Here, the data from SANS and cryo-TEM indicate that the solution contains small clusters of these cyclic polypeptoids, ranging from a single polypeptoid chain to small oligomers, while the linear counterpart shows no cluster formation. Atomistic molecular dynamics simulations reveal that the driving force for this clustering behavior is due to the interplay between the effective repulsion due to the solvation of the dipoles formed by the charged end-groups in each polypeptoid chain and the attractive forces due to dipole–dipole interactions and the solvophobic effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2021

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents computational fluid dynamics simulations of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the airstream oscillator system (AOS) in operation. Test section flow angle can be modulated using sinusoidally oscillating vanes which are located in pairs within the contraction section of the tunnel circuit. A compilation of previously published experimental data for validation has been reviewed that contains hot-wire anemometer and fast-response probe data acquired in an empty test section with the vanes oscillating in-phase and out-of-phase. These experimental data have not previously been compared with computational results. The present results provide an opportunity for validation of the simulations of the sinusoidal AOS flow field in the TDT test section. Computed tunnel centerline, lateral and cross-section contour data of flow angle and vorticity show that resolving the vane tip vortex is critical to accurately predicting the test section center flow angle. The computed results show, as do the experimental data, what appears to be a resonance at higher frequencies. It is currently not known what is the cause of that resonance. Further refinement of the mesh is expected to improve the comparisons with experimental data

Robert E Bartels↗
DOE OSTI2020

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗
DOE OSTI2015

Effect of Cu2+ Activation on Interfacial Water Structure at the Sphalerite Surface as Studied by Molecular Dynamics Simulation

In the first part of this paper, an experimental contact angle study of the fresh and Cu2+ activated sphalerite-ZnS surface as well as the covellite-CuS (001) surface is reported describing the increased hydrophobic character of the surface during Cu2+ activation. In addition to these experimental results, the fresh sphalerite-ZnS (110), copper-zinc sulfide-CuZnS2 (110), villamaninite- CuS2 (100), and covellite-CuS (001) surfaces were examined using Molecular Dynamics Simulation (MDS). Our MDS results on the behavior of interfacial water at the fresh sphalerite-ZnS (110), copper-zinc sulfide-CuZnS2 (110), villamaninite-CuS2 (100), and covellite-CuS (001) surfaces include simulated contact angles, water number density distribution, water dipole orientation, water residence time, and hydrogen-bonding considerations. The copper content at the Cu2+ activated sphalerite surface seems to account for the increased hydrophobicity as revealed by both experimental and MD simulated contact angle measurements. The relatively greater hydrophobic character developed at the Cu2+ activated sphalerite surface and at the copper-zinc sulfide surface has been described by MDS, based on the structure of interfacial water and its dynamic properties. L.X.D. acknowledges funding from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences.

Jin, Jiaqi↗
DOE OSTI2014

Atomic Movies of Laser-Induced Structural and Phase Transformations from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations of laser-materials interactions are playing an important role in investigation of complex and highly non-equilibrium processes involved in short pulse laser processing and surface modification. This role is defined by the ability of MD simulations to reveal in-depth information on the structural and phase transformations induced by the laser excitation and, at the same time, to provide clear visual representations, or “atomic movies,” of laser-induced dynamic processes. This chapter provides a brief overview of recent progress in the description of laser coupling and relaxation of photo-excited states in metals, semiconductors, insulators and molecular systems within the general framework of the classical MD technique and presents several examples of MD simulations of laser melting, generation of crystal defects, photomechanical spallation, explosive boiling and molecular entrainment in laser ablation. Possible directions of future progress in atomistic modeling of laser-materials interactions and the potential role of MD simulations in the design of an integrated multiscale computational model capable of accounting for interrelations between processes occurring at different time- and length-scales are discussed.

Classical Molecular Dynamic↗
DOE OSTI2015

The high pressure structure and equation of state of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) up to 20 GPa: X-ray diffraction measurements and first principles molecular dynamics simulations

Recent theoretical studies of 2,6-diamino-3,5-dinitropyrazine-1-oxide (C{sub 4}H{sub 4}N{sub 6}O{sub 5} Lawrence Livermore Molecule No. 105, LLM-105) report unreacted high pressure equations of state that include several structural phase transitions, between 8 and 50 GPa, while one published experimental study reports equation of state (EOS) data up to a pressure of 6 GPa with no observed transition. Here we report the results of a synchrotron-based X-ray diffraction study and also ambient temperature isobaric-isothermal atomistic molecular dynamics simulations of LLM-105 up to 20 GPa. We find that the ambient pressure phase remains stable up to 20 GPa; there is no indication of a pressure induced phase transition. We do find a prominent decrease in b-axis compressibility starting at approximately 13 GPa and attribute the stiffening to a critical length where inter-sheet distance becomes similar to the intermolecular distance within individual sheets. The ambient temperature isothermal equation of state was determined through refinements of measured X-ray diffraction patterns. The pressure-volume data were fit using various EOS models to yield bulk moduli with corresponding pressure derivatives. We find very good agreement between the experimental and theoretically derived EOS.

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

A novel energy conversion based method for velocity correction in molecular dynamics simulations

Molecular dynamics (MD) simulation has become an important tool for studying micro- or nano-scale dynamics and the statistical properties of fluids and solids. In MD simulations, there are mainly two approaches: equilibrium and non-equilibrium molecular dynamics (EMD and NEMD). In this paper, a new energy conversion based correction (ECBC) method for MD is developed. Unlike the traditional systematic correction based on macroscopic parameters, the ECBC method is developed strictly based on the physical interaction processes between the pair of molecules or atoms. The developed ECBC method can apply to EMD and NEMD directly. While using MD with this method, the difference between the EMD and NEMD is eliminated, and no macroscopic parameters such as external imposed potentials or coefficients are needed. With this method, many limits of using MD are lifted. The application scope of MD is greatly extended.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

GPU-accelerated Tersoff potentials for massively parallel Molecular Dynamics simulations

The Tersoff potential is one of the empirical many-body potentials that has been widely used in simulation studies at atomic scales. Unlike pair-wise potentials, the Tersoff potential involves three-body terms, which require much more arithmetic operations and data dependency. In this contribution, we have implemented the GPU-accelerated version of several variants of the Tersoff potential for LAMMPS, an open-source massively parallel Molecular Dynamics code. Compared to the existing MPI implementation in LAMMPS, the GPU implementation exhibits a better scalability and offers a speedup of 2.2X when run on 1000 compute nodes on the Titan supercomputer. On a single node, the speedup ranges from 2.0 to 8.0 times, depending on the number of atoms per GPU and hardware configurations. The most notable features of our GPU-accelerated version include its design for MPI/accelerator heterogeneous parallelism, its compatibility with other functionalities in LAMMPS, its ability to give deterministic results and to support both NVIDIA CUDA- and OpenCL-enabled accelerators. Our implementation is now part of the GPU package in LAMMPS and accessible for public use.

Computer Science↗