Search NASA⌕ Search

SEARCH · Search NASA

Results for “non adiabatic transitions”

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.

49 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of dilution and pressure on combustion characteristics within externally heated micro-tubes

Here this study investigates the impact of pressure and dilution on combustion behavior of of C 2H6/air/N2 mixtures within an externally heated micro-tube using both experimental and analytical techniques. Experiments are conducted for pressures up to 10 bar at multiple dilution levels and document three flame regimes, i.e. weak flames, FREI and normal flames. Characteristic temperatures for flame stabilization, as well as and ignition and extinction are obtained by applying a novel thin filament pyrometry method (TFP), which uses a thin silicon carbide (SiC) filament that enables temperature measurements for both reacting and non-reacting cases. Results show a decrease of characteristic temperatures with increasing pressure, which can be which can be partially compensated by a decrease of O2 content. Further, regime transitions occur at lower velocity as the pressure increases, whereas weak flame positions are insensitive to dilution at a given pressure. Predictions from the analytical model compare favorably with experimental results. The most prominent finding is that as dilution increases, the transition from weak flame to FREI becomes less abrupt and is completely lost for high dilution levels (i.e. low O2 content). Both experimental and analytical results indicate that the unstable FREI combustion regime does not exist at high dilution levels, i.e. weak flame and normal flame regimes are no longer clearly separated. As FREI are interpreted as a cyclic transition from an ignition zone to a propagating flame, results for highly diluted conditions are interpreted as a combustion regime where classical flame propagation ceases to exist.

42 ENGINEERING↗
DOE OSTI2026

Temporal entanglement transitions in the periodically driven Ising chain

Periodically driven quantum systems can host non-equilibrium phenomena without static analogs, including in their entanglement dynamics. Here, we discover (TET) in a Floquet spin chain, which correspond to a quantum phase transition in the spectrum of the entanglement Hamiltonian and are signaled by dynamical spontaneous symmetry breaking. We identify the symmetry principles underlying these transitions: they appear when the driven Hamiltonian preserves global symmetry (here, ℤ2), the initial state respects this symmetry, and the reduced density matrix carries weight in both subsystem-parity sectors, with TET occurring precisely when the sector weights become equal (given the previous two conditions are also satisfied). Intriguingly, we find these transitions across a broad range of driving frequencies (from adiabatic to high-frequency regime) and independently of drive details, where they manifest as periodic, sharp entanglement spectrum reorganizations marked by the Schmidt-gap closure, a vanishing entanglement echo, and symmetry-quantum-number flips, while remaining invisible to conventional local observables. At high frequencies, the entanglement Hamiltonian acquires an intrinsic timescale decoupled from the drive period, rendering the transitions genuine steady-state features. Finite-size scaling reveals universal critical behavior with correlation-length exponent 𝜈 = 1, matching equilibrium Ising universality despite its emergence from purely dynamical mechanisms decoupled from static criticality. Furthermore, our work establishes TET as novel features in Floquet quantum matter.

36 MATERIALS SCIENCE↗
DOE OSTI2026

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

thornado-hydro: A Discontinuous Galerkin Method for Supernova Hydrodynamics with Nuclear Equations of State

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge–Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters — a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states — from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 M ⊙ progenitor. Herein, we find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS — more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

Nuclear astrophysics (1129)↗
DOE OSTI2021

Discovery of multiple p-mode pulsation frequencies in the roAp star, HD 86181

ABSTRACT We report the frequency analysis of a known roAp star, HD 86181 (TIC 469246567), with new inferences from Transiting Exoplanet Survey Satellite (TESS) data. We derive the rotation frequency to be νrot = 0.48753 ± 0.00001 d−1. The pulsation frequency spectrum is rich, consisting of two doublets and one quintuplet, which we interpret to be oblique pulsation multiplets from consecutive, high-overtone dipole, quadrupole, and dipole modes. The central frequency of the quintuplet is 232.7701 d−1 (2.694 mHz). The phases of the sidelobes, the pulsation phase modulation, and a spherical harmonic decomposition all show that the quadrupole mode is distorted. Following the oblique pulsator model, we calculate the rotation inclination, i, and magnetic obliquity, β, of this star, which provide detailed information about the pulsation geometry. The i and β derived from the best fit of the pulsation amplitude and phase modulation to a theoretical model, including the magnetic field effect, slightly differ from those calculated for a pure quadrupole, indicating the contributions from ℓ = 4, 6, 8,... are small. Non-adiabatic models with different envelope convection conditions and physics configurations were considered for this star. It is shown that models with envelope convection almost fully suppressed can explain the excitation at the observed pulsation frequencies.

