Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2

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↗

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.↗

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↗

Mechanistic Insights into Photocatalyzed H 2 Dissociation on Au Clusters

Localized surface plasmon resonances (LSPRs) have attracted much recent attention for their potential in promoting chemical reactions with light. However, the mechanism of LSPR-induced chemical reactions is still not clear, even for H 2 dissociation on metal nanoparticles. Here, we investigate the mechanism for photoinduced H 2 dissociation using a simple H 2 @Au 6 model. Our time-dependent density functional theory calculations indicate that the initial excitation is largely restricted to the metal cluster, involving intraband excitation that produces hot electrons (HEs). However, diabatization via overlapping orbitals reveals two types of nested electronic states, one involving excitations of the metallic electrons, namely, the HE states, and the other concerned with charge transfer (CT) to the adsorbate antibonding σ* orbital. Dissociation of H 2 thus takes place by transitions from the former to the latter. Quantum dynamics simulations on the diabatic CT states suggest rapid dissociation of H 2 , while no such dissociation occurs on diabatic HE states. Our research provides a clear physical picture of photoinduced H 2 dissociation on Au clusters, which has important implications in plasmonic facilitated photocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiglobal diabatic potential energy matrix for the N–H photodissociation of methylamine

We constructed an analytic diabatic potential energy matrix (DPEM) that describes the N–H photodissociation of methylamine; the electronic state space includes the ground and first excited singlet states. The input for the fit was calculated by extended multi-state complete active space second-order perturbation theory. The data were diabatized using the dipole–quadrupole diabatization method in which we incorporated a coordinate-dependent weighting scheme for the contribution of the quadrupole moments. Therefore, to make the resulting potential energy surfaces semiglobal, we extended the anchor points reactive potential method, a multiscale approach that assigns the internal coordinates to categories with different levels of computational treatment. Key aspects of the adiabatic potential energy surfaces obtained by diagonalizing the DPEM agree with the available experimental and theoretical data at energies relevant for photochemical studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

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↗

Full configuration-interaction study of the ionic-neutral curve crossing in LiF

Full configuration-interaction (FCI) calculations are used to assess the relative ability of methods for truncating the n-particle expansion in describing the ionic-neutral curve crossing between the two lowest Sigma(+) states of LiF. While the FCI calculations yield a smooth dipole moment function, MRCI calculations based on CASSCF orbitals optimized for the lowest state at all r values yield a discontinuous dipole moment function. However, when the orbitals are optimized using a state-averaged CASSCF procedure, with equal weights for the ionic and neutral solutions, both the CASSCF and MRCI dipole moment functions are smooth and in reasonable agreement with the FCI. No single-reference-based method is found to work satisfactorily. Potential curves for the lowest two Sigma(+) states are determined in both the adiabatic and diabatic representations using a large atomic natural orbit Gaussian basis set and a state-averaged CASSCF/MRCI treatment of electron correlation.

Bauschlicher, Charles W., Jr.↗

Test particle propagation in magnetostatic turbulence. 3: The approach to equilibrium

The asymptotic behavior, for large time, of the quasi-linear diabatic solutions and their local approximations is considered. A time averaging procedure is introduced which yields the averages of these solutions over time intervals which contain only large time values. A discussion of the quasi-linear diabatic solutions which is limited to those solutions that are bounded from below as functions of time is given. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity the time averaged quasi-linear diabatic solutions must approach isotropy (mu-independence). The first derivative with respect to mu of these solutions is also considered. This discussion is limited to first derivatives which are bounded functions of time. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity, the time averaged first derivative must approach zero everywhere in mu except at mu = 0 where it must approach a large value which is calculated. The impact of this large derivative on the quasi-linear expansion scheme is discussed. An H-theorem for the first local approximation to the quasi-linear diabatic solutions is constructed. Without time averaging, the H-theorem is used to determine sufficient conditions for the first local approximate solutions to asymptote, with increasing time, to exactly the same final state which the time averaged quasi-linear diabatic solutions must approach as discussed above.

Klimas, A. J.↗

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↗

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↗

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↗

Mechanisms behind the Springtime North Pacific ENSO Teleconnection Bias in Climate Models

Previous studies have shown that models overestimate the strength of ENSO teleconnections to the North Pacific during springtime, but the underlying reasons for this bias remain unknown. Here, in this work, the relative contributions from basic-state and thermodynamic/dynamic forcing factors are disentangled through idealized experiments with the Community Earth System Model and a range of stationary wave modeling experiments. It is revealed that in CESM1 the diabatic heating biases over the tropical Indian Ocean and tropical central-western Pacific jointly favor a cyclonic (anticyclonic) circulation bias to occur in the North Pacific during the springtime of El Niño (La Niña) events. On one hand, the difference in the modeled and observed climatological basic state does not lead to the bias formation directly, as the diabatic heating biases are the primary cause. On the other hand, the springtime basic state is conducive to a more vigorous stationary wave response to the biased diabatic heating than the wintertime state, and this explains why the teleconnection bias occurs during springtime but not in winter. An iterative bias-correction approach is then implemented in the atmospheric model component of CESM1 to verify the linkage between the tropical diabatic heating bias and the teleconnection bias. Moreover, this explanation is shown to be relevant in other models of phase 5 of the Coupled Model Intercomparison Project (CMIP5) as a strong relationship is found between biases in ENSO-related tropical central-western Pacific/Indian Ocean precipitation and North Pacific circulation across models in spring.

54 ENVIRONMENTAL SCIENCES↗

Exploring spin symmetry-breaking effects for static field ionization of atoms: Is there an analog to the Coulson–Fischer point in bond dissociation?

Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. Here, this paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree-Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single M s (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree-Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