Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “average atom”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evidence of distinct nonpolar/polar ordering at long/short ranges in relaxor ferroelectrics

Distinct atomic nonpolar/polar ordering is determined at long, intermediate, and short ranges using synchrotron x-ray diffraction, Raman scattering, and pair distribution function data in K 0.5 Na 0.5 NbO 3 -doped BaTiO 3 -based relaxor ferroelectrics, viz., (1 - x) Ba 0.9 Sr 0.1 TiO 3 -x KNN 50. Two different polar phases with distinct symmetries are identified at short/intermediate ranges for x = 0.10, where a monoclinic phase at short ranges gradually transforms into a rhombohedral phase at intermediate ranges. Here, in contrast, a nonpolar phase with average cubic symmetry is found to be stable at long ranges for a wide temperature range (100 ≤ T ≤ 500 K). The polar behavior of short/intermediate-range ordering is quantified in terms of the amplitude of ferroelectric phonon mode $Γ^-_4$ (q = 0, 0, 0) (which corresponds to the zone center of the cubic Brillouin zone). The amplitude of the $Γ^-_4$ phonon mode increases with decreasing temperature, suggesting the enhancement of ferroelectric polarization at short/intermediate ranges on lowering temperatures. Enhanced ferroelectric polarization at low temperature is also visible and evident by ferroelectrostriction and is quantified by the increase in magnitude of spontaneous volume ferroelectrostriction. Consequently, Zero thermal expansion is observed at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In situ microstructural evolution in face-centered and body-centered cubic complex concentrated solid-solution alloys under heavy ion irradiation

This study characterizes the microstructural evolution of single-phase complex concentrated solid-solution alloy (CSA) compositions under heavy ion irradiation with the goal of evaluating mechanisms for CSA radiation tolerance in advanced fission systems. Three such alloys, Cr 18 Fe 27 Mn 27 Ni 28 , Cr 15 Fe 35 Mn 15 Ni 35 , and equimolar NbTaTiV, along with reference materials (pure Ni and E90 for the Cr-FeMnNi family and pure V for NbTaTiV) were irradiated at 50 K and 773 K with 1 MeV Kr ++ ions to various levels of displacements per atom (dpa) using in-situ transmission electron microscopy. Cryogenic irradiation resulted in small defect clusters and faulted dislocation loops as large as 12 nm in face-centered cubic (FCC) CSAs. With thermal diffusion suppressed at cryogenic temperatures, defect densities were lower in all CSAs than in their less compositionally complex reference materials indicating that point defect production is reduced during the displacement cascade stage. High temperature irradiation of the two FCC CSA resulted in the formation of interstitial dislocation loops which by 2 dpa grew to an average size of 27 nm in Cr 18 Fe 27 Mn 27 Ni 28 and 10 nm in Cr 15 Fe 35 Mn 15 Ni 35 . This difference in loop growth kinetics was attributed to the difference in Mn-content due to its effect on the nucleation rate by increasing vacancy mobility or reducing the stacking-fault energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Superconductivity in a strange metal

In a gas of charged particles with the density n, mass m, and charge e, the electrical conductivity σ is given by the Drude formula: σ = ne 2 /mΓ, where Γ is the scattering rate. Standard metals are well described by Landau's Fermi Liquid (FL) theory, in which the electric current is carried by quasi-particles, low-energy excitations of the FL that resemble electrons with some effective mass m*. The scattering rate Γ= 1/τ= v F /l, where t is the (momentum) relaxation) time, l is the mean-free path, and v F is the Fermi velocity, can be expressed ("Matthiessen's Rule") as a sum of contributions from various scattering channels: Γ = Γ 0 + Γ el-el + Γ el-ph + …, where Γ 0 describes scattering on lattice imperfections, Γ el-el the electron-electron scattering, Γ el-ph the electron-phonon scattering, etc. Of these, Γ 0 = v F /l 0 , where l 0 is the average distance between the defects, is temperature-independent. Γ el-el should scale as T 2 because of Fermi statistics; for two electrons to scatter on one another, both must come from the "Debye shell" of the width k B T/E F , where k B is the Boltzmann constant and E F is the Fermi energy. Γ el-ph typically grows as T 5 , so we expect this to overwhelm the other channels at a high enough T. However, since l cannot be shorter than the distance between the atoms, the total Γ saturates at Mott-Ioffe-Regel (MIR) limit, roughly v F /a 0 , where v F is the Fermi velocity and a 0 is the lattice constant. The resistivity should also saturate at low T, at ρ 0 = m*v F /ne 2 l 0 T. FL theory also describes other electronic properties; e.g., it predicts that in the magnetic field B, the resistivity of the metal should increase with B 2 , because σ(B) = σ(B=0)/(1 + (ω c /Γ) 2 ), where ω c = eB/m*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

