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.
We investigate the time-resolved dynamics of interatomic Coulombic decay in NeKr dimers after neon inner-valence photoionization using a synchrotron light source.
This Perspective illustrates the growing use of ultrafast laser induced Coulomb explosion imaging methods for tracking time-evolving molecular structures and advancing understanding of gas phase molecular fragmentation processes.
We have developed a local, linear theoretical model for lower hybrid drift waves that can be used for plasmas in the weakly collisional regime. Two cases with typical plasma and field parameters for the current sheet of the magnetic reconnection experiment have been studied. For a case with a low electron beta (β e =0.25, high guide field case), the quasi-electrostatic lower hybrid drift wave is unstable, while the electromagnetic lower hybrid drift wave has a positive growth rate for a high-β e case (β e =8.9, low guide field case). For both cases, including the effects of Coulomb collisions reduces the growth rate but collisional impacts on the dispersion and growth rate are limited (≲20%).
Exact moments of the Boltzmann collision operator are calculated in the irreducible Hermitian moment expansion written in terms of the random-velocity variable of each species. The formulas are presented in closed, algebraic form and can be straightforwardly implemented in computer algebra systems. They are valid for two arbitrary masses, temperatures, and flow velocities, and hence include all other existing results derived for distribution functions expanded with respect to reference states of one temperature and flow velocity. In comparison, the Landau collisional moments are good approximations for large Coulomb logarithm and small relative flow velocity, but they fail to predict the correct behavior of most collisional moments for large relative flow even for weakly coupled plasmas.
Here, the interplay of stress, disorder, and Coulomb screening dictating the mobility of doped cadmium oxide (CdO) is examined using Raman spectroscopy to identify the mechanisms driving dopant incorporation and scattering within this emerging infrared optical material. Specifically, multi-wavelength Raman and UV-vis spectroscopies are combined with electrical Hall measurements on a series of yttrium (X = Y) and indium (X = In) doped X:CdO thin-films. Hall measurements confirm n-type doping and establish carrier concentrations and mobilities. Spectral fitting along the low-frequency Raman combination bands, especially the TA+TO(X) mode, reveals that the evolution of strain and disorder within the lattice as a function of dopant concentration is strongly correlated with mobility. Coupling between the electronic and lattice environments was examined through analysis of first- and second-order longitudinal–optical phonon–plasmon coupled modes that monotonically decrease in energy and asymmetrically broaden with increasing dopant concentration. By fitting these trends to an impurity-induced Fröhlich model for the Raman scattering intensity, exciton–phonon and exciton–impurity coupling factors are quantified. These coupling factors indicate a continual decrease in the amount of ionized impurity scattering with increasing dopant concentration and are not as well correlated with mobility. This shows that lattice strain and disorder are the primary determining factors for mobility in donor-doped CdO. In aggregate, the study confirms previously postulated defect equilibrium arguments for dopant incorporation in CdO while at the same time identifying paths for its further refinement.
Positive ion desorption following electron impact dissociative ionization of ethane adsorbed on Ar, Kr, and Xe multilayers has been studied as a function of incident electron energy from threshold to 100 eV. Based on the dependence of ion yields on the identity of the rare gas, it is likely that the majority of ethane molecules undergo indirect ionization following hole transfer from the ionized underlying rare gas. Further, this has also been corroborated by density of states calculations showing the energetic alignment of the outer valence states of ethane and the condensed rare gas ionization energies. Due to the near-resonant nature of charge transfer for single-hole states, the ethane molecular ion is excited to different final ionic states on different rare gases, which leads to differences in ion desorption yields and branching ratios. The quantitative yields increase with increasing ionization energy gap between the rare gas and ethane, in the order Ar > Kr >Xe. The large increase in yields from 25 eV onwards for all rare gases is likely due to the formation and decay of two-hole states on neighboring rare gas and ethane molecules due to interatomic and intermolecular Coulomb decay (ICD) and not electron transfer mediated decay (ETMD). The ICD and ETMD pathways become accessible when the incoming electron has sufficient energy to excite the inner valence ns level of the rare gas to a Rydberg state or ionize it. The experimental findings are supported by calculations of thresholds, density of states for the final configurations of these processes, and coupling strengths for hole transfer between ethane and rare gases. The fragment ion branching ratios vary with energy from threshold to about 35 eV, showing the fragmentation pattern changes with the mode of hole transfer and availability of excess energy. Sigma C–C bonds are more likely to break than C–H bonds in the mid-20 eV range, and this effect is most pronounced for Xe, followed by Kr, and then Ar.
Site-selective probing of iodine 4d orbitals at 13.1 nm was used to characterize the photolysis of CH 2 I 2 and CH 2 BrI initiated at 202.5 nm. Time-dependent fragment ion momenta were recorded using Coulomb explosion imaging mass spectrometry and used to determine the structural dynamics of the dissociating molecules. Correlations between these fragment momenta, as well as the onset times of electron transfer reactions between them, indicate that each molecule can undergo neutral three-body photolysis. For CH 2 I 2 , the structural evolution of the neutral molecule was simultaneously characterized along the C–I and I–C–I coordinates, demonstrating the sensitivity of these measurements to nuclear motion along multiple degrees of freedom.
The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.
Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.
We have conducted an experimental study on the photo double ionization (PDI) of carbon dioxide dimers and oxygen dimers, while focusing on the dissociation dynamics upon single photon absorption. The results in terms of the kinetic energy and angular distributions of the charged particles show unambiguous experimental evidence of intermolecular Coulombic decay (ICD) in carbon dioxide dimers. In the oxygen dimer, the results show that ICD is accompanied by knock-off ionization mechanisms.
First-principles real-time time-dependent density-functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, which has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can be comparable to shift currents under experimentally accessible conditions.
Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.
We present a versatile cumulant mapping algorithm for analyzing correlated particle emission, offering insights into complex electronic and nuclear dynamics. Recently, we have demonstrated the use of cumulant mapping to extract information-rich correlations between the momenta of multiple fragments produced in Coulomb explosion imaging experiments [C. Cheng et al., Phys. Rev. Lett. 130, 093001 (2023)]. We define cumulant mapping in terms of histograms, enabling fast computation of linear (additive) observables. However, applying the same algorithm to nonlinear (nonadditive) observables poses challenges, as the computation time of conventional estimators scales nonlinearly with data size. To overcome this, we develop estimators and an accompanying algorithm to enable computationally efficient estimation of the cumulant of interest. Comparisons of computation times and signal-to-noise ratios reveal the superior performance of our approach. This method is demonstrated on the (D+, D+, C+, O+) dissociation channel of CD 2 O 4+ produced in a strong-field ionization experiment. Additionally, Poisson statistics are used to simulate the two methods and provide insights into the efficiency of our algorithm. The proposed methodology unlocks efficient computation of cumulant mapping for a broader range of complex systems and observables, such as the laser pulse dependence of ionization dynamics.
Here, we consider a screened Coulomb interaction between electrons in graphene and determine their dynamic response functions, such as a longitudinal and a transverse electric conductivity and a polarization function and compare them to the corresponding quantities in the short-range interaction model. The calculations are performed to all orders for short-range interaction by taking into account the self-energy renormalization of the electron velocity and using a ladder approximation to account for the vertex corrections, ensuring that the Ward identity (charge conservation law) is satisfied. Our findings predict a resonant response of interacting electron-hole pairs at a particular frequency below the threshold $\textit{qv = ω}$ and further predict an instability for sufficiently strong interactions.
We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the 4 He and 16 O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for 4He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and experimental data available in the literature, obtaining a fair agreement. This is a first and necessary step towards initiating a program for computing neutrino-nucleus interactions from first principles in coupled-cluster theory and supporting the experimental long-baseline neutrino program with a state-of-the-art theory that can reach medium-mass nuclei.
In this work, the shape and collectivity of 106 Cd was investigated via a sub-barrier-energy Coulomb excitation experiment performed at the NSCL ReA3 facility using the JANUS setup. Transition matrix elements between low-lying states were found to agree with adopted values, and information on the shape and collectivity of higher-lying states was extracted for the first time. Locally-optimized large-scale shell-model calculations were found to describe well the B(E2) transition strengths but failed to reproduce the spectroscopic quadrupole moments Q s . An analysis of the E2 rotational invariants and the normalized quadrupole moment q s indicates that this may be due to a significant degree of triaxiality in 106 Cd which is not captured by the present shell-model calculations. Analogous calculations for the Fe isotopes (two protons below the Z = 28 magic number) reveal the critical role of high-j neutron configurations for the description of quadrupole moments in the heavy Fe and Cd isotopes (two protons below magic Z = 50), but this effect is insufficient to explain the shape of 106 Cd, posing a puzzle for the understanding of nuclear structure towards N = 50.
Background: Electric-quadrupole ( E 2 ) strengths relate to the underlying quadrupole deformation of a nucleus and present a challenge for many nuclear theories. Additionally, mirror nuclei in the vicinity of the line of N = Z represent a convenient laboratory for testing deficiencies in such models, making use of the isospin symmetry of the systems. Purpose: Uncertainties associated with literature E 2 strengths in Mg 23 are some of the largest in T z = | 1 2 | nuclei in the s d shell. The purpose of the present paper is to improve the precision with which these values are known, to enable better comparison with theoretical models. Methods: Coulomb-excitation measurements of Mg 23 and Na 23 were performed at the TRIUMF-ISAC facility using the TIGRESS spectrometer. They were used to determine the E 2 matrix elements of mixed E 2 / M 1 transitions. Results: Reduced E 2 transition strengths, B ( E 2 ) , were extracted for Mg 23 and Na 23 . Their precision was improved by factors of approximately 6 for both isotopes, while agreeing within uncertainties with previous measurements. Conclusions: A comparison was made with both shell-model and ab initio valence-space in-medium similarity renormalization group calculations. Valence-space in-medium similarity renormalization group calculations were found to underpredict the absolute E 2 strength, in agreement with previous studies.
Here, the Coulomb excitation of 102 Ru was performed with beams of 12 C and 16 O ions. The beam particles scattered at forward angles were momentum analyzed with a magnetic spectrograph. The resolution achieved enabled the populations of the $2^{+}_{1}$ state, the unresolved $2^{+}_{2}$/ $4^{+}_{1}$, and $2^{+}_{4}$/$3^{-}_{1}$, doublets of states, and the $3^{-}_{2}$ state to be determined as a function of the scattering angle. These populations are compared with gosia calculations, yielding B(E2; $2^{+}_{1}$→$0^{+}_{1}$) = 41.5 ± 2.3 W.u., B(E2; $2^{+}_{2}$→$0^{+}_{1}$) = 1.75 ± 0.11 W.u., B(E3; $3^{-}_{1}$→$0^{+}_{1}$) = 31.5 ± 3.5 W.u., and B(E3; $3^{-}_{2}$→$0^{+}_{1}$) = 6.8 ± 0.5 W.u. The B(E3; $3^{-}_{1}$→$0^{+}_{1}$) value is significantly larger than previously measured. The weakly populated $2^{+}_{3}$ state, presumed to be a member of the band built on the $0^{+}_{2}$ state, was observed clearly for a single angle only, and a fit to its population results in B(E2; $2^{+}_{3}$→$0^{+}_{1}$) = 0.053 ± 0.011 W.u. Using the known γ -ray branching ratios for the $2^{+}_{3}$ level, the B(E2; $2^{+}_{3}$→$0^{+}_{2}$) value is calculated to be 18 ± 4 W.u., substantially less than the B(E2; $2^{+}_{1}$→$0^{+}_{1}$) . This suggests that the deformation of the $0^{+}_{2}$ state is lower than that of the $0^{+}_{1}$ state. The results are compared with beyond-mean-field calculations with the Gogny-D1S interaction using the symmetry-conserving configuration-mixing method.