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

BEYONDPLANCK XIV. Polarized foreground emission between 30 and 70 GHz

We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequency channel. This analysis represents the first implementation of true multi-resolution component separation applied to CMB observations for both amplitude and spectral energy distribution (SED) parameters. For synchrotron emission, we approximate the SED as a power-law in frequency and find that the low signal-to-noise ratio of the current data strongly limits the number of free parameters that may be robustly constrained. We partition the sky into four large disjoint regions (High Latitude; Galactic Spur; Galactic Plane; and Galactic Center), each associated with its own power-law index. We find that the High Latitude region is prior-dominated, while the Galactic Center region is contaminated by residual instrumental systematics. The two remaining regions appear to be signal-dominated, and for these we derive spectral indices of $β$ S Spur = -3.17 ± 0.06 and $β$ S Plane = -3.03 ± 0.07, in good agreement with previous results. For thermal dust emission we assume a modified blackbody model and we fit a single power-law index across the full sky. We find $β$ d = 1.64 ± 0.03, which is slightly steeper than reported from Planck HFI data, but still statistically consistent at the 2σ confidence level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

From BEYONDPLANCK to COSMOGLOBE: Preliminary WMAP Q -band analysis

We present the first application of the COSMOGLOBE analysis framework by analyzing nine-year WMAP time-ordered observations that uses similar machinery to that of BEYONDPLANCK for the Planck Low Frequency Instrument (LFI). We analyzed only the Q-band (41 GHz) data and report on the low-level analysis process based on uncalibrated time-ordered data to calibrated maps. Most of the existing BEYONDPLANCK pipeline may be reused for WMAP analysis with minimal changes to the existing codebase. The main modification is the implementation of the same preconditioned biconjugate gradient mapmaker used by the WMAP team. Producing a single WMAP Q1-band sample requires 22 CPU-hrs, which is slightly more than the cost of a Planck 44 GHz sample of 17 CPU-hrs; this demonstrates that a full end-to-end Bayesian processing of the WMAP data is computationally feasible. In general, our recovered maps are very similar to the maps released by the WMAP team, although with two notable differences. In terms of temperature, we find a ~2 μK quadrupole difference that most likely is caused by different gain modeling, while in polarization we find a distinct 2.5 μK signal that has been previously referred to as poorly measured modes by the WMAP team. In the COSMOGLOBE processing, this pattern arises from temperature-to-polarization leakage from the coupling between the CMB Solar dipole, transmission imbalance, and sidelobes. No traces of this pattern are found in either the frequency map or TOD residual map, suggesting that the current processing has succeeded in modeling these poorly measured modes within the assumed parametric model by using Planck information to break the sky-synchronous degeneracies inherent in the WMAP scanning strategy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spatially resolved gas and stellar kinematics in compact starburst galaxies

