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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric C Geistfeld↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric Geistfeld↗
NASA NTRS1983

Voyager 2 observations of energetic particle variations in the Ganymede wake region - A possible acceleration mechanism

Voyager 2's passage through the downstream corotation region of Ganymede found disturbances in the field and particle environment. Large fluctuations in the intensities and energy spectra of ions in the 0.1-to 4.0-MeV interval were also observed with the low energy charged particle experiment. All ion species were evidently affected, up through medium (C, N, O) nuclei. Electrons intensities, while disturbed, did not behave similarly to the ions. These effects are probably associated with the wake produced by the absorption of Jovian magnetospheric corotating plasma by Ganymede. A two-dimensional theoretical model for perturbed field in the wake region is proposed and it is shown by direct numerical simulation of exact particle trajectories that such perturbations are capable of modulating the charged particle energies to about the observed amount for the observed magnetic and plasma parameters. Analysis shows that the energization process is dependent on the pitch angles of the energetic ions and the process itself is nonadiabatic (the Larmor radius is approximately the size of the disturbance region).

Tariq, G.F.↗
DOE OSTI2024

Quantum stabilization and flat hydrogen-based bands of nitrogen-doped lutetium hydride

We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.

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

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

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

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

Theory of tunneling between two-dimensional electron layers driven by spin pumping: Adiabatic regime and beyond

Tunneling spectroscopy between parallel two-dimensional (2D) electronic systems provides a powerful method to probe the underlying electronic properties by measuring tunneling conductance. Here, in this work, we present a theoretical framework for spin transport in 2D-to-2D tunneling systems, driven by spin pumping. This theory applies to a vertical heterostructure where two layers of metallic 2D electron systems are separated by an insulating barrier, with one layer exchange coupled to a magnetic layer driven at resonance. Utilizing a nonperturbative Floquet-Keldysh formalism, we derive general expressions for the tunneling spin and charge currents across a broad range of driving frequencies, extending beyond the traditional adiabatic pumping regime. At low frequencies, we obtain analytical results that recover the known behaviors in the adiabatic regime. However, at higher frequencies, our numerical findings reveal significant deviations in the dependence of spin and charge currents on both frequency and precession angle. This work offers fresh insights into the role of magnetization dynamics in tunneling transport, opening up new avenues for exploring nonadiabatic spin pumping phenomena.

Green's function methods↗
DOE OSTI2023

Polaron absorption in aligned conjugated polymer films: breakdown of adiabatic treatments and going beyond the conventional mid-gap state model

This study provides the first experimental polarized intermolecular and intramolecular optical absorption components of field-induced polarons in regioregular poly(3-hexylthiophene-2,5-diyl), rr-P3HT, a polymer semiconductor. Highly aligned rr-P3HT thin films were prepared by a high temperature shear-alignment process that orients polymer backbones along the shearing direction. rr-P3HT in-plane molecular orientation was measured by electron diffraction, and out-of-plane orientation was measured through series of synchrotron X-ray scattering techniques. Then, with molecular orientation quantified, polarized charge modulation spectroscopy was used to probe mid-IR polaron absorption in the ℏω = 0.075 — 0.75 eV range and unambiguously assign intermolecular and intramolecular optical absorption components of hole polarons in rr-P3HT. This data represents the first experimental quantification of these polarized components and allowed long-standing theoretical predictions to be compared to experimental results. The experimental data is discrepant with predictions of polaron absorption based on an adiabatic framework that works under the Born–Oppenheimer approximation, but the data is entirely consistent with a more recent nonadiabatic treatment of absorption based on a modified Holstein Hamiltonian. This nonadiabatic treatment was used to show that both intermolecular and intramolecular polaron coherence break down at length scales significantly smaller than estimated structural coherence in either direction. Importantly, this strongly suggests that polaron delocalization is fundamentally limited by energetic disorder in rr-P3HT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Charge carrier nonadiabatic dynamics in non-metal doped graphitic carbon nitride

