Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ISM”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

BEYONDPLANCK X. Planck Low Frequency Instrument frequency maps with sample-based error propagation

We present Planck Low Frequency Instrument (LFI) frequency sky maps derived within the BEYONDPLANCK framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical, and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples, each of which corresponds to one possible realization of the various modeled instrumental systematic corrections, including correlated noise, time-variable gain, as well as far sidelobe and bandpass corrections. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products, including astrophysical component maps, angular power spectra, and cosmological parameters. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support – for the first time – full Bayesian error propagation for these effects at full angular resolution. We compared our posterior mean maps with traditional frequency maps delivered by the Planck Collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of Δg 0 /g 0 = 5 x 10 -5 , which is nominally 40 times smaller than that reported by Planck 2018. However, we also note that the original Planck 2018 estimate has a nontrivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced from the posterior samples directly, without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modeling and error propagation, and may play an important role for future Cosmic Microwave Background (CMB) B-mode experiments aiming at nanokelvin precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BeyondPlanck: I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω | d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, we identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-1/f correlated noise in the 30 and 44 GHz channels. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ 2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The need for multicomponent dust attenuation in modeling nebular emission: Constraints from SDSS-IV MaNGA

A fundamental assumption adopted in nearly every extragalactic study that analyzes optical emission lines is that the attenuation of different emission lines can be described by a single attenuation curve, scaled by a single reddening parameter, usually E(B - V). Here we show this assumption fails in many cases with important implications for derived results. We developed a new method to measure the differential nebular attenuation among three kinds of transitions: the Balmer lines of hydrogen; high-ionization transitions (> 13.6 eV) including [Ne III], [O III], and [S III]; and low-ionization transitions (≲13.6 eV) including [O II], [N II], and [S II]. This method bins the observed data in a multidimensional space spanned by attenuation-insensitive line ratios. Within each small bin, the variations in nebular parameters such as the metallicity and ionization parameter are negligible compared to the variation in the nebular attenuation. This allowed us to measure the nebular attenuation using both forbidden lines and Balmer lines. We applied this method to a sample of 2.4 million star-forming (SF) spaxels from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We found that the attenuation of high ionization lines and Balmer lines can be well described by a single Fitzpatrick (1999, PASP, 111, 63) extinction curve with R V = 3.1. However, no single attenuation curve can simultaneously account for these transitions and the derived attenuation of low-ionization lines. This strongly suggests that different lines have different effective attenuations, likely because spectroscopy at hundreds of parsecs to kiloparsecs of resolution mixes multiple physical regions that exhibit different intrinsic line ratios and different levels of attenuation. As a result, the assumption that different lines follow the same attenuation curve breaks down. Using a single attenuation curve determined by Balmer lines to correct attenuation-sensitive forbidden line ratios could bias the nebular parameters derived by 0.06–0.25 dex at A V = 1, depending on the details of the dust attenuation model. Observations of a statistically large sample of H II regions with high spatial resolutions and large spectral coverage are vital for improved modeling and deriving accurate corrections for this effect.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Limited impact of jet-induced feedback in the multi-phase nuclear interstellar medium of 4C12.50

4C12.50 (IRAS 13451+1232) at z = 0.122 is an ultraluminous infared radio galaxy that has often been proposed as a prime candidate for the link between ultraluminous infared galaxies and young radio galaxies. It is also an interesting target to investigate whether and how radio-induced feedback affects the evolution of galaxies in the early phases of radio activity. We study, in detail for the first time, the hot (≥1500 K) molecular gas in 4C12.50. The potential impact of the radio jet on this gas phase, as well as on the star formation activity, are investigated. We also study the ionised (including coronal) gas as traced by the near-infrared lines. Using near-infrared long slit spectroscopy obtained with EMIR on GTC and X-shooter on VLT, we analyse the emission line spectrum of the ionised, coronal, and, specially, the hot molecular gas in the western nucleus hosting the compact radio jet. Based on high spatial resolution ALMA CO(2–1) data, we also revise the location of 4C12.50 in the Kennicutt-Schmidt diagram in order to investigate whether star formation is suppressed. 4C12.50 hosts (2.1 ± 0.4) × 10 4 M ⊙ of hot molecular gas. An unusually high rotational temperature T rot = 3020 ± 160 K is inferred. The molecular gas mass obeys a power-law temperature distribution, dM H2 /dT ∝ T –5 , from T ~300 K and up to ~3000 K. Both results support the idea that shocks (probably induced by the radio jet) contribute to the heating and excitation of the hot molecular gas. A molecular outflow is not detected. The coupling of the outflowing ionised and neutral outflows with the hot molecular gas is poor. Contrary to other studies, we claim that there is no evidence for star formation suppression in this object. If radio-induced feedback can regulate the star formation activity in galaxies, 4C12.50 is a promising candidate to reveal this phenomenon in action. However, we find no solid evidence for a current or past impact of this mechanism on the evolution of this system, neither by clearing out the dusty central cocoon efficiently, nor by suppressing the star formation activity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multiple emission components in the Cygnus cocoon detected from Fermi -LAT observations

Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The impact of environmental effects on active galactic nuclei: A decline in the incidence of ionized outflows

