Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum chemical 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

High field magnetometry with hyperpolarized nuclear spins

Abstract Quantum sensors have attracted broad interest in the quest towards sub-micronscale NMR spectroscopy. Such sensors predominantly operate at low magnetic fields. Instead, however, for high resolution spectroscopy, the high-field regime is naturally advantageous because it allows high absolute chemical shift discrimination. Here we demonstrate a high-field spin magnetometer constructed from an ensemble of hyperpolarized 13 C nuclear spins in diamond. They are initialized by Nitrogen Vacancy (NV) centers and protected along a transverse Bloch sphere axis for minute-long periods. When exposed to a time-varying (AC) magnetic field, they undergo secondary precessions that carry an imprint of its frequency and amplitude. For quantum sensing at 7T, we demonstrate detection bandwidth up to 7 kHz, a spectral resolution < 100mHz, and single-shot sensitivity of 410pT $$/\sqrt{{{{{{{{\rm{Hz}}}}}}}}}$$ / Hz . This work anticipates opportunities for microscale NMR chemical sensors constructed from hyperpolarized nanodiamonds and suggests applications of dynamic nuclear polarization (DNP) in quantum sensing.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Programmable simulations of molecules and materials with reconfigurable quantum processors

Simulations of quantum chemistry and quantum materials are believed to be among the most important applications of quantum information processors. However, realizing practical quantum advantage for such problems is challenging because of the prohibitive computational cost of programming typical problems into quantum hardware. Here we introduce a simulation framework for strongly correlated quantum systems represented by model spin Hamiltonians that uses reconfigurable qubit architectures to simulate real-time dynamics in a programmable way. Our approach also introduces an algorithm for extracting chemically relevant spectral properties via classical co-processing of quantum measurement results. We develop a digital–analogue simulation toolbox for efficient Hamiltonian time evolution using digital Floquet engineering and hardware-optimized multi-qubit operations to accurately realize complex spin–spin interactions. As an example, we propose an implementation based on Rydberg atom arrays. In addition, we show how detailed spectral information can be extracted from the dynamics through snapshot measurements and single-ancilla control, enabling the evaluation of excitation energies and finite-temperature susceptibilities from a single dataset. To illustrate the approach, we show how to use the method to compute key properties of a polynuclear transition-metal catalyst and two-dimensional magnetic materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

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

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

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

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Novel transition dynamics of topological solitons

Continuous phase transitions can be classified into ones characterized by local-order parameters and others that need additional topological constraints. The critical dynamics near the former transitions have been extensively studied, but the latter is less understood. We fill this gap in knowledge by studying the transition dynamics to a parity-breaking topological ground state called the chiral soliton lattice in quantum chromodynamics at finite temperature, baryon chemical potential, and an external magnetic field. We find a slowing down of the soliton’s translational motion as the critical magnetic field approaches, while the local relaxation rate to the solitonic state remains finite. The time required to converge to the stationary state strongly depends on the symmetry of the initial configuration, determining whether translational motion occurs during the dynamic process. Published by the American Physical Society 2024

Nishimura, Kentaro↗
DOE OSTI2020

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked TATB high explosive crystal. By scale-bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers. While shear bands quickly cool,they remain chemically activated and support increased reaction rates without the local heating typically evoked by the hot spot paradigm. We describe this phenomenon as chemical activation through shear banding

Chemistry - Inorganic↗
DOE OSTI2022

Dynamical electron correlation and the chemical bond. I. Covalent bonds in AH and AF (A = B–F)

Dynamical electron correlation has a major impact on the computed values of molecular properties and the energetics of molecular processes. This study focused on the effect of dynamical electron correlation on the spectroscopic constants ($R$ e , $ω$ e , $D$ e ), and potential energy curves, Δ$E$($R$), of the covalently bound AH and AF molecules, A = B–F. The changes in the spectroscopic constants (Δ$R$ e , Δ$ω$ e , Δ$D$ e ) caused by dynamical correlation are erratic and, at times, even surprising. These changes can be understood based on the dependence of the dynamical electron correlation energies of the AH and AF molecules as a function of the bond distance, i.e., Δ$E$ DEC ($R$). At large $R$, the magnitude of Δ$E$ DEC ($R$) increases nearly exponentially with decreasing $R$, but this increase slows as $R$ continues to decrease and, in many cases, even reverses at very short $R$. The changes in Δ$E$ DEC ($R$) in the region around $R$ e were as unexpected as they were surprising, e.g., distinct minima and maxima were found in the curves of Δ$E$ DEC ($R$) for the most polar molecules. The variations in Δ$E$ DEC ($R$) for $R$ ≲ $R$ e are directly correlated with major changes in the electronic structure of the molecules as revealed by a detailed analysis of the spin-coupled generalized valence bond wave function. In conclusion, the results reported here indicate that we have much to learn about the nature of dynamical electron correlation and its effect on chemical bonds and molecular properties and processes.

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

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Influence of pump laser fluence on ultrafast myoglobin structural dynamics

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump–probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore. As multiphoton absorption may force the protein response into non-physiological pathways, it is of great concern whether this experimental approach allows valid conclusions to be drawn vis-à-vis biologically relevant single-photon-induced reactions. Here we describe ultrafast pump–probe SFX experiments on the photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe–CO bond distance (predicted by recent quantum wavepacket dynamics) are seen to depend strongly on pump laser energy, in line with quantum chemical analysis. Our results confirm both the feasibility and necessity of performing ultrafast TR-SFX pump–probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing both the design and the interpretation of ultrafast TR-SFX pump–probe experiments such that mechanistically relevant insight emerges.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Shake-Down Spectroscopy as State- and Site-Specific Probe of Ultrafast Chemical Dynamics

