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 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 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↗
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 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↗
DOE OSTI2022

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin–orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3 (A) leading to the NH(X 3 Σ – , a 1 Δ) + H 2 channel

Here we simulate the photodissociation of NH 3 originating from its first excited singlet state S 1 into the NH 2 + H (radical) and NH + H 2 (molecular) channels. The states considered are the ground singlet state S 0 , the first excited singlet state S 1 and the lowest-lying triplet state T 1 , which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S 0 , S 1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T 1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H 2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH 2 + H channel. The NH(X 3 Σ_) radical produced through the intersystem crossing from S 0 to T 1 is rare (~0.2%) compared to NH(a 1 Δ) due to the process being spin forbidden.

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

The Low-Lying Electronic States of NO 2 : Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum

Nitrogen dioxide, NO 2 , is a free radical composed of the two most abundant elements in Earth’s atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner–Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N( 4 S u ) and O( 3 P g ) atoms. In this work, three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four ($\tilde X$, $\tilde A$, $\tilde B$, and $\tilde C$) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground $\tilde X$ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm –1 above the zero-point energy) have a root-mean-squared error of 16.5 cm –1 . In addition, dipole moments for each of the lowest four electronic states—and the transition dipoles between them—were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ~1.5% with no empirical adjustments.

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

Shuttling Majorana zero modes in disordered and noisy topological superconductors

The braiding of Majorana zero modes (MZMs) forms the fundamental building block for topological quantum computation. Braiding protocols which involve the physical exchange of MZMs are typically envisioned on a network of topological superconducting wires. An important component of these protocols is the transport of MZMs, which can be performed by using electric gates to locally tune sections of the wire between topologically trivial and nontrivial phases. In this work, we numerically simulate this transport by tuning a single section of a superconducting wire which contains either disorder (uncorrelated and correlated) or noise. We focus on the impact of these additional effects on the diabatic error, which describes unwanted transitions between the ground state and excited states. We show that the behavior of the average diabatic error is predominantly controlled by the statistics of the minimum bulk energy gap which is suppressed in the presence of disorder. The increase in diabatic error can be several orders of magnitude and is most deleterious when the disorder correlation length is a finite fraction of the transport distance and negligible when these lengths are far apart. In the presence of noise, the diabatic error is significantly enhanced due to optical transitions which depend on the minimum bulk energy gap as well as the frequency modes present in the noise. The results presented here serve to further characterize the diabatic error in disordered and noisy settings, which are important considerations in practical implementations of physical braiding schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

PotLib 2023: New version of a potential energy surface library for chemical systems

POTLIB is a library of global and semiglobal potential energy surface subprograms. The library currently features 410 entries, including both single-state entries and multi-state entries. When one calls the routine of a single-state entry, it returns the ground-electronic-state adiabatic potential energy surface at the input geometry. In addition, some entries also return the gradient of the surface. Further, when one calls a multi-state entry, it returns a diabatic potential energy matrix (DPEM). If the entry also has the gradient of the DPEM, one can compute adiabatic surfaces, their gradients, and the nonadiabatic coupling vectors (NACs) from the DPEM and its gradient by diagonalization. Some but not all the routines conform to one of a set of standard interfaces. The goal is to facilitate chemical dynamics research by collecting and disseminating a comprehensive collection of state-of-the-art potential energy routines (developed by a wide, international group of researchers) with systematic and well-defined interfaces for use with chemical dynamics programs. Systems in the library include CHArO 2 , CHN 2 O + , CH 2 O, CH 2 O 2 , CH 3 N 2 , CH 3 O, CH 4 , CH 4 Br, CH 4 Cl, CH 4 F, CH 4 O, CH 4 OCl, CH 4 OF, CH 4 OH, CH 5 , CH$^{+}_{5}$, CH 5 N, CH 5 O 2 , CH 8 O 2 , C 2 H 2 N 2 O, C 2 H 2 O, C 2 H 4 N, C 2 H 4 O 4 , C 2 H 6 Cl, C 2 H 6 F, C 2 H 6 H, C 2 H 6 O, C 2 H 6 OH, C 2 H 7 , C 2 O 2 , C 3 H 4 O 2 , C 3 H 7 NO, C 6 H 6 O, C 6 H 6 S, C 7 H8S, HBrCl, HCl 2 , HF 2 , HI 2 , HLiF, HNaF, HO 2 , HO 3 , HOBr, HSiO, H 2 Br, H 2 ClO, H 2 F, H 2 F 2 , H 2 FO, H 2 Na, H 2 O, H 2 O 2 , H 2 OBr, H 3 , H 3 Cl, H 3 ClN, H 3 ClO, H 3 N, H 3 O, H 3 O 2 , H$^{+}_{3}$, H 3 S, H 4 ClSi, H 4 N, H 4 NO, H 4 O 2 , H 5 GeO, H 5 Si, H$^{+}_{7}$, Al m H n , Al n , Ar N O, K 2 Rb 2 , NO 2 , N 2 O, N 2 O 2 , N 3 , N 4 , O 3 , and O 4 .

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

Air-sea interactions yielding rapid Beaufort Sea ice losses during the 2021 ONR THINICE Pilot Field Campaign

