Mapping the Magnetic Interstellar Medium in Three Dimensions over the Full Sky with Neutral Hydrogen

Recent analyses of 21 cm neutral hydrogen (H I) emission have demonstrated that H I gas is organized into linear filamentary structures that are preferentially aligned with the local magnetic field, and that the coherence of these structures in velocity space traces line-of-sight magnetic field tangling. On this basis, we introduce a paradigm for modeling the properties of magnetized, dusty regions of the interstellar medium (ISM), using the orientation of H I structure at different velocities to map "magnetically coherent" regions of space. We construct three-dimensional (position-position-velocity) Stokes parameter maps using H I4PI full-sky spectroscopic H I data. We compare these maps, integrated over the velocity dimension, to Planck maps of the polarized dust emission at 353 GHz. Without any free parameters governing the relation between H I intensity and dust emission, we find that our Q and U maps are highly correlated (r > 0.75) with the 353 GHz Q and U maps of polarized dust emission observed by Planck and reproduce many of its large-scale features. The E/B ratio of the dust emission maps agrees well with the H I-derived maps at large angular scales (l less than or similar to 120), supporting the interpretation that this asymmetry arises from the coupling of linear density structures to the Galactic magnetic field. We demonstrate that our 3D Stokes parameter maps constrain the 3D structure of the Galactic ISM and the orientation of the interstellar magnetic field.