Search NASA⌕ Search

SEARCH · Search NASA

Results for “electronically nonadiabatic”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Cluster Theory for Nonadiabatic Dynamics: Nuclear Gradients and Nonadiabatic Couplings in Similarity Constrained Coupled Cluster Theory

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801–4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1972

Trapping boundary and field-line motion during geomagnetic storms.

Observation that the high-latitude trapping boundary for 20-keV electrons and 100-keV protons became very thin in the early morning hours during two intense substorms. The gradients were too steep to be maintained by drifting particles, so they must have been produced locally over the nightside of the earth. The flux gradient is seen to move at speeds in excess of 100 km/sec. Plasma appears to move away from the tail and around the earth at these high speeds during the sudden expansion phases of the substorms. The rapid plasma motion requires the presence of fluctuating electric fields that sometimes exceed 50 to 100 mV/m at a geomagnetic latitude of 30 deg on the L = 5 field line. These observations fit best into a model that contains two field-aligned sheet currents. The high electric fields that accompany the rapid plasma flow can produce nonadiabatic acceleration of 0.1- to 1-MeV electrons and protons.

Kaufmann, R. L.↗
NASA NTRS1995

A study of weak anisotropy in electron pressure in the tail current sheet

We adopt a magnetotail model with stretched field lines where ion motions are generally nonadiabatic and where it is assumed that the pressure anisotropy resides only in the electron pressure tensor. We show that the magnetic field lines with p(perpendicular) greater than p(parallel) are less stretched than the corresponding field lines in the isotropic model. For p(parallel) greater than p(perpendicular), the magnetic field lines become more and more stretched as the anisotropy approaches the marginal firehose limit, p(parallel) = p(perpendicular) + B(exp 2)/mu(sub 0). We also show that the tail current density is highly enhanced at the firehose limit, a situation that might be subject to a microscopic instability. However, we emphasize that the enhancement in the current density is notable only near the center of the tail current sheet (z = 0). Thus it remains unclear whether any microscopic instability can significantly alter the global magnetic field configuration of the tail. By comparing the radius of the field-line curvature at z = 0 with the particle's gyroradius, we suspect that even the conventional adiabatic description of electrons may become questionable very close to the marginal firehose limit.

Lee, D.-Y.↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
NASA NTRS2012

Solar Wind Electron Interaction with the Dayside Lunar Surface and Crustal Magnetic Fields: Evidence for Precursor Effects

Electron distributions measured by Lunar Prospector above the dayside lunar surface in the solar wind often have an energy dependent loss cone, inconsistent with adiabatic magnetic reflection. Energy dependent reflection suggests the presence of downward parallel electric fields below the spacecraft, possibly indicating the presence of a standing electrostatic structure. Many electron distributions contain apparent low energy (<100 eV) upwardgoing conics (58% of the time) and beams (12% of the time), primarily in regions with non-zero crustal magnetic fields, implying the presence of parallel electric fields and/or wave-particle interactions below the spacecraft. Some, but not all, of the observed energy dependence comes from the energy gained during reflection from a moving obstacle; correctly characterizing electron reflection requires the use of the proper reference frame. Nonadiabatic reflection may also play a role, but cannot fully explain observations. In cases with upward-going beams, we observe partial isotropization of incoming solar wind electrons, possibly indicating streaming and/or whistler instabilities. The Moon may therefore influence solar wind plasma well upstream from its surface. Magnetic anomaly interactions and/or non-monotonic near surface potentials provide the most likely candidates to produce the observed precursor effects, which may help ensure quasi-neutrality upstream from the Moon.

Solar↗
DOE OSTI2020

Photoinduced Dynamics with Constrained Vibrational Motion: FrozeNM Algorithm

Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. In this paper, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.

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

How to calculate the rate constants for nonradiative transitions between the M S components of spin multiplets?

Predicting the rates of spin-dependent processes characterized by nonradiative transitions between electronic states with different spin multiplicities is important for understanding the mechanisms of many photochemical and catalytic reactions. To calculate these rates, it is necessary to define the spin state representation and the couplings between these states that drives the interstate transitions. In this work, we describe three different approaches to calculating the spin-orbit coupling (SOC), transition probabilities, and rate constants between the MS components of the electronic states with different spin multiplicities. We implemented these approaches in our nonadiabatic statistical theory (NAST) software package, which predicts the transition probabilities and rate constants of spin-dependent processes using information obtained from electronic structure calculations. Here, we discuss the advantage and drawbacks of each approach and, as an example, calculate the rate constants for transitions between the spin states in the active site model of the protein rubredoxin.

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

An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics

