Search NASA⌕ Search

SEARCH · Search NASA

Results for “Quantum molecular 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Solvent effects on triplet yields in BODIPY-based photosensitizers

Here, we employ molecular dynamics simulations and quantum rate theories to elucidate the complex condensed-phase dynamics underpinning triplet-state formation in organic photosensitizers. Using models informed by first-principles calculations complete with a molecular representation of solvents of different polarities, we elucidate the interplay of the internal and environmental interactions underlying triplet yield. We find that triplet yields depend sensitively on the dielectric stabilization of the charge transfer intermediate that facilitates a transition into the triplet manifold. Our results illustrate the importance of molecularly detailed models in understanding the excited-state internal charge-transfer dynamics of photochemically relevant organic molecules.

Coello Escalante, Leonardo [University of Californ↗
DOE OSTI2020

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

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

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of Molecular Structure of Quinones and Carbon Electrode Surfaces on the Interfacial Electron Transfer Process

Quinones can undergo thermodynamically reversible proton-coupled electron transfer reactions and are being applied as electroactive compounds in aqueous organic batteries. However, the electrochemical reversibility of these compounds is affected not only by their molecular structure but also by the properties of a carbon-based electrode surface. This study combines experimental and theoretical approaches to understand this dependence. Specifically, we study the electron transfer kinetics of two synthesized quinone derivatives and two commercially available ones with a glassy carbon, a highly ordered pyrolytic graphite, and a high-edge-density graphite electrode (HEDGE). The electrochemical reversibility is notably improved on the HEDGE, which shows a higher density of defects and presents oxygenated functional groups at its surface. The electron transfer kinetics are controlled by adsorbed species onto the HEDGE. Molecular dynamics simulation and quantum mechanics calculations suggest defects with oxygen-containing functional groups, such as C-O and C$=$O, on HEDGE surfaces drive the interaction with the functional groups of the molecules, during physisorption from van der Waals forces. The presence of sulfonic acid side groups and a greater number of aromatic rings in the molecular structure may contribute to a higher stabilization of quinone derivatives on HEDGEs. We propose that high-performance carbon-based electrodes can be obtained without catalysts for organic batteries, by the engineering of carbon-based surfaces with edge-like defects and oxygenated functional groups.

25 ENERGY STORAGE↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2023

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. Furthermore, this approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Understanding Trigger Linkage Dynamics in Energetic Materials Using Mixed Picramide Nitrate Ester Explosives

The ability to predict the handling sensitivity of new organic energetic materials has been a longstanding goal. We report the synthesis and characterization of six new nitropicramide energetic materials with mixed functional groups that mimic known explosives such as nitroglycerin, erythritol tetranitrate (ETN), and pentaerythritol tetranitrate (PETN). The molecules have been studied theoretically using quantum molecular dynamics (QMD) simulations and density functional theory (DFT) calculations to identify the weakest bond in the reactants - the trigger-linkages - which control handling sensitivity, and to quantify their specific enthalpies of explosion. In good accord with the drop weight impact sensitivity data, our calculations predict that the sensitivities of the molecules are very similar owing to the small variations of the energy output and rates of trigger linkage rupture. In addition, both the QMD and DFT calculations point to the nitropicramide N–NO 2 bonds as the trigger linkages rather than the more typical O–NO 2 bonds. We propose that the switch of the trigger linkage from the nitrate esters to the nitramine groups arises from the strongly electron withdrawing character of the adjacent trinitrobenzene groups.

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

Tellurene Photodetector with High Gain and Wide Bandwidth

Two-dimensional (2D) semiconductors have been extensively explored as a new class of materials with great potential. In particular, black phosphorus (BP) has been considered to be a strong candidate for applications such as high-performance infrared photodetectors. However, the scalability of BP thin film is still a challenge, and its poor stability in the air has hampered the progress of the commercialization of BP devices. Herein, we report the use of hydrothermal-synthesized and air-stable 2D tellurene nanoflakes for broadband and ultrasensitive photodetection. The tellurene nanoflakes show high hole mobilities up to 458 cm 2 /V·s at ambient conditions, and the tellurene photodetector presents peak extrinsic responsivity of 383 A/W, 19.2 mA/W, and 18.9 mA/W at 520 nm, 1.55 μm, and 3.39 μm light wavelength, respectively. Because of the photogating effect, high gains up to 1.9 × 10 3 and 3.15 × 10 4 are obtained at 520 nm and 3.39 μm wavelength, respectively. At the communication wavelength of 1.55 μm, the tellurene photodetector exhibits an exceptionally high anisotropic behavior, and a large bandwidth of 37 MHz is obtained. The photodetection performance at different wavelength is further supported by the corresponding quantum molecular dynamics (QMD) simulations. Our approach has demonstrated the air-stable tellurene photodetectors that fully cover the short-wave infrared band with ultrafast photoresponse.

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

Hit Expansion of a Noncovalent SARS-CoV-2 Main Protease Inhibitor

Inhibition of the SARS-CoV-2 main protease (M pro ) is a major focus of drug discovery efforts against COVID-19. Here we report a hit expansion of non-covalent inhibitors of M pro . Starting from a recently discovered scaffold (The COVID Moonshot Consortium. Open Science Discovery of Oral Non-Covalent SARS-CoV-2 Main Protease Inhibitor Therapeutics. bioRxiv 2020.10.29.339317) represented by an isoquinoline series, we searched a database of over a billion compounds using a cheminformatics molecular fingerprinting approach. We identified and tested 48 compounds in enzyme inhibition assays, of which 21 exhibited inhibitory activity above 50% at 20 μM. Among these, four compounds with IC 50 values around 1 μM were found. Interestingly, despite the large search space, the isoquinolone motif was conserved in each of these four strongest binders. Room-temperature X-ray structures of co-crystallized protein–inhibitor complexes were determined up to 1.9 Å resolution for two of these compounds as well as one of the stronger inhibitors in the original isoquinoline series, revealing essential interactions with the binding site and water molecules. Molecular dynamics simulations and quantum chemical calculations further elucidate the binding interactions as well as electrostatic effects on ligand binding. The results help explain the strength of this new non-covalent scaffold for M pro inhibition and inform lead optimization efforts for this series, while demonstrating the effectiveness of a high-throughput computational approach to expanding a pharmacophore library.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

IC w22_phadiagractox Highlight: Composition-dependent melting temperature of a compound [Poster]

Figure 1 displays, ambient melting points of the Au-Cu system as a function of Cu atomic percent: present model vs. experimental data from four different sources. Figure 2 displays, ambient melting points of the U-O system as a function of O atomic percent: present model vs. both experiment and quantum molecular dynamics (QMD) simulations on the stoichiometric U-O compounds.

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

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2020

Simple and efficient algorithms for training machine learning potentials to force data

Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information than the energy alone. It may appear that training a model to this large quantity of force data would introduce significant computational costs. Actually, training to all available force data should only be a few times more expensive than training to energies alone. Here, we present a new algorithm for efficient force training, and benchmark its accuracy by training to forces from real-world datasets for organic chemistry and bulk aluminum.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