Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-adiabaticity”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-throughput and data-driven search for stable optoelectronic AMSe 3 materials

The rapid advancement in emerging optoelectronic technologies demands highly efficient, affordable, and ecofriendly materials. In this context, ternary chalcogenides, especially ternary selenides, show early promise as a material class due to their stability and remarkable electronic, optical, and transport properties. In this work, we integrate first-principles-based high-throughput computations with machine learning (ML) techniques to predict the thermodynamic stability and optoelectronic properties of 920 valency-satisfied selenide compounds. Through investigating polymorphism, our study reveals the edge-sharing orthorhombic Pnma phase (NH 4 CdCl 3 -type) as the most stable structure for most ternary selenides. High-fidelity supervised ML models are trained and tested to accelerate stability and band gap predictions. These data-driven models pin down the most influential features that dominantly control key material characteristics. The multistep high-throughput computations identify the ternary selenides with optimal direct band gaps, light carrier masses, and strong optical absorption edges. The extensive materials screening considering phase stability, toxicity, and defect tolerance, finally identifies the seven most suitable candidates for photovoltaic applications. Two of these final compounds, SrZrSe 3 and SrHfSe 3 , have already been synthesized in a single-phase form, with the latter showing an optically suitable band gap, aligning well with our findings. The non-adiabatic molecular dynamics reveal sufficiently long photoexcited charge carrier lifetimes (on the order of nanoseconds) in some of these selected selenide materials, indicating their exciting characteristics. Overall, our study suggests a robust in silico framework that can be extended to screen large datasets of various material classes for identifying promising photoactive candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
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↗
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 OSTI2024

Measuring the molecular Migdal effect with neutron scattering on diatomic gases

The Migdal effect is a key inelastic signal channel which could be used to detect low-mass dark matter, but it has never been observed experimentally using Standard Model probes. Here we propose a conceptual design for an experiment which could detect the Migdal effect in diatomic molecules through low-energy neutron scattering, and we provide the requirements on the beam spectrum and reducible backgrounds such that a detection may be achieved. The enhancement of the Migdal rate through non-adiabatic couplings, which are absent in isolated atoms, combined with the distinctive photon energies of electronic transitions in CO, suggest that a positive detection of the molecular Migdal effect may be possible with modest beam times at existing neutron facilities. Published by the American Physical Society 2024

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Spin-Controllable Dynamics in Defect-Engineered Carbon Nanotubes as Single Photon Emitters: Data-Driven Modeling and Computations

Quantum technologies, such as quantum computing and sensing, require efficient single-photon emission (SPE) sources that operate at room temperature in telecom wavelengths. While several materials can serve as SPE sources, no single platform meets all the criteria for efficiency, ambient operation, and scalability. Single-walled carbon nanotubes (SWCNTs) with covalently attached molecules offer a promising solution. Their SPE can be easily tuned via modifications of the SWCNT's diameter, chirality, and bonded molecules, enabling emission across near-IR to telecom wavelengths at ambient conditions. However, to fully realize the potential of SWCNTs and unlock their quantum capabilities, a deeper understanding of how structural defects from molecular adducts affect their emission and competing photoexcited processes is essential. To address this gap in our knowledge, this project combined quantum chemistry calculations with data-driven methods of cheminformatics (QSAR) and machine learning (ML). The developed computational approaches have provided several design strategies for covalent functionalization of SWCNTs to improve their optical response. The collaboration with Los Alamos National Lab (LANL) enabled direct comparison of computational and experimental data, facilitating method validation. This partnership was enhanced through access to LANL's Center for Integrated Nanotechnologies (CINT) utilizing User Facility Program and summer internships, which provided three NDSU graduate students with hands-on experience at LANL. The outcomes of this project included (1) Advancing the current stage of computational methods in accurate modeling of non-adiabatic spin-dependent photoexcited dynamics and its applicability to nanosystems consisting of thousands of atoms, realized as open-access codes linked to existing DFT-based software; (2) Establishing the relationship between the structure of adducts and SWCNTs and intrinsic excitonic and spin properties of defect states for guiding novel synthetic strategies and experimental probes of chemically functionalized SWCNTs as near-IR emitting materials; (3) Generating virtual libraries of hypothetical functionalized SWCNTs for virtual screening of their chemical structures and optical properties, leveraging new functionalities of SWCNTs; (4) Offering a unique experience for NDSU graduate students that prepared them for future scientific careers related to materials modeling and big data processing. These results were summarized in 12 published journal papers and 3 recently submitted papers. One of a key finding is that the position of defect sites on the SWCNT surface primarily drives the emission redshift (up to 100 meV), while the polarity of the defect-inducing molecules has a much smaller effect (~10 meV). However, the electron-donating or withdrawing properties of a molecule influence selecting reactivity of defect sites. These insights important for optimizing synthetic protocols for desired emissions in SWCNTs. We also revealed that the interaction between two defects at various positions on the SWCNT enhances the redshift and optical activity of states, favoring strong near-IR emission. This suggests that manipulations in defect concentrations is a promising strategy for controlling efficient emission. Mostly important, the defect position was found controllable by the spin states of photoexcited intermediates: Excited aromatic molecules form ortho defects with SWCNTs at their singlet states in the presence of oxygen, while oxygen-free conditions favor para defects via the triplet-state mechanism. Additionally, a heat-activated [2+2] cycloaddition reaction facilitates divalent defect formation with fewer bonding positions that narrows emission bands. These groundbreaking findings have been experimentally validated and significantly advance our understanding of defect chemistry in SWCNTs. Using a novel encoding technique and 3D-MoRSE descriptors, we developed highly accurate ML/QSAR models to predict both the 3D structure and optical properties of SWCNTs with chemical defects. This model enabled the creation of a virtual library of 125,556 structures, providing new insights into the relationship between SWCNT-defect structure and emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗
NASA NTRS1962

