Faraday tomography of broadband radio polarization surveys enables us to study magnetic fields and their interaction with the interstellar medium (ISM). Such surveys include the Global Magneto-Ionic Medium Survey (GMIMS), which covers the northern and southern hemispheres at sim 300--1800 MHz. In this work, we used the GMIMS High Band South (1328--1768 MHz), also named the Southern Twenty-centimeter All-sky Polarization Survey (STAPS), which observes the southern sky at a resolution of 18arcmin . To extract the key parameters of the magnetized ISM from STAPS, we computed the Faraday moments of the tomographic data cubes. These moments include the total polarized intensity, the mean Faraday depth weighted by the polarized intensity, the weighted dispersion of the Faraday spectrum, and its skewness. We compared the Faraday moments to those calculated over the same frequency range in the northern sky (using the Dominion Radio Astrophysical Observatory, DRAO), in a strip of $360 that overlaps with STAPS coverage. We find that the total polarized intensity is generally dominated by diffuse emission that decreases at longitudes of $l The Faraday moments reveal a variety of polarization structures. Low-intensity regions at high latitudes usually have a single Faraday depth component. Due to its insufficiently large frequency coverage, STAPS cannot detect Faraday thick structures. Comparing the Faraday depths from STAPS to total rotation measures from extragalactic sources suggests that STAPS frequencies are high enough that the intervening ISM causes depolarization to background emission at intermediate and high Galactic latitudes. Where they overlap, the STAPS and DRAO surveys exhibit broad correspondence but differ in polarized intensity by a factor of sim 1.8.

Aims. We investigate the geometry of the magnetic field toward the Radcliffe wave, a coherent part of the nearby Local Arm of 3 kpc in length recently discovered via three-dimensional dust mapping. Methods. We used archival stellar polarization in the optical and new measurements in the near-infrared to trace the magnetic field as projected on the plane of the sky. Our new observations cover the portion of the structure that is closest to the Sun, between Galactic longitudes of 122 degrees and 188 degrees. Results. The polarization angles of stars immediately behind the Radcliffe wave appear to be aligned with the structure as projected on the plane of the sky. The observed magnetic field configuration is inclined with respect to the Galactic disk at an angle of 18 degrees. This departure from a geometry parallel to the plane of the Galaxy is contrary to previous constraints from more distant stars and polarized dust emission. We confirm that the polarization angle of stars at larger distances shows a mean orientation parallel to the Galactic disk. Conclusions. We discuss the implications of the observed morphology of the magnetic field for models of the large-scale Galactic magnetic field, as well as formation scenarios for the Radcliffe wave itself.

Utilizing Planck polarized dust emission maps at 353 GHz and new large-area maps of the neutral hydrogen (HI) 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 degrees) sky. We find that the correlation between p(353) and f(CNM) is qualitatively distinct from the p(353)-HI column density (NHI) relationship. At low column densities (N-HI < 4 x 10(20) cm(-2)) where p(353) and N-HI 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) similar to 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 (UNM). Our results constitute a new large-area constraint on the average relative disorder of magnetic fields between the neutral phases of the ISM, and are consistent with the physical picture of a more magnetically aligned CNM column forming out of a disordered WNM.

We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of the polarization of starlight in a region of the sky of about four square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the three-dimensional (3D) volume occupied by the observed stars. We used this method to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates ( l , b ) ≈ (103.3°, 22.3°), we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes and which we identified as being an intersection of the wall of the Local Bubble and the Cepheus Flare. Farther clouds, with a distance of up to 2 kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in H I spectra in this region of the sky. We found that the orientation of the plane-of-sky component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the plane-of-sky component of the magnetic field in dusty regions. This demonstrates that the inversion of large data volumes, as expected from the P ASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization.

Context . The calibration of optical polarimeters relies on the use of stars with negligible polarization (i.e., unpolarized standard stars) for determining the instrumental polarization zero point. For wide-field polarimeters, calibration is often done by imaging the same star over multiple positions in the field of view (FoV), which is a time-consuming process. A more effective technique is to target fields containing multiple standard stars. While this method has been used for fields with highly polarized stars, there are no such sky regions with well measured unpolarized standard stars. Aims . We aim to identify sky regions with tens of stars exhibiting negligible polarization that are suitable for a zero-point calibration of wide-field polarimeters. Methods . We selected stars in regions with extremely low reddening, located at high Galactic latitudes. We targeted four ~40′ × 40′ fields in the northern and eight in the southern equatorial hemispheres. Observations were carried out at the Skinakas Observatory and the South African Astronomical Observatory. Results . We found two fields in the north and seven in the south characterized by a mean polarization lower than p < 0.1%. Conclusions . At least 9 out of the 12 fields can be used for a zero-point calibration of wide-field polarimeters.

