Search NASA⌕ Search

SEARCH · Search NASA

Results for “diabatic states”

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 OSTI2020

On the nonadiabatic collisional quenching of OH(A) by H 2 : a four coupled quasi-diabatic state description

A four-state diabatic potential energy matrix (DPEM), H d , for the description of the nonadiabatic quenching of OH(A 2 Σ + ) by collisions with H 2 is reported. The DPEM is constructed as a fit to adiabatic energies, energy gradients, and derivative couplings obtained exclusively from multireference configuration interaction wave functions. A four-adiabatic-electronic-state representation is used in order to describe all energetically accessible regions of the nuclear coordinate space. Partial permutation-inversion symmetry is incorporated into the representation. Here, the fit is based on electronic structure data at 42 882 points, described by over 1.6 million least squares equations with a root mean square (mean unsigned) error of 178(83) cm –1 . Comparison of ab initio and H d determined minima, saddle points, and energy minimized points on C 2v , C s , C ∞v , and C 1 (noncoplanar) portions of two conical intersection seams are used to establish the accuracy of the H d .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗
DOE OSTI2022

Diabatic Valence-Hole States in the C 2 Molecule: “Putting Humpty Dumpty Together Again”

Despite the long history of spectroscopic studies of the C$_2$ molecule, fundamental questions about its chemical bonding are still being hotly debated. The complex electronic structure of C$_2$ is a consequence of its dense manifold of near-degenerate, low-lying electronic states. A global multi-state diabatic model is proposed here to disentangle the numerous configuration interactions within four symmetry manifolds of C$_2$ ($^{1}\Pi_g$, $^{3}\Pi_g$, $^{1}\Sigma_u^+$, and $^{3}\Sigma_u^+$). The key concept of our model is the existence of two "valence-hole" configurations, $2\sigma_g^22\sigma_u^11\pi_{u}^33\sigma_g^2$ for $^{1,3}\Pi_g$ states and $2\sigma_g^22\sigma_u^11\pi_{u}^43\sigma_g^1$ for $^{1,3}\Sigma_u^+$ states that derive from $3\sigma_g\leftarrow2\sigma_u$ electron promotion. The lowest-energy state from each of the four C$_2$ symmetry species is dominated by this type of valence-hole configuration at its equilibrium internuclear separation. Finally, as a result of their large binding energy (nominal bond order of 3) and correlation with the 2s$^2$2p$^2$+2s2p$^3$ separated-atom configurations, the presence of these valence-hole configurations has a profound impact on the $global$ electronic structure and unimolecular dynamics of C$_2$.

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

PyCDFT: A Python package for constrained density functional theory

In this paper, we present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.

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

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

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 OSTI2022

Simulation of the photodetachment spectra of the nitrate anion (NO 3 ₋ ) in the B ~ 2 E' energy range and non-adiabatic electronic population dynamics of NO 3

The photodetachment spectrum of the nitrate anion (NO 3 ₋ ) in the energy range of the NO 3 second excited state is simulated from first principles using quantum wave packet dynamics. The prediction at 10 K and 435 K relies on the use of an accurate full-dimensional fully coupled five state diabatic potential model utilizing an artificial neural network. The ability of this model to reproduce experimental spectra was demonstrated recently for the lower energy range. Analysis of the spectra indicates a weaker Jahn–Teller coupling compared to the first excited state. The detailed non-adiabatic dynamics is studied by computing the population dynamics. An ultra-fast non-statistical radiationless decay is found only among the Jahn–Teller components, which is followed by a slow statistical non-radiative decay among the different state manifolds. The latter is reproduced perfectly by a simple first order kinetics model. The dynamics in the second excited state is not affected by the presence of a conical intersection with the first excited state manifold.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Non-adiabatic quantum interference and complex formation in ultracold collisions of Rb with KRb

Ultracold elastic collisions of 87 Rb with 40 K 87 Rb in its ground vibrational and rotational state are investigated using a first principles based theoretical methodology. Full-dimensional ab initio computed potential energy surfaces are reported that include the two lowest-lying electronic states, their conical intersection, non-adiabatic couplings and an accurate long-range behavior. A numerically exact time-independent quantum dynamics method in hyperspherical coordinates is used to compute the elastic scattering cross sections, rate coefficients and collision lifetime spectrum. The quantum scattering calculations include all degrees of freedom and treat both electronic states and their non-adiabatic couplings using a two-state diabatic representation. The theoretically computed elastic rate coefficient is in good agreement with the recently reported experimental value. Significant non-adiabatic quantum interference effects are shown to originate from the unique properties of ultracold collisions and the geometric phase associated with the conical intersection. A high-resolution collision energy grid is used to investigate the origin of the experimentally reported long-lived 3-body collision complexes.

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