Thermodynamic and Transport Property Correlation Formulas for Equilibrium Air from 1,000 Degrees Kelvin to 15,000 Degrees Kelvin

The thermodynamic properties, density and temperature, as well as transport property parameters involving viscosity, Prandtl number (including diffusion effects), and gaseous radiation absorption coefficients have been correlated as a function of enthalpy at four pressure levels (10 (sup -1), 10 (sup 0), 10, and 10 (sup 2) atmospheres). The correlation formulas are written in a generalized form for which coefficients for a particular property and pressure level are tabulated. The correlation formulas are useful in digital computer programs for non-adiabatic viscous flow problems.

AIR↗
NASA NTRS1972

Very high temperature laminar flow of a gas through the entrance region of a cooled tube - Numerical calculations and experimental results.

The laminar flow equations in differential form are solved numerically on a digital computer for flow of a very high temperature gas through the entrance region of an externally cooled tube. The solution method is described and calculations are carried out in conjunction with experimental measurements. The agreement with experiment is good, with the result indicating relatively large energy and momentum losses in the highly cooled flows considered where the pressure is nearly uniform along the flow and the core flow becomes non-adiabatic a few diameters downstream of the inlet. The effects of a large range of Reynolds number and Mach number (viscous dissipation) are also investigated.

Back, L. H.↗
NASA NTRS1976

Retrodictive determinism

With respect to irreversible, non-homeomorphic maps, contravariant and covariant tensor fields have distinctly natural covariance and transformational behavior. For thermodynamic processes which are non-adiabatic, the fact that the process cannot be represented by a homeomorphic map emphasizes the logical arrow of time, an idea which encompasses a principle of retrodictive determinism for covariant tensor fields.

Kiehn, R. M.↗
NASA NTRS1978

Direct heating surface combustor

The combustor utilizes a non-adiabatic flame to provide low-emission combustion for gas turbines. A fuel-air mixture is directed through a porous wall, the other side of which serves as a combustion surface. A radiant heat sink disposed adjacent to and spaced from the combustion surface controls the combustor flame temperature in order to prevent the formation of oxides of nitrogen. A secondary air flow cools the heat sink. Additionally, up to 100% of secondary air flow is mixed with the combustion products at the direct heating surface combustor to dilute such products thereby reducing exit temperature. However, if less than 100% secondary air is mixed to the combustor, the remainder may be added to the combustion products further downstream.

Beremand, D. G.↗
NASA NTRS1980

Theoretical studies of solar oscillations