Active galactic nuclei (AGNs) have generally been considered to be less frequent in denser environments due to a lower number of galaxy-galaxy interactions and/or the removal of their gas-rich reservoirs by the dense intergalactic medium. However, recent observational and theoretical works suggest that the effect of ram-pressure stripping acting on galaxies in dense environments might reduce the angular momentum of their gas, causing it to infall towards the super massive black hole at their centre, activating the AGN phase. In this work we explore the connection between environment and nuclear activity by evaluating the variation in the incidence of ionized outflows, a common phenomenon associated with nuclear activity, in AGNs across different environments. We select a sample of approximately 3300 optical AGNs from the Sloan Digital Sky Survey Data Release 13, which we match with a group catalogue to identify galaxies in isolation or residing in groups. We further probe their environments through the projected distance to the central galaxy of the group or cluster and the projected surface density to the fifth neighbour (δ 5 ). The presence of ionized outflows is determined through the modelling of the [O III] λ5007 emission line. We find that at lower masses (< 10 10.3 M ⊙ ), the fraction of ionized outflows is significantly lower in satellite AGNs (~7%) than in isolated (~22%) AGNs, probably due to their different AGN luminosity, L [O III] , in this stellar mass range. The fraction of outflows decreases at distances closer to the central galaxy of the group or cluster for all satellite AGNs; however, only the lower-mass ones (10 9 - 10 10.3 M ⊙ ) display a significant decline with δ 5 . Although this study does not include AGNs in the densest regions of galaxy clusters, our findings suggest that AGNs in dense environments accrete less gas than those in the field, potentially due to the removal of the gas reservoirs via stripping or starvation, leading to a negative connection between environment and AGN activity. Based on our results, we propose that the observed change in the incidence of outflows when moving towards denser regions of groups and clusters could contribute to the higher gas metallicities of cluster galaxies compared to field galaxies, especially at lower masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Star formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

The stellar mass (M $\star$ ) and the star formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, their history, and future evolution. Infrared (IR) photometry is widely used to measure the M $\star$ and SFR of galaxies because the near-IR traces the continuum emission of the majority of their stellar populations (SPs), and the mid/far-IR traces the dust emission powered by star-forming activity. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their SPs. It aims to provide accurate and precise IR-photometry SFR and M $\star$ calibrations that account for SP age and extinction while providing quantification of their scatter. We used the CIGALE spectral energy distribution (SED) fitting code to create model SEDs of galaxies with a wide range of star formation histories, dust content, and interstellar medium properties. We fit the relations between M $\star$ and SFR with IR and optical photometry of the model-galaxy SEDs with the Markov chain Monte Carlo (MCMC) method. As an independent confirmation of the MCMC fitting method, we performed a machine-learning random forest (RF) analysis on the same data set. The RF model yields similar results to the MCMC fits, thus validating the latter. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST NIR-F200W and MIRI-F2100W. It also provides mass-to-light ratio calibrations based on the WISE band-1, the JWST NIR-F200W, and the optical u - r or g - r colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. The proposed calibrations show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the M $\star$ calibration yields significantly better results for dusty and high-SFR galaxies where dust emission can otherwise bias the estimations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr, down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The MeerKAT Fornax Survey: II. The rapid removal of H I from dwarf galaxies in the Fornax cluster