Atmospheric Diabatic Heating in Different Weather States and the General Circulation

Analysis of multiple global satellite products identifies distinctive weather states of the atmosphere from the mesoscale pattern of cloud properties and quantifies the associated diabatic heating/cooling by radiative flux divergence, precipitation, and surface sensible heat flux. The results show that the forcing for the atmospheric general circulation is a very dynamic process, varying strongly at weather space-time scales, comprising relatively infrequent, strong heating events by ''stormy'' weather and more nearly continuous, weak cooling by ''fair'' weather. Such behavior undercuts the value of analyses of time-averaged energy exchanges in observations or numerical models. It is proposed that an analysis of the joint time-related variations of the global weather states and the general circulation on weather space-time scales might be used to establish useful ''feedback like'' relationships between cloud processes and the large-scale circulation.

Cooling↗
DOE OSTI2022

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. Here, we detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

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

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

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

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2021

Clock transitions guard against spin decoherence in singlet fission

Short coherence times present a primary obstacle in quantum computing and sensing applications. In atomic systems, clock transitions (CTs), formed from avoided crossings in an applied Zeeman field, can substantially increase coherence times. We show how CTs can dampen intrinsic and extrinsic sources of quantum noise in molecules. Conical intersections between two periodic potentials form CTs in electron paramagnetic resonance experiments of the spin-polarized singlet fission photoproduct. We report on a pair of CTs for a two-chromophore molecule in terms of the Zeeman field strength, molecular orientation relative to the field, and molecular geometry.

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

Collisional quenching of O(1D) by O(3P)

Metastable O(1D) atoms may be quenched in collisions with ground state O(3P) atoms by transitions in the avoided crossing regions of the three lowest 3Pi(g) states of O2 of which the lowest separates to O(3P) + O(3P) and the two upper to O(3P) + O(1D). Quantal calculations of the adiabatic potential energy curves of the 3Pi(g) states are carried out with particular attention to an avoided crossing region in the lowest two states around a nuclear separation of 3.2a(0). Diabatic potential matrix elements are constructed from the adiabatic curves by imposing the requirement that they be smooth everywhere. A multi-state diabatic formulation is used to describe the scattering and the cross-sections for the collision-induced quenching of O(1D) atoms are calculated.

Yee, J.-H.↗
NASA NTRS1990

Theoretical study of the spectroscopy of NO(+)

The spectroscopy of the NO(+) molecule has been studied using multireference configuration-interaction calculations. The peak at 21.7 eV in the photoelectron spectrum of NO is assigned to the (2) 3Pi state (c 3Pi) and the vibrational progression centered at 23.1 eV is assigned to the diabatic (3) 1Pi state (B 1Pi). Several additional bound states with binding energies in excess of 1 eV are characterized in the present work. Radiative lifetimes are presented for many of the states.

Partridge, Harry↗
NASA NTRS2003

Diabatic Initialization of Mesoscale Models in the Southeastern United States: Can 0 to 12h Warm Season QPF be Improved?

It is well known that numerical warm season quantitative precipitation forecasts lack significant skill for numerous reasons. Some are related to the model--it may lack physical processes required to realistically simulate convection or the numerical algorithms and dynamics employed may not be adequate. Others are related to initialization-mesoscale features play an important role in convective initialization and atmospheric observation systems are incapable of properly depicting the three-dimensional stability structure at the mesoscale. The purpose of this study is to determine if a mesoscale model initialized with a diabatic initialization scheme can improve short-term (0 to 12h) warm season quantitative precipitation forecasts in the Southeastern United States. The Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory is used to diabatically initialize the Pennsylvania State University/National center for Atmospheric Research (PSUNCAR) Mesoscale Model version 5 (MM5). The SPORT Center runs LAPS operationally on an hourly cycle to produce analyses on a 15 km covering the eastern 2/3 of the United States. The 20 km National Centers for Environmental Prediction (NCEP) Rapid Update Cycle analyses are used for the background fields. Standard observational data are acquired from MADIS with GOES/CRAFT Nexrad data acquired from in-house feeds. The MM5 is configured on a 140 x 140 12 km grid centered on Huntsville Alabama. Preliminary results indicate that MM5 runs initialized with LAPS produce improved 6 and 12h QPF threat scores compared with those initialized with the NCEP RUC.

Lapenta, William M.↗