On the Difference Between Pulsar Radio Emission Beams from the Two Poles
Pith reviewed 2026-05-07 07:33 UTC · model grok-4.3
The pith
Pulsar emission beams from the two magnetic poles differ in size, with none of the studied cases showing similar azimuth widths.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The emission beams from the two magnetic poles of a pulsar may be generally dissimilar in size. Among eight double-pole pulsars, none shows similar azimuth width between poles, only two show potentially similar beam radii, and six show comparable emission intensities in specific parameter spaces. The results imply that physical conditions for pair production and acceleration differ between poles and that the polar cap contains randomly distributed active regions due to local magnetic field or surface variations.
What carries the argument
Rotating vector model with aberration and retardation effects, used to determine intrinsic emission regions from observed pulse windows and position angle swings, combined with conal and fan beam models to compare radius, azimuth width, and intensity between main pulse and interpulse.
If this is right
- Pair production and particle acceleration must operate under different physical conditions at the two poles rather than symmetrically.
- The polar cap contains inhomogeneous active emission regions that vary randomly rather than filling uniformly.
- Pulsar beam models need to treat the two poles as potentially independent instead of mirror images.
- Differences in local magnetic field structure or surface properties can produce observable asymmetries in radio emission.
Where Pith is reading between the lines
- Models of neutron star magnetospheres may need to incorporate pole-to-pole variations in magnetic field curvature or temperature to explain the observed beam differences.
- Population studies that assume symmetric beams when estimating pulsar birth rates or beaming fractions could systematically miscalculate the fraction of detectable objects.
- Higher-resolution polarization data on additional interpulse sources could test whether the asymmetry correlates with spin period, age, or magnetic field strength.
Load-bearing premise
The rotating vector model with aberration and retardation effects accurately recovers the true emission regions inside the observed pulse windows, and the conal and fan beam models provide an unbiased way to compare beam parameters between the two poles.
What would settle it
Discovery of multiple additional double-pole pulsars in which the main pulse and interpulse have matching beam radii, azimuth widths, and intensities across the same model parameter ranges.
Figures
read the original abstract
The long-standing assumption of symmetric radio emission beams from the two magnetic poles of pulsars is challenged by observational evidence of asymmetry and underfill. Direct testing of this symmetry remains difficult for most pulsars. As an indirect test, we collected polarization profiles of 11 interpulse pulsars observed with the Five-hundred-meter Aperture Spherical radio Telescope, MeerKAT, and Parkes. We developed a rotating vector model incorporating aberration and retardation effects to fit the position angle swings of selected pulsars, thereby determining the intrinsic emission region corresponding to the observed pulse windows. Based on both the conal and fan beam models, we compared three key parameters-beam radius, magnetic azimuth width, and emission intensity-between the intrinsic emission regions of the main pulse and interpulse. Among the eight pulsars with a confirmed double-pole geometry, none exhibits similarity in the azimuth width. Only two show potentially similar beam radii, while six demonstrate comparable emission intensities within specific parameter spaces. These results indicate that the emission beams from the two magnetic poles of a pulsar may be generally dissimilar in size, suggesting that the physical conditions governing pair production and particle acceleration differ between the two poles. The random distribution of active emission regions further implies inhomogeneity within the polar cap, which may originate from the differences in local magnetic field structure or surface properties.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes polarization profiles of 11 interpulse pulsars observed with FAST, MeerKAT, and Parkes. A rotating-vector model incorporating aberration and retardation is fitted to the position-angle swings to map observed pulse windows onto intrinsic emission regions in magnetic coordinates. For the eight pulsars with confirmed double-pole geometry, beam radius, magnetic azimuth width, and emission intensity are compared between main-pulse and interpulse regions under both conal and fan-beam models. The results indicate that none of the eight show similar azimuth widths, only two show similar radii, and six show comparable intensities within certain parameter spaces, leading to the conclusion that emission beams from the two poles are generally dissimilar in size and that the polar caps are inhomogeneous.
Significance. If the geometric reconstructions prove robust, the work supplies direct observational evidence against the long-standing assumption of symmetric radio beams from both magnetic poles. It implies that pair-production and acceleration conditions can differ between poles and that active emission regions are randomly distributed within the polar cap, with possible origins in local magnetic-field or surface inhomogeneities. These findings would constrain emission models and magnetospheric physics. The multi-telescope dataset and application to a sample of double-pole sources are strengths; however, the conclusions rest on the fidelity of the RVM mapping and the quantitative criteria used for similarity.
major comments (2)
- [§3.2] §3.2 (RVM fitting procedure): No posterior sampling, Monte-Carlo error propagation, or sensitivity analysis over emission height r_em is reported. Because aberration/retardation shifts scale linearly with r_em (typically 100–1000 km) and the PA swing amplitude depends on both α and β, the derived intrinsic beam radii and azimuth widths inherit a degeneracy. The statements that “none exhibits similarity in the azimuth width” and “only two show potentially similar beam radii” therefore depend on the specific height chosen for each pulsar; modest changes in r_em can alter which pairs are classified as similar. A marginalization over plausible height ranges or explicit uncertainty bands on the mapped parameters is required before the dissimilarity claim can be considered robust.
- [Results] Results section (parameter comparison tables/figures): The criteria used to decide “similar” or “comparable” are not defined quantitatively (e.g., ratio within 1.2, difference < 3σ, or overlap of 1σ intervals). Without tabulated best-fit values, uncertainties, or a clear similarity metric for beam radius, azimuth width, and intensity, it is impossible to assess whether the reported differences are statistically significant or could be consistent with similarity once modeling uncertainties are included. This directly undermines the central claim that the beams are “generally dissimilar in size.”
minor comments (3)
- [§2] The selection criteria for the 11 pulsars and the subset of eight double-pole cases should be stated explicitly, including any cuts on signal-to-noise or PA-swing quality.
- [Tables] Tables listing the fitted geometric parameters (α, β, r_em, beam radius, azimuth width) together with their uncertainties would improve reproducibility and allow readers to judge the similarity statements directly.
- [§4] Notation for “magnetic azimuth width” in the fan-beam model should be defined mathematically (e.g., Δφ_m or equivalent) and distinguished from the conal-beam equivalent to avoid ambiguity when comparing the two models.
Simulated Author's Rebuttal
We are grateful to the referee for their thorough and constructive review of our manuscript. The comments correctly identify areas where additional analysis and clarification will strengthen the robustness of our conclusions on the dissimilarity of pulsar emission beams. We address each major comment below and will incorporate the suggested improvements in the revised version.
read point-by-point responses
-
Referee: [§3.2] §3.2 (RVM fitting procedure): No posterior sampling, Monte-Carlo error propagation, or sensitivity analysis over emission height r_em is reported. Because aberration/retardation shifts scale linearly with r_em (typically 100–1000 km) and the PA swing amplitude depends on both α and β, the derived intrinsic beam radii and azimuth widths inherit a degeneracy. The statements that “none exhibits similarity in the azimuth width” and “only two show potentially similar beam radii” therefore depend on the specific height chosen for each pulsar; modest changes in r_em can alter which pairs are classified as similar. A marginalization over plausible height ranges or explicit uncertainty bands on the mapped parameters is required before the dissimilarity claim can be considered robust.
Authors: We agree that the original analysis lacked a systematic sensitivity study over emission height and did not propagate uncertainties from r_em. The fitting relied on χ² minimization with adopted r_em values drawn from the literature. To address the degeneracy, the revised manuscript will include an explicit sensitivity analysis in which r_em is varied from 100 to 1000 km for each of the eight double-pole pulsars. We will demonstrate that the azimuth widths differ by factors of at least 1.8 in all cases across this range, while beam radii remain similar in only the same two pulsars. Uncertainty bands will be shown as the envelope of models with χ² within 10% of the minimum. Although a full Bayesian posterior sampling was not performed in the original work owing to computational cost, the sensitivity analysis provides a practical quantification of the height dependence and will be added to the paper. revision: yes
-
Referee: [Results] Results section (parameter comparison tables/figures): The criteria used to decide “similar” or “comparable” are not defined quantitatively (e.g., ratio within 1.2, difference < 3σ, or overlap of 1σ intervals). Without tabulated best-fit values, uncertainties, or a clear similarity metric for beam radius, azimuth width, and intensity, it is impossible to assess whether the reported differences are statistically significant or could be consistent with similarity once modeling uncertainties are included. This directly undermines the central claim that the beams are “generally dissimilar in size.”
Authors: The referee is correct that the original text did not supply quantitative similarity thresholds or tabulated best-fit values with uncertainties. This made independent assessment of the differences difficult. In the revised manuscript we will add a new table listing the best-fit beam radius, magnetic azimuth width, and relative intensity (with 1σ uncertainties from the fit and r_em sensitivity) for the main pulse and interpulse of each pulsar. We will define the metrics explicitly: beams are similar in radius or azimuth width if their ratio lies between 2/3 and 3/2 and the absolute difference is less than twice the combined uncertainty; intensities are comparable if they agree within a factor of three. Application of these criteria leaves the conclusions unchanged—none of the eight pulsars show similar azimuth widths, only two show similar radii, and six show comparable intensities—while allowing readers to judge statistical significance directly. revision: yes
Circularity Check
No significant circularity; derivation uses external observations and standard models
full rationale
The paper's chain begins with collected polarization profiles from independent telescope observations (FAST, MeerKAT, Parkes). It applies the established rotating-vector model augmented by known aberration and retardation effects to fit position-angle swings and map observed windows onto magnetic coordinates. Beam parameters are then extracted and compared using the standard conal and fan-beam frameworks. None of these steps reduces by the paper's own equations to quantities defined solely in terms of its fitted outputs or to self-citations; the dissimilarity conclusion is an empirical comparison of independently reconstructed regions. The derivation remains self-contained against external pulsar-geometry benchmarks.
Axiom & Free-Parameter Ledger
free parameters (3)
- beam radius
- magnetic azimuth width
- emission intensity
axioms (3)
- domain assumption The rotating vector model with aberration and retardation effects accurately maps observed position angle swings to intrinsic emission regions.
- domain assumption Conal and fan beam models are appropriate representations for comparing emission geometry parameters between the two poles.
- domain assumption The subsample of eight pulsars with confirmed double-pole geometry is representative for generalizing about beam dissimilarity.
Reference graph
Works this paper leans on
-
[1]
Hogg, D. W., & O’Neil, M. 2015, IEEE Transactions on Pattern Analysis and Machine Intelligence, 38, 252, doi: 10.1109/TPAMI.2015.2448083
-
[2]
Arzamasskiy, L. I., Beskin, V. S., & Pirov, K. K. 2017, MNRAS, 466, 2325, doi: 10.1093/mnras/stw3139
-
[3]
1996, ApJ, 470, 1111, doi: 10.1086/177935
Wolszczan, A. 1996, ApJ, 470, 1111, doi: 10.1086/177935
-
[4]
Ashton, G., et al. 2019, Astrophys. J. Suppl., 241, 27, doi: 10.3847/1538-4365/ab06fc
-
[5]
2010, ApJ, 715, 1282, doi: 10.1088/0004-637X/715/2/1282 Baub¨ ock, M., Psaltis, D., &¨Ozel, F
Bai, X.-N., & Spitkovsky, A. 2010, ApJ, 715, 1282, doi: 10.1088/0004-637X/715/2/1282
-
[6]
Blaskiewicz, M., Cordes, J. M., & Wasserman, I. 1991, ApJ, 370, 643, doi: 10.1086/169850
-
[7]
R., Karastergiou, A., & Johnston, S
Brook, P. R., Karastergiou, A., & Johnston, S. 2019, MNRAS, 488, 5702, doi: 10.1093/mnras/stz2092
-
[8]
Desvignes, G., Kramer, M., Lee, K., et al. 2019, Science, 365, 1013, doi: 10.1126/science.aav7272
-
[9]
2003, ApJ, 598, 1201, doi: 10.1086/379052
Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201, doi: 10.1086/379052
-
[10]
Dyks, J., Rudak, B., & Demorest, P. 2010, MNRAS, 401, 1781, doi: 10.1111/j.1365-2966.2009.15679.x
-
[11]
Gangadhara, R. T., & Gupta, Y. 2001, ApJ, 555, 31, doi: 10.1086/321439
-
[12]
2006, ApJ, 650, 1048, doi: 10.1086/506982
Gil, J., Melikidze, G., & Zhang, B. 2006, ApJ, 650, 1048, doi: 10.1086/506982
-
[13]
Gil, J. A., & Lyne, A. G. 1995, MNRAS, 276, L55, doi: 10.1093/mnras/276.1.L55
-
[14]
Goldreich, P., & Julian, W. H. 1969, ApJ, 157, 869, doi: 10.1086/150119
-
[15]
E., Lupsasca, A., & Philippov, A
Gralla, S. E., Lupsasca, A., & Philippov, A. 2017, ApJ, 851, 137, doi: 10.3847/1538-4357/aa978d
-
[16]
Han, J. L., Wang, C., Wang, P. F., et al. 2021, Research in Astronomy and Astrophysics, 21, 107, doi: 10.1088/1674-4527/21/5/107 15
-
[17]
Han, J. L., Zhou, D. J., Wang, C., et al. 2025, Research in Astronomy and Astrophysics, 25, 014001, doi: 10.1088/1674-4527/ada3b7
-
[18]
Statistics and Computing 29(5), 891 (2019) https: //doi.org/10.1007/s11222-018-9844-0
Higson, E., Handley, W., Hobson, M., & Lasenby, A. 2019, Statistics and Computing, 29, 891, doi: 10.1007/s11222-018-9844-0
-
[19]
Hulse, R. A., & Taylor, J. H. 1975, ApJL, 195, L51, doi: 10.1086/181708
-
[20]
2019, MNRAS, 485, 640, doi: 10.1093/mnras/stz400
Johnston, S., & Karastergiou, A. 2019, MNRAS, 485, 640, doi: 10.1093/mnras/stz400
-
[21]
2018, MNRAS, 474, 4629, doi: 10.1093/mnras/stx3095
Johnston, S., & Kerr, M. 2018, MNRAS, 474, 4629, doi: 10.1093/mnras/stx3095
-
[22]
2019, MNRAS, 490, 4565, doi: 10.1093/mnras/stz2865
Johnston, S., & Kramer, M. 2019, MNRAS, 490, 4565, doi: 10.1093/mnras/stz2865
-
[23]
Application of the rotating vector model
Johnston, S., Kramer, M., Karastergiou, A., et al. 2023, MNRAS, 520, 4801, doi: 10.1093/mnras/stac3636
-
[24]
2024, Monthly Notices of the Royal Astronomical Society, 530, 4839, doi: 10.1093/mnras/stae1175
Karastergiou, A. 2024, MNRAS, 530, 4839, doi: 10.1093/mnras/stae1175
-
[25]
The Thousand-Pulsar-Array programme on MeerKAT – I. Science objectives and first results
Johnston, S., Karastergiou, A., Keith, M. J., et al. 2020, MNRAS, 493, 3608, doi: 10.1093/mnras/staa516
-
[26]
Kaspi, V. M., Lyne, A. G., Manchester, R. N., et al. 2000, ApJ, 543, 321, doi: 10.1086/317103
-
[27]
Keith, M. J., Johnston, S., Weltevrede, P., & Kramer, M. 2010, MNRAS, 402, 745, doi: 10.1111/j.1365-2966.2009.15926.x
-
[28]
Komesaroff, M. M. 1970, Nature, 225, 612, doi: 10.1038/225612a0
-
[29]
Kramer, M., & Johnston, S. 2008, MNRAS, 390, 87, doi: 10.1111/j.1365-2966.2008.13780.x —. 2026, MNRAS, 547, staf2258, doi: 10.1093/mnras/staf2258
-
[30]
Lee, K. J., Cui, X. H., Wang, H. G., Qiao, G. J., & Xu, R. X. 2009, ApJ, 703, 507, doi: 10.1088/0004-637X/703/1/507
-
[31]
V., Melchior, A.-L., & Zolotukhin, I
Lee, K. J., Du, Y. J., Wang, H. G., et al. 2010, MNRAS, 405, 2103, doi: 10.1111/j.1365-2966.2010.16600.x
-
[32]
2014, MNRAS, 441, 690, doi: 10.1093/mnras/stu564
Lomiashvili, D., & Lyutikov, M. 2014, MNRAS, 441, 690, doi: 10.1093/mnras/stu564
-
[33]
Lorimer, D. R., Stairs, I. H., Freire, P. C., et al. 2006, ApJ, 640, 428, doi: 10.1086/499918
-
[34]
Lyne, A. G., & Manchester, R. N. 1988, MNRAS, 234, 477, doi: 10.1093/mnras/234.3.477
-
[35]
Lyne, A. G., Burgay, M., Kramer, M., et al. 2004, Science, 303, 1153, doi: 10.1126/science.1094645
-
[36]
Maciesiak, K., Gil, J., & Ribeiro, V. A. R. M. 2011, MNRAS, 414, 1314, doi: 10.1111/j.1365-2966.2011.18471.x
-
[37]
Manchester, R. N. 1995, Journal of Astrophysics and Astronomy, 16, 107, doi: 10.1007/BF02714828
-
[38]
Manchester, R. N., Kramer, M., Stairs, I. H., et al. 2010, ApJ, 710, 1694, doi: 10.1088/0004-637X/710/2/1694
-
[39]
2024, ApJ, 966, 46, doi: 10.3847/1538-4357/ad381c
Meng, L., Zhu, W., Kramer, M., et al. 2024, ApJ, 966, 46, doi: 10.3847/1538-4357/ad381c
-
[40]
2016, MNRAS, 460, 3063, doi: 10.1093/mnras/stw1186
Mitra, D., Rankin, J., & Arjunwadkar, M. 2016, MNRAS, 460, 3063, doi: 10.1093/mnras/stw1186
-
[41]
Muslimov, A. G., & Harding, A. K. 2003, ApJ, 588, 430, doi: 10.1086/368162
-
[42]
Perera, B. B. P., Kim, C., McLaughlin, M. A., et al. 2014, ApJ, 787, 51, doi: 10.1088/0004-637X/787/1/51
-
[43]
McLaughlin, M. A., & Lyutikov, M. 2012, ApJ, 750, 130, doi: 10.1088/0004-637X/750/2/130
-
[44]
2021, MNRAS, 508, 4249, doi: 10.1093/mnras/stab2775
Posselt, B., Karastergiou, A., Johnston, S., et al. 2021, MNRAS, 508, 4249, doi: 10.1093/mnras/stab2775
-
[45]
2022, The Thousand-Pulsar-Array program on MeerKAT – IX
Posselt, B., Karastergiou, A., Johnston, S., et al. 2022, The Thousand-Pulsar-Array program on MeerKAT – IX. The time-averaged properties of the observed pulsar population: data set, Zenodo, doi: 10.5281/ZENODO.7272361
-
[46]
The time-averaged properties of the observed pulsar population
Posselt, B., Karastergiou, A., Johnston, S., et al. 2023, MNRAS, 520, 4582, doi: 10.1093/mnras/stac3383
-
[47]
Qiao, G. J., Lee, K. J., Wang, H. G., Xu, R. X., & Han, J. L. 2004, ApJL, 606, L49, doi: 10.1086/421048
-
[48]
Qiao, G. J., Liu, J. F., Zhang, B., & Han, J. L. 2001, A&A, 377, 964, doi: 10.1051/0004-6361:20011188
-
[49]
Qiu, J. L., Tong, H., & Wang, H. G. 2023, ApJ, 958, 78, doi: 10.3847/1538-4357/ad003f
-
[50]
Radhakrishnan, V., & Cooke, D. J. 1969, Astrophys. Lett., 3, 225
1969
-
[51]
Weisberg, J. M., & Devine, K. E. 2004, ApJ, 606, 1167, doi: 10.1086/383179
-
[52]
Rankin, J. M. 1983, ApJ, 274, 333, doi: 10.1086/161450
-
[53]
Rankin, J. M., & Rathnasree, N. 1995, Journal of Astrophysics and Astronomy, 16, 327, doi: 10.1007/BF02715608
-
[54]
Romani, R. W., & Yadigaroglu, I.-A. 1995, ApJ, 438, 314, doi: 10.1086/175076
-
[55]
C., Weltevrede, P., & Johnston, S
Rookyard, S. C., Weltevrede, P., & Johnston, S. 2015, MNRAS, 446, 3367, doi: 10.1093/mnras/stu2236
-
[56]
Ruderman, M. A., & Sutherland, P. G. 1975, ApJ, 196, 51, doi: 10.1086/153393
-
[57]
2017, MNRAS, 467, 2529, doi: 10.1093/mnras/stx204
Saha, L., & Dyks, J. 2017, MNRAS, 467, 2529, doi: 10.1093/mnras/stx204
-
[58]
2021, MNRAS, 505, 4483, doi: 10.1093/mnras/staa2811 16
Serylak, M., Johnston, S., Kramer, M., et al. 2021, MNRAS, 505, 4483, doi: 10.1093/mnras/staa2811 16
-
[59]
Skilling, J. 2004, in American Institute of Physics Conference Series, Vol. 735, Bayesian Inference and Maximum Entropy Methods in Science and Engineering: 24th International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, ed. R. Fischer, R. Preuss, & U. V. Toussaint (AIP), 395–405, doi: 10.1063/1.1835238
-
[60]
Bayesian Analysis 1(4), 833–859 (2006) https://doi.org/10.1214/06-BA127
Skilling, J. 2006, Bayesian Analysis, 1, 833 , doi: 10.1214/06-BA127
-
[61]
Speagle, J. S. 2020, MNRAS, 493, 3132, doi: 10.1093/mnras/staa278
-
[62]
Stovall, K., Freire, P. C. C., Chatterjee, S., et al. 2018, ApJL, 854, L22, doi: 10.3847/2041-8213/aaad06
-
[63]
Sun, S. N., Wang, N., Yan, W. M., & Wang, S. Q. 2025, ApJ, 983, 179, doi: 10.3847/1538-4357/adc254
-
[64]
2013, arXiv e-prints, arXiv:1304.4203, doi: 10.48550/arXiv.1304.4203
Szary, A. 2013, arXiv e-prints, arXiv:1304.4203, doi: 10.48550/arXiv.1304.4203
-
[65]
Wang, H. G., Pi, F. P., Zheng, X. P., et al. 2014, ApJ, 789, 73, doi: 10.1088/0004-637X/789/1/73
-
[66]
Wang, P. F., Han, J. L., Xu, J., et al. 2023, Research in Astronomy and Astrophysics, 23, 104002, doi: 10.1088/1674-4527/acea1f
-
[67]
Weisberg, J. M., & Taylor, J. H. 2002, ApJ, 576, 942, doi: 10.1086/341803
-
[68]
Weltevrede, P., & Johnston, S. 2008, MNRAS, 387, 1755, doi: 10.1111/j.1365-2966.2008.13382.x
-
[69]
Weltevrede, P., & Wright, G. 2009, MNRAS, 395, 2117, doi: 10.1111/j.1365-2966.2009.14643.x
-
[70]
1991, Nature, 350, 688, doi: 10.1038/350688a0
Wolszczan, A. 1991, Nature, 350, 688, doi: 10.1038/350688a0
-
[71]
Xu, Z.-H., Wang, W.-Y., Cao, S.-S., & Xu, R.-X. 2026, Research in Astronomy and Astrophysics, 26, 035014, doi: 10.1088/1674-4527/ae2b59 17 APPENDIX A.THE INNER OPEN FIELD LINES AND COORDINATE TRANSFORMATION This appendix derives the analytic equation for inner field lines whose footprints are uniformly distributed on the stellar surface and provides the t...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.