Search NASA⌕ Search

SEARCH · Search NASA

Results for “Quantum molecular dynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Equation of state of dense plasmas: Orbital-free molecular dynamics as the limit of quantum molecular dynamics for high-Z elements

The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on a semiclassical approximation and possibly on a gradient correction, is a simulation method available at high temperature. For a high-Z element such as lutetium, we examine how orbital-free molecular dynamics applied to the equation of state of a dense plasma can be regarded as the limit of quantum molecular dynamics at high temperature. For the normal mass density and twice the normal mass density, we show that the pressures calculated with the quantum approach converge monotonically towards those calculated with the orbital-free approach; we observe a faster convergence when the orbital-free approach includes the gradient correction. We propose a method to obtain an equation of state reproducing quantum molecular dynamics results up to high temperatures where this approach cannot be directly implemented. With the results already obtained for low-Z plasmas, the present study opens the way for reproducing the quantum molecular dynamics pressure for all elements up to high temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Calculations of the thermal conductivity of National Ignition Facility target materials at temperatures near 10 eV and densities near 10 g/cc using finite-temperature quantum molecular dynamics

Using quantum finite-temperature density functional theory molecular dynamics (QMD), we performed simulations of several important materials in the Inertial Confinement Fusion-National Ignition Facility nominal target designs, comprising various mixtures of proposed ablator materials (Be or CH) with the DT fuel. Simulations were done over a range of temperatures between 5 eV and 20 eV, at densities between 7.5 and {approx}12.5 g/cc. From the QMD trajectories, we calculated the electrical and thermal conductivity. We estimated the number of free electrons per atom by fitting the frequency-dependent electrical conductivity to the Drude formula. The thermal conductivity of the fuel increases with density but that of the ablator material is insensitive to modest density variations. We find that the thermal conductivity depends strongly on the ablator/fuel mix fraction but a Faber-Ziman interpolation scheme provides a reasonable approximation. We also compare our QMD results to the Hubbard and Lee-More models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Equation of state of dense plasmas by ab initio simulations: Bridging the gap between quantum molecular dynamics and orbital-free molecular dynamics at high temperature

The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on the Thomas-Fermi semiclassical approximation and possibly on a gradient correction, is the only simulation method currently available at high temperature. We show in the case of a dense boron plasma that the two approaches give pressures differing by a few percent even at temperatures as high as a few tens of electron-volts. We indicate how the pressures obtained by orbital-free molecular dynamics can be corrected in order to appear as a limit of the quantum molecular dynamics results as temperature increases. We thus obtain a method to calculate the equation of state of a dense plasma up to high temperatures where quantum molecular dynamics cannot be directly implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Quantum molecular dynamics simulation of shock-wave experiments in aluminum

We present quantum molecular dynamics calculations of principal, porous, and double shock Hugoniots, release isentropes, and sound velocity behind the shock front for aluminum. A comprehensive analysis of available shock-wave data is performed; the agreement and discrepancies of simulation results with measurements are discussed. Special attention is paid to the melting region of aluminum along the principal Hugoniot; the boundaries of the melting zone are estimated using the self-diffusion coefficient. Also, we make a comparison with a high-quality multiphase equation of state for aluminum. Independent semiempirical and first-principle models are very close to each other in caloric variables (pressure, density, particle velocity, etc.) but the equation of state gives higher temperature on the principal Hugoniot and release isentropes than ab initio calculations. Thus, the quantum molecular dynamics method can be used for calibration of semiempirical equations of state in case of lack of experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

L2 Milestone: Simulations as Data. Quantum Molecular Dynamics Contribution

Milestone description: Direct Numerical Simulations and Quantum Molecular Dynamics simulations have been used as one would use physical experiments for cases where the physical data are not available. The goal of this milestone is to assess these simulations in the manner that they are being used and the priority of the specific configurations being planned.

Condensed Matter Physics, Superconductivity & Supe↗
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 OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

O(N) Quantum Molecular Dynamics (QMD) Algorithms, Code Requirements, Release Notes, and Library Deployment

The Exascale Computing Project (ECP) Co-design Center for Particle Applications (CoPA) is addressing the challenges for particle applications to run on upcoming exascale computing architectures. This includes short- and long-range molecular dynamics (MD), particle-in-cell (PIC), and quantum molecular dynamics (QMD). This document addresses $\mathcal{O}(N)$ QMD. The goal of this CoPA team is to develop $\mathcal{O}(N)$ linear scaling methods for QMD optimized for new and next generation architectures. This will provide the capability to simulate larger systems and run for longer timeframes than is possible today. Focus is initially on the density matrix construction. This requires efficient matrix formats and linear algebra operations for distributed architectures. A significant reduction in the wallclock time for an MD timestep is necessary to achieve practical large scale QMD. Equally challenging is to develop scalable distributed implementations. This requires the exploration of efficient data formats (dense, sparse, etc.), data decomposition and distribution across compute nodes using block- and graph-based methods. Traditional programming models such as MPI/OpenMP, hybrid MPI/OpenMP/Accelerator, as well as asynchronous task-based models will be explored to achieve maximum parallelism. Parallelism at the level of a matrix operation, solver, and/or application needs to be explored to gain the most benefit on new and emerging architectures. This report serves as a “living document” of the additions and enhancements throughout the ECP CoPA project for $\mathcal{O}(N)$ QMD. The following sections address the applications we are working with, the quantum chemistry libraries (BML and PROGRESS), and our co-design exploration through the use of a proxy application.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Quantum molecular dynamics of warm dense iron and a five-phase equation of state