Shi, Fangfei (ORCID:0000000323998216)↗
DOE OSTI2021

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2 p and 3 p 1 h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

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

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electron transfer in nonpolar media

The Marcus theory of electron transfer views fluctuating orientations of permanent dipoles as the nuclear mode bringing the donor and acceptor into the tunneling resonance. Electronic polarization of the solvent is excluded as the fast mode adiabatically following the electronic density. This view, valid for solids, does not apply to molecular liquids where molecular translations (density fluctuations) modulate the induction interaction of the donor-acceptor complex with the solvent. This mechanism of promoting radiation-less electronic transitions is considered here in the framework of the perturbation liquid-state theory. The reorganization energy of electron transfer in nonpolar solvents is nonzero and reaches the values of 0.1–0.3 eV for typical molecular sizes and solvents used in applications. The reorganization energy scales quadratically with the molecular polarizability of the solvent and decays as the inverse fifth power with the size of the donor and acceptor. The combination of the entropic character of the density fluctuations, driven by re-packing of molecular cores, with the short range of induction solute–solvent interactions leads to the violation of the fluctuation–dissipation theorem for the variance of the donor-acceptor energy gap. An explicit, approximately hyperbolic, dependence of the reorganization energy on temperature is predicted. It leads to a non-Arrhenius kinetic law for the rate constant of electron transfer in nonpolar liquid solvents.

14 SOLAR ENERGY↗
DOE OSTI2023

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

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

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

On the origin of the DIII-D L-H power threshold isotope effect

Abstract The increased low to high confinement mode (L to H-mode) power threshold P L H in DIII-D low collisionality hydrogen plasmas (compared to deuterium) is shown to result from lower impurity (carbon) content, consistent with reduced (mass-dependent) physical and chemical sputtering of graphite. Trapped gyro-Landau fluid (TGLF) quasilinear calculations and local non-linear gyrokinetic CGYRO simulations confirm stabilization of ion temperature gradient (ITG) driven turbulence by increased carbon ion dilution as the most important isotope effect. In the plasma edge, electron non-adiabaticity is also predicted to contribute to the isotope dependence of thermal transport and P L H , however its effect is subdominant compared to changes from impurity isotopic behavior. This L-H power threshold reduction with increasing carbon content at low collisionality is in stark contrast to high collisionality results, where additional impurity content appears to increase the power necessary for H-mode access.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Properties of the Prompt Optical Counterpart Arising from the Cooling of Electrons in Gamma-Ray Bursts

This work extends a contemporaneous effort to study the properties of the lower-energy counterpart synchrotron emission produced by the cooling of relativistic gamma-ray burst (GRB) electrons through radiation (synchrotron and self-Compton) emission and adiabatic losses. We derive the major characteristics (pulse duration, lag time after burst, and brightness relative to the burst) of the prompt optical counterpart (POC) occurring during or after the GRB. Depending on the magnetic field lifetime, duration of electron injection, and electron transit time Δt o from hard X-ray (GRB) to optical-emitting energies, a POC may appear during the GRB pulse (of duration δt γ ) or after (delayed OC). The signature of counterparts arising from the cooling of GRB electrons is that POC pulses (Δto < δtγ) last as long as the corresponding GRB pulse (δt o ≃ δt γ ), while delayed OC pulses (Δt o > δt γ ) last as long as the transit time (δt o ≃ Δt o ). If OC variability can be measured, then another signature for this OC mechanism is that the GRB variability is passed on to POCs but not to delayed OCs. Within the GRB electron cooling model for counterparts, POCs should be on average dimmer than delayed ones (consistent with the data), and harder GRB low-energy slopes β LE should be associated more often with the dimmer POCs. The latter sets an observational bias against detecting POCs from (the cooling of electrons in) GRBs with a hard slope β LE , making it more likely that the detected POCs of such bursts arise from another mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