Structure of dynamical condensation fronts in the interstellar medium

We describe a technique that for the first time measures star formation rates (SFRs) in damped Lyα systems (DLAs) directly. We assume that massive stars form in DLAs and that the far-ultraviolet (FUV) radiation they emit heats the gas by the grain photoelectric mechanism. We infer the heating rate by equating it to the cooling rate measured by the strength of C II* λ1335.7 absorption. Since the heating rate is proportional to the product of the dust-to-gas ratio, the grain photoelectric heating efficiency, and the SFR per unit area, , we can deduce for DLAs in which the cooling rate and dust-to-gas ratio have been measured. We consider models in which the dust consists of carbonaceous grains and silicate grains. We present two-phase models in which the cold neutral medium (CNM) and warm neutral medium (WNM) are in pressure equilibrium. In the CNM model the line of sight passes through CNM and WNM gas, while in the WNM model the line of sight passes only through WNM gas. Since the grain photoelectric heating efficiency is at least an order of magnitude higher in the CNM than in the WNM, most of the C II* absorption arises in the CNM in the CNM model. We use the measured C II* absorption lines to derive for a sample of ≈30 DLAs in which κ has been inferred from element depletion patterns. We show that the inferred corresponds to an average over the star-forming volume of the DLA rather than to local star formation along the line of sight. We obtain the average and show that = 10-2.2 M☉ yr-1 kpc-2 for the CNM solution and = 10-1.3 M☉ yr-1 kpc-2 for the WNM solution. Interestingly, the SFR per unit area in the CNM solution is similar to that measured in the Milky Way interstellar medium.

Utilizing Planck polarized dust emission maps at 353 GHz and large-area maps of the neutral hydrogen (H i) cold neutral medium (CNM) fraction (f CNM), we investigate the relationship between dust polarization fraction (p 353) and f CNM in the diffuse high latitude ( b>30° ) sky. We find that the correlation between p 353 and f CNM is qualitatively distinct from the p 353–H i column density (N H i ) relationship. At low column densities (N H i < 4 × 1020 cm−2) where p 353 and N H i are uncorrelated, there is a strong positive p 353–f CNM correlation. We fit the p 353–f CNM correlation with data-driven models to constrain the degree of magnetic field disorder between phases along the line of sight. We argue that an increased magnetic field disorder in the warm neutral medium (WNM) relative to the CNM best explains the positive p 353–f CNM correlation in diffuse regions. Modeling the CNM-associated dust column as being maximally polarized, with a polarization fraction p CNM ∼ 0.2, we find that the best-fit mean polarization fraction in the WNM-associated dust column is 0.22p CNM. The model further suggests that a significant f CNM-correlated fraction of the non-CNM column (an additional 18.4% of the H i mass on average) is also more magnetically ordered, and we speculate that the additional column is associated with the unstable medium. Our results constitute a new large-area constraint on the average relative disorder of magnetic fields between the neutral phases of the interstellar medium, and are consistent with the physical picture of a more magnetically aligned CNM column forming out of a disordered WNM.

The ubiquity of very thin and lengthy cold neutral medium (CNM) has been reported by multiple authors in the H i community. Yet, the reason of how the CNM can be so long and lengthy is still in debate. In this paper, we recognize a new type of instability due to the attractive nature of the pressure force in the unstable phase. We provide a new estimation of the average CNM filament aspect ratio with the consideration of force balances at the phase boundary, which is roughly 5–20 in common CNM environment. We show that most of the cold filaments are less filamentary than what usually predicted via MHD turbulence theory or inferred from observations: The average length of CNM filament is roughly 1/2 of that in isothermal MHD turbulence with similar turbulence conditions. This suggests that the ‘cold filaments’ that are identified in observations might not be in pressure equilibrium or generated via other mechanisms.