Through quantum molecular dynamics (QMD), utilizing both Kohn-Sham (orbital-based) and orbital-free density functional theory, we calculate the equation of state of warm dense iron in the density range 7–30/ cm 3 and temperatures from 1 to 100 eV. A critical examination of the iron pseudopotential is made, from which we find a significant improvement at high pressure to the previous QMD calculations of Wang et al. [Phys. Rev. E 89, 023101 (2014)]. Here, our results also significantly extend the ranges of density and temperature that were attempted in that prior work. We calculate the shock Hugoniot and find very good agreement with experimental results to pressures over 20 TPa. Finally, these results are then incorporated with previous studies to generate a five-phase equation of state for iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes

Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lower temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics

Ferroelectric materials exhibit a rich range of complex polar topologies, but their study under far-from-equilibrium optical excitation has been largely unexplored because of the difficulty in modeling the multiple spatiotemporal scales involved quantum-mechanically. To study optical excitation at spatiotemporal scales where these topologies emerge, we have performed multiscale excited-state neural network quantum molecular dynamics simulations that integrate quantum-mechanical description of electronic excitation and billion-atom machine learning molecular dynamics to describe ultrafast polarization control in an archetypal ferroelectric oxide, lead titanate. Far-from-equilibrium quantum simulations reveal a marked photo-induced change in the electronic energy landscape and resulting cross-over from ferroelectric to octahedral tilting topological dynamics within picoseconds. The coupling and frustration of these dynamics, in turn, create topological defects in the form of polar strings. The demonstrated nexus of multiscale quantum simulation and machine learning will boost not only the emerging field of ferroelectric topotronics but also broader optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Equation of state and transport properties of warm dense helium via quantum molecular dynamics simulations

The equation of state, self-diffusion, and viscosity coefficients of helium have been investigated by quantum molecular dynamics (QMD) simulations in the warm dense matter regime. Our simulations are validated through the comparison with the reliable experimental data. The calculated principal and reshock Hugoniots of liquid helium are in good agreement with the gas-gun data. On this basis, we revisit the issue for helium, i.e., the possibility of the instabilities predicted by chemical models at around 2000 GPa and 10 g/cm{sup 3} along the pressure isotherms of 6309, 15 849, and 31 623 K. Our calculations show no indications of instability in this pressure-temperature region, which reconfirm the predictions of previous QMD simulations. The self-diffusion and viscosity coefficients of warm dense helium have been systematically investigated by the QMD simulations. We carefully test the finite-size effects and convergences of statistics, and obtain numerically converged self-diffusion and viscosity coefficients by using the Kubo-Green formulas. The present results have been used to evaluate the existing one component plasma models. Finally, the validation of the Stokes-Einstein relationship for helium in the warm dense regime is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2014

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

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

Ab Initio Nonadiabatic Quantum Molecular Dynamics

Here, the Born–Oppenheimer approximation underlies much of chemical simulation and provides the framework defining the potential energy surfaces that are used for much of our pictorial understanding of chemical phenomena. However, this approximation breaks down when the dynamics of molecules in excited electronic states are considered. Describing dynamics when the Born–Oppenheimer approximation breaks down requires a quantum mechanical description of the nuclei. Chemical reaction dynamics on excited electronic states is critical for many applications in renewable energy, chemical synthesis, and bioimaging. Furthermore, it is necessary in order to connect with many ultrafast pump–probe spectroscopic experiments. In this review, we provide an overview of methods that can describe nonadiabatic dynamics, with emphasis on those that are able to simultaneously address the quantum mechanics of both electrons and nuclei. Such ab initio quantum molecular dynamics methods solve the electronic Schrödinger equation alongside the nuclear dynamics and thereby avoid the need for precalculation of potential energy surfaces and nonadiabatic coupling matrix elements. Two main families of methods are commonly employed to simulate nonadiabatic dynamics in molecules: full quantum dynamics, such as the multiconfigurational time-dependent Hartree method, and classical trajectory-based approaches, such as trajectory surface hopping. In this review, we describe a third class of methods that is intermediate between the two: Gaussian basis set expansions built around trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