Seasonal sea ice coverage has declined in the Pacific Arctic region. Losses of the ice cover have modified the region’s surface heat and momentum fluxes providing weather systems more direct interaction with the upper ocean. For example, recent cyclones that have passed over these waters have initiated multiple air-sea interactions and have impacted mixed layer depth, wave properties, and sea ice mass budgets due to diabatic and/or advective processes. As a United States Office of Naval Research (ONR) Departmental Research initiative, one key component of the ONR THINICE program was to advance knowledge on Arctic cyclone coupling with the ocean surface, its ice cover, and understand how such coupling can lead to rapid ice loss. Leading up to the primary field campaign in August 2022, a pilot field campaign from 19 August to 13 September 2021 focused in Utqiagvik, Alaska and Longyearbyen, Svalbard was carried out. This period was marked by frequent Arctic cyclones, periods of strong low-level winds, and stints of extreme short-term Beaufort Sea ice area losses. Cyclones traversed these waters much more frequently than usual during this 26-day period, but did not exhibit extreme central pressures indicative of high intensity and impactful storms. Using ERA5 atmospheric reanalysis fields, CERES and AIRS surface flux retrievals, and in situ sensor data collected from ice-tethered profilers and ice mass balance buoys, we highlight air-sea-ice interactions in the Beaufort Sea during and between these cyclones leading up to extreme ice loss events. Our case study analysis is aimed at understanding how Arctic cyclones, including air-sea interactions prompted by their passage over the marginal ice zone, influence the end-of-summer ice conditions and ice loss events. This study provides an example where frequent, yet weak storms drive rather large marginal sea-scale losses albeit with heterogenous sea ice effects at local scales.

Ballinger, Thomas [International Arctic Research C↗
DOE OSTI2024

Diabatic Valence-Hole Concept

A global diabatization scheme, based on the “valence-hole” concept, has been previously applied to model webs of avoided crossings that exist in four electronic-state symmetry manifolds of C 2 ( 1 Π g , 3 Π g , 1 Σ u + , and 3 Σ u + ). Here, this model is extended to the electronically excited states of four more molecules: CN ( 2 Σ + ), N 2 ( 3 Π u ), SiC ( 3 Π), and Si 2 ( 3 Π g ). Many strangenesses in the spectroscopic observations (e.g., energy level structure, predissociation linewidths, and radiative lifetimes) for all four electronic state systems discussed here are accounted for by this unified model. The key concept of the model is valence-hole electron configurations: 3σ 2 4σ 1 1π 4 5σ 2 in CN ( 2 Σ + ), 2σ g 2 2σ u 1 1π u 4 3σ g 2 1π g 1 in N 2 ( 3 Π u ), 5σ 2 6σ 1 7σ 2 2π 3 in SiC ( 3 Π), and 4σ g 2 4σ u 1 5σ g 2 2π u 3 in Si 2 ( 3 Π g ), all of which have a triply occupied “valence-core” (i.e., 2σ g 2 2σ u 1 or the equivalent). These valence-hole configurations have a nominal bond order of three or higher and correlate with high-energy separated-atom limits with an np ← ns (n = 2, 3) promotion in one of the atomic constituents. On its way to dissociation, the strongly bound diabatic valence-hole state crosses multiple weakly bound or repulsive states, which belong to electron configurations with a completely filled valence-core. These curve crossings between diabatic potentials result in a network of many avoided crossings among multiple electronic states, analogous to the well-studied electronic structure landscape of ionic-covalent crossings in strongly ionic molecules. In conclusion, considering the unique role of valence-hole states in shaping the global electronic structure, the valence-hole concept should be added to our intuitive framework of chemical bonding.

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

Are North Atlantic Tropical Cyclones Modulated by the Madden–Julian Oscillation in HighResMIP AGCMs?

This study assesses the representation of the observed relationship between Atlantic tropical cyclones (TCs) and the Madden–Julian oscillation (MJO) across nine models participating in CMIP6 High Resolution Model Intercomparison Project (HighResMIP). Most models struggle to faithfully reproduce the observed impacts of the MJO on Atlantic TCs, with the primary issue being the underestimated TC activity over the Atlantic main development region (MDR). The negative biases in genesis frequency within the MDR can be further attributed to weaker-than-observed African easterly wave (AEW) activity south of ∼12°N. Errors in the diabatic heating profile within the Atlantic intertropical convergence zone lead to insufficient potential vorticity production in the lower troposphere and constrain the amplification of AEWs. In addition to the biased TC climatology, the eastward-propagating power of the MJO is consistently underestimated across all models. Nevertheless, in models with a higher eastward–westward power ratio, the simulated MJO demonstrates a stronger capacity to modulate subseasonal TC activity. Models with relatively realistic eastward propagation of the MJO also exhibit greater variance in tropical intraseasonal convection. Stronger contrasts in convective heating over the North Atlantic between phases 2–3 and phases 6–7 drive larger fluctuations in MDR shear and AEW activity over the Gulf of Mexico and West Africa, resulting in a more pronounced TC response to the MJO. Overall, our findings suggest that improved MDR TC climatology and MJO propagation are essential for models to accurately capture the observed modulations of Atlantic TCs by the MJO.

54 ENVIRONMENTAL SCIENCES↗