Graphitic carbon nitride (GCN) has attracted significant attention due to its excellent performance in photocatalytic applications. Non-metal doping of GCN has been widely used to improve the efficiency of the material as a photocatalyst. Using a combination of time-domain density functional theory with nonadiabatic molecular dynamics, we study the charge carrier dynamics in oxygen and boron doped GCN systems. The reported simulations provide a detailed time-domain mechanistic description of the charge separation and recombination processes that are of fundamental importance while evaluating the photovoltaic and photocatalytic performance of the material. The appearance of smaller energy gaps due to the presence of dopant states improves the visible light absorption range of the doped systems. At the same time, the nonradiative lifetimes are shortened in the doped systems as compared to the pristine GCN. In the case of boron doped at a carbon (B–C–GCN), the charge recombination time is very long as compared to the other two doped systems owing to the smaller electron–phonon coupling strength between the valence band maximum and the trap state. The results suggest B–C–GCN as the most suitable candidate among three doped systems studied in this work for applications in photocatalysis. Finally, this work sheds light into the influence of dopants on quantum dynamics processes that govern GCN performance and, thus, guides toward building high-performance devices in photocatalysis.

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 OSTI2020

Pb dimerization greatly accelerates charge losses in MAPbI 3 : Time-domain ab initio analysis

Metal halide perovskites constitute a new type of semiconducting materials with long charge carrier lifetimes and efficient light-harvesting. The performance of perovskite solar cells and related devices is limited by nonradiative charge and energy losses, facilitated by defects. Combining nonadiabatic molecular dynamics and time-domain density functional theory, we demonstrate that charge losses depend strongly on the defect chemical state. By considering an extra Pb atom in CH 3 NH 3 PbI 3 , which is a common defect in lead halide perovskites, we investigate its influence on charge trapping and recombination. In a chemically inert form as a Pb interstitial, the extra Pb atom has only a mild influence on charge recombination. However, if the extra Pb atom binds to a native Pb atom to form a dimer, the charge trapping and recombination are greatly accelerated because the Pb-dimer creates a localized midgap trap state that couples strongly to the perovskite valence band edge. Holes disappear from the valence band two orders of magnitude faster than in the pristine perovskite and recombine with conduction band electrons one order of magnitude faster. The simulations identify the phonon modes involved in the nonradiative charge trapping and recombination and highlight the importance of rapid decoherence within the electronic subsystem for long carrier lifetimes. The detailed atomistic analysis of the charge trapping and recombination mechanisms enriches the understanding of defect properties and provides theoretical guidance for improving perovskite performance.

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

Suppressed superexchange interactions in the cuprates by bond-stretching oxygen phonons

We study a multiorbital Hubbard–Su-Schrieffer-Heeger model for the one-dimensional (1D) corner-shared cuprates in the adiabatic and nonadiabatic limits using exact diagonalization and determinant quantum Monte Carlo. At half filling and in the adiabatic limit, lattice dimerization can be achieved only over a narrow range of couplings slightly below a critical coupling g c . Beyond this value, the sign of the effective hopping changes, and the lattice becomes unstable. Strong lattice fluctuations replace the dimerization state in the nonadiabatic case. We also examine the model's temperature-dependent uniform and dynamical magnetic susceptibilities and compare them to the results of an effective spin-1/2 Heisenberg model. In doing so, we demonstrate that lattice fluctuations induced by the e-ph interaction suppress the effective superexchange interaction when g < g c . Furthermore, our results elucidate the effect of bond-stretching phonons in the parent cuprate compounds in general and are particularly relevant to 1D cuprates, where strong e-ph interactions have recently been inferred.

36 MATERIALS SCIENCE↗
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 OSTI2025

Excited-State Densities from Time-Dependent Density Functional Response Theory

While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.

approximation↗
DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
NASA NTRS1997

An Accurate Potential Energy Surface for H2O

We have carried out extensive high quality ab initio electronic structure calculations of the ground state potential energy surface (PES) and dipole moment function (DMF) for H2O. A small adjustment is made to the PES to improve the agreement of line positions from theory and experiment. The theoretical line positions are obtained from variational ro-vibrational calculations using the exact kinetic energy operator. For the lines being fitted, the root-mean-square error was reduced from 6.9 to 0.08 /cm. We were then able to match 30,092 of the 30,117 lines from the HITRAN 96 data base to theoretical lines, and 80% of the line positions differed less than 0.1 /cm. About 3% of the line positions in the experimental data base appear to be incorrect. Theory predicts the existence of many additional weak lines with intensities above the cutoff used in the data base. To obtain results of similar accuracy for HDO, a mass dependent correction to the PH is introduced and is parameterized by simultaneously fitting line positions for HDO and D2O. The mass dependent PH has good predictive value for T2O and HTO. Nonadiabatic effects are not explicitly included. Line strengths for vibrational bands summed over rotational levels usually agree well between theory and experiment, but individual line strengths can differ greatly. A high temperature line list containing about 380 million lines has been generated using the present PES and DMF

Schwenke, David W.↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