We present MeerKAT Fornax Survey atomic hydrogen (H I) observations of the dwarf galaxies located in the central ~2.5 × 4 deg 2 of the Fornax galaxy cluster (R vir ~2°). The H I images presented in this work have a 3σ column density sensitivity between 2.7 and 50 × 10 18 cm –2 over 25 km s –1 for spatial resolution between 4 and 1 kpc. We are able to detect an impressive M HI = 5 × 10 5 M ⊙ 3σ point source with a line width of 50 km s –1 at a distance of 20 Mpc. We detected H I in 17 out of the 304 dwarfs in our field, with 14 out of the 36 late-type dwarfs (LTDs) and three out of the 268 early-type dwarfs (ETDs). The H I-detected LTDs have likely just joined the cluster and are on their first infall as they are located at large clustocentric radii, with comparable MHI and mean stellar surface brightness at fixed luminosity as blue, star-forming LTDs in the field. By contrast, the H I-detected ETDs have likely been in the cluster longer than the LTDs and acquired their H I through a recent merger or accretion from nearby H I. Eight of the H I-detected LTDs host irregular or asymmetric H I emission and disturbed or lopsided stellar emission. There are two clear cases of ram pressure shaping the H I, with the LTDs displaying compressed H I on the side closest to the cluster centre and a one-sided, starless tail pointing away from the cluster centre. The H I-detected dwarfs avoid the most massive potentials (i.e. cluster centre and massive galaxies), consistent with massive galaxies playing an active role in the removal of H I. We created a simple toy model to quantify the timescale of H I stripping in the cluster by reproducing the observed M r' –M HI relation. We find that a M HI = 10 8 M ⊙ dwarf is stripped in ~240 Myr. The model is consistent with our observations, where low-mass LTDs are directly stripped of their H I from a single encounter and more massive LTDs can harbour a disturbed H I morphology due to longer times or multiple encounters being required to fully strip their H I. This is the first time dwarf galaxies with M HI ≲ 1 × 10 6 M ⊙ have been detected and resolved beyond the local group and in a galaxy cluster.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Searching for Milky Way twins: Radial abundance distribution as a strict criterion

We search for Milky Way-like galaxies among a sample of approximately 500 galaxies. The characteristics we considered of the candidate galaxies are the following: stellar mass M $\star$ , optical radius R 25 , rotation velocity V rot , central oxygen abundance (O/H) 0 , and abundance at the optical radius (O/H) R 25 . If the values of R 25 and M $\star$ of the galaxy were close to that of the Milky Way, then the galaxy was referred to as a structural Milky Way analogue (sMWA). The oxygen abundance at a given radius of a galaxy is defined by the evolution of that region (astration level, that is, the fraction of gas converted into stars, as well as gas exchange with the surroundings), and we then assumed that the similarity of (O/H) 0 and (O/H) R 25 in two galaxies suggests a similarity in their (chemical) evolution. If the values of (O/H) 0 and (O/H) R 25 in the galaxy were close to that of the Milky Way, then the galaxy was referred to as an evolutionary Milky Way analogue (eMWA). If the galaxy was simultaneously an eMWA and sMWA, then the galaxy was considered a Milky Way twin. We find that the position of the Milky Way on the (O/H) 0 –(O/H) R 25 diagram shows a large deviation from the general trend in the sense that the (O/H) R 25 in the Milky Way is appreciably lower than in other galaxies of similar (O/H) 0 . This feature of the Milky Way evidences that its (chemical) evolution is not typical. We identify four galaxies (NGC 3521, NGC 4651, NGC 2903, and MaNGA galaxy M-8341-09101) that are simultaneously sMWA and eMWA and can therefore be considered as Milky Way twins. In previous studies, Milky Way-like galaxies were selected using structural and morphological characteristics, that is, sMWAs were selected. We find that the abundances at the centre and at the optical radius (evolutionary characteristics) provide a stricter criterion for selecting real Milky Way twins.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Structure of iso-density sets in supersonic isothermal turbulence

Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. The relations between these two-point statistics and the geometrical features of iso-density sets such as the volume fraction, the surface density, the curvature, and fractal characteristics are provided. An exact scale-by-scale budget equation is further derived revealing the role of the turbulent cascade and dilation on the structural evolution of the density field. Although applicable to many flow situations, this tool is here first invoked for characterising supersonic isothermal turbulence, using data from the currently best-resolved numerical simulation. Results: We show that iso-density sets are surface fractals rather than mass fractals, with dimensions that markedly differ between dilute, neutral, and dense regions. The surface–size relation is established for different iso-density values. We further find that the turbulent cascade of iso-density sets is directed from large towards smaller scales, in agreement with the classical picture that turbulence acts to concentrate more surface into smaller volumes. Intriguingly, there is no range of scales that complies with a constant transfer rate in the cascade, challenging our fundamental understanding of interstellar turbulence. Finally, we recast the virial theorem in a new formulation drawing an explicit relation between the aforementioned geometrical measures and the dynamics of iso-density sets.

