Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonadiabatic 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

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I – (aq), but this initial step in the genesis of e – (aq) has never been simulated directly using ab initio molecular dynamics. Here, we report the first such simulations, combining ground- and excited-state simulations of I – (aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e – (aq) starting from the CTTS state: I – (aq) + hν → I –* (aq) → I • (aq) + e – (aq). Notably, the equilibrated ground state of e – (aq) evolves from I –* (aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Carter-Fenk, Kevin↗
DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

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

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Charge-Transfer Luminescence in a Molecular Donor–Acceptor Complex: Computational Insights

Donor–acceptor molecular complexes are a popular class of materials utilizing charge-transfer states for practical applications. A recent class of donor–acceptor dyads based on the fluorescent BODIPY functionalized with triphenylamine (TPA) shows the peculiar property of dual fluorescence. It is hypothesized that instead of the sensitized charge-transfer state being optically dark, it provides an additional bright radiative pathway. In this study, we use time-dependent density functional theory to characterize the energetic alignment of excitonic and charge-transfer states in a BODIPY-TPA molecular complex. We observe that using a long-range exchange corrected functional in combination with state-specific solvation scheme gives a qualitatively correct alignment of the exciton and charge-transfer states and an enhancement in oscillator strength for the equilibrium solvated charge-transfer state, in agreement with experiment. Furthermore, this work provides rationalization of charge-transfer state emission and provides a foundation to explore charge-transfer using ab initio excited-state nonadiabatic dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Computational Framework for Simulations of Dissipative Nonadiabatic Dynamics on Hybrid Oscillator-Qubit Quantum Devices

Here, we introduce a computational framework for simulating nonadiabatic vibronic dynamics on circuit quantum electrodynamics (cQED) platforms. Our approach leverages hybrid oscillator-qubit quantum hardware with midcircuit measurements and resets, enabling the incorporation of environmental effects such as dissipation and dephasing. To demonstrate its capabilities, we simulate energy transfer dynamics in a triad model of photosynthetic chromophores inspired by natural antenna systems. We specifically investigate the role of dissipation during the relaxation dynamics following photoexcitation, where electronic transitions are coupled to the evolution of quantum vibrational modes. Our results indicate that hybrid oscillator-qubit devices, operating with noise levels below the intrinsic dissipation rates of typical molecular antenna systems, can achieve the simulation fidelity required for practical computations on near-term and early fault-tolerant quantum computing platforms.

Hamiltonians↗
DOE OSTI2023

Probing C–I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C 4 H 3 IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Furthermore, gradual absorption changes up to ~0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n ⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C–I bond, enabling a nonadiabatic transition to a dissociative π→σ C-I * state. For the subsequent fragmentation, a relatively narrow bond-length region along the C–I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

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

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗
DOE OSTI2024

Inter-Oligomer Interaction Influence on Photoluminescence in Cis-Polyacetylene Semiconductor Materials

Semiconducting conjugated polymers (CPs) are pivotal in advancing organic electronics, offering tunable properties for solar cells and field-effect transistors. Here, we carry out first-principle calculations to study individual cis-polyacetylene (cis-PA) oligomers and their ensembles. The ground electronic structures are obtained using density functional theory (DFT), and excited state dynamics are explored by computing nonadiabatic couplings (NACs) between electronic and nuclear degrees of freedom. We compute the nonradiative relaxation of charge carriers and photoluminescence (PL) using the Redfield theory. Our findings show that electrons relax faster than holes. The ensemble of oligomers shows faster relaxation compared to the single oligomer. The calculated PL spectra show features from both interband and intraband transitions. The ensemble shows broader line widths, redshift of transition energies, and lower intensities compared to the single oligomer. This comparative study suggests that the dispersion forces and orbital hybridizations between chains are the leading contributors to the variation in PL. It provides insights into the fundamental behaviors of CPs and the molecular-level understanding for the design of more efficient optoelectronic devices.

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

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carrier Localization and Spontaneous Formation of Two-Dimensional Polarization Domain in Halide Perovskites

Halide perovskites are known for their rich phase diagram and superior performance in diverse optoelectronics applications. The latter property is often attributed to the long electron-hole recombination time, whose underlying physical mechanism has been a long-standing controversy. In this Letter, we investigate the transport and localization properties of electron and hole carriers in a prototypical halide perovskite (CsPbBr 3 ), through ab initio tight-binding nonadiabatic dynamics approach for large-scale (tens of nm size) supercell calculations, to simulate electron and ion dynamics on the same footing. Here, we found distinct structural, lattice polarization, and electron-phonon coupling properties at low (below 100 K) and high temperatures, consistent with experimental observations. In particular, at low temperature we find spontaneous formation of polar grain boundaries in the nonpolar bulk systems, which result in two-dimensional polarization patterns that serve to localize and separate electrons and holes. We reveal phonon-assisted variable-range hopping mostly responsible for low-temperature transport, and their characteristic frequency correlates with temperature-dependent phonon power spectrum and energy oscillation frequency in nonadiabatic dynamics. We answer the critical questions of long electron-hole recombination lifetime at low temperature and offer the correlation among polarization domains, electron-phonon couplings, and photocarrier dynamics.

Carrier generation & recombination↗
DOE OSTI2021

Vibronic Coupling and Exciton Chirality: Electronic and Structural Rearrangement between Helical to Zero Momentum Molecular Exciton States

Helical porphyrin aggregates are attractive macromolecular nanostructures for biomimetic light-harvesting and energy transfer in solar energy technology due to their large light absorbing cross sections and promising energy conservation benefits from chiral-induced spin selectivity. However, with these soft materials, it is imperative that we understand how the nuclear degrees of freedom impact the excitonic energy landscape and dynamics, particularly pertaining to how it influences the chiral angular momentum of the excitons. To this end, we have measured time-resolved depolarization ratios using femtosecond stimulated Raman spectroscopy to uncover the vibronic coupling guided by molecular vibrations between excitons of different angular momentum in helical tetra(sulfonatophenyl)porphyrin aggregates. We find that while transient absorption anisotropy reveals rapid (1 ps lifetime) exciton rotation from helical to achiral states in the Q-band, time-resolved vibrational depolarization ratios evolve on remarkably faster time scales (<500 fs lifetime) for the totally symmetric 1530 cm –1 vibration and the nontotally symmetric 1540 cm –1 vibration. The time-resolved depolarization ratios of the 1540 cm –1 vibration, in particular, show strong evidence of vibronic coupling and nonadiabatic exciton behavior that mediates the nonradiative transition between chiral and achiral excitons. Furthermore, our study shows that for rational design of helical molecular aggregates for exciton transport, the vibronic coupling between excitons of different angular momenta driven by molecular vibrations must be considered especially in the case where chiral-induced spin selectivity is desired.

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

Exploring the Unusual Reactivity of the Hydrated Electron with CO 2

Many questions remain about the reactions of the hydrated electron despite decades of study. Of particular note is the fact that they do not appear to follow the Marcus theory of electron transfer reactions, a feature that has yet to be explained. To investigate these issues, we use ab initio molecular dynamics (AIMD) simulations to investigate a one of the better studied reactions, the hydrated electron reduction of CO 2 . The rate constant for the hydrated electron-CO 2 reaction complex to react to form CO 2 - is, for the first time, estimated from AIMD simulations. Results at 298 and 373 K show the rate constant is insensitive to temperature, consistent with the low measured activation energy for the reaction, and the implications of this behavior are examined. The sampling provided by the simulations yields insight into the reaction mechanism. The reaction is found to involve both solvent reorganization and changes in the carbon dioxide structure. The latter lead to significant vibrational excitation of the bending and symmetric stretch vibrations in the CO 2 - product, indicating the reaction is vibrationally nonadiabatic. The former is estimated from calculation of an approximate collective solvent coordinate and the free energy in this coordinate is determined. Furthermore, these results indicate that AIMD simulations can reasonably estimate hydrated electron reaction activation energies and provide new insight into the mechanism that can help illuminate the features of this unusual chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nonadiabatically Driven Quantum Interference Effects in the Ultracold K + KRb → Rb + K 2 Chemical Reaction

The K + KRb → Rb + K 2 chemical reaction is the first ultracold atom–diatom chemical reaction for which experimental results have been reported for temperatures below 1 μK more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding nonadiabatic effects, has reported a rate coefficient that is about 35% below the experimental value. Here, we report the first nonadiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

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

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

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

Photoselective isotope fractionation dynamics of N 2 with cosmo and atmospheric chemistry perspectives

Stable isotope ratio measurements provide valuable insights into a broad range of natural processes, from planetary atmospheres and climate to interstellar chemistry. Nitrogen, which has two stable isotopes, exhibits varying isotope ratios across the solar system. To model these observations, the isotope fraction as a function of energy is essential. At the Advanced Light Source (ALS), we measured the photodissociation of molecular nitrogen (N 2 ) with vacuum UV photons where a single photon is sufficiently energetic to dissociate the strong bond. The nitrogen atoms produced are scavenged with H 2 to form ammonia, whose isotopic makeup is determined. Blending the experiments with dynamical computations that include the shielding of light, we examine the isotopic composition and electronic atomic states produced. The measured photodissociation of N 2 at a natural isotopic composition with a frequency broad light beam exceptionally strongly favors the formation of the heavier nitrogen isotope, 15 N. Computations concur and suggest that the maximum in the quantum yield reflects significant variations in the specific electronic quantum states of the product N atoms that have quite different reactivities. Our quantum computations show that at similar energies, photodissociation of 14 N 14 N and 15 N 14 N can lead to different product channels. The computed dynamics include extensive state-selective spin–orbit and nonadiabatic couplings affecting the light absorption and dissociation pathways that proceed via the triplet manifold of states. Our results are relevant for future exploration missions, both in situ and sample-return and for other molecules such as O 2 and CO.

Komarova, Ksenia [Hebrew Univ. of Jerusalem (Israe↗
DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