The primary source of systematic uncertainty in the quest for the B-mode polarization of the Cosmic Microwave Background (CMB) introduced by primordial gravitational waves is polarized thermal emission from Galactic dust. Therefore, accurate characterization and separation of the polarized thermal dust emission is an essential step in distinguishing such a faint CMB B-mode signal. We provide a modelling framework to simulate polarized thermal dust emission based on the model described in Ghosh et al. (2017), making use of both the Planck dust and Effelsberg-Bonn HI surveys over the northern Galactic cap. Our seven-parameter dust model, incorporating both HI gas in three different column density templates as a proxy for spatially variable dust intensity and a phenomenological model of Galactic magnetic field, is able to reproduce both 1- and 2-point statistics of the observed dust polarization maps seen by Planck at 353 GHz over a selected low-column density region in the northern Galactic cap. This work has important applications in assessing the accuracy of component separation methods and in quantifying the confidence level of separating polarized Galactic emission and the CMB B-mode signal, as is needed for ongoing and future CMB missions.

Context. There are significant amounts of H2 in the Milky Way. Due to its symmetry H2 does not radiate at radio frequencies. CO is thought to be a tracer for H2, however CO is formed at significantly higher opacities than H2. Thus, toward high Galactic latitudes significant amounts of H2 are hidden and called CO-dark. Aims. We demonstrate that the dust-to-gas ratio is a tool to identify locations and column densities of CO-dark H2. Methods. We adopt the hypothesis of a constant E(B-V)/NH ratio, independent of phase transitions from HI to H2. We investigate the Doppler temperatures TD, from a Gaussian decomposition of HI4PI data, to study temperature dependencies of E(B-V)/NHI. Results. The E(B-V)/NHI ratio in the cold HI gas phase is high in comparison to the warmer one. We consider this as evidence that cold HI gas toward high Galactic latitudes is associated with H2. Beyond CO-bright regions we find for TD < 1165 K a correlation (NHI + 2NH2 )/NHI prop -log T_D. In combination with a factor XCO = 4.0 10 20 cm^-2 (K km s^-1 )-1 this yields for the full-sky NH /E(B-V) sim 5.1 to 6.7 10^21 cm^-2 mag^-1, compatible with X-ray scattering and UV absorption line observations. Conclusions. Cold HI with T_D < 1165 K contains on average 46% CO-dark H2. Prominent filaments have TD < 220 K and typical excitation temperatures Tex sim 50 K. With a molecular gas fraction of > 61% they are dominated dynamically by H2.

Context. The interstellar medium is affected by turbulence and observed H i structures in channel maps are shaped by turbulent motions. It is taken for granted by a few theoreticians that observed H i structures do not represent real density enhancement but velocity caustics, caused by velocity crowding. This interpretation was questioned and objections by Clark et al. led to violent debates. Aims. To settle the discussion we verify theoretical key parameters by using Effelsberg Bonn H i Survey (EBHIS) observations. Methods. We apply unsharp masking to determine filamentary H i structures at high spatial frequencies. In addition we use Gaussian parameters to distinguish the cold neutral medium (CNM) from observed H i column densities. We compare power spectra and spatial distributions of dust and H i column densities, distinguishing CNM and multiphase column densities at various velocity widths. Results. Observations contradict the Velocity Channel Analysis (VCA) postulate that the spectral index should steepen with the width of the velocity window. We rather find that the thin slice spectral index depends strongly on the H i phase composition. Multiphase power spectra are steeper for regions with cold gas. VCA denies such H i phase dependencies on the power distribution. Separating the CNM we find that the power spectra are significantly flatter than those for the multiphase H i composite. We observe excess CNM power for small scale structures originating from cold dust bearing filaments that are embedded in the CNM. Spectral indices for narrow channel widths depend on the Doppler temperature of the H i gas. In presence of enhanced small scale H i structure the far infrared emission from dust is also enhanced. Conclusions. Small scale cold filamentary H i structures are predominantly caused by density enhancements due to phase transitions rather than by velocity caustics.

