Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ab initio quantum dynamics of plasmonic charge carriers

Surface plasmons respond strongly to electric fields and generate energetic (hot) charge carriers that can be used in optoelectronic devices. However, utilization of hot carriers has to outcompete fast carrier relaxation in metals. Often, the reported efficiencies of nanoscale devices based on plasmon excitations are low and the mechanisms of device operation remain inconsistent within the field. Further developments hinge on fundamental understanding of the nature and kinetics of plasmon decay processes. In this work, we provide an overview of time domain ab initio modeling of hot carrier dynamics, demonstrating that atomistic details of the surface structure of plasmonic materials, and their chemical interaction with semiconducting and molecular charge acceptors and substrates, play key roles in governing the dynamics pathways. Considered case studies represent various dynamics channels and illustrate solar energy and optoelectronic applications.

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

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗
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 OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

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

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

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

Constraining nucleon effective masses with flow and stopping observables from the S π RIT experiment

Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow (υ 1 ), elliptic flow (υ 2 ) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the SπRIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, S 0 , and slope of the symmetry energy, L, at saturation density, isoscalar effective mass, $m^{⁎}_{S}$/$m_{N}$ , isovector effective mass, $m^{⁎}_{υ}$/$m_{N}$ and the in-medium cross-section enhancement factor η. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the SπRIT data suggest a narrow range of preferred values for $m^{⁎}_{S}$/$m_{N}$, $m^{⁎}_{υ}$/$m_{N}$ and η.

Astronomy & Astrophysics↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials (Final Technical Report)

The goal of this project was to conduct dynamic compression experiments and predictive simulations to reveal the fundamental high-energy-density (HED) physics of amorphous carbon. These results are essential for assessing amorphous carbon as a potential ablator material for next generation of inertial confinement fusion (ICF) capsules. We made significant progress in exploring the HED properties of amorphous carbon through experiments at Omega EP Laser and the European XFEL, in addition to billion-atom, quantum-accurate molecular dynamics (MD) simulations. Through our joint experimental and simulation program, we mapped the phase diagram of amorphous carbon, uncovering its range of metastability and identifying phase transitions to diamond and liquid carbon along the Hugoniot and at higher pressures using double shock compression pathways. Our findings indicate that amorphous carbon melts at significantly lower shock pressures than high-density carbon (diamond). However, nanocrystalline diamond nucleates across a broad range of pressures and temperatures. This emergence of the nanocrystalline microstructure during compression can negatively impact the planarity of the shock front and potentially trigger ablator/fuel mixing during Inertial Fusion Energy (IFE) applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Geophysical Impacts and Spectroscopic Identification of a Hydrous Iron Sulfate on Icy Worlds

Over geologic time-scales, large volumes of exogenic sulfur ions from Io's plasma torus have been supplied to the surface of Europa and Ganymede, which, combined with recent interpretations of orbiter images, dynamical modeling, and surface-subsurface exchange, suggests further sulfur transport into the interior of the icy worlds. These observations motivate mixed-phase spectral modeling for interpreting orbiter spectroscopy data and determination of hydration states of candidate surface materials including hydrous sulfates. In this work, we present a combined experimental and theoretical study of the low temperature and high pressure vibrational spectral signature of the iron-sulfate monohydrate endmember, szomolnokite (FeSO 4 ·H 2 O). By employing synchrotron Fourier-transform infrared spectroscopy (FTIR) in the diamond anvil cell up to 23 GPa and down to 20 K, we explore the extreme range of pressure-temperature domains relevant to icy environments throughout our solar system and beyond. Combined with our density-functional theory quantum-mechanics molecular dynamics results, we demonstrate that experimentally observed infrared features in the O-H stretching region commonly associated with nH 2 O (n > 1) hydration states can be attributed to a pure monohydrate without the need for pressure-induced exsolved ice, other coexisting hydrous iron sulfates, or strong overtone and combination modes. We further discuss the possibility of lateral variations in density and shear properties on icy worlds associated with temperature variations and the high-pressure phases of kieserite group monohydrated sulfates.

Geosciences↗
DOE OSTI2021

A Multireference Approach to Electron and Electron–Nuclear Dynamics in Nanomaterials (Final Report)

Many important chemical and physical phenomena involve dynamics on large number of electronic states. Thus, there is a critical need to develop methods to simulate dynamics in dense manifolds of states. Towards this end, we have: a) developed the multiple cloning in dense manifolds of states (MCDMS) method, which is capable of accurately modeling the quantum mechanical coherence between populations on a large number of electronic states, b) implemented MCDMS into the free, open-source PySpawn software package, c) developed graphics processing unit-accelerated algorithms modeling electron dynamics in light fields via Floquet time-dependent configuration interaction (F-TDCI), and d) critically compared different orbital bases in order to achieve an accurate and efficient F-TDCI expansion. This grant ended in August 2020, when our group moved to from Michigan State University to Stony Brook University, where this project continues under grant number DE-SC0021643.

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

Cooperative diffusion in body-centered cubic iron in Earth and super-Earths’ inner core conditions

Abstract The physical chemistry of iron at the inner-core conditions is key to understanding the evolution and habitability of Earth and super-Earth planets. Based on full first-principles simulations, we report cooperative diffusion along the longitudinally fast ⟨ 111 ⟩ directions of body-centered cubic (bcc) iron in temperature ranges of up to 2000–4000 K below melting and pressures of ∼300–4000 GPa. The diffusion is due to the low energy barrier in the corresponding direction and is accompanied by mechanical and dynamical stability, as well as strong elastic anisotropy of bcc iron. These findings provide a possible explanation for seismological signatures of the Earth’s inner core, particularly the positive correlation between P wave velocity and attenuation. The diffusion can also change the detailed mechanism of core convection by increasing the diffusivity and electrical conductivity and lowering the viscosity. The results need to be considered in future geophysical and planetary models and should motivate future studies of materials under extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. Furthermore, these backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si–O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coordination and thermophysical properties of select trivalent lanthanides in LiCl–KCl

The coordination chemistry of various fission and decay products, such as actinides and lanthanides, are crucial to the commercial deployment of molten salt reactors as they can affect the thermophysical properties. In this report we examined the structure, coordination environment, and physical properties such as the density and the vibrational density of states for three lanthanide species, namely Ce, Eu, and Sm in the LiCl–KCl eutectic system using a combination of quantum mechanics simulations and spectroscopic experiments. Quantum mechanics molecular dynamics (QM-MD) modelling was employed to determine the physical properties of each system resulting in accurate local coordination of each species. Then, the vibrational density of states (DOS) was determined using a two-phase thermodynamic modelling which was then compared to the experimentally obtained Raman spectra of the species in molten LiCl–KCl having the eutectic composition. We find that Ce 3+ , Eu 3+ and Sm 3+ all adopt octahedral local coordination environments in the eutectic salt composition in good agreement with experimental results. Ce3+ is found to fluctuate between an octahedral six-coordinated and a seven-coordinated structure due to the increased local proximity of Cl in the eutectic salt, resulting in a lower fluidicity/diffusivity than the other trivalent lanthanides studied. The thermophysical properties of the eutectic composition with trivalent lanthanides were not significantly different from the pure eutectic salt composition, but several changes were noted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