Using 3D dust maps and Planck polarized dust emission data, we investigate the influence of the 3D geometry of the nearby interstellar medium (ISM) on the statistics of the dust polarization on large ( 80 ' ) scales. We test recent models that assume that the magnetic field probed by the polarized dust emission is preferentially tangential to the Local Bubble wall, but we do not find an imprint of the Local Bubble geometry on the dust polarization fraction. We also test the hypothesis that the complexity of the 3D dust distribution drives some of the measured variation of the dust polarization fraction. We compare sight lines with similar total column densities and find that, on average, the dust polarization fraction decreases when the dust column is substantially distributed among multiple components at different distances. Conversely, the dust polarization fraction is higher for sight lines where the dust is more concentrated in 3D space. This finding is statistically significant for the dust within 1.25 kpc, but the effect disappears if we only consider dust within 270 pc. In conclusion, we find that the extended 3D dust distribution, rather than solely the dust associated with the Local Bubble, plays a role in determining the observed dust polarization fraction at 80 '. This conclusion is consistent with a simple analytical prediction and remains robust under various modifications to the analysis. These results illuminate the relationship between the 3D geometry of the ISM and tracers of the interstellar magnetic field. We discuss implications for our understanding of the polarized dust foreground to the cosmic microwave background.

The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring formation and evolution scenarios. We include structures from galactic-scale filaments to tenth-of-a-parsec scale filaments, comprising both molecular and atomic structures, from both observational and theoretical perspectives. In addition to the literature overview, we assemble a large amount of catalogue data from different surveys and provide the most comprehensive census of filamentary structures to date. Our census consists of 22 803 filamentary structures, facilitating a holistic perspective and new insights. We use our census to conduct a meta-analysis, leading to a description of filament properties over four orders of magnitudes in length and eight in mass. Our analysis emphasises the hierarchical and dynamical nature of filamentary structures. Filaments do not live in isolation, nor they generally resemble static structures close to equilibrium. We propose that accretion during filament formation and evolution sets some of the key scaling properties of filaments. This highlights the role of accretion during filament formation and evolution and also in setting the initial conditions for star formation. Overall, the study of filamentary structures during the past decade has been observationally driven. While great progress has been made on measuring the basic properties of filaments, our understanding of their formation and evolution is clearly lacking. In this context, we identify a number of directions and questions we consider most pressing for the field.

ABSTRACT We measure the mean-squared polarization fraction of a sample of 6282 Galactic cold clumps at 353 GHz, consisting of Planck Galactic cold clump (PGCC) catalogue category 1 objects [flux densities measured with signal-to-noise ratio (S/N) > 4]. At 353 GHz, we find the mean-squared polarization fraction, which we define as the mean-squared polarization divided by the mean-squared intensity, to be (4.79 ± 0.44) × 10−4 equation to an $11\, \sigma$ detection of polarization. We test if the polarization fraction depends on the clumps’ physical properties, including flux density, luminosity, Galactic latitude, and physical distance. We see a trend towards increasing polarization fraction with increasing Galactic latitude, but find no evidence that polarization depends on the other tested properties. The Simons Observatory, with angular resolution of order 1 arcmin and noise levels between 22 and $54\, \mu$K−arcmin at high frequencies, will substantially enhance our ability to determine the magnetic field structure in Galactic cold clumps. At $\ge 5\, \sigma$ significance, we predict the Simons Observatory will detect at least ∼12 000 cold clumps in intensity and ∼430 cold clumps in polarization. This number of polarization detections would represent a two orders of magnitude increase over the current Planck results. We also release software that can be used to mask these Galactic cold clumps in other analyses.

LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is planned to orbit the Sun-Earth Lagrangian point L2, where it will map the cosmic microwave background (CMB) polarization over the entire sky for three years, with three telescopes in 15 frequency bands between 34 and 448 GHz, to achieve an unprecedented total sensitivity of 2.2$\mu$K-arcmin, with a typical angular resolution of 0.5$^\circ$ at 100 GHz. The primary scientific objective of LiteBIRD is to search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. We provide an overview of the LiteBIRD project, including scientific objectives, mission and system requirements, operation concept, spacecraft and payload module design, expected scientific outcomes, potential design extensions and synergies with other projects.

ABSTRACT We characterize the kinematic and magnetic properties of H i filaments located in a high Galactic latitude region (165° < α < 195° and 12° < δ < 24°). We extract three-dimensional filamentary structures using fil3d from the Galactic Arecibo L-Band Feed Array H i (GALFA-H i) survey 21-cm emission data. Our algorithm identifies coherent emission structures in neighbouring velocity channels. Based on the mean velocity, we identify a population of local and intermediate velocity cloud (IVC) filaments. We find the orientations of the local (but not the IVC) H i filaments are aligned with the magnetic field orientations inferred from Planck 353 GHz polarized dust emission. We analyse position–velocity diagrams of the velocity-coherent filaments, and find that only 15 per cent of filaments demonstrate significant major-axis velocity gradients with a median magnitude of 0.5 km s−1 pc−1, assuming a fiducial filament distance of 100 pc. We conclude that the typical diffuse H i filament does not exhibit a simple velocity gradient. The reported filament properties constrain future theoretical models of filament formation.

LiteBIRD is a planned JAXA-led 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, $\delta r$, down to $\delta 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 SED, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compare the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the HFT (High-Frequency Telescope) frequency range is shifted logarithmically towards higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measure 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 better ability to discriminate between models through higher $\chi^2$ sensitivity. (abridged)

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

We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (H i) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, H i is strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on primordial gravitational waves. We separate the H i emission into distinct velocity components and detect dust polarization correlated with the local Galactic H i but not with the H i associated with Magellanic Stream i. We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary H i morphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the Hi morphology template correlates in B modes at a ∼10%–65% level over the multipole range 20 < ℓ < 200 with the BICEP/Keck maps, which contain contributions from dust, CMB, and noise components. We measure the spectral index of the filamentary dust component spectral energy distribution to be β = 1.54 ± 0.13. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity H i. Finally, we explore the morphological parameter space in the H i-based filamentary model.

We present evidence for scale-independent misalignment of interstellar dust filaments and magnetic fields. We estimate the misalignment by comparing millimeter-wave dust-polarization measurements from Planck with filamentary structures identified in neutral-hydrogen (H i ) measurements from H i 4PI. We find that the misalignment angle displays a scale independence (harmonic coherence) for features larger than the H i 4PI beamwidth (16.′2). We additionally find a spatial coherence on angular scales of  ( 1 ° ) . We present several misalignment estimators formed from the auto- and cross-spectra of dust-polarization and H i -based maps, and we also introduce a map-space estimator. Applied to large regions of the high-Galactic-latitude sky, we find a global misalignment angle of ∼2°, which is robust to a variety of masking choices. By dividing the sky into small regions, we show that the misalignment angle correlates with the parity-violating TB cross-spectrum measured in the Planck dust maps. The misalignment paradigm also predicts a dust EB signal, which is of relevance in the search for cosmic birefringence but as yet undetected; the measurements of EB are noisier than those of TB , and our correlations of EB with misalignment angle are found to be weaker and less robust to masking choices. We also introduce an H i -based dust-polarization template constructed from the Hessian matrix of the H i intensity, which is found to correlate more strongly than previous templates with Planck dust B modes.

Recent observations from the MUSTANG2 instrument on the Green Bank Telescope have revealed evidence of enhanced long-wavelength emission in the dust spectral energy distribution (SED) in the Orion Molecular Cloud (OMC) 2/3 filament on 25″ (0.1 pc) scales. Here we present a measurement of the SED on larger spatial scales (map size 0.°5–3° or 3–20 pc), at somewhat lower resolution (120″, corresponding to 0.25 pc at 400 pc) using data from the Herschel satellite and Atacama Cosmology Telescope (ACT). We then extend the 120″-scale investigation to other regions covered in the Herschel Gould Belt Survey (HGBS), specifically the dense filaments in the southerly regions of Orion A, Orion B, and Serpens-S. Our data set in aggregate covers approximately 10 deg 2 , with continuum photometry spanning from 160 μ m to 3 mm. These OMC 2/3 data display excess emission at 3 mm, though less (10.9% excess) than what is seen at higher resolution. Strikingly, we find that the enhancement is present even more strongly in the other filaments we targeted, with an average excess of 42.4% and 30/46 slices showing an inconsistency with the modified blackbody to at least 4 σ . Applying this analysis to the other targeted regions, we lay the groundwork for future high-resolution analyses. Additionally, we also consider a two-component dust model motivated by Planck results and an amorphous grain dust model. While both of these have been proposed to explain deviations in emission from a generic modified blackbody, we find that they have significant drawbacks, requiring many spectral points or lacking experimental data coverage.

Feedhorn- and orthomode transducer- (OMT) coupled transition edge sensor (TES) bolometers have been designed and micro-fabricated to meet the optical specifications of the LiteBIRD high frequency telescope (HFT) focal plane. We discuss the design and optical characterization of two LiteBIRD HFT detector types: dual-polarization, dual-frequency-band pixels with 195/280 GHz and 235/337 GHz band centers. Results show well-matched passbands between orthogonal polarization channels and frequency centers within 3% of the design values. The optical efficiency of each frequency channel is conservatively reported to be within the range 0.64 $$-$$ 0.72, determined from the response to a cryogenic, temperature-controlled thermal source. These values are in good agreement with expectations and either exceed or are within 10% of the values used in the LiteBIRD sensitivity forecast. Lastly, we report a measurement of loss in Nb/SiN $$_x$$ /Nb microstrip at 100 mK and over the frequency range 200–350 GHz, which is comparable to values previously reported in the literature.

Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments.

Recent measurements of Galactic polarized dust emission have found a nonzero TB signal, a correlation between the total intensity and the B-mode polarization component. We present evidence that this parity-odd signal is driven by the relative geometry of the magnetic field and the filamentary interstellar medium in projection. Using neutral hydrogen morphology and Planck polarization data, we find that the angle between intensity structures and the plane-of-sky magnetic field orientation is predictive of the signs of Galactic TB and EB. Our results suggest that magnetically misaligned filamentary dust structures introduce nonzero TB and EB correlations in the dust polarization, and that the intrinsic dust EB can be predicted from measurements of dust TB and TE over the same sky mask. We predict correlations between TE, TB, EB, and EE/BB, and confirm our predictions using synthetic dust polarization maps from magnetohydrodynamic simulations. We introduce and measure a scale-dependent effective magnetic misalignment angle, ψℓdust∼5° for 100 ≲ ℓ ≲ 500, and predict a positive intrinsic dust EB with amplitude DℓEB≲2.5μKCMB2 for the same multipole range at 353 GHz over our sky mask. Both the sign and amplitude of the Galactic EB signal can change with the sky area considered. Our results imply that searches for parity violation in the cosmic microwave background must account for the nonzero Galactic EB and TB signals, necessitating revision of existing analyses of the evidence for cosmic birefringence.

If a single line of sight (LOS) intercepts multiple dust clouds with different spectral energy distributions and magnetic field orientations, then the frequency scaling of each of the Stokes Q and U parameters of the thermal dust emission may be different, a phenomenon we refer to as LOS frequency decorrelation. We present first evidence for LOS frequency decorrelation in Planck data using independent measurements of neutral-hydrogen (H I ) emission to probe the 3D structure of the magnetized interstellar medium (ISM). We use H I -based measurements of the number of clouds per LOS and the magnetic field orientation in each cloud to select two sets of sightlines: (i) a target sample of pixels that are likely to exhibit LOS frequency decorrelation and (ii) a control sample of pixels that lack complex LOS structure. We test the null hypothesis that LOS frequency decorrelation is not detectable in Planck 353 and 217 GHz polarization data at high Galactic latitudes. We reject the null hypothesis at high significance based on data that show that the combined effect of polarization angle variation with frequency and depolarization are detected in the target sample. This detection is robust against the choice of cosmic microwave background (CMB) map and map-making pipeline. The observed change in polarization angle due to LOS frequency decorrelation is detectable above the Planck noise level. The probability that the detected effect is due to noise alone ranges from 5 × 10 −2 to 4 × 10 −7 , depending on the CMB subtraction algorithm and treatment of residual systematic errors; correcting for residual systematic errors consistently increases the significance of the effect. Within the target sample, the LOS decorrelation effect is stronger for sightlines with more misaligned magnetic fields, as expected. With our sample, we estimate that an intrinsic variation of ~15% in the ratio of 353 to 217 GHz polarized emission between clouds is sufficient to reproduce the measured effect. Our finding underlines the importance of ongoing studies to map the three-dimensional structure of the magnetized and dusty ISM that could ultimately help component separation methods to account for frequency decorrelation effects in CMB polarization studies.

The spatial distribution of neutral hydrogen (H i ) emission is a powerful probe of interstellar medium physics. The small-scale structure in H i channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H i observations and 50,985 quasar spectra from the Sloan Digital Sky Survey. We measure the equivalent widths of interstellar Na i D 1 and Na i D 2 absorption, and robustly conclude that together they depend nearly four times as strongly on the column density of small-scale structure in H i than on either the large-scale H i structure or the total H i column. This is inconsistent with the hypothesis that small-scale channel map structure is driven by velocity crowding. Instead, the data favor the interpretation that this emission structure predominantly originates in cold, dense interstellar material, consistent with a clumpy cold neutral medium.