The kinematics of galaxies provide valuable insights into their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. We aim to identify what governs the kinematics in a sample of SDSS-selected nearby starburst galaxies, by obtaining spatially resolved measurements of the gas and stellar kinematics. We obtained near-infrared integral-field K-band spectroscopy with VLT/SINFONI for 15 compact starburst galaxies. We derived the integrated as well as spatially resolved stellar and gas kinematics. The stellar kinematics were derived from the CO absorption bands, and Paα and Brγ emission lines were used for the gas kinematics. Based on the integrated spectra, we find that the majority of galaxies have gas and stellar velocity dispersion that are comparable. A spatially resolved comparison shows that the six galaxies that deviate show evidence for a bulge or stellar feedback. Two galaxies are identified as mergers based on their double-peaked emission lines. In our sample, we find a negative correlation between the ratio of the rotational velocity over the velocity dispersion (v rot /σ) and the star formation rate surface density. We propose a scenario where the global kinematics of the galaxies are determined by gravitational instabilities that affect both the stars and gas. This process could be driven by mergers or accretion events. Effects of stellar feedback on the ionised gas are more localised and detected only in the spatially resolved analysis. The mass derived from the velocity dispersion provides a reliable mass even if the galaxy cannot be spatially resolved. The technique used in this paper is applicable to galaxies at low and high redshift with the next generation of infrared-focussed telescopes (JWST and ELT).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling

Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions (SEDs) into the main luminosity contributors (e.g. stellar populations of distinct age and metallicity or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years, with unknown ramifications to the current understanding of galaxy evolution. The objective of this work is to estimate the relations between stellar properties (e.g. total mass, mean age, and mean metallicity) of SF galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through the more common purely stellar spectral synthesis approach. The main galaxy sample from SDSS DR7 was analysed with two distinct population synthesis codes: FADO, which estimates self-consistently both the stellar and nebular contributions to the SED, and the original version of STARLIGHT, as representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as ~0.06 dex for the overall population of galaxies and ~0.12 dex for SF galaxies (galaxies with EW(Hα) > 3 Å), with the most prominent difference between both codes in the two populations being in the light-weighted mean stellar age. FADO presents a broader range of mean stellar ages and metallicities for SF galaxies than STARLIGHT, with the latter code preferring metallicity solutions around the solar value (Z ⊙ = 0.02). A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(Hα) > 75 Å) reveals that the light contributions of simple stellar populations (SSPs) younger than ≤10 7 (10 9 ) years in STARLIGHT are higher by ~5.41% (9.11%) compared to FADO. Moreover, FADO presents higher light contributions from SSPs with metallicity ≤ Z ⊙ /200 (Z ⊙ /50) of around 8.05% (13.51%) when compared with STARLIGHT. This suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to ~0.17 dex and overestimating the light-weighted metallicity by up to ~0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(Hα) and decreasing total stellar mass. Moreover, comparing SF results from FADO in a purely stellar mode with the previous results qualitatively suggests that differences between codes are primarily due to mathematical and statistical differences and secondarily due to the impact of the nebular continuum modelling approach (or lack thereof). However, it is challenging to adequately quantify the relative role of each factor since they are likely interconnected. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties (e.g. mean age and mean metallicity) of galaxies with relatively high star formation rates (e.g. dwarf spirals, ‘green peas’, and starburst galaxies). In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity, which are factors worth bearing in mind in light of future high-resolution spectroscopic surveys at higher redshifts (e.g. MOONS and 4MOST-4HS).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Disturbed, diffuse, or just missing? A global study of the H I content of Hickson compact groups

Hickson compact groups (HCGs) are dense configurations of four to ten galaxies, whose H I morphology appears to follow an evolutionary sequence of three phases, with gas initially confined to galaxies, then significant amounts spread throughout the intra-group medium, and finally with almost no gas remaining in the galaxies themselves. It has also been suggested that several groups may harbour a diffuse H I component that is resolved out by interferometric observations. The H I deficiency of HCGs is expected to increase as the H I morphological phase progresses along the evolutionary sequence. If this is the case, H I deficiency would be a rough proxy for the age and evolutionary state of a HCG. We aim to test this hypothesis for the first time using a large sample of HCGs and to investigate the evidence for diffuse H I in HCGs. We performed a uniform reduction of all publicly available VLA H I observations (38 HCGs) with a purpose-built pipeline that also maximises the reproducibility of this study. The resulting H I data cubes were then analysed with the latest software tools to perform a manual separation of emission features into those belonging to galaxies and those extending into the intra-group medium. We thereby classified the H I morphological phase of each group as well as quantified their H I deficiency compared to galaxies in isolation. We find little evidence that H I deficiency can be used as a proxy for the evolutionary phase of a compact group in either of the first two phases, with the distribution of H I deficiency being consistent in both. However, for the final phase, the distribution clearly shifts to high H I deficiencies, with more than 90% of the expected H I content typically missing. Across all HCGs studied, we identify a few cases where there is strong evidence for a diffuse gas component in the intra-group medium, which might be detectable with improved observations. We also classify a new sub-phase where groups contain a lone H I -bearing galaxy, but are otherwise devoid of gas. The new morphological phase we have identified is likely the result of an evolved, gas-poor group acquiring a new, gas-rich member. The large spread of H I deficiencies in the first two morphological phases suggests that there is a broad range of initial H I content in HCGs, which is perhaps influenced by large-scale environment, and that the timescale for morphological changes is, in general, considerably shorter than the timescale for the destruction or consumption of neutral gas in these systems.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