Optical depth variations in the Galactic neutral interstellar medium (ISM) with spatial scales from hundreds to thousands of astronomical units have been observed through H i absorption against pulsars and continuum sources, while extremely small structures with spatial scales of tens of astronomical units remain largely unexplored. The nature and formation of such tiny-scale atomic structures (TSAS) need to be better understood. We report a tentative detection of TSAS with a signal-to-noise ratio of 3.2 toward PSR B1557−50 in the second epoch of two Parkes sessions just 0.36 yr apart, which are the closest-spaced spectral observations toward this pulsar. One absorption component showing marginal variations has been identified. Based on the pulsar’s proper motion of 14 mas yr−1 and the component’s kinematic distance of 3.3 kpc, the corresponding characteristic spatial scale is 17 au, which is among the smallest sizes of known TSAS. Assuming a similar line-of-sight (LOS) depth, the tentative TSAS cloud detected here is overdense and overpressured relative to the cold neutral medium (CNM), and can radiatively cool fast enough to be in thermal equilibrium with the ambient environment. We find that turbulence is not sufficient to confine the overpressured TSAS. We explore the LOS elongation that would be required for the tentative TSAS to be at the canonical CNM pressure, and find that it is ∼5000—much larger than filaments observed in the ISM. We see some evidence of line width and temperature variations in the CNM components observed at the two epochs, as predicted by models of TSAS-like cloud formation colliding warm neutral medium flows.

Cold neutral hydrogen (H i) is a crucial precursor for molecular gas formation and can be studied via H i absorption. This study investigates H i absorption in low column density regions of the Small and Large Magellanic Clouds (SMC and LMC) using the Galactic-ASKAP H i (GASKAP-H i) survey, conducted by the Australian Square Kilometer Array Pathfinder (ASKAP). We select 10 SMC directions in the outer regions and 18 LMC directions, with four in the outskirts and 14 within the main disk. Using the radiative transfer method, we decompose the emission and absorption spectra into individual cold neutral medium (CNM) and warm neutral medium (WNM) components. In the SMC, we find H i peak optical depths of 0.09–1.16, spin temperatures of ∼20–50 K, and CNM fractions of 1%–11%. In the LMC, optical depths range from 0.03–3.55, spin temperatures from ∼10–100 K, and CNM fractions from 1%–100%. The SMC’s low CNM fractions likely result from its low metallicity and large LOS depth. Additionally, the SMC’s outskirts show lower CNM fractions than the main body, potentially due to increased CNM evaporation influenced by the hot Magellanic Corona. Shell motions dominate the kinematics of the majority of CNM clouds in this study and likely supply cold H i to the Magellanic Stream. In the LMC, high CNM fraction clouds are found near supergiant shells, where thermal instability induced by stellar feedback promotes WNM-to-CNM transition. Although no carbon monoxide has been detected, enhanced dust shielding in these areas helps maintain the cold H i.

Obscuration of quasars by accreted gas and dust, or dusty intervening galaxies, can cause active galactic nuclei (AGN) to be missed in optically selected surveys. Radio observations can overcome this dust bias. In particular, radio surveys searching for H i absorption inform us on how the AGN can impact on the cold neutral gas medium within the host galaxy, or the population of intervening galaxies through the observed line of sight gas kinematics. We present the results of an H i absorption line survey at 0.4 &amp;lt; z &amp;lt; 1 towards 34 obscured quasars with the Australian SKA Pathfinder (ASKAP) commissioning array. We detect three H i absorption lines, with one of these systems previously unknown. Through optical follow-up for two sources, we find that in all detections the H i gas is associated with the AGN, and hence that these AGN are obscured by material within their host galaxies. Most of our sample are compact, and in addition, are either gigahertz peaked spectrum (GPS), or steep spectrum (CSS) sources, both thought to represent young or recently re-triggered radio AGN. The radio spectral energy distribution classifications for our sample agree with galaxy evolution models in which the obscured AGN has only recently become active. Our associated H i detection rate for GPS and compact SS sources matches those of other surveys towards such sources. We also find shallow and asymmetric H i absorption features, which agrees with previous findings that the cold neutral medium in compact radio galaxies is typically kinematically disturbed by the AGN.

