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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Growth of bi- and tri-layered graphene on silicon carbide substrate via molecular dynamics simulation

Molecular dynamics (MD) simulation with simulated annealing method is used to study the growth process of bi- and tri-layered graphene on a 6H-SiC (0001) substrate via molecular dynamics simulation. Tersoff-Albe-Erhart (TEA) potential is used to describe the inter-atomic interactions among the atoms in the system. The formation temperature, averaged carbon-carbon bond length, pair correlation function, binding energy and the distance between the graphene formed and the SiC substrate are quantified. The growth mechanism, graphitization of graphene on the SiC substrate and characteristics of the surface morphology of the graphene sheet obtained in our MD simulation compare well to that observed in epitaxially grown graphene experiments and other simulation works.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Characterization of radiation damage in TiO 2 using molecular dynamics simulations

Molecular dynamics simulations are carried out to characterize irradiation effects in TiO 2 rutile, for wide ranges of temperatures (300–900 K) and primary knock-on atom (PKA) energies (1–10 keV). The number of residual defects decreases with increased temperature and decreased PKA energy, but is independent of PKA type. In the ballistic phase, more oxygen than titanium defects are produced, however, the primary residual defects are titanium vacancies and interstitials. Defect clustering depends on the PKA energy, temperature, and defect production. For some 10 keV PKAs, the largest cluster of vacancies at the peak of the ballistic phase and after annealing has up to ≈1200 and 100 vacancies, respectively. For the 10 keV PKAs at 300 K, the energy storage, primarily in residual Ti vacancies and interstitials, is estimated at 140–310 eV. It decreases with increased temperature to as little as 5–180 eV at 900 K. In conclusion, selected area electron diffraction patterns and radial distribution functions confirm that although localized amorphous regions form during the ballistic phase, TiO 2 regains full crystallinity after annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Bergman-type medium range order in amorphous Zr77Rh23 alloy studied by ab initio molecular dynamics simulations

In recent years, some arguments about the existence of medium-range order (MRO) in the Zr-Rh system have been put forward. However, research on the structural features of the Zr-Rh binary alloy at the atomic level is still lacking. This study uses ab initio molecular dynamics simulations to systematically study the local structures of Zr77Rh23 from the liquid to the glassy states. Pair correlation function (PCF), coordination number (CN), Honeycutt–Anderson(HA) index, bond-angle distribution functions, and the Voronoi tessellation method are used to reveal a clear icosahedral-like configuration in the amorphous Zr77Rh23 alloy. It is noteworthy that the splitting in the second peak of the partial PDF implies the existence of a medium range order (MRO) in the Zr77Rh23 system. We obtain the local order in three-dimensional atomic density distributions by using a new atomistic cluster alignment (ACA) method. Interestingly, a Bergman-type MRO is observed in the glassy Zr77Rh23. Furthermore, the spatial distribution and connections of icosahedral-like clusters are shown to further demonstrate the MRO network. Our findings shed light on the nature of atomic local structures of amorphous Zr77Rh23 alloy and have important implications to understanding the formation of various MROs in metallic glasses.

Guo, Y. R.↗
DOE OSTI2020

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗
DOE OSTI2021

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2015

Thermal Interface Conductance Between Aluminum and Silicon by Molecular Dynamics Simulations

The thermal interface conductance between Al and Si was simulated by a non-equilibrium molecular dynamics method. In the simulations, the coupling between electrons and phonons in Al are considered by using a stochastic force. The findings demonstrate the size dependence of the interface thermal conductance and the effect of electron–phonon coupling on the interface thermal conductance. To understand the mechanism of interface resistance, the vibration power spectra are calculated. We postulate that the atomic level disorder near the interface is an important aspect of interfacial phonon transport, which leads to a modification of the phonon states near the interface. There, the vibrational spectrum near the interface greatly differs from the bulk. This change in the vibrational spectrum affects the results predicted by AMM and DMM theories and indicates new physics is involved with phonon transport across interfaces.

42 ENGINEERING↗
DOE OSTI2019

Effects of electron-phonon coupling and electronic thermal conductivity in high energy molecular dynamics simulations of irradiation cascades in nickel

