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.
This erratum corrects a factor 2 error in previous work on current-dependent TDDFT.
When compared with the National Ignition Facility’s (NIF) original soft x-ray opacity spectrometer, which used a convex cylindrical design, an elliptically shaped design has helped to increase the signal-to-noise ratio and eliminated nearly all reflections from alternate crystal planes. The success of the elliptical geometry in the opacity experiments has driven a new elliptical geometry crystal with a spectral range covering 520–1100 eV. When coupled with the primary elliptical geometry, which spans 1000–2100 eV, the new sub-keV elliptical geometry helps to cover the full iron L-shell and major oxygen transitions important to solar opacity experimentation. The new design has been built and tested by using a Henke x-ray source and shows the desired spectral coverage. Additional plans are underway to expand these opacity measurements into a mode of time-resolved detection, ~1 ns gated, but considerations for the detector size and photometrics mean a crystal geometry redesign. Here, the new low-energy geometry, including preliminary results from the NIF opacity experiments, is presented along with the expansion plans into a time-resolved platform.
It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.
Spatial distributions of electrons ionized and scattered from ultra-low-pressure gases are proposed and experimentally demonstrated as a method to directly measure the intensity of an ultra-high-intensity laser pulse. Analytic models relating the peak scattered electron energy to the peak laser intensity are derived and compared to paraxial Runge–Kutta simulations highlighting two models suitable for describing electrons scattered from weakly paraxial beams (f#>5) for intensities in the range of 1018−1021 W cm−2. Scattering energies are shown to be dependent on gas species, emphasizing the need for specific gases for given intensity ranges. Direct measurements of the laser intensity at full power of two laser systems are demonstrated, both showing a good agreement between indirect methods of intensity measurement and the proposed method. One experiment exhibited the role of spatial aberrations in the scattered electron distribution, motivating a qualitative study on the effect. We propose the use of convolutional neural networks as a method for extracting quantitative information on the spatial structure of the laser at full power. We believe the presented technique to be a powerful tool that can be immediately implemented in many high-power laser facilities worldwide.
Computational optimal transport is used to analyze the difference between pairs of continuous molecular spectra. It is demonstrated that transport distances which are derived from this approach may be a more appropriate measure of the difference between two continuous spectra than more familiar measures of distance under many common circumstances. Associated with the transport distances is the transport map which provides a detailed analysis of the difference between two molecular spectra and is a key component of our study of quantitative differences between two continuous spectra. The use of optimal transport for comparing molecular spectra is developed in detail here with a set of model spectra, so that the discussion is self-contained. The difference between the transport distance and more common definitions of distance is elucidated for some well-chosen examples and it is shown where transport distances may be very useful alternatives to standard definitions of distance. The transport distance between a theoretical and experimental electronic absorption spectrum for SO 2 is studied and it is shown how the theoretical spectrum can be modified to fit the experimental spectrum better adjusting the theoretical band origin and the resolution of the theoretical spectrum. In conclusion, this analysis includes the calculation of transport maps between the theoretical and experimental spectra suggesting future applications of the methodology.
We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.
Explore the source record for details and available documents.
We have epitaxially stabilized a series of monoclinic (Al x Ga 1-x-y In y ) 2 O 3 alloys by careful choice of molecular beam epitaxy growth conditions, which balance alloy growth with suboxide desorption. The films are pseudomorphic to (010) β-Ga 2 O 3 substrates at thicknesses up to 150 nm with compositions ranging from (Al 0.01 Ga 0.83 In 0.16 ) 2 O 3 to (Al 0.24 Ga 0.75 In 0.03 ) 2 O 3 . The absorption edge shifts from approximately 4.62-5.14 eV with coincidently increasing Al and decreasing In mole fractions. J-V measurements reveal an increase in resistivity over four orders of magnitude with a maximum value of 4.2 x 10 5 Ω-cm for (Al 0.17 Ga 0.76 In 0.07 ) 2 O 3 , which has nearly identical lattice parameters (both in-plane and out-of-plane) to the underlying β-Ga 2 O 3 . Scanning transmission electron microscopy of this sample reveals a mostly uniform and single crystalline film, though we identify areas of non-uniform In incorporation and some γ-phase inclusions. This work demonstrates the feasibility of thick layers lattice-matched to β-Ga 2 O 3 with increased bandgap compared to phase-separation limited (Al,Ga) 2 O 3 . These alloys can enable higher bandgap epitaxial dielectrics and high sheet charge density transistors by increasing the conduction band offset with respect to β-Ga 2 O 3 .
Molecular dynamics (MD) simulations are a powerful tool for modeling warm and hot dense matter. Density functional theory (DFT) MD simulations are often preferred in dense plasmas in order to accurately model quantum electronic structure. However, DFT-MD simulations neglect interaction effects due to fluctuations in excited states. In this work, we present an MD approach that uses excited state method pseudoatoms to run dense plasma simulations with many different core-electron configurations at classical MD speeds. We also allow for transitions between different configurations in our simulations and find that these fluctuations are especially important for highly excited states. Our results suggest that finite configuration lifetimes that are comparable to the inverse ion plasma frequency need to be accounted for in order to accurately model ion distributions in dense plasma simulations. We also demonstrate that excited state fluctuations have a direct impact on ion plasma microfields, generate different plasma microfields for different excitation levels, and thereby induce absorption–emission line shape asymmetries even in steady-state plasmas.
In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.
A statistical approach to treating energy transfer in low temperature scattering was found to capture some quantum effects.
We simulate vibrational strong coupling (VSC) and vibrational ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O–H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Furthermore, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H 2 O molecules especially under V-USC, which could play a role in ground-state chemistry.
Predicting the rates of spin-dependent processes characterized by nonradiative transitions between electronic states with different spin multiplicities is important for understanding the mechanisms of many photochemical and catalytic reactions. To calculate these rates, it is necessary to define the spin state representation and the couplings between these states that drives the interstate transitions. In this work, we describe three different approaches to calculating the spin-orbit coupling (SOC), transition probabilities, and rate constants between the MS components of the electronic states with different spin multiplicities. We implemented these approaches in our nonadiabatic statistical theory (NAST) software package, which predicts the transition probabilities and rate constants of spin-dependent processes using information obtained from electronic structure calculations. Here, we discuss the advantage and drawbacks of each approach and, as an example, calculate the rate constants for transitions between the spin states in the active site model of the protein rubredoxin.
In this paper, we present time-resolved ultraviolet-pump x-ray probe Auger spectra of 2-thiouracil. An ultraviolet induced shift towards higher kinetic energies is observed in the sulfur 2p Auger decay. The difference Auger spectra of pumped and unpumped molecules exhibit ultrafast dynamics in the shift amplitude, in which three phases can be recognized. In the first 100 fs, a shift towards higher kinetic energies is observed, followed by a 400 fs shift back to lower kinetic energies and a 1 ps shift again to higher kinetic energies. We use a simple Coulomb-model, aided by quantum chemical calculations of potential energy states, to deduce a C-S bond expansion within the first 100 fs. The bond elongation triggers internal conversion from the photoexcited S 2 to the S 1 state. Based on timescales, the subsequent dynamics can be interpreted in terms of S 1 nuclear relaxation and S 1 -triplet internal conversion.
Here, endofullerene photoionization of open-shell N@C 60 is investigated using a coupled-channel, correlated R-matrix approach with additional modifications for a finite spherical-shell fullerene potential. Particular focus is on the initial-state dependence of the photoionization cross section. Comparisons show significant differences between the 1s 2 2s 2 2p 3 ( 4 S) initial state and the metastable 2 D and 2 P states, as well as strong confinement oscillations.
The prediction of electronic structure for strongly correlated molecules represents a promising application for near-term quantum computers. Significant attention has been paid to ground state wavefunctions, but excited states of molecules are relatively unexplored. In this work, we consider the adaptive, problem-tailored (ADAPT)-variational quantum eigensolver (VQE) algorithm, a single-reference approach for obtaining ground states, and its state-averaged generalization for computing multiple states at once. We demonstrate for both rectangular and linear H4, as well as for BeH2, that this approach, which we call multistate-objective, Ritz-eigenspectral (MORE)-ADAPT-VQE, can make better use of small excitation manifolds than an analogous method based on a single-reference ADAPT-VQE calculation, q-sc-EOM. In particular, MORE-ADAPT-VQE is able to accurately describe both avoided crossings and crossings between states of different symmetries. In addition to more accurate excited state energies, MORE-ADAPT-VQE can recover accurate transition dipole moments in situations where traditional ADAPT-VQE and q-sc-EOM struggle. These improvements suggest a promising direction toward the use of quantum computers for difficult excited state problems.
High-Z atoms are excellent laboratories to study the combination of relativistic and many-electron correlation effects in electronic structure and dynamics. In the present work, the relativistic-random-phase approximation at different levels of truncation is employed to explore the final-state correlation effects in the photoelectron dynamics of the heaviest known element in the Periodic Table: oganesson (Og) with Z = 118. Here, the focus of this work is to illustrate the relativistic effects resulting from coupling different photoionization channels arising from spin-orbit split subshells, termed the spin-orbit-interaction-activated interchannel coupling effect. Comparison with the photoelectron dynamics of Rn on a qualitative level is also carried out, since Og is a homologue of Rn. The photoelectron dynamics of 7p, 7s, and 6d subshells are investigated and differences between Og and Rn, due to enhanced relativistic effects, are investigated.
Here, a study of Cooper minima (CM) arising from the photoionization of 6$\textit{s}$, 5$\textit{s}$, and 4$\textit{s}$ subshells of high-$\textit{Z}$ atoms has been performed using Dirac-Fock (DF), two-channel relativistic-random-phase approximation (RRPA), and fully coupled RRPA. The results show huge splittings between $ns → εp_{3/2}$ and $ns → εp_{1/2}$ CM which increase with $\textit{Z}$ owing primarily to the relativistic interactions (spin orbit) that are attractive for the $εp_{1/2}$ final state but repulsive for the corresponding $εp_{3/2}$. In addition, it was found that correlation in the form of interchannel coupling (essentially configuration interaction in the final continuum states) plays a huge role in determining the location of the CM. For 6$\textit{s}$ photoionization, the 6$s → εp_{3/2}$ and 6$s → εp_{1/2}$ CM behave completely different as a function of $\textit{Z}$ . It was also found that for 5$\textit{s}$ and 4$\textit{s}$ photoionization, the CM move below the threshold, with increasing $\textit{Z}$ , and, at high enough $\textit{Z}$ , the 5s → εp3/2 and 4s → εp3/2 CM re-emerge into the continuum. The calculations have been carried out for the $\textit{ns}$ subshells of Hg ($\textit{Z}$ = 80), Rn ($\textit{Z}$ = 86), Ra ($\textit{Z}$ = 88), No ($\textit{Z}$ = 102), Cn ($\textit{Z}$ = 112), and Og ($\textit{Z}$ = 118).