We study the correlation between the non-thermal velocity dispersion ( $\sigma_{nth}$ ) and the length scale (L) in the neutral interstellar medium (ISM) using a large number of Hi gas components taken from various published Hi surveys and previous Hi studies. We notice that above the length-scale (L) of 0.40 pc, there is a power-law relationship between $\sigma_{nth}$ and L. However, below 0.40 pc, there is a break in the power law, where $\sigma_{nth}$ is not significantly correlated with L. It has been observed from the Markov chain Monte Carlo (MCMC) method that for the dataset of L $\gt$ 0.40 pc, the most probable values of intensity (A) and power-law index (p) are 1.14 and 0.55, respectively. Result of p suggests that the power law is steeper than the standard Kolmogorov law of turbulence. This is due to the dominance of clouds in the cold neutral medium. This is even more clear when we separate the clouds into two categories: one for L is $\gt$ 0.40 pc and the kinetic temperature ( $T_{k}$ ) is $\lt$ 250 K, which are in the cold neutral medium (CNM) and for other one where L is $\gt$ 0.40 pc and $T_{k}$ is between 250 and 5 000 K, which are in the thermally unstable phase (UNM). Most probable values of A and p are 1.14 and 0.67, respectively, in the CNM phase and 1.01 and 0.52, respectively, in the UNM phase. A greater number of data points is effective for the UNM phase in constructing a more accurate estimate of A and p, since most of the clouds in the UNM phase lie below 500 K. However, from the value of p in the CNM phase, it appears that there is a significant difference from the Kolmogorov scaling, which can be attributed to a shock-dominated medium.

We performed one-dimensional hydrodynamic simulations with detailed cooling, heating and chemical processes to examine the thermal stability of shocked gas in cold neutral medium (CNM) and molecular clouds. We find that both CNM and molecular clouds can be thermally unstable in the cooling layer behind the shock wave. The characteristic wavelength of the thermal instability ranges from 10^-5 pc to 0.1 pc in the CNM, and from 10^-7 pc to 0.1 pc in the molecular clouds. This coincides with the size of observed tiny scale structures in the CNM and molecular clouds, indicating that the thermal instability in the shocked gas could be a formation mechanism of these tiny structures in the interstellar medium. We have also calculated the e-folding number of the thermal instability to estimate the amplification of the density fluctuation in the shocked gas. Density perturbations in the CNM grow by a factor of exp(5)~150, whereas the perturbations in the molecular clouds grow only by a factor of a few behind a high Mach number shock. The amplification factor is larger at lower densities and higher velocities. Formation of very small-scale structures by thermal instability in shocked gas is more effective in lower densities.

Neutral hydrogen (H I) emission exhibits complex morphology that encodes rich information about the physics of the interstellar medium. We apply the scattering transform (ST) to characterize the H I emission structure via a set of compact and interpretable coefficients, and find a connection between the H I emission morphology and H I cold neutral medium (CNM) phase content. Where H I absorption measurements are unavailable, the H I phase structure is typically estimated from the emission via spectral line decomposition. Here, we present a new probe of the CNM content using measures that are solely derived from H I emission spatial information. We apply the ST to GALFA-H I data at high Galactic latitudes (), and compare the resulting coefficients to CNM fraction measurements derived from archival H I emission and absorption spectra. We quantify the correlation between the ST coefficients and the measured CNM fraction (f CNM), finding that the H I emission morphology encodes substantial f CNM-correlating information and that ST-based metrics for small-scale linearity are particularly predictive of f CNM. This is further corroborated by the enhancement of the I 857/N HI ratio with larger ST measures of small-scale linearity. These results are consistent with the picture of regions with higher CNM content being more populated by small-scale filamentary H I structures. Our work illustrates a physical connection between the H I morphology and phase content, and suggests that future phase decomposition methods can be improved by making use of both H I spectral and spatial information.