First experimental determination of the ⁴⁰Ar(𝑛, 2𝑛)³⁹Ar reaction cross section and ³⁹Ar production in Earth’s atmosphere

The cosmogenic ³⁹Ar(t1∕2 = 268 years) isotope of argon is used for geophysical dating and tracing of underground and ocean water, as well as ice owing to its appropriate half-life and chemical inertness as a noble gas; ³⁹Ar serves also in nuclear weapon test monitoring. We measured for the first time the total cross section of the main ³⁹Ar cosmogenic production reaction in the atmosphere, namely ⁴⁰Ar(𝑛, 2𝑛)39 Ar, using 14.8 ± 0.3 MeV neutrons. The neutrons, produced by a deuterium-tritium generator, impinged on a stainless steel sphere filled with Ar gas highly enriched in the ⁴⁰Ar isotope and were monitored by a stack of fast-neutron activation foils. The reaction yield was measured by atom counting of long-lived ³⁹Ar with noble gas accelerator mass spectrometry and, independently, by decay counting relative to atmospheric argon (³⁹Ar/Ar= 8.12 × 10−16 ). A total ⁴⁰Ar(𝑛, 2𝑛)39 Ar cross section of 610 ± 100 mb was determined at 14.8 ± 0.3 MeV incident neutron energy. This result serves as a benchmark for recent theoretical calculations and evaluations, found to reproduce well the experimental total cross section. We use these energy-dependent theoretical cross sections together with experimental spectra of cosmogenic neutrons at different altitudes to calculate the global average rate of neutron-induced ³⁹Ar atmospheric production, resulting in 770 ± 240 39 Ar atoms/cm²/day. The secular equilibrium between the ³⁹Ar calculated production rate and radioactive decay rate leads to a partial isotopic abundance ³⁹Ar/Ar = (5.9 ± 1.8) × 10−16 , showing that ≈73% of atmospheric ³⁹Ar is produced by cosmogenic neutrons, the remaining part believed to be induced by muons and high-energy 𝛾 rays. The ⁴⁰Ar(𝑛, 2𝑛)³⁹Ar cross section at 14 MeV is also a key parameter for quantifying the anthropogenic contribution to atmospheric 39 Ar produced during the thermonuclear tests of the 1960s. We estimate that anthropogenic ³⁹Ar accounts for roughly 20% of the present atmospheric inventory.

39 Ar atmospheric production↗
DOE OSTI2025

Reduced thermal resistance of Al-rich AlGaN HEMTs via top-side diamond integration

We report back-end-of-line growth of nanocrystalline diamond (NCD) on ultrawide bandgap (UWBG) high Al content aluminum gallium nitride (AlGaN) channel high electron mobility transistors for thermal management. A thin (∼15 nm) silicon nitride (SiN x ) interlayer was deposited to protect the device surface before performing a low temperature (500 °C) NCD growth process in an attempt to protect the gates on these fully fabricated devices. Notably, atomic force microscopy showed that the maximum lateral grain size exceeded 300 nm even though the film thickness was ∼250 nm. Comparing electrical (DC) performance before and after NCD growth, the gate leakage increased by ∼10 2 after NCD growth. Despite the lower NCD growth temperature, intermixing of the Ni and Au was observed in the Schottky gate metal stack; however, we believe there is another mechanism, possibly hydrogen-related, that is responsible for the measured increase in gate leakage. Regarding thermal management, the device-level thermal resistance (quantified using the average gate temperature rise measured by thermoreflectance imaging) was reduced by 29% through the incorporation of the top-side diamond film. Using time-domain thermoreflectance, the thermal conductivity of the ≈250 nm thick NCD film was measured to be 45 ± 25 W m −1 K −1 . This is expected to be at least 5× greater than the thermal conductivity of the thin disordered AlGaN alloy. There could also be a coupled electrothermal component contributing to the reduced temperature rise from electric field spreading and consequent heat spreading. This study demonstrates a promising first step toward device-level thermal management of high power UWBG Al-rich AlGaN devices.