Possible sources for the excitation of the solar 5 minute oscillations were investigated and a linear non-adiabatic stability code was applied to a preliminary study of the solar g-modes with periods near 160 minutes. Although no definitive conclusions concerning the excitation of these modes were reached, the excitation of the 5 minute oscillations by turbulent stresses in the convection zone remains a viable possibility. Theoretical calculations do not offer much support for the identification of the 160 minute global solar oscillation (reported by several independent observers) as a solar g-mode. A significant advance was made in attempting to reconcile mixing-length theory with the results of the calculations of linearly unstable normal modes. Calculations show that in a convective envelope prepared according to mixing length theory, the only linearly unstable modes are those which correspond to the turbulent eddies which are the basic element of the heuristic mixing length theory.

Goldreich, P.↗
NASA NTRS1980

An artificial energy method for calculating flows with shocks

The artificial-viscosity method, first proposed by von Neumann and Richtmyer, introduces an artificial viscous pressure term in regions of compression such that an increase in entropy occurs in shock transition zones. The paper describes how dissipative flows can be induced by reducing the total energy available for adiabatic processes in shock zones. A class of inviscid fluid flows, called semiflows, is described in which the flows exhibit thermodynamic differences. Induced dissipative flows modify the pressure in regions of compression in a manner analogous to the artificial-viscosity method and for a gas, the effect is equivalent to suitably modifying the gas constant in the equation of state. By employing MacCormack's method and the usual non-adiabatic equations, numerical solutions of a Riemann problem are compared with the modified artificial energy method, showing that the dissipation effect predicted by the analytical formulation is reflected in the numerical method as well.

Rose, M. E.↗
NASA NTRS1985

Observational and theoretical studies of the evolving structure of baroclinic waves: Attractor dynamics of global weather systems

Significant progress was made in developing a six-component dry model of mid-latitude baroclinic wave evolution, and a ten-component moist model. These models include representations of airflow over mountians and non-adiabatic processes. As a first step, the properties of a reduced, three component, baroclinic system are studied. The classical baroclinic stability criteria emerge as the properties of one of the three equilibria admitted, with the remaining two equilibria forming the attractor region for the observed weather activity. Passive tracer evolution in a baroclinic wave was studied. A model in which only constituents trapped within low atmospheric levels are considered, has been extended to include the possibility for fluxes into the upper levels of the atmosphere. The reported results for the lower levels achieved previously are shown to be qualitatively similar to those obtained by these new calculations in which the vertical flux constraint is removed. An extensive study of an 8-year record of global outgoing longwave radiation for the Northern Hemisphere reveals that blocking events exhibit only a weak signature of blocking highs, as measured by relatively low values of the ratio of the standard deviation to the mean value of the observed long wave outgoing radiation. Though present in many cases, the signature is not a strongly distinctive feature of the blocking episodes.

Saltzman, B.↗
NASA NTRS1988

Experimental verification of the thermodynamic properties for a jet-A fuel

Thermodynamic properties for a Jet-A fuel were determined by Shell Development Company in 1970 under a contract for NASA Lewis Research Center. The polynomial fit necessary to include Jet-A fuel (liquid and gaseous phases) in the library of thermodynamic properties of the NASA Lewis Chemical Equilibrium Program is calculated. To verify the thermodynamic data, the temperatures of mixtures of liquid Jet-A injected into a hot nitrogen stream were experimentally measured and compared to those calculated by the program. Iso-octane, a fuel for which the thermodynamic properties are well known, was used as a standard to calibrate the apparatus. The measured temperatures for the iso-octane/nitrogen mixtures reproduced the calculated temperatures except for a small loss due to the non-adiabatic behavior of the apparatus. The measurements for Jet-A were corrected for this heat loss and showed excellent agreement with the calculated temperatures. These experiments show that this process can be adequately described by the thermodynamic properties fitted for the Chemical Equilibrium Program.

Graciasalcedo, Carmen M.↗
NASA NTRS1990

Current carriers in the near-earth cross-tail current sheet during substorm growth phase

Throughout most of the growth phase of a substorm, the cross-tail current at x about -10 Re can be supplied by the curvature drift of a bi-directional field aligned distribution of 1 keV electrons. Just prior to its local disruption after substorm onset, the cross-tail current in the now thin (about 400 km) current sheet is carried by the cross-tail serpentine motion of non-adiabatic ions (Speiser, 1965). The instability of this latter current leads to the local disruption of the near-earth current sheet.

Mitchell, D. G.↗