79 ASTRONOMY AND ASTROPHYSICS↗
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 + , S2, 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 OSTI2024

Spectrum and extension of the inverse-Compton emission of the Crab Nebula from a combined Fermi -LAT and H.E.S.S. analysis

The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy γ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula’s γ-ray emission between 1 GeV and ∼100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package GAMMAPY, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing γ-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Theoretical study of the HCS + –H 2 van der Waals complex: potential energy surface, rovibrational bound states, and rotationally inelastic collisional cross sections

An accurate ground-state intermolecular potential energy surface (PES) was calculated for the HCS + –H 2 complex. The surface was constructed from 3023 ab initio energies, computed with explicitly correlated coupled-cluster theory, CCSD(T)-F12b, with extrapolation to the complete basis set limit (VTZ-F12/VQZ-F12). The new 4D PES was used to compute rovibrational energies, rota- tional constants, and close-coupling quantum scattering calculations at low collision energies. A symmetry-adapted Lanczos algorithm was used to variationally compute the low-lying rovibrational bound states. Rotational constants for states of two isomers were determined from the energy levels and reported in hopes of motivating experiments. The rotationally inelastic state-to-state cross sections of HCS + by collision with para- and ortho-H 2 were computed and compared, and trends and propensities are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An AMUSING look at the host of the periodic nuclear transient ASASSN-14ko reveals a second AGN

ABSTRACT We present Multi-Unit Spectroscopic Explorer (MUSE) integral-field spectroscopy of ESO 253−G003, which hosts a known active galactic nucleus (AGN) and the periodic nuclear transient ASASSN-14ko, observed as part of the All-weather MUse Supernova Integral-field of Nearby Galaxies survey. The MUSE observations reveal that the inner region hosts two AGN separated by $1.4\pm 0.1~\rm {kpc}$ (≈1${_{.}^{\prime\prime}}$7). The brighter nucleus has asymmetric broad permitted emission-line profiles and is associated with the archival AGN designation. The fainter nucleus does not have a broad emission-line component but exhibits other AGN characteristics, including $\hbox{$v_{\rm {FWHM}}$} \approx 700~\hbox{km~s$^{-1}$}$ forbidden line emission, $\rm{\log _{10}(\rm{[O\,\small {III}]}/\rm{H\beta})} \approx 1.1$, and high-excitation potential emission lines, such as [Fe vii] λ6086 and He ii λ4686. The host galaxy exhibits a disturbed morphology with large kpc-scale tidal features, potential outflows from both nuclei, and a likely superbubble. A circular relativistic disc model cannot reproduce the asymmetric broad emission-line profiles in the brighter nucleus, but two non-axisymmetric disc models provide good fits to the broad emission-line profiles: an elliptical disc model and a circular disc + spiral arm model. Implications for the periodic nuclear transient ASASSN-14ko are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

WALLABY pre-pilot survey: two dark clouds in the vicinity of NGC 1395

ABSTRACT We present the Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pre-pilot observations of two ‘dark’ H i sources (with H i masses of a few times 108 $\rm {M}_\odot$ and no known stellar counterpart) that reside within 363 kpc of NGC 1395, the most massive early-type galaxy in the Eridanus group of galaxies. We investigate whether these ‘dark’ H i sources have resulted from past tidal interactions or whether they are an extreme class of low surface brightness galaxies. Our results suggest that both scenarios are possible, and not mutually exclusive. The two ‘dark’ H i sources are compact, reside in relative isolation, and are more than 159 kpc away from their nearest H i-rich galaxy neighbour. Regardless of origin, the H i sizes and masses of both ‘dark’ H i sources are consistent with the H i size–mass relationship that is found in nearby low-mass galaxies, supporting the possibility that these H i sources are an extreme class of low surface brightness galaxies. We identified three analogues of candidate primordial ‘dark’ H i galaxies within the TNG100 cosmological, hydrodynamic simulation. All three model analogues are dark matter dominated, have assembled most of their mass 12–13 Gyr ago, and have not experienced much evolution until cluster infall 1–2 Gyr ago. Our WALLABY pre-pilot science results suggest that the upcoming large-area H i surveys will have a significant impact on our understanding of low surface brightness galaxies and the physical processes that shape them.

79 ASTRONOMY AND ASTROPHYSICS↗