Lundh, James Spencer [U.S. Naval Research Laborato↗
DOE OSTI2022

Electronic structure of α - RuCl 3 by fixed-node and fixed-phase diffusion Monte Carlo methods

Layered material α-RuCl 3 has caught wide attention due to its possible realization of Kitaev's spin liquid and its electronic structure that involves the interplay of electron-electron correlations and spin-orbit effects. Several DFT+U studies have suggested that both electron-electron correlations and spin-orbit effects are crucial for accurately describing the band gap. This work studies the importance of these two effects using fixed-node and fixed-phase diffusion Monte Carlo calculations both in spin-averaged and explicit spin-orbit formalisms. In the latter, the Slater-Jastrow trial function is constructed from two-component spin orbitals using our recent quantum Monte Carlo (QMC) developments and thoroughly tested effective core potentials. Our results show that the gap in the ideal crystal is already accurately described by the spin-averaged case, with the dominant role being played by the magnetic ground state with significant exchange and electron correlation effects. We find qualitative agreement between hybrid DFT, DFT+U, and QMC. In addition, QMC results agree very well with available experiments, and we identify the values of exact Fock exchange mixing that provide comparable gaps. Explicit spin-orbit QMC calculations reveal that the effect of spin-orbit coupling on the gap is minor, of the order of 0.2 eV, which corresponds to the strength of the spin orbit of the Ru atom.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into radiation resistance of titanium alloys from displacement cascade simulations

Radiation damage in beam window materials limits the use of high power proton beams in high energy physics research. The alloy Ti-6Al-4V is presently used as a beam window material but a prospective alternative, Ti-15V-3Cr-3Sn-3Al has been proposed. Since both these alloys contain dual phases at room temperature, we compare the radiation damage in the α and ß-phases of these two materials via primary knock-on atom (PKA) cascade simulations in the 10-40keV energy range. At PKA energies 30 and 40keV, the number of Frenkel pairs in the ballistic stage is higher in the ß-phase of Ti-6Al-4V than that in the ß-phase of Ti-15V-3Cr-3Sn-3Al almost by a factor of 2. The α-phase, of both these alloys, by far outperforms the ß-phases of the two alloys, both in terms of damage during the ballistic stage and in terms of the surviving defects. The average displacement threshold energy (E d ) in the α-phase of both alloys was found to be 66eV while that in the ß-phase of Ti-15-3 was 55 eV and in the ß-phase of Ti-6-4 was 46. Further, while the number of surviving defects is almost equal in both alloys, the vacancy and interstitial clustering mechanisms differ notably, which can impact the degree of radiation hardening and loss of ductility. Our simulations show larger vacancy and interstitial clusters form in Ti-6Al-4V as compared to that in Ti-15-3-3-3 alloy. These results indicate that Ti-15-3-3-3 alloys may be a promising candidate for next generation beam window material with a higher radiation tolerance than the existing Ti-6-4 alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

First-principles study of the structures and redox mechanisms of Ni-rich lithium nickel manganese cobalt oxides

To reduce the cobalt (Co) content in lithium-ion batteries, Ni-rich (high-Ni) lithium nickel manganese cobalt oxides (NMC) are pursued as one of the next-generation cathode materials. However, there is still debate on the crystal and electronic structures of the baseline, LiNiO 2 . Density Functional Theory (DFT) calculations were performed to provide a theoretical understanding of Ni-rich NMC. First, it was found that the commonly used $R\bar{3}m$ structure for LiNiO 2 is metallic, contrary to the experimentally reported mix-conducting behavior. Among the four different space groups, $R\bar{3}m$, C2/m, P2 1 /c, and P2/c, P2/c with charge disproportionation of Ni 2+ and Ni 4+ is the most energetically stable and semiconducting structure of LiNiO 2 . Therefore, the atomic structures of representative Ni-rich NMC were built by partially replacing Ni with Co or Mn in the P2/c LiNiO 2 to form Li x Ni y Mn z Co 1-y-z O 2 . In the fully lithiated (x=1.0) high Ni content NMC (y>0.5), the oxidation state of all Mn ions becomes 4+, while Co ions still maintain 3+, and part of the Ni ions become 3+ to compensate for the charge. Upon delithiation, the local environment shows more variation of the charge states on the transition metal (TM) ions. The average oxidation on each TM follows a sequence of losing electrons that starts from Ni 2+ to Ni 3+ , then oxidizing Ni 3+ and Co 3+ , while Mn 4+ remains electrochemically inactive till x=0. Finally, a general relationship for the oxidation state change in each TM as a function of x is derived and shows agreement with both modeling and experimental data.

25 ENERGY STORAGE↗
DOE OSTI2022

Optimizing the Internal Cell Structure of the Radon Monitor Based on Electrostatic Collection Method

Abstract As Rn-222 decays, an alpha particle is emitted and the residual polonium nucleus recoils in the opposite direction. At the end of the recoil path, 88% of the polonium atoms have a positive charge and 12% are neutral. The electric potential distribution in the 60 ml hemispherical internal cell of the radon monitor based on electrostatic collection is studied for reducing the combined probability of the positively charged Po-218 and the OH− produced by the ionization of water vapour in the air. The COMSOL software is used to simulate the electric potential distribution in the internal cell of the radon monitor based on the electrostatic collection method at 27°C, a pressure of 0.1 Mpa. For improving the collection efficiency of Po-218 ions, the average collection time along vertical and oblique lines is calculated when the upper surface of the internal cell is plastic, uncharged metal and charged metal, respectively. Assuming that the gas in the internal cell is uniformly distributed, the results show that if the upper surface of the hemispherical internal cell is plastic or uncharged metal, the electric potential formed in the internal cell is more uniform, and it is beneficial to reduce the total collection time of the positively charged Po-218 ions, thereby improving collection efficiency. The simulation results can be used as an effective reference for optimizing the design of the internal cell structure of the radon monitor based on electrostatic collection method.

Fan, Zhongkai (ORCID:0000000181247459)↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detecting rare neutral atomic-carbon absorbers with a deep neural network

ABSTRACT C i absorbers play an important role as indicators for exploring the presence of cold gas in the interstellar medium of galaxies. However, the current data base of C i absorbers is very limited due to their weak absorption feature and rarity. Here, we report results from a search of C i λλ1560, 1656 absorption lines using Mg ii absorbers as signposts with modified deep learning algorithms, which provides a very quick way to search for weak C i absorber candidates. A total of 107 C i absorbers were detected, which nearly doubles the size of previously known samples. In addition, we found 17 C i absorbers to be associated with 2175 Å dust absorbers (2DAs), i.e. about 16 per cent C i absorbers are associated with 2DAs. Comparing the average dust depletion patterns of C i absorbers with those of damped Lyman α absorbers (DLAs), Mg ii absorbers, Ca ii absorbers, and 2175 Å dust absorbers (2DAs) shows that C i absorbers generally have environments with more dust than DLAs, Mg ii, and Ca ii absorbers, but similar to dust in 2DAs. Similarity between the dust depletion pattern of C i absorbers to that of the warm disc in the Milky Way indicates that C i absorption clouds are possibly associated with disc components in distant galaxies. Therefore, C i absorbers are confirmed to be excellent probes to trace cold gas and dust in the Universe.

Ge, Jian↗
DOE OSTI2022

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