SEARCH · Engineering Papers
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.
Adjoint sensitivity equations are presented, which can be solved simultaneously (i.e., forward in time) with the dynamics of a nonlinear neural network. These equations provide the foundations for a new methodology which enables the implementation of temporal learning algorithms in a highly efficient manner.
Abstract Structural imaging of transient excited-state species is a key goal of molecular physics, promising to unveil rich information about the dynamics underpinning photochemical transformations. However, separating the electronic and nuclear contributions to the spectroscopic observables is challenging, and typically requires the application of high-level theory. Here, we employ site-selective ionisation via ultrashort soft X-ray pulses and time-resolved Coulomb explosion imaging to interrogate structural dynamics of the ultraviolet photochemistry of carbon disulfide. This prototypical system exhibits the complex motifs of polyatomic photochemistry, including strong non-adiabatic couplings, vibrational mode couplings, and intersystem crossing. Immediately following photoexcitation, we observe Coulomb explosion signatures of highly bent and stretched excited-state geometries involved in the photodissociation. Aided by a model to interpret such changes, we build a comprehensive picture of the photoinduced nuclear dynamics that follows initial bending and stretching motions, as the reaction proceeds towards photodissociation.
Electromagnetic ion cyclotron (EMIC) waves are an important magnetospheric emission, which is excited near the magnetic equator with frequencies below the proton gyro-frequency. The source of bee energy for wave growth is provided by temperature anisotropy of ring current (RC) ions, which develops naturally during inward convection from the plasma sheet These waves strongly affect the dynamic s of resonant RC ions, thermal electrons and ions, and the outer radiation belt relativistic electrons, leading to non-adiabatic particle heating and/or pitch-angle scattering and loss to the atmosphere. The rate of ion and electron scattering/heating is strongly controlled by the Wave power spectral and spatial distributions, but unfortunately, the currently available observational information regarding EMIC wave power spectral density is poor. So combinations of reliable data and theoretical models should be utilized in order to obtain the power spectral density of EMIC waves over the entire magnetosphere throughout the different storm phases. In this study, we present the simulation results, which are based on two coupled RC models that our group has developed. The first model deals with the large-scale magnetosphere-ionosphere electrodynamic coupling, and provides a self-consistent description of RC ions/electrons and the magnetospheric electric field. The second model is based on a coupled system of two kinetic equations, one equation describes the RC ion dynamics and another equation describes the power spectral density evolution of EMIC waves, and self-consistently treats a micro-scale electrodynamic coupling of RC and EMIC waves. So far, these two models have been applied independently. However, the large-scale magnetosphere-ionosphere electrodynamics controls the convective patterns of both the RC ions and plasmasphere altering conditions for EMIC wave-particle interaction. In turn, the wave induced RC precipitation Changes the local field-aligned current distributions and the ionospheric conductances, which are crucial for a large-scale electrodynamics. The initial results from this new self-consistent model of the magnetospheric electric field, RC and EMIC waves will be shown in this presentation.
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..
The origin of the charge density wave phases in the kagome metal compound AV 3 Sb 5 is still under great scrutiny. Here, we combine diffuse and inelastic x-ray scattering to identify a 3-dimensional precursor of the charge order at the L point that condenses into a CDW through a first order phase transition. The quasi-elastic critical scattering indicates that the dominant contribution to the diffuse precursor is the elastic central peak without phonon softening. However, the inelastic spectra show a small broadening of the Einstein-type phonon mode on approaching T CDW . Our results point to the situation where the Fermi surface instability at the L point is of order-disorder type with critical growth of quasi-static domains. The experimental data indicate that the CDW consists on an alternating Star of David and trihexagonal distortions and its dynamics goes beyond the classical weak-coupling scenario and is discussed within strong-electron phonon coupling and non-adiabatic models.
The DMC meeting, since 1965, has a great history of providing a unique platform and focal point for the gathering of experimentalists and theoreticians in the field of chemical dynamics. At this biannual conference which spans five days, the latest developments in bimolecular reactions, unimolecular reactions, reactions in clusters and solutions, reactions at interfaces, reactions in combustion, atmospheres, and interstellar media, as well as cold chemistry, new spectroscopic techniques and non-adiabatic reactions will be discussed by experts in these fields. The 2021 meeting was postponed due to Covid but is back in 2023 including the traditional awarding of two Herschbach Medals, honoring important contributions to the fields of theoretical and experimental dynamics respectively. The 2023 meeting will be organized by Richard Dawes (Chair) and Amy Mullin (Vice Chair). The DMC format encourages participation of young scientists through invited and contributed talks as well as posters. This proposal seeks support for students, post-docs and invited scientists to attend the meeting.
Mechanistic insights into photodissociation dynamics of transition metal carbonyls, like Fe(CO) 5 , are fundamental for understanding active catalytic intermediates. Although extensively studied, the structural dynamics of these systems remain elusive. Using ultrafast X-ray scattering, we uncover the photochemistry of Fe(CO) 5 in real space and time, observing synchronous oscillations in atomic pair distances, followed by a prompt rotating CO release preferentially in the axial direction. This behavior aligns with simulations, reflecting the interplay between the axial Fe-C distances’ potential energy landscape and non-adiabatic transitions between metal-to-ligand charge-transfer states. Additionally, we characterize a secondary delayed CO release associated with a reduction of Fe-C steady state distances and structural dynamics of the formed Fe(CO) 4 . Our results quantify energy redistribution across vibration, rotation, and translation degrees of freedom, offering a microscopic view of complex structural dynamics, enhancing our grasp on Fe(CO) 5 photodissociation, and advancing our understanding of transition metal catalytic systems.
Time-dependent density functional theory continues to draw a large number of users in a wide range of fields exploring myriad applications involving electronic spectra and dynamics. Although in principle exact, the predictivity of the calculations is limited by the available approximations for the exchange-correlation functional. In particular, it is known that the exact exchange-correlation functional has memory-dependence, but in practise adiabatic approximations are used which ignore this. Here we review the development of non-adiabatic functional approximations, their impact on calculations, and challenges in developing practical and accurate memory-dependent functionals for general purposes.
Here, we explore the dynamics of energy transfer between two nanocrystal quantum dots placed within an optical microcavity. By adjusting the coupling strength between the cavity photon mode and the quantum dots, we have the capacity to fine-tune the effective coupling between the donor and acceptor. Introducing a non-adiabatic parameter, γ, governed by the coupling to the cavity mode, we demonstrate the system’s capability to shift from the overdamped Förster regime (γ $\ll$ 1) to an underdamped coherent regime (γ $\gg$ 1). In the latter regime, characterized by swift energy transfer rates, the dynamics are influenced by decoherence time. To illustrate this, we study the exciton energy transfer dynamics between two closely positioned CdSe/CdS core/shell quantum dots with sizes and separations relevant to experimental conditions. Employing an atomistic approach, we calculate the excitonic level arrangement, exciton–phonon interactions, and transition dipole moments of the quantum dots within the microcavity. These parameters are then utilized to define a model Hamiltonian. Subsequently, we apply a generalized non-Markovian quantum Redfield equation to delineate the dynamics within the polaritonic framework.
We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.
Here, we present an investigation of the relaxation dynamics of deuterated water molecules after direct photo-double ionization at 61 eV. We focus on the very rare D + + O + + D reaction channel in which the sequential fragmentation mechanisms were found to dominate the dynamics. Aided by theory, the state-selective formation and breakup of the transient OD + (a 1 Δ, b 1 Σ + ) is traced, and the most likely dissociation path—OD + : a 1 Δ or b 1 Σ + → A 3 Π → X 3 Σ – → B 3 Σ – —involving a combination of spin–orbit and non-adiabatic charge transfer transitions is determined. The multi-step transition probability of this complex transition sequence in the intermediate fragment ion is directly evaluated as a function of the energy of the transient OD + above its lowest dissociation limit from the measured ratio of the D + + O + + D and competing D + + D + + O sequential fragmentation channels, which are measured simultaneously. Our coupled-channel time-dependent dynamics calculations reproduce the general trends of these multi-state relative transition rates toward the three-body fragmentation channels.
Our simulation demonstrates the robustness of the Coulomb explosion imaging technique in studying methyl iodide photodissociation, and shows that it can be effectively used for imaging non-adiabatic transitions in coordinate space.
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.
Design-specific control over the transitions between excited electronic states with different spin multiplicities is of the utmost importance in molecular and materials chemistry. Previous studies have indicated that the combination of spin-orbit and vibronic effects, collectively termed the spin-vibronic effect, can accelerate quantum-mechanically forbidden transitions at non-adiabatic crossings. However, it has been difficult to identify precise experimental manifestations of the spin-vibronic mechanism. Here we present coherence spectroscopy experiments that reveal the interplay between the spin, electronic and vibrational degrees of freedom that drive efficient singlet-triplet conversion in four structurally related dinuclear Pt(II) metal-metal-to-ligand charge-transfer (MMLCT) complexes. Photoexcitation activates the formation of a Pt-Pt bond, launching a stretching vibrational wavepacket. The molecular-structure-dependent decoherence and recoherence dynamics of this wavepacket resolve the spin-vibronic mechanism. Further, we find that vectorial motion along the Pt-Pt stretching coordinates tunes the singlet and intermediate-state energy gap irreversibly towards the conical intersection and subsequently drives formation of the lowest stable triplet state in a ratcheting fashion. This work demonstrates the viability of using vibronic coherences as probes to clarify the interplay among spin, electronic and nuclear dynamics in spin-conversion processes, and this could inspire new modular designs to tailor the properties of excited states.
A recently published analysis of current sheets has updated the classic Harris 1D static solution by considering multiple classes of charged particle trajectories in a generalized and dynamic current sheet. It uses a 1D PIC simulation to describe dynamic pinching and bifurcation of the current sheet. These 1D results strongly suggest that properties of the inflowing plasma, including the plasma beta, have an important effect on the equilibrium thickness of the pinched current sheet. Since 1D studies cannot describe magnetic reconnection, the time appears right to carry such 1D studies over to 2D or 3D simulations to explore current sheet thickness effects on reconnection. The Magnetospheric Multiscale Mission (MMS), with its well-resolved multipoint measurements of collisionless plasma and fields, has found that collisionless reconnection is accompanied by non-adiabatic motions of electrons that only occur in magnetic structures with a narrow scale comparable to electron inertial lengths (d e ). The recent 1D studies suggest that a plasma pinch to such scales may only occur for inflowing magnetized plasmas with relatively low plasma beta. We conclude that a parametric exploration of simulated and observed reconnection inflow conditions, particularly plasma beta, should shed light on the enablement of reconnection in collisionless plasmas.
We investigate the stability of superconducting strings as bound states of strings and fermion zero modes at both the classical and quantum levels. The dynamics of these superconducting strings can result in a stable configuration, known as a vorton . We mainly focus on global strings, but the majority of the discussion can be applied to local strings. Using lattice simulations, we study the classical dynamics of superconducting strings and confirm that they relax to the vorton configuration through Nambu-Goldstone boson radiation, with no evidence of over-shooting that would destabilize the vorton. We explore the tunneling of fermion zero modes out of the strings. Both our classical analysis and quantum calculations yield consistent results: the maximum energy of the zero mode significantly exceeds the fermion mass, in contrast to previous literature. Additionally, we introduce a world-sheet formalism to evaluate the decay rate of zero modes into other particles, which constitute the dominant decay channel. We also identify additional processes that trigger zero-mode decay due to non-adiabatic changes of the string configuration. In these decay processes, the rates are suppressed by the curvature of string loops, with exponential suppression for large masses of the final states. We further study the scattering with light charged particles surrounding the string core produced by the zero-mode current and find that a wide zero-mode wavefunction can enhance vorton stability.
Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.
Turbulent and compressed sheath regions preceding interplanetary coronal mass ejections strongly impact electron dynamics in the outer radiation belt. Changes in electron flux can occur on timescales of tens of minutes, which are unlikely to be captured by a two-satellite mission. The recently released Global Positioning System (GPS) data set generally has shorter revisit times (at L ~ 4–8) owing to the large number of satellites in the constellation equipped with energetic particle detectors. Investigating electron fluxes at energies from 140 keV to 4 MeV and sheaths observed in 2012–2018, we show that the flux response to sheaths on a timescale of 6 hr, previously reported from Van Allen Probes (RBSP) data, is reproduced by GPS measurements. Furthermore, GPS data enables derivation of the response on a timescale of 30 min, which further confirms that the energy and L-shell dependent changes in electron flux are associated with the impact of the sheath. Sheath-driven loss is underestimated over longer timescales as the electrons recover during the ejecta. We additionally show the response of electron phase space density (PSD), which is a key quantity in identifying non-adiabatic loss from the system and electron energization through wave-particle interactions. The PSD response is calculated from both RBSP and GPS data for the 6 hr timescale, as well as from GPS data for the 30 min timescale. The response is divided based on the geoeffectiveness of the sheaths revealing that electrons are effectively accelerated only during geoeffective sheaths, while loss commonly occurs during all sheaths.