Search NASA⌕ Search

SEARCH · Search NASA

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

UQ4QM: Uncertainty Quantification for Quantum Materials (Chemical Dynamics Initiative Final Report)

The Uncertainty Quantification for Quantum Materials (UQ4QM) LDRD project focused on developing solutions for strategic two project areas the Chemical Dynamics Initiative (CDi) Use Case 3. This included data-driven approaches toward understanding heterogeneous data, such as multi-fidelity data or data with unknown spatial or temporal perturbations. The applications were varied but the project aims pursued solutions which were amenable to common UQ approaches in order to maximize impact.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗
DOE OSTI2024

Polariton-assisted incoherent to coherent excitation energy transfer between colloidal nanocrystal quantum dots

Here, we explore the dynamics of energy transfer between two nanocrystal quantum dots placed within an optical microcavity. By adjusting the coupling strength between the cavity photon mode and the quantum dots, we have the capacity to fine-tune the effective coupling between the donor and acceptor. Introducing a non-adiabatic parameter, γ, governed by the coupling to the cavity mode, we demonstrate the system’s capability to shift from the overdamped Förster regime (γ $\ll$ 1) to an underdamped coherent regime (γ $\gg$ 1). In the latter regime, characterized by swift energy transfer rates, the dynamics are influenced by decoherence time. To illustrate this, we study the exciton energy transfer dynamics between two closely positioned CdSe/CdS core/shell quantum dots with sizes and separations relevant to experimental conditions. Employing an atomistic approach, we calculate the excitonic level arrangement, exciton–phonon interactions, and transition dipole moments of the quantum dots within the microcavity. These parameters are then utilized to define a model Hamiltonian. Subsequently, we apply a generalized non-Markovian quantum Redfield equation to delineate the dynamics within the polaritonic framework.

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

Geometric phase in coupled cluster theory

It has been well-established that the topography around conical intersections between excited electronic states is incorrectly described by coupled cluster and many other single reference theories (the intersections are “defective”). Despite this, we show both analytically and numerically that the geometric phase effect (GPE) is correctly reproduced upon traversing a path around a defective excited-state conical intersection (CI) in coupled cluster theory. The theoretical analysis is carried out by using a non-Hermitian generalization of the linear vibronic coupling approach. Interestingly, the approach qualitatively explains the characteristic (incorrect) shape of the defective CIs and CI seams. Moreover, the validity of the approach and the presence of the GPE indicate that defective CIs are local (and not global) artifacts. This implies that a sufficiently accurate coupled cluster method could predict nuclear dynamics, including geometric phase effects, as long as the nuclear wavepacket never gets too close to the conical intersections.

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

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

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

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li( 2 S) + NaLi( a 3 Σ + ) → Na( 2 S) + Li 2 ( a 3 Σ u + )

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. Here we perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a 3 Σ + ) and Li( 2 S) to form Li 2 (a 3 Σ u + ) and Na( 2 S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li 2 Na interaction at short range, which we attribute to a relatively small number of open Li 2 (a 3 Σ u + ) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗
DOE OSTI2024

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

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

ZnO quantum dot–molecule conjugates: Chemical interactions, charge dynamics, and spin polarization

Conjugates between molecules and quantum dots (QDs) have been explored for a range of potential applications from photocatalysis and photovoltaics to quantum information science technologies. A particularly ubiquitous material in many of these applications are ZnO QDs since they can accept and transport electrons and can also act as hosts for unique spin states. Conjugates between molecular light absorbers and ZnO QDs have been explored for decades as components in dye-sensitized solar cells. Recently, these materials have also attracted interest for their ability to produce spin-polarized states upon photoexcitation. The current paper employs a series of light absorbing perylene molecules with different ZnO QD sizes to explore key features of these QD–molecule conjugates: (1) chemical interactions, (2) charge dynamics, and (3) spin polarization. The chemical interactions between the molecules and QDs are determined with binding equilibria and reveal dramatic impact of ligand size. The charge transfer dynamics from photoexcited perylenes to ZnO QDs were found to depend exponentially on the linker length. Finally, time-resolved electron paramagnetic resonance experiments reveal that these conjugates generate spin-polarized states in the form of radical pairs and triplets. These spin states hold promise as potential qubits and also offer an avenue to efficiently sensitize molecular triplets.

Hernandez, Frida S. [Amherst College, MA (United S↗
DOE OSTI2024

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

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

A fractional calculus framework for open quantum dynamics: From Liouville to Lindblad to memory kernels

Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini–Kossakowski–Sudarshan–Lindblad equation uniquely characterizes Markovian completely positive and trace-preserving (CPTP) evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here, we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner–Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima–Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for modeling non-Markovian quantum dynamics in chemical physics and physical chemistry, providing a CPTP-preserving, computationally efficient surrogate for structured condensed-phase environments where long-time memory and dissipation play a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Fast methods for multisite charge transfer processes. I. Constrained, state averaged CASSCF(1,n) and CASSCF(2n − 1,n) simulations

We design a dynamically weighted state-averaged constrained complete active space self-consistent field (DW-SA-cCASSCF) algorithm to treat electrons or holes moving between n molecular fragments (where n can be larger than 2). Within such a so-called eDSCn/hDSCn approach, we consider configurations that are mutually single excitations of each other, and we apply a generalized set of constraints to tailor the method for studying charge transfer problems. The constrained optimization problem is efficiently solved using a DIIS-SQP algorithm, thus maintaining computational efficiency. We demonstrate the method for a finite Su–Schrieffer–Heeger chain, successfully reproducing the expected exponential decay of diabatic couplings with distance. When combined with a gradient, the current extension immediately enables efficient nonadiabatic dynamics simulations of complex multi-state charge transfer processes.

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

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2024

Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S 2 ). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. Herein we use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S 2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