Tracking the multifarious ultrafast electronic and structural changes occurring in a molecule during a photochemical transformation is a challenging endeavor that benefits from recent experimental and computational progress in time-resolved techniques. Measurements of valence electronic states, which provide a global picture of the bonding structure of the molecule, and core electronic states, which provide insight into the local environment, traditionally require different approaches and are often studied separately. Here, we demonstrate that X-ray pulses from a seeded free-electron laser (FEL) enable the measurement of high-resolution, time-resolved X-ray photoelectron spectra (XPS) that capture weak satellite states resulting from shake-down processes in a valence-excited molecule. This approach effectively combines the advantages of both valence- and core-state investigations. We applied this method to investigate photoexcited CS2 molecules, where the role of internal conversion (IC) and intersystem crossing (ISC) in determining the predissociation dynamics is controversial. We present XPS spectra from photoexcited CS2, obtained at the FERMI FEL. High-resolution measurements, compared to the corresponding spectra obtained from accurate multireference quantum chemical calculations, reveal that shake-down satellite channels are highly sensitive to both valence electronic and geometric changes. Previous studies of the predissociation dynamics have led to uncertain assignments of the branching between singlet and triplet excited states. We derive a propensity rule that demonstrates the spin-selectivity of the shake-downs. This selectivity allows us to unequivocally assign contributions from the bright and dark singlet excited states, with populations tracked along the predissociation dynamic pathway.

Thompson, Henry J. [Univ. of Southampton (United K↗
DOE OSTI2023

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca 2+ and Mg 2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca 2+ and Mg 2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca 2+ and Mg 2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca 2+ and Mg 2+ . Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Here, our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

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

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗
DOE OSTI2021

Experimental Studies of the Dynamics of Adsorbed and Condensed Phase Iodine Species Formed in Severe Reactor Accidents (Final Project Report)

We have studied the chemical transformations of semi-volatile organic iodide species (RI) in the condensed phase as well as adsorbed onto the surfaces of environmentally relevant particles, e.g., soot (carbon) and dust (silica) under real-world conditions. The combination of chosen experimental methods provide a molecular-level insight into the chemical dynamics of organic iodides. Throughout this project, we have developed new experimental, quantum chemistry, and data analysis capabilities. All experimental work was performed with 127 I.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Digital Twin for Chemical Science (DTCS) v0.01

Directly visualizing the trajectories of chemistry can unravel novel insights into the behavior of catalysts, gas phase reactions, photo-induced dynamics, and building blocks for quantum information processing. The ability of explicitly identifying, tracking, and tagging the exchange of matter, hence the annihilation and creation of new chemical species, can be best realized through a close coupling of theory and experiment. While the synchrotron-based characterization facilities propelled rapidly in its hardware, providing higher brightness, better resolution, and more precision, the software infrastructure is lagging. We developed DTCS (Digital Twin for Chemical Science) v.01, a central platform that faithfully mimics advanced instrumentations in Scientific User Facilities, by solving a variety of technical challenges in data acquisition, analysis, and model-driven interpretation. Rooted in physics and accelerated by AI, we validated this concept by direct comparison with precise experimental X-ray Photoelectron Spectroscopy (XPS) observations using a ubiquitous metal-water interfacial scenario, i.e., Ag/H2O as our main narrative. The DTCS v.01 input mirrors how the bench chemists work, with the output directly linked to the end station computer, thereby providing a user-friendly, knowledge-driven, and accessible user experience with mechanistic insights standardized in a way that are ready to be published, versioned, and transferred flexibly.

Qian, Jin↗
DOE OSTI2021

Large-Scale, GPU-Enhanced DFTB Approaches for Probing Multi-Component Alloys

The major goals of the project are to develop and apply GPU-enhanced density functional theory tight binding (DFTB) calculations for efficient simulations of structural materials of alloy systems. The figure below depicts the GPU algorithmic developments on the left that will be carried out, and the solid-state structures on the right are representative of the complex alloys that can be computed with the GPU-enhanced DFTB approach. The methodology and tools developed in this project encompass accurate intermolecular potentials and GPU enhancements to the density functional tight binding (DFTB) approach for high-throughput ab initio molecular dynamics calculations of multi-component alloys at elevated temperatures. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of large, multi-component alloys at the predictive quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum nature of chemical systems; however, these methods cannot tackle the large sizes relevant to these multi-component systems. The DFTB-based ab initio molecular dynamics approach (coupled with our in-house GPU capabilities for enhanced speed) used in this project provides a viable approach for probing these large systems at a quantum level of detail that is significantly faster than current first-principle methods Collectively, the capabilities developed in this project directly respond to DOE HBCU-OMI, AOI 2-1 initiatives by (1) enabling accurate and efficient predictions and (2) bringing a fundamental understanding of structural interactions in these complex systems at elevated temperatures.

20 FOSSIL-FUELED POWER PLANTS↗