The neutral interstellar medium (ISM) of the Milky Way is magnetized, but the only direct measurements of magnetic field strength — Zeeman splitting of the 21-cm hydrogen line — rest on a tiny sample: the Arecibo Millennium Survey measured the line-of-sight field toward only about 42 radio sources, far too few to establish the field-strength distribution in the cold neutral medium (CNM) or its variation across Galactic structures. An Astro2020 white paper by Clark, Heiles & Robishaw (2019) laid out the follow-up program: (1) a statistically meaningful, unbiased absorption-line Zeeman survey over hundreds of sightlines to determine the distribution and fluctuations of within and between structures such as the Orion/Eridanus superbubble; (2) densely sampled emission-line Zeeman mapping of shocked superbubble walls to test the reported factor-of-two field enhancement and an apparent field reversal; and (3) confirmation or refutation — and physical explanation — of the authors' own serendipitous detection of linear polarization of the 21-cm line at the 0.14%–0.35% level toward 13 of 18 Millennium sources, unexpected for a collisionally populated hyperfine transition. All three threads are open as of August 2026: the Millennium Survey remains the largest absorption Zeeman dataset, only three new sightlines from the FAST (Five-hundred-meter Aperture Spherical radio Telescope) have been added, emission-Zeeman work has been dormant since the mid-1990s, and no published confirmation or refutation of the linear-polarization claim exists.
Background
The interstellar medium (ISM) — the gas between stars — is threaded by a magnetic field of a few microgauss that competes with gravity, turbulence, and cosmic rays in governing gas dynamics and star formation. The diffuse atomic phase is conventionally divided into the warm neutral medium (WNM, roughly 5000–10000 K) and the cold neutral medium (CNM, roughly 20–200 K), both traced by the 21-cm hyperfine line of atomic hydrogen (HI). Magnetic fields are usually mapped indirectly: the linear polarization of starlight and of far-infrared dust emission, produced by elongated dust grains aligned with the field, reveals only the plane-of-sky field orientation , while Zeeman splitting is the only direct probe of field strength. In the Zeeman effect, the two circularly polarized components of a spectral line are split in frequency by an amount proportional to the line-of-sight field ; for the 21-cm line the splitting is a few Hz per microgauss, orders of magnitude smaller than the km/s-scale line widths, so the signal must be recovered statistically from Stokes V (circular polarization) spectra — a measurement limited by both sensitivity and instrumental systematics.
The existing census of in the CNM rests almost entirely on the Arecibo "Millennium" survey of 21-cm absorption against compact background radio sources (Heiles & Troland 2003, 2005), which observed only about 42 sightlines — roughly 7.5% of the suitably bright (above 1 Jy) sources accessible to Arecibo — so that major Galactic structures (the Perseus region, the Orion/Eridanus superbubble, the North Polar Spur superbubble, and the Radio Loop 3 superbubble) are each sampled by only a handful of lines of sight. In emission, where one may target arbitrary sky positions, Zeeman mapping by Heiles (1989) sampled 73 positions in shocked superbubble walls and found field strengths about twice those of randomly selected sightlines, together with an apparent large-scale field reversal across the Orion/Eridanus superbubble — but with sampling too coarse for any statistically firm statement. Separately, high-dynamic-range 21-cm surveys revealed that the CNM is threaded by long, thin "HI fibers" aligned with the field orientation traced by starlight polarization and polarized dust emission, so that HI morphology itself encodes field-orientation information; combining such orientation data with distances (e.g., from Gaia and dust-reddening maps) opens a path to three-dimensional magnetic maps — a path blocked precisely where field strength, not orientation, is required, and one that also matters for modeling polarized Galactic foregrounds in cosmic-microwave-background experiments.
Against this background, the Astro2020 decadal-survey white paper of Clark, Heiles & Robishaw (2019) reported a serendipitous discovery made while reanalyzing the Millennium data for calibration purposes: 13 of 18 test sources expected to be unpolarized show detectable linear polarization (in Stokes Q and U) at levels of 0.14%–0.35%. This was unexpected because the 21-cm hyperfine transition is commonly idealized as a two-level system populated by collisions, which should produce no linear polarization. The authors noted that Lyman-alpha photons are also an important populating agent for the 21-cm levels, and that anisotropic radiation could differentially populate the magnetic sublevels — a mechanism related to the "atomic alignment" polarization predicted by Yan & Lazarian (2007) — and speculated about a connection to the anisotropic HI fibers. The white paper therefore proposed an integrated program: an unbiased absorption-Zeeman survey of the CNM over hundreds of sightlines; densely sampled emission-Zeeman mapping of shocked superbubble walls; confirmation or refutation of the linear-polarization claim; the single-dish instrumentation these require; and integration with dust and starlight polarization into three-dimensional maps of magnetic field strength and orientation.
Problem Statement
Determine the strength, sign, and spatial distribution of the magnetic field in the diffuse neutral interstellar medium by analysing the existing published measurements and specifying the observational program set out in Clark, Heiles & Robishaw (2019), and resolving the unexplained linear polarization of the 21-cm line that paper reported. The program comprises three requested results, all of which are required:
-
(Census) A statistically meaningful, unbiased determination of in the CNM from the existing 21-cm absorption-Zeeman sample, together with the design of the unbiased absorption-line Zeeman survey that would extend it: hundreds of sightlines (the white paper notes 566 suitably bright sources in the Arecibo sky alone) selected without regard to field outcome, with the sensitivity and calibration required, to determine the distribution and fluctuations of magnetic-field strength within and between major Galactic structures (the Perseus region, the Orion/Eridanus superbubble, the North Polar Spur superbubble, and Radio Loop 3) and replacing the roughly 42-sightline Millennium sample on which all current statistics rest.
-
(Structure-linked strength) An assessment, from the Heiles (1989) emission measurements and their sampling limits, of (a) the reported factor-of-two field-strength enhancement of compressed wall and filament gas relative to randomly selected sightlines and (b) whether the apparent large-scale field reversal across the Orion/Eridanus superbubble is real and coherent, together with the design of the densely sampled Zeeman mapping of 21-cm emission in shocked superbubble walls and filaments — dense enough relative to the HI structural detail, unlike the 73-position sampling of Heiles (1989) — that would verify or refute them.
-
(The polarization anomaly) An adjudication of the reported linear polarization of the 21-cm line at the 0.14%–0.35% level seen toward 13 of 18 Millennium sources, based on the published Stokes Q and U measurements, together with the design of the observational follow-up — well-calibrated full-Stokes (I, Q, U, V) 21-cm measurements that include the off-source emission spectra needed to separate instrumental from astrophysical polarization — that would confirm or refute it; and, if the signal is astrophysical, identification of its physical mechanism accounting for the observed polarization fraction and its spectral behavior.
Three-dimensional maps of magnetic-field strength and orientation in the local diffuse ISM, combining these results with polarized dust emission and starlight polarization (including modern high-accuracy stellar polarimetry and distance information), are the motivating context of the program: as the white paper notes, HI structure carries information about the local field orientation, and a statistical sample of Zeeman splitting observations would supply the crucial missing component, the field strength. The development of sensitive single-dish full-Stokes instrumentation, polarization calibration, and weak-signal extraction algorithms is the enabling means toward these measurements. The three results are linked: the census supplies the field-strength statistics, the emission mapping ties strength to identified structures, and resolution of the polarization anomaly either removes a systematic concern for all 21-cm polarimetry or establishes a new magnetic diagnostic.
References
- Clark, S. E., Heiles, C., & Robishaw, T. (2019). Astro2020 Science White Paper: Magnetic Fields and Polarization in the Diffuse Interstellar Medium. arXiv:1903.07671, DOI: 10.48550/arXiv.1903.07671. https://arxiv.org/abs/1903.07671
- Heiles, C., & Troland, T. H. (2003). The Millennium Arecibo 21 Centimeter Absorption-Line Survey. I. Techniques and Gaussian Fits. The Astrophysical Journal Supplement Series 145, 329. DOI: 10.1086/367785. https://doi.org/10.1086/367785
- Heiles, C., & Troland, T. H. (2005). The Millennium Arecibo 21 Centimeter Absorption-Line Survey. IV. Statistics of Magnetic Field, Column Density, and Turbulence. The Astrophysical Journal 624, 773. DOI: 10.1086/428896. https://doi.org/10.1086/428896
- Heiles, C. (1989). Magnetic fields, pressures, and thermally unstable gas in prominent H I shells. The Astrophysical Journal 336, 808. DOI: 10.1086/167051. https://doi.org/10.1086/167051
- Yan, H., & Lazarian, A. (2007). Polarization from Aligned Atoms as a Diagnostic of Circumstellar, Active Galactic Nuclei, and Interstellar Magnetic Fields. II. Atoms with Hyperfine Structure. The Astrophysical Journal 657, 618. DOI: 10.1086/510847. https://doi.org/10.1086/510847
- Nowotka, M., Clark, S. E., Burkhart, B., Fissel, L., Ching, T.-C., Robishaw, T., & Heiles, C. (2025). Joint Analysis of HI Absorption Zeeman Measurements and the Morphology of Filamentary HI Emission. arXiv:2508.20065, DOI: 10.48550/arXiv.2508.20065. https://arxiv.org/abs/2508.20065
- Robishaw, T., Nowotka, M., Ching, T.-C., Green, J. A., Richards, A. M. S., Etoka, S., Gray, M., Clark, S. E., Bourke, T. L., & Fish, V. (2026). Measuring Magnetic Fields Near and Far with the SKA via the Zeeman Effect. arXiv:2606.24810, DOI: 10.48550/arXiv.2606.24810. https://arxiv.org/abs/2606.24810
- Chen, X., Ching, T.-C., Li, D., Heiles, C., Robishaw, T., Du, X., Krco, M., Jiang, P., Yang, Q., & Lu, J. (2025). Beam Measurements of Full Stokes Parameters for the FAST L-band 19-beam Receiver. The Astronomical Journal 169, 158. DOI: 10.3847/1538-3881/ada8aa. https://doi.org/10.3847/1538-3881/ada8aa
- Ching, T.-C., Li, D., Heiles, C., Li, Z.-Y., Qian, L., Yue, Y. L., Tang, J., & Jiao, S. H. (2022). An early transition to magnetic supercriticality in star formation. Nature 601, 49. DOI: 10.1038/s41586-021-04159-x. https://doi.org/10.1038/s41586-021-04159-x
- Islam, M. N., Robishaw, T., Veidt, B., Brar, P., & Du, X. (2024). Finite element model verification of the John A. Galt 26-m telescope through photogrammetry. Proceedings of SPIE 13094 (Ground-based and Airborne Telescopes X), 130943G. DOI: 10.1117/12.3020829. https://doi.org/10.1117/12.3020829
- Clark, S. E., & Hensley, B. S. (2019). Mapping the Magnetic Interstellar Medium in Three Dimensions over the Full Sky with Neutral Hydrogen. The Astrophysical Journal 887, 136. DOI: 10.3847/1538-4357/ab5803. https://doi.org/10.3847/1538-4357/ab5803
- Murray, C. E., Stanimirović, S., Heiles, C., Dickey, J. M., McClure-Griffiths, N. M., Lee, M.-Y., Goss, W. M., & Killerby-Smith, N. (2021). The MACH HI Absorption Survey. I. Physical Conditions of Cold Atomic Gas outside of the Galactic Plane. The Astrophysical Journal Supplement Series 256(2), 37. DOI: 10.3847/1538-4365/ac0f0b. https://doi.org/10.3847/1538-4365/ac0f0b
- Roshi, D. A., Aponte, N., Araya, E., Arce, H., Baker, L. A., Baan, W., Becker, T. M., Breakall, J. K., Brown, R. G., Brum, C. G. M., Busch, M., Campbell, D. B., Cohen, T., Cordova, F., Deneva, J. S., Devogele, M., Dolch, T., Fernandez-Rodriguez, F. O., Ghosh, T., Goldsmith, P. F., Gurvits, L. I., Haynes, M., Heiles, C., Hessel, J. W. T., Hickson, D., Isham, B., Kerr, R. B., Kelly, J., Kiriazes, J. J., Lautenbach, J., Lebron, M., Lewandowska, N., Magnani, L., Manoharan, P. K., Margot, J. L., Marshall, S. E., McGilvray, A. K., Mendez, A., Minchin, R., Negron, V., Nolan, M. C., Olmi, L., Paganelli, F., Palliyaguru, N. T., Pantoja, C. A., Paragi, Z., Parshley, S. C., Peek, J. E. G., Perera, B. B. P., Perillat, P., Pinilla-Alonso, N., Quintero, L., Radovan, H., Raizada, S., Robishaw, T., Route, M., Salter, C. J., Santoni, A., Santos, P., Sau, S., Selvaraj, D., Smith, A. J., Sulzer, M., Vaddi, S., Vargas, F., Venditti, F. C. F., Venkataraman, A., Verkouter, H., Virkki, A. K., Vishwas, A., Weinreb, S., Werthimer, D., Wolszczan, A., & Zambrano-Marin, L. F. (2021). The Future Of The Arecibo Observatory: The Next Generation Arecibo Telescope. arXiv:2103.01367, DOI: 10.48550/arXiv.2103.01367. https://arxiv.org/abs/2103.01367
- Shane, B., Burkhart, B., Fissel, L., Clark, S. E., Mocz, P., & Foley, M. M. (2024). Tracing 3-D Magnetic Field Structure Using Dust Polarization and the Zeeman Effect. arXiv:2411.10286, DOI: 10.48550/arXiv.2411.10286. https://arxiv.org/abs/2411.10286
- Booth, R. A., Ordog, A., Brown, J.-A., Landecker, T. L., Hill, A. S., West, J. L., Lei, M., Clark, S. E., Bracco, A., Dickey, J. M., & Carretti, E. (2025). A three-dimensional model for the reversal in the local large-scale interstellar magnetic field. arXiv:2512.03332, DOI: 10.48550/arXiv.2512.03332. https://arxiv.org/abs/2512.03332
- Li, B., Tan, J., & Mao, Y. (2021). Linear Polarization of the 21 cm Line from the Epoch of Reionization. The Astrophysical Journal 918, 14. DOI: 10.3847/1538-4357/ac09e4. https://doi.org/10.3847/1538-4357/ac09e4