Context: The interstellar medium (ISM) on all scales is full of structures that can be used as tracers of processes that feed turbulence. Aims: We used HI survey data to derive global properties of the angular power distribution of the local ISM. Methods: HI4PI observations on an nside = 1024 HEALPix grid and Gaussian components representing three phases, the cold, warm, and unstable lukewarm neutral medium (CNM, WNM, and LNM), were used for velocities $|v_{\mathrm{LSR}}| \leq 25$ \kms. For high latitudes $|b| > 20\deg$ we generated apodized maps. After beam deconvolution we fitted angular power spectra. Results: Power spectra for observed column densities are exceptionally well defined and straight in log-log presentation with 3D power law indices $\gamma \geq -3$ for the local gas. For intermediate velocity clouds (IVCs) we derive $\gamma = -2.6$ and for high velocity clouds (HVCs) $\gamma = -2.0$. Single-phase power distributions for the CNM, LNM, and WNM are highly correlated and shallow with $ \gamma \sim -2.5$ for multipoles $l \leq 100$. Excess power from cold filamentary structures is observed at larger multipoles. The steepest single-channel power spectra for the CNM are found at velocities with large CNM and low WNM phase fractions. Conclusions: The phase space distribution in the local ISM is configured by phase transitions and needs to be described with three distinct different phases, being highly correlated but having distributions with different properties. Phase transitions cause locally hierarchical structures in phase space. The CNM is structured on small scales and is restricted in position-velocity space. The LNM as an interface to the WNM envelops the CNM. It extends to larger scales than the CNM and covers a wider range of velocities. Correlations between the phases are self-similar in velocity.

Context. A large fraction of the interstellar medium can be characterized as a multiphase medium. The neutral hydrogen gas is bistable with a cold and warm neutral medium (CNM and WNM respectively) but there is evidence for an additional phase at intermediate temperatures, a lukewarm neutral medium (LNM) that is thermally unstable. Aims. We use all sky data from the HI4PI survey to separate these neutral HI phases with the aim to determine their distribution and phase fractions in the local interstellar medium. Methods. HI4PI observations, gridded on an nside = 1024 HEALPix grid, were decomposed into Gaussian components. From the frequency distribution of the velocity dispersions we infer three separate linewidth regimes. Accordingly we extract the HI line emission corresponding to the CNM, LNM, and WNM. We generated all-sky maps of these phases in the local HI gas with -8 < v_LSR < 8 km/s. Results. Each of the HI phases shows distinct structures on all scales. The LNM never exists as a single phase but contributes on average 41% of the HI. The CNM is prominent only for 22% of the sky, contributes there on average 34% but locally up to 60% of the HI and is associated with dust at temperatures T_dust ~ 18.6 K. Embedded cold filaments show a clear anti-correlation between CNM and LNM. Also the smoothly distributed WNM is anti-correlated with the CNM. It contributes for the rest of the sky 39% with dust associated at temperatures T_dust ~ 19.4 K. Conclusions. The CNM in filaments exists on small scales. Here the observed anti-correlation between LNM and CNM implies that both, filaments and the surrounding more extended LNM, must have a common origin.

LOFAR detected toward 3C 196 linear polarization structures which were found subsequently to be closely correlated with cold filamentary HI structures. The derived direction-dependent HI power spectra revealed marked anisotropies for narrow ranges in velocity, sharing the orientation of the magnetic field as expected for magneto hydrodynamical turbulence. Using the Galactic portion of the Effelsberg-Bonn HI Survey we continue our study of such anisotropies in the HI distribution in direction of two WSRT fields, Horologium and Auriga; both are well known for their prominent radio-polarimetric depolarization canals. At 349 MHz the observed pattern in total intensity is insignificant but polarized intensity and polarization angle show prominent ubiquitous structures with so far unknown origin. Apodizing the HI survey data by applying a rotational symmetric 50 percent Tukey window, we derive average and position angle dependent power spectra. We fit power laws and characterize anisotropies in the power distribution. We use a Gaussian analysis to determine relative abundances for the cold and warm neutral medium. For the analyzed radio-polarimetric targets significant anisotropies are detected in the HI power spectra; their position angles are aligned to the prominent depolarization canals, initially detected by WSRT. HI anisotropies are associated with steep power spectra. Steep power spectra, associated with cold gas, are detected also in other fields. Radio-polarimetric depolarization canals are associated with filamentary HI structures that belong to the cold neutral medium (CNM). Anisotropies in the CNM are in this case linked to a steepening of the power-spectrum spectral index, indicating that phase transitions in a turbulent medium occur on all scales. Filamentary HI structures, driven by thermal instabilities, and radio-polarimetric filaments are associated with each other.