The two-temperature model has been applied to investigate the effects of the electronic subsystem on 150 keV Ni ion cascades in nickel using molecular dynamics simulation. We explore the effects of the magnitude of the electron-phonon coupling and the electronic thermal conductivity on defect production and cluster formation. It has been found that stronger electron-phonon coupling allows larger and more rapid energy feedback to the atomic subsystem, leading to reduction of number of point defects and suppression of the formation of larger defect clusters. It was observed that larger electronic thermal conductivity results in slightly increased number of point defects and larger size vacancy clusters. The latter takes place because of suppression of point defects recombination in faster cooling areas of initial damage.

Zarkadoula, Eva↗
DOE OSTI2019

Dynamical properties of a room temperature ionic liquid: Using molecular dynamics simulations to implement a dynamic ion cage model [plus supplemental information]

The transport behavior of ionic liquids (ILs) is pivotal for a variety of applications, especially when ILs are used as electrolytes. Many aspects of the transport dynamics of ILs remain to be understood. Here, a common ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BmimNTf 2 ), was studied with molecular dynamics simulations. The results show that BmimNTf 2 displays typical structural relaxation, subdiffusive behavior, and a breakdown of the Stokes-Einstein diffusion relation as in glass-forming liquids. In addition, the simulations show that the translational dynamics, reorientation dynamics, and structural relaxation dynamics are well described by the Vogel-Fulcher-Tammann equation like fragile glass forming liquids. Building on previous work that employed ion cage models, it was found that the diffusion dynamics of the cations and anions were well described by a hopping process random walk where the step time is the ion cage lifetime obtained from the cage correlation function. Detailed analysis of the ion cage structures indicated that the electrostatic potential energy of the ion cage dominates the diffusion dynamics of the caged ion. The ion orientational relaxation dynamics showed that ion reorientation is a necessary step for ion cage restructuring. The dynamic ion cage model description of ion diffusion presented here may have implications for designing ILs to control their transport behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting the Dissolution Rate of Borosilicate Glasses using QSPR analysis based on molecular dynamics simulations

Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition-property relationships of glass and other materials with complex structures. A series of 20 borosilicate and boroaluminosilicate glasses have been modeled by using MD simulations with recently developed effective potentials. Short- and medium-range structures of these glasses were analyzed and, based on these structural information, QSPR analysis of the initial dissolution rate (r0) was made and compared with measured r0 at 90°C and pH 9 using various structural descriptors such as percentage of bridging oxygen species, network connectivity and average ring size. The structural descriptors, Fnet, containing energetic information such as single bond strength and other structural information were also used. It was found that overall network connectivity, average ring size and Fnet give reasonable predictions of the r0 of studied glasses, given the conditions that the glasses are homogeneous and dissolve congruently. Modifying glass compositions to account preferential release of modifiers gives a better prediction for incongruently dissolving glasses. The results were compared with our recent work of predicting glass dissolution behavior from compositions using the topological-constraints-based models.

Du, Jincheng↗
DOE OSTI2017

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Dynamics of water confined in lyotropic liquid crystals: Molecular dynamics simulations of the dynamic structure factor

The properties of water under confinement are of practical and fundamental interest. Here in this work we study the properties of water in the self-assembled lyotropic phases of gemini surfactants with a focus on testing the standard analysis of quasi-elastic neutron scattering (QENS) experiments. In QENS experiments the dynamic structure factor is measured and fit to models to extract the translational diffusion constant, D T , and rotational relaxation time, τ R . We test this procedure by using simulation results for the dynamic structure factor, extracting the dynamic parameters from the fit as is typically done in experiments, and comparing the values to those directly measured in the simulations. We find that the decoupling approximation, where the intermediate scattering function is assumed to be a product of translational and rotational contributions, is quite accurate. The jump-diffusion and isotropic rotation models, however, are not accurate when the degree of confinement is high. In particular, the exponential approximations for the intermediate scattering function fail for highly confined water and the values of D T and τ R can differ from the measured value by as much as a factor of two. Other models have more fit parameters, however, and with the range of energies and wave-vectors accessible to QENS, the typical analysis appears to be the best choice. In the most confined lamellar phase, the dynamics are sufficiently slow that QENS does not access a large enough time scale and neutron spin echo measurements would be a valuable technique in addition to QENS.

36 MATERIALS SCIENCE↗
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↗