In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.

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

Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born–Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner–Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ 2 , where κ = (m/M)1/4 is the usual Born–Oppenheimer parameter. The lowest order term is the usual Born–Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born–Oppenheimer theory. Here, we show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born–Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

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

Unveiling the spatial distribution of molecular coherences at conical intersections by covariance X-ray diffraction signals

Significance Monitoring the real-space motion of the electrons directly involved in ultrafast molecular dynamics is a challenging goal, hindered by the strong contributions of the remaining inactive electrons. Molecular quantum coherences can provide a direct window into elementary electronic dynamics, but they could not be separately imaged so far. Here, we demonstrate a time- and frequency-resolved X-ray diffraction technique, which can employ existing stochastic X-ray pulses to separately image the coherence contributions. The technique can selectively follow the relevant electrons involved in the ultrafast dynamics and identify the confined molecular regions covered by their motion. It enables readily feasible measurements and can be straightforwardly extended to full three-dimensional imaging via X-ray and electron diffraction.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Ultrafast Adsorbate Excitation Probed with Subpicosecond-Resolution X-Ray Absorption Spectroscopy

We use a pump-probe scheme to measure the time evolution of the C K-edge x-ray absorption spectrum from CO/Ru(0001) after excitation by an ultrashort high-intensity optical laser pulse. Because of the short duration of the x-ray probe pulse and precise control of the pulse delay, the excitation-induced dynamics during the first picosecond after the pump can be resolved with unprecedented time resolution. By comparing with density functional theory spectrum calculations, we find high excitation of the internal stretch and frustrated rotation modes occurring within 200 fs of laser excitation, as well as thermalization of the system in the picosecond regime. The ~100 fs initial excitation of these CO vibrational modes is not readily rationalized by traditional theories of nonadiabatic coupling of adsorbates to metal surfaces, e.g., electronic frictions based on first order electron-phonon coupling or transient population of adsorbate resonances. Finally, we suggest that coupling of the adsorbate to nonthermalized electron-hole pairs is responsible for the ultrafast initial excitation of the modes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

In Silico Discovery of Multistep Chemistry Initiated by a Conical Intersection: The Challenging Case of Donor–Acceptor Stenhouse Adducts

Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ~2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Finally, our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.

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

The non-adiabatic nanoreactor: towards the automated discovery of photochemistry

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S 0 /S 1 and S 1 /S 2 conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

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

Adiabatic criteria for outer-sphere bimolecular electron-transfer reactions

A model is presented for outer-sphere bimolecular electron-transfer reactions which is correct in the adiabatic, nonadiabatic, and intermediate dynamical regimes for an overdamped solvent coordinate. From this model, the conditions for the transfer to be adiabatic or nonadiabatic are deduced. The time-scale separations needed to adequately describe the process as an average over (distant dependent) unimolecular rates are described.

Onuchic, Jose Nelson↗
DOE OSTI2020

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

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

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

Nonadiabatic coupling of the dynamical structure to the superconductivity in YSr 2 Cu 2.75 Mo 0.25 O 7.54 and Sr 2 CuO 3.3

A crucial issue in cuprates is the extent and mechanism of the coupling of the lattice to the electrons and the superconductivity. Here we report Cu K edge extended X-ray absorption fine structure measurements elucidating the internal quantum tunneling polaron (iqtp) component of the dynamical structure in two heavily overdoped superconducting cuprate compounds, tetragonal YSr 2 Cu 2.75 Mo 0.25 O 7.54 with superconducting critical temperature, T c = 84 K and hole density p = 0.3 to 0.5 per planar Cu, and the tetragonal phase of Sr 2 CuO 3.3 with T c = 95 K and p = 0.6. In YSr 2 Cu 2.75 Mo 0.25 O 7.54 changes in the Cu-apical O two-site distribution reflect a sequential renormalization of the double-well potential of this site beginning at T c , with the energy difference between the two minima increasing by ~6 meV between T c and 52 K. Sr 2 CuO 3.3 undergoes a radically larger transformation at T c , >1-Å displacements of the apical O atoms. The principal feature of the dynamical structure underlying these transformations is the strongly anharmonic oscillation of the apical O atoms in a double-well potential that results in the observation of two distinct O sites whose Cu–O distances indicate different bonding modes and valence-charge distributions. The coupling of the superconductivity to the iqtp that originates in this nonadiabatic coupling between the electrons and lattice demonstrates an important role for the dynamical structure whereby pairing occurs even in a system where displacements of the atoms that are part of the transition are sufficiently large to alter the Fermi surface. The synchronization and dynamic coherence of the iqtps resulting from the strong interactions within a crystal would be expected to influence this process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