The Planck survey has quantified polarized Galactic foregrounds and established that they are a main limiting factor in the quest for the cosmic microwave background (CMB) B-mode signal induced by primordial gravitational waves during cosmic inflation. Accurate separation of the Galactic foregrounds therefore binds this quest to our understanding of the magnetized interstellar medium (ISM). The two most relevant empirical results from analysis of Planck data are line of sight depolarization arising from fluctuations of the Galactic magnetic field orientation and alignment of filamentary dust structures with the magnetic field at high Galactic latitude. Furthermore, Planck and HI emission data in combination indicate that most of the filamentary dust structures are in the cold neutral medium. The goal of this paper is to test whether these salient observational results, taken together, can account fully for the statistical properties of the dust polarization over a selected low column density region comprising 34 % of the southern Galactic cap ($b \le -30\deg$). To do this, we construct a dust model that incorporates HI column density maps as tracers of the dust intensity structures and a phenomenological description of the Galactic magnetic field. By adjusting the parameters of the dust model, we were able to reproduce the Planck dust observations at 353 GHz in the selected region. Realistic simulations of the polarized dust emission enabled by such a dust model are useful for testing the accuracy of component separation methods, studying non-Gaussianity, and constraining the amount of decorrelation with frequency.

Measurement of the Galactic neutral atomic hydrogen (HI) column density, NHI, and brightness temperatures, Tb, is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy HI datasets has been the Leiden/Argentine/Bonn Survey (LAB). We release the HI 4$\pi$ survey (HI4PI), an all-sky database of Galactic HI, which supersedes the LAB survey. The HI4PI survey is based on data from the recently completed first coverage of the Effelsberg-Bonn HI Survey (EBHIS) and from the third revision of the Galactic All-Sky Survey (GASS). EBHIS and GASS share similar angular resolution and match well in sensitivity. Combined, they are ideally suited to be a successor to LAB. The new HI4PI survey outperforms the LAB in angular resolution (16.2', FWHM) and sensitivity (RMS: 43 mK). Moreover, it has full spatial sampling and thus overcomes a major drawback of LAB, which severely undersamples the sky. We publish all-sky column density maps of the neutral atomic hydrogen in the Milky Way, along with full spectroscopic data, in several map projections including HEALPix.

We investigate data from the Galactic Effelsberg--Bonn HI Survey (EBHIS), supplemented with data from the third release of the Galactic All Sky Survey (GASS III) observed at Parkes. We explore the all sky distribution of the local Galactic HI gas with $|v_{\rm LSR}| < 25 $ kms$^{-1}$ on angular scales of 11' to 16'. Unsharp masking (USM) is applied to extract small scale features. We find cold filaments that are aligned with polarized dust emission and conclude that the cold neutral medium (CNM) is mostly organized in sheets that are, because of projection effects, observed as filaments. These filaments are associated with dust ridges, aligned with the magnetic field measured on the structures by Planck at 353 GHz. The CNM above latitudes $|b|>20^\circ$ is described by a log-normal distribution, with a median Doppler temperature $T_{\rm D} = 223$ K, derived from observed line widths that include turbulent contributions. The median neutral hydrogen (HI) column density is $N_{\rm HI} \simeq 10^{19.1}\,{\rm cm^{-2}}$. These CNM structures are embedded within a warm neutral medium (WNM) with $N_{\rm HI} \simeq 10^{20} {\rm cm^{-2}}$. Assuming an average distance of 100 pc, we derive for the CNM sheets a thickness of $< 0.3$ pc. Adopting a magnetic field strength of $B_{\rm tot} = (6.0 \pm 1.8)\mu$G, proposed by Heiles & Troland 2005, and assuming that the CNM filaments are confined by magnetic pressure, we estimate a thickness of 0.09 pc. Correspondingly the median volume density is in the range $ 14 < n < 47 {\rm cm^{-3}}$.

The Galactic All-Sky Survey is a survey of Galactic atomic hydrogen emission in the southern sky observed with the Parkes 64-m Radio Telescope. The first data release (GASS I) concerned survey goals and observing techniques, the second release (GASS II) focused on stray radiation and instrumental corrections. We seek to remove the remaining instrumental effects and present a third data release. We use the HEALPix tessellation concept to grid the data on the sphere. Individual telescope records are compared with averages on the nearest grid position for significant deviations. All averages are also decomposed into Gaussian components with the aim of segregating unacceptable solutions. Improved priors are used for an iterative baseline fitting and cleaning. In the last step we generate 3-D FITS data cubes and examine them for remaining problems. We have removed weak, but systematic baseline offsets with an improved baseline fitting algorithm. We have unraveled correlator failures that cause time dependent oscillations; errors cause stripes in the scanning direction. The remaining problems from radio frequency interference (RFI) are spotted. Classifying the severeness of instrumental errors for each individual telescope record (dump) allows us to exclude bad data from averages. We derive parameters that allow us to discard dumps without compromising the noise of the resulting data products too much. All steps are reiterated several times: in each case, we check the Gaussian parameters for remaining problems and inspect 3-D FITS data cubes visually. We find that in total ~1.5% of the telescope dumps need to be discarded in addition to ~0.5% of the spectral channels that were excluded in GASS II.The new data release facilitates data products with improved quality. A new web interface, compatible with the previous version, is available for download of GASS III FITS cubes and spectra.

The question, whether the stellar populations in the Milky Way take part in flaring of the scale heights as observed for the HI gas is a matter of debate. Standard mass models for the Milky Way assume a constant scale height for each of the different stellar distributions. However, there is mounting evidence that at least some of the stellar distributions reach at large galactocentric distances high altitudes that are incompatible with a constant scale height. We discuss recent observational evidence for stellar flaring and compare it with HI data from the Leiden/Argentine/Bonn (LAB) survey. Within the systemic and statistical uncertainties we find a good agreement between both.

Context. In the HI line profiles in the Leiden-Argentina-Bonn (LAB) all-sky database, we have found a population of very cold HI clouds. So far, the role of these clouds in the interstellar medium (ISM) has remained unclear. Aims. In this paper, we attempt to confirm the existence of the narrow-line HI emission (NHIE) clouds by using the data from the Parkes Galactic all-sky survey (GASS) and try to find their place among other coldest constituents of the ISM. Methods. We repeat the search of NHIE with the GASS data and derive or compile some preliminary estimates for the distribution, temperatures, distances, linear sizes, column and number densities, masses, and the composition of NHIE clouds, and compare these data with corresponding estimates for HI self-absorption (HISA) features, the Planck cold clumps (CC), and infrared dark clouds (IRDC). Results. We demonstrate that from LAB and GASS we can separate comparable NHIE complexes, and the properties of the obtained NHIE clouds are very similar to those of HISA features, but both of these types of clouds are somewhat warmer and more extended and have lower densities than the cores in the Planck CC and IRDC. Conclusions. We conclude that NHIE may be the same type of clouds as HISA, but in different observing conditions, in the same way as the Planck CC and IRDC are most likely similar ISM structures in different observing conditions and probably in slightly different evolutionary stages. Both NHIE and HISA may be an intermediate phase between the diffuse cold neutral medium and star-forming molecular clumps represented by the Planck CC and IRDC.

In the previous papers of this series, we have decomposed into Gaussian components all the HI 21-cm line profiles of the Leiden-Argentina-Bonn (LAB) database, and studied statistical distributions of the obtained Gaussians. Now we are interested in separation from the general database of the components the "clouds" of closely spaced similar Gaussians. In this paper we describe the new cloud-finding algorithm. To separate the clouds of similar Gaussians, we start with the single-link hierarchical clustering procedure in five-dimensional (longitude, latitude, velocity, Gaussian width and height) space, but make some modifications to accommodate it to the large number of components. We also use the requirement that each cloud may be represented at any observed sky position by only one Gaussian and take into account the similarity of global properties of the merging clouds. As a test, we apply the algorithm for finding the clouds of the narrowest HI 21-cm line components. Using the full sky search for cold clouds, we easily detect the coldest known HI clouds and demonstrate that actually they are a part of a long narrow ribbon of cold clouds. We model these clouds as a part of a planar gas ring, deduce their spatial placement, and discuss their relation to supernova shells in the solar neighborhood. We conclude that the proposed algorithm satisfactorily solves the posed task. We guess that the study of the narrowest HI 21-cm line components may be a useful tool for finding the structure of neutral gas in solar neighborhood.

Traditionally IVC and HVC were defined to be concentrations of HI gas, with line-of-sight velocities that are inconsistent with data on the differential rotation of the Galaxy. We demonstrate that IVCs and HVCs can be identified from certain density enhancements in (V_C, FWHM) distribution of Gaussians, representing the Galactic HI 21 cm radio lines. We study the Gaussians, which parameters fall into the regions of the phase-space density enhancements about $(V_C, FWHM) = (-131, 27), (164, 26) and (-49 km/s, 23 km/s). The sky distribution of the Gaussians, corresponding to the first two concentrations, very well represents the sky distribution of HVCs, as obtained on the basis of the traditional definition of these objects. The Gaussians of the last concentration correspond to IVCs. Based on this identification, the division line between IVCs and HVCs can be drawn at about |V_C| = 74 km/s, and IVCs can be identified down to velocities of about |V_C| = 24 km/s. Traces of both IVCs and HVCs can also be seen in the sky distribution of Gaussians with FWHM = 7.3 km/s. In HVCs, these cold cores have small angular dimensions and low observed brightness temperatures T_b. In IVCs, the cores are both larger and brighter. This definition of IVCs and HVCs is less dependant than the traditional one, on the differential rotation model of the Galaxy. The consideration of line-width information may enable IVCs and HVCs to be better distinguished from each other, and from the ordinary Galactic HI

Context. Gas within a galaxy is forced to establish pressure balance against gravitational forces. The shape of an unperturbed gaseous disk can be used to constrain dark matter models. Aims. We derive the 3-D HI volume density distribution for the Milky Way out to a galactocentric radius of 40 kpc and a height of 20 kpc to constrain the Galactic mass distribution. Methods. We used the Leiden/Argentine/Bonn all sky 21-cm line survey. The transformation from brightness temperatures to densities depends on the rotation curve. We explored several models, reflecting different dark matter distributions. Each of these models was set up to solve the combined Poisson-Boltzmann equation in a self-consistent way and optimized to reproduce the observed flaring. Results. Besides a massive extended halo of M ~ 1.8 10^{12} Msun, we find a self-gravitating dark matter disk with M=2 to 3 10^{11} Msun, including a dark matter ring at 13 < R < 18.5 kpc with M = 2.2 to 2.8 10^{10} Msun. The existence of the ring was previously postulated from EGRET data and coincides with a giant stellar structure that surrounds the Galaxy. The resulting Milky Way rotation curve is flat up to R~27 kpc and slowly decreases outwards. The \hi gas layer is strongly flaring. The HWHM scale height is 60 pc at R = 4 kpc and increases to ~2700$ pc at R=40 kpc. Spiral arms cause a noticeable imprint on the gravitational field, at least out to R = 30 kpc. Conclusions. Our mass model supports previous proposals that the giant stellar ring structure is due to a merging dwarf galaxy. The fact that the majority of the dark matter in the Milky Way for $R \la 40$ kpc can be successfully modeled by a self-gravitating isothermal disk raises the question of whether this massive disk may have been caused by similar merger events in the past.

Context: To investigate the properties of the 21-cm radio-lines of galactic neutral hydrogen, the profiles of "The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI" are decomposed into Gaussian components. Aims: The width distribution of the obtained components is analysed and compared with similar studies by other authors. Methods: The study is based on an automatic profile decomposition algorithm. As the Gaussians obtained for the complex HI profiles near the galactic plane cannot be directly interpreted in terms of the properties of gas clouds, we mainly study the selected simplest profiles in a limited velocity range. The selection criteria are described and their influence on the results discussed. Results: Considering only the simplest HI profiles, we demonstrate that for Gaussians with relatively small LSR velocities (-9 <= V_C <= 4 km/s) it is possible to distinguish three or four groups of preferred line-widths with mean widths FWHM = 3.9 +- 0.6, 11.8 +- 0.5, 24.1 +- 0.6, and 42 +- 5 km/s. Verschuur previously proposed similar line-width regimes, but with somewhat larger widths. He used a human-assisted decomposition for a nearly 50 times smaller database and we discuss systematic differences in analysis and results. Conclusions: The line-widths of about 3.9 and 24.1 km/s are well understood in the framework of traditional models of the two-phase interstellar medium. The components with the widths around 11.8 km/s indicate that a considerable fraction (up to about 40%) of the HI gas is thermally unstable. The reality and the origin of the broad lines with the widths of about 42 km/s is more obscure. These, however, contain only about 4% of the total observed column densities.

We have analyzed the Gaussian decomposition of the Leiden/Dwingeloo Survey (LDS) of galactic neutral hydrogen for the presence of Gaussians probably not directly related to galactic HI emission. It is demonstrated that at least three classes of such components can be distinguished. The narrowest Gaussians mostly represent stronger random noise peaks in profiles and some still uncorrected radio-interferences. Many of slightly wider weak Gaussians are caused by increased uncertainties near the profile edges and with the still increasing width the baseline problems become dominating among weak components. Statistical criteria are given for separation of the parameter space regions, most likely populated with the problematic components from those where the Gaussians are with higher probability describing the actual Milky Way HI emission. The same analysis is applied to the Leiden/Argentina/Bonn survey (LAB). It is demonstrated that the selection criteria for dividing the parameter space are to a great extent independent of the particular survey in use. The presence of the baseline problems in the LDS is indicated by the peculiarities of the distribution of the widest Gaussians in the sky. A similar plot for the northern part of the LAB demonstrates considerably lower numbers of spurious components, but there are still problems with the southern part of the LAB. The strange characteristics of the observational noise in the southern part of the LAB are pointed out.

We study the properties of all major HVC complexes from a sample compiled 1991 by Wakker & van Woerden (WvW). We use the Leiden/Argentine/Bonn all sky 21-cm line survey and decompose the profiles into Gaussian components. We find a well defined multi-component structure for most of the HVC complexes. The cold HVC phase has lines with typical velocity dispersions of sigma = 3 km/s and exists only within more extended broad line regions, typically with sigma = 12 km/s. The motions of the cores relative to the envelopes are characterized by Mach numbers M = 1.5. The center velocities of the cores within a HVC complex have typical dispersions of 20 km/s. Remarkable is the well defined two-component structure for some prominent HVC complexes in the outskirts of the Milky Way: Complex H, the Magellanic Stream and the Leading Arm. There might be some indications for an interaction between HVCs and disk gas at intermediate velocities. This is possible for complex H, M, C, WB, WD, WE, WC, R, G, GCP, and OA, but not for complex A, MS, ACVHV, EN, WA, and P. Conclusions: The line widths, determined by us, imply that estimates of HVC masses, as far as derived from the WvW database, need to be scaled up by a factor 1.4. Correspondingly, guesses for the external pressure of a confining coronal gas need to be revised upward by a factor of 2. The HVC multi-phase structure implies in general that currently the halo pressure is significantly underestimated. In consequence, the HVC multi-phase structure may indicate that most of the complexes are circum-galactic. HVCs have turbulent energy densities which are an order of magnitude larger than that of comparable clumps in the Galactic disk.

We study the hypothesis that high-velocity clouds (HVCs) may originate from instabilities within the gaseous phase of the Galactic halo. According to the hydrostatic equilibrium model of Kalberla & Kerp (1998), we determine the probability distribution of occurrence of instabilities within the Galactic halo. Instabilities may produce condensations within the Galactic halo beyond a z-distances of 4 kpc, which are accelerated by gravity and decelerated by drag-forces. We present evidence that the HI high-velocity dispersion component, discovered by Kalberla et al. (1998) is associated with high-velocity-halo gas. The physical properties of this high-velocity-halo gas are similar to the recently detected highly-ionized HVCs by Sembach et al. (1995, 1998). Accordingly, the HI high-velocity dispersion component may be the neutral tracer of the turbulent gas motions within the Galactic halo. Our calculations demonstrate, that the sky-averaged signal of in-falling condensations does not differ significantly from the HI emission of the turbulent Galactic halo gas.