CHANG-ES XXXIX. Magnetic field structure in edge-on galaxies: Stacking Stokes parameters
Pith reviewed 2026-05-08 19:04 UTC · model grok-4.3
The pith
Stacking Stokes Q and U cubes from 27 edge-on galaxies recovers faint polarized halo emission and shows an X-shaped pattern.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Stacking aligned, scaled, and convolved Stokes Q and U cubes from 27 star-forming late-type edge-on galaxies, followed by RM synthesis, yields maps with a clear X-shaped pattern in the polarization plane, polarized emission detected up to 9 kpc above the disk, stronger PI on the approaching side, and an approximately 60 percent decrease near the minor axis. Synthetic data tests establish that the method is valid when the underlying distributions of PI, χ0, and RM are tightly constrained, although it introduces a systematic RM uncertainty of 90 rad m^{-2} and underestimates the recovered polarized intensity. No global RM pattern is recovered from the stack.
What carries the argument
Stacking of aligned Stokes Q and U spectra across the sample followed by RM synthesis on the combined cubes.
If this is right
- Faint polarized emission in galactic halos can be recovered out to at least 9 kpc when data from multiple galaxies are stacked.
- The X-shaped polarization pattern appears consistently across the sample and matches earlier single-galaxy results.
- Polarized intensity is stronger on the approaching side of the galaxies in the stack.
- Polarized intensity drops by about 60 percent near the minor axis in the halo region.
- No coherent global rotation measure pattern emerges from the stacked data.
Where Pith is reading between the lines
- Broader frequency coverage in future observations could reduce the reported 90 rad m^{-2} systematic RM uncertainty and improve Faraday depth resolution.
- The stacking technique could be applied to larger or more distant samples if galaxies are pre-selected to ensure comparable underlying PI, χ0, and RM distributions.
- Combining stacked results with high-resolution maps of individual galaxies would help separate common halo features from galaxy-to-galaxy variations.
Load-bearing premise
The 27 galaxies share sufficiently similar and tightly constrained distributions of polarized intensity, intrinsic polarization angle, and rotation measure that stacking does not introduce major unaccounted biases.
What would settle it
A test in which individual galaxies with independently measured differing RM distributions are stacked and the output X-pattern or average RM deviates systematically by more than 90 rad m^{-2} from the expected range would falsify the validity of the stacking approach.
Figures
read the original abstract
Galactic magnetic fields regulate star formation and cosmic-ray (CR) transport, and understanding their three-dimensional structure, particularly in star-forming late-type galaxies, is key to constraining galactic CR transport. We explore the validity of stacking Stokes $Q$ and $U$ spectra, to infer about the intrinsic polarisation characteristics of star-forming galaxies. To prepare the stacking experiment, we align, scale, convolve, and reproject $C$-band (6 GHz) Stokes $Q$ and Stokes $U$ cubes of 27 star-forming late-type edge-on galaxies. On the stacked cubes, we perform RM-synthesis and discuss the derived polarised intensity (PI), polarisation angle ($\chi_0$), and RM maps. Synthetic data tests demonstrate that stacking Stokes $Q$ and $U$ spectra is valid for tightly constrained underlying distributions of PI, $\chi_0$, and RM. For underlying PI, $\chi_0$, and RM distributions that represent star-forming galaxies, stacking introduces a systematic uncertainty of $\delta_\mathrm{RM}^\mathrm{sys}=90 \mathrm{rad m^{-2}}$ and significantly underestimates the recovered PI. Stacking results reveal a clear X-shaped pattern in the polarisation plane, consistent with prior findings, detecting polarised emission up to 9 kpc above the galactic disc. We find stronger PI on the approaching side of galaxies. Furthermore, we find a decrease in PI in the galactic halo of $\sim 60$% near the galaxy's minor axis. A global RM pattern, as reported in a previous study, cannot be confirmed. Based on our analysis, we present stacking of Stokes $Q$ and Stokes $U$ cubes as an effective tool to recover faint polarised emission in the halo of nearby galaxies, if the underlying distributions of PI, $\chi_0$, and RM are tightly constrained. Our findings motivate future studies using broader-band data to increase the resolution in Faraday depth.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates stacking of aligned, scaled, convolved, and reprojected Stokes Q and U cubes from 27 star-forming late-type edge-on galaxies at C-band (6 GHz). RM-synthesis is applied to the stacked cubes to derive maps of polarised intensity (PI), intrinsic polarisation angle (χ₀), and rotation measure (RM). Synthetic tests are used to assess stacking validity, showing that the method recovers faint halo emission only for tightly constrained underlying distributions of PI, χ₀, and RM; for representative star-forming galaxy distributions, it underestimates PI and introduces δ_RM^sys = 90 rad m^{-2}. Observational results include an X-shaped polarisation pattern extending to 9 kpc above the disc, stronger PI on the approaching side, a ~60% PI decrease near the minor axis in the halo, and no confirmed global RM pattern. The authors conclude that Stokes Q/U stacking is an effective tool for faint halo emission if the distributions are tightly constrained.
Significance. If the validity condition is met by the sample, the work offers a practical method to enhance detection of faint polarised emission in galactic halos, revealing magnetic field structures relevant to cosmic-ray transport models. The synthetic tests provide a clear framework for quantifying biases (PI underestimation and systematic RM uncertainty), which is a methodological strength.
major comments (2)
- [Abstract and synthetic tests section] Abstract and synthetic tests section: The validity of the stacking method is explicitly conditional on 'tightly constrained' distributions of PI, χ₀, and RM. However, the manuscript does not report the measured standard deviations or variances of PI, χ₀, and RM from the 27 individual aligned cubes, nor does it quantitatively compare these to the parameter ranges of the 'tightly constrained' versus 'representative' synthetic cases that produce the noted PI underestimation and δ_RM^sys = 90 rad m^{-2}. This comparison is required to determine whether the real ensemble satisfies the condition under which the method is declared valid.
- [Results and discussion] Results and discussion: The non-confirmation of a global RM pattern (as reported in a prior study) must account for the impact of the 90 rad m^{-2} systematic uncertainty; without this, the comparison to previous work is incomplete and the interpretation of the stacked RM map is weakened.
minor comments (2)
- [Methods] Methods: Provide more detail on the precise alignment parameters, scaling factors, and selection criteria for the 27 galaxies to allow reproducibility and assessment of potential selection biases.
- [Figures] Figures: Label the vertical scale heights explicitly on the stacked maps and indicate the noise levels or detection thresholds used to claim emission out to 9 kpc.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important points regarding the justification of the stacking method and the interpretation of the RM results. We address each major comment below and will revise the manuscript to incorporate the suggested clarifications.
read point-by-point responses
-
Referee: The validity of the stacking method is explicitly conditional on 'tightly constrained' distributions of PI, χ₀, and RM. However, the manuscript does not report the measured standard deviations or variances of PI, χ₀, and RM from the 27 individual aligned cubes, nor does it quantitatively compare these to the parameter ranges of the 'tightly constrained' versus 'representative' synthetic cases that produce the noted PI underestimation and δ_RM^sys = 90 rad m^{-2}. This comparison is required to determine whether the real ensemble satisfies the condition under which the method is declared valid.
Authors: We agree that explicit reporting of the observed dispersions is necessary to substantiate the applicability of the stacking results. In the revised manuscript, we will add a new subsection (or table) in the synthetic tests section that tabulates the standard deviations of PI, χ₀, and RM measured directly from the 27 aligned individual cubes. These values will be compared quantitatively to the input ranges of both the 'tightly constrained' and 'representative' synthetic ensembles, allowing readers to assess whether the real data satisfy the validity criterion. revision: yes
-
Referee: The non-confirmation of a global RM pattern (as reported in a prior study) must account for the impact of the 90 rad m^{-2} systematic uncertainty; without this, the comparison to previous work is incomplete and the interpretation of the stacked RM map is weakened.
Authors: We concur that the systematic RM uncertainty must be folded into the comparison with earlier studies. In the revised discussion section, we will explicitly state that the 90 rad m^{-2} uncertainty could mask or dilute a global RM pattern of the amplitude reported previously. We will add a quantitative statement on the minimum detectable RM gradient given this uncertainty and will qualify the non-confirmation accordingly, while retaining the observational result that no coherent large-scale pattern is detected above the noise and uncertainty floor. revision: yes
Circularity Check
No significant circularity; validity claim rests on synthetic tests and data application rather than self-referential reduction.
full rationale
The paper's core procedure—aligning, stacking Stokes Q/U cubes, performing RM-synthesis, and interpreting results—relies on observational processing and separate synthetic validation tests rather than any equation or parameter that reduces to its own fitted inputs by construction. Self-citations to prior CHANG-ES work appear only for context on X-shaped patterns and are not load-bearing for the stacking validity argument. The conditional nature of the conclusion (effective only if distributions are tightly constrained) is stated explicitly but does not create a definitional loop; the paper does not derive the constraint from the stacked data itself. No fitted-input-called-prediction, ansatz smuggling, or renaming of known results occurs in the derivation chain.
Axiom & Free-Parameter Ledger
free parameters (1)
- δ_RM^sys =
90 rad m^{-2}
axioms (2)
- domain assumption The 27 galaxies have sufficiently similar distributions of PI, χ0, and RM for stacking to be valid.
- standard math Standard assumptions of Faraday rotation measure synthesis apply to the stacked Stokes Q and U data.
Reference graph
Works this paper leans on
-
[1]
Arshakian, T. G. & Beck, R. 2011, MNRAS, 418, 2336 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167 Astropy Collaboration, Price-Whelan, A. M., Sip˝ocz, B. M., et al. 2018, AJ, 156, 123 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33
work page 2011
- [2]
- [3]
- [4]
- [5]
-
[6]
Brentjens, M. A. & de Bruyn, A. G. 2005, A&A, 441, 1217
work page 2005
-
[7]
Cabral, B. & Leedom, L. C. 1993, in Proceedings of the 20th annual conference on Computer graphics and interactive techniques, 263–270 CASA Team, Bean, B., Bhatnagar, S., et al. 2022, PASP, 134, 114501
work page 1993
-
[8]
Chiu, H. H. S., Ruszkowski, M., Thomas, T., Werhahn, M., & Pfrommer, C. 2024, ApJ, 976, 136
work page 2024
-
[9]
2021, CARTA: The Cube Analysis and Rendering Tool for Astronomy de Bruyn, A
Comrie, A., Wang, K.-S., Hsu, S.-C., et al. 2021, CARTA: The Cube Analysis and Rendering Tool for Astronomy de Bruyn, A. G. 1977, A&A, 58, 221
work page 2021
-
[10]
Dettmar, R. J. 1990, A&A, 232, L15
work page 1990
-
[11]
English, J., Richardson, M. L. A., Ferrand, G., & Deg, N. 2024, CosmosCanvas: Useful color maps for different astrophysical properties, Astrophysics Source Code Library, record ascl:2401.005 Ferrière, K. & Terral, P. 2014, A&A, 561, A100
work page 2024
-
[12]
Fletcher, A., Beck, R., Shukurov, A., Berkhuijsen, E. M., & Horellou, C. 2011, MNRAS, 412, 2396
work page 2011
- [13]
-
[14]
Heald, G. 2009, in IAU Symposium, V ol. 259, Cosmic Magnetic Fields: From Planets, to Stars and Galaxies, ed. K. G. Strassmeier, A. G. Kosovichev, & J. E. Beckman, 591–602
work page 2009
-
[15]
Heald, G. H., Heesen, V ., Sridhar, S. S., et al. 2022, MNRAS, 509, 658
work page 2022
- [16]
-
[17]
Heesen, V ., O’Sullivan, S. P., Brüggen, M., et al. 2023, A&A, 670, L23
work page 2023
-
[18]
Heesen, V ., Stein, M., Pourjafari, N., et al. 2025, A&A, 699, A243
work page 2025
-
[19]
Heesen, V ., Wiegert, T., Irwin, J., et al. 2024, A&A, 691, A273
work page 2024
-
[20]
Henriksen, R. N. 2017, MNRAS, 469, 4806
work page 2017
-
[21]
Henriksen, R. N. 2022, A&A, 658, A101
work page 2022
-
[22]
Horellou, C., Beck, R., Berkhuijsen, E. M., Krause, M., & Klein, U. 1992, A&A, 265, 417
work page 1992
- [23]
-
[24]
Hutschenreuter, S., Anderson, C. S., Betti, S., et al. 2022, A&A, 657, A43
work page 2022
- [25]
- [26]
-
[27]
A., Schmidt, P., Damas-Segovia, A., et al
Irwin, J. A., Schmidt, P., Damas-Segovia, A., et al. 2017, MNRAS, 464, 1333
work page 2017
- [28]
-
[29]
Kamphuis, P., Holwerda, B. W., Allen, R. J., Peletier, R. F., & van der Kruit, P. C. 2007, A&A, 471, L1
work page 2007
-
[30]
Kleimann, J., Schorlepp, T., Merten, L., & Becker Tjus, J. 2019, ApJ, 877, 76
work page 2019
- [31]
- [32]
-
[33]
P., Reinhardt, M., & Simard-Normandin, M
Kronberg, P. P., Reinhardt, M., & Simard-Normandin, M. 1977, A&A, 61, 771
work page 1977
-
[34]
B., Armillotta, L., Ostriker, E
Linzer, N. B., Armillotta, L., Ostriker, E. C., & Quataert, E. 2025, ApJ, 988, 214
work page 2025
-
[35]
Makarov, D., Prugniel, P., Terekhova, N., Courtois, H., & Vauglin, I. 2014, A&A, 570, A13
work page 2014
-
[36]
Mohan, N. & Rafferty, D. 2015, PyBDSF: Python Blob Detection and Source
work page 2015
-
[37]
C., Krause, M., Basu, A., et al
Mora-Partiarroyo, S. C., Krause, M., Basu, A., et al. 2019, A&A, 632, A11
work page 2019
- [38]
-
[39]
Myserlis, I. & Contopoulos, I. 2021, A&A, 649, A94 O’Sullivan, S. P., Brüggen, M., Vazza, F., et al. 2020, MNRAS, 495, 2607
work page 2021
-
[40]
Pakmor, R., van de V oort, F., Bieri, R., et al. 2020, MNRAS, 498, 3125
work page 2020
-
[41]
Pandhi, A., Gaensler, B. M., Pleunis, Z., et al. 2025, ApJ, 982, 146
work page 2025
-
[42]
Purcell, C. R., Van Eck, C. L., West, J., Sun, X. H., & Gaensler, B. M. 2020, RM- Tools: Rotation measure (RM) synthesis and Stokes QU-fitting, Astrophysics Source Code Library, record ascl:2005.003
work page 2020
-
[43]
Qazi, Y ., Shukurov, A., Tharakkal, D., Gent, F. A., & Bendre, A. B. 2024, MN- RAS, 527, 7994
work page 2024
- [44]
-
[45]
Robitaille, T. & Bressert, E. 2012, APLpy: Astronomical Plotting Library in
work page 2012
-
[46]
Ruzmaikin, A. A., Sokolov, D. D., & Shukurov, A. M. 1988, Magnetic Fields of
work page 1988
-
[47]
Shukurov, A., Rodrigues, L. F. S., Bushby, P. J., Hollins, J., & Rachen, J. P. 2019, A&A, 623, A113
work page 2019
-
[48]
Sike, B., Thomas, T., Ruszkowski, M., Pfrommer, C., & Weber, M. 2025, ApJ, 987, 204
work page 2025
-
[49]
2025, Master’s thesis, Queen’s University
Skeggs, N. 2025, Master’s thesis, Queen’s University
work page 2025
-
[50]
Sokoloff, D. D., Bykov, A. A., Shukurov, A., et al. 1998, MNRAS, 299, 189
work page 1998
- [51]
-
[52]
Stein, M., Kleimann, J., Adebahr, B., et al. 2025, A&A, 696, A112
work page 2025
-
[53]
Stein, Y ., Dettmar, R. J., Beck, R., et al. 2020, A&A, 639, A111
work page 2020
-
[54]
Stil, J. M., Keller, B. W., George, S. J., & Taylor, A. R. 2014, ApJ, 787, 99
work page 2014
-
[55]
Sun, X. H. & Reich, W. 2009, A&A, 507, 1087
work page 2009
-
[56]
Tabatabaei, F. S., Minguez, P., Prieto, M. A., & Fernández-Ontiveros, J. A. 2018, Nature Astronomy, 2, 83
work page 2018
-
[57]
Unger, M. & Farrar, G. R. 2024, ApJ, 970, 95 Van Eck, C. L., Purcell, C. R., Baidoo, L., et al. 2026, arXiv e-prints, arXiv:2601.20092
-
[58]
2023, Science Advances, 9, eade7233
Vernstrom, T., West, J., Vazza, F., et al. 2023, Science Advances, 9, eade7233
work page 2023
- [59]
-
[60]
Werhahn, M., Pfrommer, C., & Girichidis, P. 2021, MNRAS, 508, 4072
work page 2021
-
[61]
Whitmore, B. C., Lucas, R. A., McElroy, D. B., et al. 1990, AJ, 100, 1489
work page 1990
- [62]
- [63]
-
[64]
2014, AJ, 148, 127 Article number, page 15 of 29 A&A proofs:manuscript no
Yim, K., Wong, T., Xue, R., et al. 2014, AJ, 148, 127 Article number, page 15 of 29 A&A proofs:manuscript no. aa58397-25 Appendix A: Properties of the galaxy sample Table A.1: Fundamental information about the initial galaxy sample. Galaxy RA Dec PA Q AS D d ang dphy Bmaj Bmaj Nd Beam RMMW [H:M:S] [D:M:S] [deg] [Mpc] [ ′] [kpc] [ ′′] [kpc] [rad m −2] NGC ...
work page 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.