We present the first localized detections of the cold neutral medium (CNM) in IC 10, offering a rare view of dense atomic gas in a low-metallicity ( Z / Z ⊙ ∼ 0.27) dwarf galaxy. As a low-metallicity starburst, IC 10’s interstellar medium conditions could reflect small scale physics conditions that mirror those of early galaxies, providing a unique window into the heating and cooling processes that shaped the interstellar medium in early-Universe environments. Leveraging the high angular (&lt;5″ ∼ 15 pc) and spectral (0.4 km s −1 ) resolution of the Local Group L -band Survey, we searched for H I absorption against nine continuum radio sources and detected absorption along three sightlines corresponding to internal radio emission sources within IC 10. Using Gaussian decomposition and radiative transfer, we characterize the CNM, deriving spin temperatures of ∼30–55 K, column densities of (0.6–3.0) × 10 21 cm −2 , cold H I fractions of ∼21%–37%, and line widths of ∼5.6–13.6 km s −1 . For each individual detection of H I absorption, we find corresponding molecular emission from 12 CO ( J = 1–0), HCO + ( J = 1–0), and HCN ( J = 1–0) at similar velocities and with comparable line widths, indicating a well-mixed cold atomic and molecular medium. In IC 10, the CNM shows a clear kinematic connection to the high-density ISM, implying a stronger dynamical coupling with molecular gas than in the Milky Way, in line with expectations for low-metallicity environments. At the ∼15 pc scales probed by slightly extended H II regions in IC 10, unresolved CNM clouds likely contribute to line blending, so the observed broad H I line widths may partly reflect spatial and kinematic averaging.

We present 21-cm absorption measurements towards 12 radio continuum sources with previously identified thermally-unstable warm neutral medium (WNM). These observations were obtained with the Expanded Very Large Array (EVLA) and were complemented with the HI emission spectra obtained with the Arecibo Observatory. Out of 12 sources, HI absorption was detected along 5 lines of sight (seven new absorption features in total), resulting in a detection rate of ~42%. While our observations are sensitive to the WNM with a spin temperature T_s<3000 K, we detected only two wide absorption lines with T_s=400-900 K. These temperatures lie above the range allowed for the cold neutral medium (CNM) by the thermal equilbrium models and signify the thermally unstable WNM. Several absorption features have an optical depth of only a few x10^{-3}. While this is close or lower than what is theoretically expected for the CNM, we show that these weak lines are important for constraining the fraction of the thermally unstable WNM. Our observations demonstrate that, for the first time, high bandpass stability can be achieved with the VLA, allowing detection of absorption lines with a peak optical depth of ~10^{-3}.

We discuss the detectability of high-redshift galaxies via [C II] 158 μm line emission by coupling an analytic model with cosmological smoothed particle hydrodynamics (SPH) simulations that are based on the concordance Λ cold dark matter (CDM) model. Our analytic model describes a multiphase interstellar medium (ISM) irradiated by the far-ultraviolet (FUV) radiation from local star-forming regions, and it calculates thermal and ionization equilibrium between cooling and heating. The model allows us to predict the mass fraction of a cold neutral medium (CNM) embedded in a warm neutral medium (WNM). Our cosmological SPH simulations include a treatment of radiative cooling/heating, star formation, and feedback effects from supernovae and galactic winds. Using our method, we make predictions for the [C II] luminosity from high-redshift galaxies that can be directly compared with upcoming observations by the Atacama Large Millimeter Array (ALMA) and the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We find that the number density of high-redshift galaxies detectable by ALMA and SPICA via [C II] emission depends significantly on the amount of neutral gas, which is highly uncertain. Our calculations suggest that, in a CDM universe, most [C II] sources at z = 3 are faint objects with Sν < 0.01 mJy. Lyman break galaxies (LBGs) brighter than RAB = 23.5 mag are expected to have flux densities Sν = 1-3 mJy depending on the strength of galactic wind feedback. The recommended observing strategy for ALMA and SPICA is to aim at very bright LBGs or star-forming DRG/BzK galaxies.

We investigated H i absorption toward a single pulsar, PSR J1644−4559, and its variability over timescales from days to years, using Murriyang, CSIRO’s Parkes Radio Telescope. Our 19 epochs of spectral observations, spanning 1.2 yr with intervals as short as 1 day, provide the most comprehensive cadence coverage for monitoring H i absorption to date. We identified two significant detections of tiny-scale atomic structure (TSAS) with spatial scales ranging from a lower limit of ∼11 au to an upper limit of 165 au, both exhibiting integrated signal-to-noise ratios exceeding 5.0. We find a relationship between linear size and optical depth variation in the cold neutral medium (CNM) component hosting the TSAS, described by a power-law relationship, Δτ int = Δτ 0(ΔL)(α−2)/2, with α = 4.1 ± 0.4. This is the first observational evidence explicitly connecting TSAS to turbulence in CNM. This power-law index is significantly steeper than previously reported values for the CNM, where α ranges from 2.3 to 2.9, but similar to those observed in the warm ionized gas. Additionally, we observe no significant variation in α across the entire range of spatial scales traced in our study, indicating that turbulence may be cascading down and dissipating at smaller scales. While there is no precise proper-motion measurement for this pulsar, our estimates for the turbulence dissipation in the CNM place the lower and upper limits at less than 0.03 au and 0.4 au, respectively.

Radiative processes play a pivotal role in shaping the thermal and chemical states of gas across diverse astrophysical environments, from the interstellar medium (ISM) to the intergalactic medium. We present a hybrid cooling model for cosmological simulations that incorporates a comprehensive treatment of radiative processes, including parametrizations of the interstellar radiation field, cosmic ray rates, and dust physics. The model uses the chimes chemical network and combines on-the-fly non-equilibrium calculations with quasi-equilibrium cooling rates. The quasi-equilibrium rates account for the time-dependent free electron fractions of elements tracked in non-equilibrium, balancing computational efficiency with physical accuracy. We evaluate the performance under various conditions, including the thermal evolution of primordial gas at the cosmic mean density, the properties of the warm and cold neutral media in Milky Way-like galaxies, and the atomic-to-molecular hydrogen transition. We demonstrate that thermal equilibrium predictions for the neutral phases of the ISM underestimate the median gas pressures in simulations of isolated galaxies by up to 0.5 dex. Finally we find that the atomic-to-molecular hydrogen transition is shifted to lower densities by up to 1 dex if oxygen is not included in the chemical network. Our work provides a robust framework for studying the multiphase ISM and its role in galaxy formation and evolution.

Faraday tomography of polarimetric observations at low frequency in the radio is a unique tool for studying the structure of the magneto-ionic diffuse interstellar medium (ISM) based on Faraday depth. LOFAR data below 200 MHz have revealed a plethora of features in polarization, whose origin remains unknown. Previous studies have highlighted the remarkable association of such features with tracers of the magnetized-neutral ISM, such as interstellar dust and atomic hydrogen (HI). However, the physical conditions responsible for the correlation between magneto-ionic and neutral media have not been clarified yet. In this Letter we further investigate the correlation between LOFAR data and the HI spectroscopic observations at 21 cm from the Effelsberg-Bonn HI Survey (EBHIS). We focus on the multiphase properties of the HI gas. We present the first statistical study on the morphological correlation between LOFAR tomographic data and the cold (CNM), lukewarm (LNM), and warm (WNM) neutral medium HI phases. We use the Regularized Optimization for Hyper-Spectral Analysis approach to decompose the HI phases based on a Gaussian decomposition of the HI spectra. We study four fields of view – Fields 3C196, A, B, and C – and find, in at least the first two, a significant correlation between the LOFAR and EBHIS data using the histograms of oriented gradients (HOG) feature. The absence of a correlation in Fields B and C is caused by a low signal-to-noise ratio in polarization. The observed HOG correlation in Fields 3C196 and A is associated with all HI phases and it is surprisingly dominant in the CNM and LNM phases. We discuss possible mechanisms that would explain the correlation between CNM, LNM, and WNM with polarized emission at Faraday depths up to 10 rad m−2. Our results show how the complex structure of the ionic medium seen by the LOFAR tomographic data is tightly related to phase transition in the diffuse and magnetized neutral ISM traced by HI spectroscopic data.