arxiv: 2604.15501 · v1 · submitted 2026-04-16 · 🌌 astro-ph.HE

Recognition: unknown

Refined Constraints on the Hard X-ray Polarization of the Crab Pulsar and Nebula Derived from an Extended XL-Calibur Dataset

Matthew G. Baring , Jacob Casey , Sohee Chun , Ephraim Gau , Tomohiro Hakamata , Kun Hu , Daiki Ishi , Fabian Kislat

show 21 more authors

M\'ozsi Kiss Merlin Kole Henric Krawczynski Haruki Kuramoto Lindsey Lisalda Bingkun Liu Yoshitomo Maeda Hironori Matsumoto Shravan Vengalil Menon Takuya Miyazawa Kaito Murakami Takashi Okajima Mark Pearce Brian Rauch Kentaro Shirahama Sean Spooner Hiromitsu Takahashi Sayana Takatsuka Yuusuke Uchida Varun Andrew Thomas West

Authors on Pith no claims yet

Pith reviewed 2026-05-10 09:42 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords Crab pulsarX-ray polarizationnebulaphase recoverysynchrotron emissionballoon-bornetiming reconstructionGPS

0 comments

The pith

By treating the Crab pulsar as its own clock, XL-Calibur recovers timing to refine nebular hard X-ray polarization to 27.7 percent.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper addresses incomplete timing data from a balloon-borne X-ray polarimeter due to GPS issues during observations of the Crab pulsar and nebula. By using the pulsar's regular pulses as a timing reference and applying a Markov-Chain Monte-Carlo fit to recover phase information, nearly all the affected data becomes usable for analysis. This yields a polarization degree for the nebular emission of about 28 percent at an angle aligned with the Crab's spin axis. A sympathetic reader would care because it strengthens the evidence that hard X-ray light comes from synchrotron processes in the inner nebula regions, and shows a practical way to salvage timing in space observations. Phase-resolved results further support emission models distinguishing pulsar and nebula contributions.

Core claim

The paper establishes that a Markov-Chain Monte-Carlo joint fit of phase offsets and frequency derivatives during GPS-off intervals reconstructs timing with sufficient accuracy to include 95 percent of the affected data in the polarization analysis. This extended dataset confirms the nebular emission polarization degree of (27.7 ± 4.9)% at a polarization angle of 127.2° ± 5.1°, which remains aligned with the Crab's spin axis and is consistent with synchrotron emission from the inner nebula. Phase-resolved measurements indicate strong polarization in the off-pulse and bridge intervals, with the pulsar peaks weakly constrained but agreeing with lower-energy observations.

What carries the argument

A Markov-Chain Monte-Carlo framework that uses the Crab pulsar's 33 ms period to jointly fit phase offsets and frequency derivatives for timing reconstruction during GPS-off periods.

If this is right

The nebular hard X-ray polarization is confirmed at 27.7 percent and aligned with the spin axis, supporting synchrotron emission from the inner nebula.
Off-pulse and bridge phases exhibit strong polarization while the pulsar peaks are weakly constrained.
The findings reinforce that hard X-ray emission arises primarily in the nebular torus and wind regions.
The phase recovery method increases the usable dataset to nearly 100 percent of observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This timing recovery technique may prove useful for other balloon-borne or space-based instruments facing intermittent clock signal losses.
Future studies could test whether applying the method to different energy bands or sources yields similar improvements in polarization precision.
Combining the phase-resolved hard X-ray data with soft X-ray polarization measurements could better map the emission geometry across the Crab system.

Load-bearing premise

The Markov-Chain Monte-Carlo reconstruction accurately determines the phase tags without introducing systematic biases into the polarization measurements.

What would settle it

A significant discrepancy between the polarization results obtained with and without the recovered GPS-off data, or an independent verification showing that the fitted phases do not match the actual pulsar arrival times within the required precision.

Figures

Figures reproduced from arXiv: 2604.15501 by Andrew Thomas West, Bingkun Liu, Brian Rauch, Daiki Ishi, Ephraim Gau, Fabian Kislat, Haruki Kuramoto, Henric Krawczynski, Hiromitsu Takahashi, Hironori Matsumoto, Jacob Casey, Kaito Murakami, Kentaro Shirahama, Kun Hu, Lindsey Lisalda, Mark Pearce, Matthew G. Baring, Merlin Kole, M\'ozsi Kiss, Sayana Takatsuka, Sean Spooner, Shravan Vengalil Menon, Sohee Chun, Takashi Okajima, Takuya Miyazawa, Tomohiro Hakamata, Varun, Yoshitomo Maeda, Yuusuke Uchida.

**Figure 1.** Figure 1: On-source (red) and off-source (blue) count rates during the three Crab observations (Days 2–4). Each cluster of points corresponds to a single observation day. Observations of Cyg X−1 occur during gaps between observations when the Crab is out of view, see H. Awaki et al. (2025b) for more details on these observations. The modulation of the on-source rate is caused by changes in pointing elevation, with t… view at source ↗

**Figure 2.** Figure 2: Time segmentation of the Crab observation showing GPS-on and GPS-off intervals. The x-axis is truncated to only depict observations of the Crab. While the GPS problem persisted during the Cyg X−1 observations, the same data loss did not occur due to the lower required precision for absolute timing for the source. Source events with valid GPS timestamps are shown on the top panel in blue, while those lackin… view at source ↗

**Figure 3.** Figure 3: Folded light curve of the Crab pulsar derived from GPS-on data shown in blue, overlaid with the best fit 18-term Fourier template in red. The pulse profile is shown over two full rotations for clarity. The agreement between the data and the smooth model confirms the adequacy of the harmonic expansion and supports it’s use as a reference template for phase recovery in GPS-off intervals. P1 in the GPS-on dat… view at source ↗

**Figure 4.** Figure 4: Propagation of phase uncertainties from the MCMC posterior distributions into total phase error across the Day 4 interval. Contributions from each parameter offset: ϕ5 (orange), νe (green), and ˙νe (red), are shown individually, along with the combined total phase uncertainty (blue). The maximum error at any point remains well below the width of phase intervals used for phase-resolved polarization analys… view at source ↗

**Figure 5.** Figure 5: The top panel depicts the time- and phase-resolved pulse profile, with phase on the horizontal and time on the vertical axis. The color scale corresponds to the number of events in a given bin. Brighter regions depict bins with higher count rates. Clear vertical lines are visible at ϕ ∼ 1 and ∼ 1.3 due to the pulsar peaks. The middle panel shows the reconstructed, time-integrated pulse profile for the Da… view at source ↗

**Figure 6.** Figure 6: Similar to [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

**Figure 9.** Figure 9: Energy-resolved polarization measurements of the Crab Nebula emission from IXPE and OSO−8, reported in W. Wei et al. (2025) and M. C. Weisskopf et al. (1978) respectively, and compared with the expanded XL-Calibur results. Data from XL-Calibur follows the trend where as energy increases, the nebular polarization approaches the pulsar’s spin axis. The arrow length is proportional to the PD in each energy … view at source ↗

**Figure 10.** Figure 10: Phase-resolved polarization degree (left) and polarization angle (right) for the Crab pulsar in the 19 – 64 keV band, as measured by XL Calibur (red stars), using updated phase tagging when compared to [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗

read the original abstract

We present updated hard X-ray polarization measurements of the Crab pulsar and nebula obtained with the balloon-borne polarimeter XL-Calibur in the ~19-64 keV energy range. During the flight, intermittent GPS-failure resulted in poorly constrained timing for ~38% of the Crab dataset. By implementing a new phase-recovery method that reconstructs timing during extended GPS-off intervals, phase tag data is recovered for ~95% of the GPS-off dataset, increasing the precision of the phase-resolved analysis. Phase-information for the data is recovered by using the Crab pulsar, with its 33 ms period, as an external timing source. Using a Markov-Chain Monte-Carlo framework to jointly fit phase offsets and frequency derivatives, sufficient phase accuracy is achieved, across multiple periods without GPS for a phase-resolved analysis. This enables inclusion of nearly the full dataset in the polarization study. The polarization degree of the nebular emission is found to be (27.7${\pm}$4.9)% at a polarization angle of 127.2{\deg}${\pm}$5.1{\deg} confirming previous XL-Calibur results and remaining aligned with the Crab's spin axis, consistent with synchrotron emission from the inner nebula. Phase-resolved measurements show that the off-pulse and bridge intervals exhibit a strong polarization, while the pulsar peaks, although weakly constrained, remain in agreement with the softer-energy trends of IXPE. These findings reinforce a scenario in which hard X-ray emission arises primarily in the nebular torus and wind regions. The successful recovery of precise phase tagging from GPS-off data demonstrates the capacity to use the pulsar as an external clock even in the case of sparsely populated data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The advance is a phase-recovery trick that salvages most of the GPS-off data by using the pulsar itself as a clock, which tightens the nebular polarization numbers without changing the overall picture.

read the letter

The paper's main new piece is the MCMC method that reconstructs phase tags for the 38 % of the flight where GPS dropped out. They fit offsets and frequency derivatives against the known 33 ms pulsar period and recover timing for 95 % of those intervals, so nearly the full dataset goes into the polarization analysis. That is a practical step for balloon data with timing gaps, and it directly produces the reported nebular polarization of 27.7 ± 4.9 % at 127.2 ± 5.1 degrees, which sits on top of the earlier XL-Calibur values and lines up with the spin axis and IXPE trends at softer energies. The phase-resolved results for off-pulse and bridge intervals also look consistent with synchrotron expectations from the inner nebula and torus. Those are the concrete gains: more exposure and a modest tightening of the error bars on an already-measured quantity. The data are real observations, so the numbers carry weight even if they are incremental. The soft spot is the phase-recovery step itself. The text states that sufficient accuracy is reached, but it does not show recovery tests on simulated light curves, residual comparisons between GPS-on and recovered segments, or checks on how prior choices affect the polarization output. Without those, an undetected systematic offset could still mix pulsed and nebular components in the bins. That is the one place where the chain from raw counts to final PD and PA is least anchored. This paper is for readers who track hard X-ray polarimetry of the Crab or who need timing solutions for gapped balloon or satellite data. A specialist in synchrotron modeling or in planning follow-up observations would find the numbers and the method useful. It is solid enough to send to referees; the measurement is grounded and the technique is worth documenting, even if the validation of the timing solution needs one more round of explicit tests.

Referee Report

1 major / 2 minor

Summary. The manuscript reports refined hard X-ray polarization measurements of the Crab pulsar and nebula from an extended XL-Calibur balloon flight dataset in the ~19-64 keV band. A new MCMC-based phase-recovery technique is presented that uses the 33 ms pulsar as an external clock to reconstruct timing during GPS-off intervals (which affected ~38% of the exposure), recovering phase tags for ~95% of that subset and enabling inclusion of nearly the full dataset. The central results are a nebular polarization degree of (27.7 ± 4.9)% at a polarization angle of 127.2° ± 5.1°, consistent with prior XL-Calibur values and aligned with the Crab spin axis, plus phase-resolved measurements that show strong polarization in off-pulse and bridge intervals while remaining compatible with IXPE trends at lower energies.

Significance. If the phase-recovery accuracy holds, the work tightens constraints on the emission geometry at hard X-ray energies, reinforcing that the nebular torus and wind dominate the polarized flux and that the geometry is consistent with synchrotron processes. The timing-recovery method itself is a practical advance for balloon or space-borne instruments that encounter intermittent GPS or clock failures, potentially increasing the usable exposure in future observations.

major comments (1)

[MCMC phase-recovery method] In the description of the MCMC phase-recovery method: the statement that 'sufficient phase accuracy is achieved' for the phase-resolved polarization analysis lacks quantitative support. No recovery-bias tests on simulated light curves, residual-phase jitter comparisons between GPS-on segments and recovered GPS-off segments, or sensitivity checks to prior choices are reported. Because this step determines whether ~38% of the total exposure can be safely included without mixing pulsed and nebular components, an undetected systematic offset would directly propagate into the reported nebular PD, PA, and the phase-binned results.

minor comments (2)

[Abstract] The abstract contains LaTeX formatting artifacts (e.g., '127.2$ {deg} $ {pm} $5.1$ {deg} $') that should be rendered cleanly in the final version.
[Data analysis section] A brief expansion of the background-subtraction and systematic-error budget (currently summarized only at a high level) would improve reproducibility, even if the central results remain unchanged.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address the single major comment below and agree that additional quantitative validation of the MCMC phase-recovery method is warranted to support the inclusion of the GPS-off data. We will revise the manuscript to incorporate the requested tests and supporting details.

read point-by-point responses

Referee: [MCMC phase-recovery method] In the description of the MCMC phase-recovery method: the statement that 'sufficient phase accuracy is achieved' for the phase-resolved polarization analysis lacks quantitative support. No recovery-bias tests on simulated light curves, residual-phase jitter comparisons between GPS-on segments and recovered GPS-off segments, or sensitivity checks to prior choices are reported. Because this step determines whether ~38% of the total exposure can be safely included without mixing pulsed and nebular components, an undetected systematic offset would directly propagate into the reported nebular PD, PA, and the phase-binned results.

Authors: We thank the referee for identifying this gap. The manuscript currently asserts that 'sufficient phase accuracy is achieved' via the MCMC joint fit of phase offsets and frequency derivatives but does not present the quantitative validation steps requested. We agree this is a substantive omission given the impact on the ~38% GPS-off exposure. In the revised manuscript we will add a dedicated subsection (or appendix) that includes: (i) recovery-bias tests on simulated light curves with known injected offsets and derivatives to quantify recovered-phase accuracy; (ii) direct residual-phase jitter comparisons between GPS-on segments and the recovered GPS-off segments; and (iii) sensitivity checks to the choice of priors in the MCMC. These additions will provide the missing quantitative evidence that phase errors remain small enough to avoid component mixing and bias in the reported polarization results. We view this as a straightforward and necessary improvement. revision: yes

Circularity Check

0 steps flagged

No significant circularity in measurement-based claims

full rationale

The paper reports direct empirical polarization measurements (PD = 27.7 ± 4.9% at PA = 127.2° ± 5.1°) from balloon-borne XL-Calibur data in the 19-64 keV band. The phase-recovery technique employs MCMC fitting of offsets and derivatives during GPS-off intervals to enable inclusion of additional data, but the polarization values are obtained from the processed photon events and are not quantities that reduce to the timing-fit parameters by construction. Prior XL-Calibur results are cited only for confirmation, not as load-bearing justification for the central claims. No self-definitional loops, fitted inputs renamed as predictions, or ansatz smuggling appear in the derivation chain.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard assumptions of X-ray polarimetry and synchrotron emission models plus the validity of the new timing recovery; no new physical entities are introduced.

free parameters (2)

phase offsets and frequency derivatives
Jointly fitted via MCMC during GPS-off intervals to recover timing; these are data-driven parameters required for the phase-resolved polarization analysis.
polarization degree and angle
Fitted parameters for the nebular emission and phase bins; reported with uncertainties.

axioms (2)

domain assumption The Crab pulsar's 33 ms period is stable enough to serve as an external clock for phase reconstruction over multiple periods without GPS.
Invoked when using the pulsar to recover timing in GPS-off data.
domain assumption Hard X-ray emission in the Crab arises primarily from synchrotron processes in the nebular torus and wind regions.
Used to interpret the measured polarization alignment with the spin axis.

pith-pipeline@v0.9.0 · 5741 in / 1424 out tokens · 39414 ms · 2026-05-10T09:42:40.475818+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages

[1]

2021, Astroparticle Physics, 126, 102529, doi: 10.1016/j.astropartphys.2020.102529

Abarr, Q., Awaki, H., Baring, M., et al. 2021, Astroparticle Physics, 126, 102529, doi: 10.1016/j.astropartphys.2020.102529

work page doi:10.1016/j.astropartphys.2020.102529 2021
[2]

G., Chun, S., et al

Aoyagi, M., Bose, R. G., Chun, S., et al. 2024, Astroparticle Physics, 158, 102944, doi: 10.1016/j.astropartphys.2024.102944 16

work page doi:10.1016/j.astropartphys.2024.102944 2024
[3]

G., Bose, R., et al

Awaki, H., Baring, M. G., Bose, R., et al. 2025a, MNRAS, 540, L34, doi: 10.1093/mnrasl/slaf026

work page doi:10.1093/mnrasl/slaf026
[4]

G., Bose, R., et al

Awaki, H., Baring, M. G., Bose, R., et al. 2025b, ApJ, 994, 37, doi: 10.3847/1538-4357/ae0f1d

work page doi:10.3847/1538-4357/ae0f1d
[5]

2026, arXiv e-prints, arXiv:2602.09886, doi: 10.48550/arXiv.2602.09886

Bouchet, T., Laurent, P., Cangemi, F., & Rodriguez, J. 2026, arXiv e-prints, arXiv:2602.09886, doi: 10.48550/arXiv.2602.09886

work page doi:10.48550/arxiv.2602.09886 2026
[6]

2023, Nature Astronomy, 7, 602, doi: 10.1038/s41550-023-01936-8

Bucciantini, N., Ferrazzoli, R., Bachetti, M., et al. 2023, Nature Astronomy, 7, 602, doi: 10.1038/s41550-023-01936-8

work page doi:10.1038/s41550-023-01936-8 2023
[8]

2017, Scientific Reports, 7, 7816, doi: 10.1038/s41598-017-07390-7

Chauvin, M., Flor´ en, H.-G., Friis, M., et al. 2017, Scientific Reports, 7, 7816, doi: 10.1038/s41598-017-07390-7

work page doi:10.1038/s41598-017-07390-7 2017
[10]

2018, MNRAS, 477, L45, doi: 10.1093/mnrasl/sly027 de Jager, O

Chauvin, M., Flor´ en, H.-G., Friis, M., et al. 2018, MNRAS, 477, L45, doi: 10.1093/mnrasl/sly027 de Jager, O. C., Raubenheimer, B. C., & Swanepoel, J. W. H. 1989, A&A, 221, 180

work page doi:10.1093/mnrasl/sly027 2018
[11]

W., Lang, D., & Goodman, J

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306, doi: 10.1086/670067 Gonz´ alez-Caniulef, D., Heyl, J., Fabiani, S., et al. 2025, A&A, 693, A152, doi: 10.1051/0004-6361/202451815

work page doi:10.1086/670067 2013
[12]

K., & Kalapotharakos, C

Harding, A. K., & Kalapotharakos, C. 2017, ApJ, 840, 73, doi: 10.3847/1538-4357/aa6ead

work page doi:10.3847/1538-4357/aa6ead 2017
[13]

E., Lacey , C

Hobbs, G., Lyne, A. G., & Kramer, M. 2010, MNRAS, 402, 1027, doi: 10.1111/j.1365-2966.2009.15938.x

work page doi:10.1111/j.1365-2966.2009.15938.x 2010
[14]

2012, The Astrophysical Journal, 749, 2, doi: 10.1088/0004-637X/749/1/2

Contopoulos, I. 2012, The Astrophysical Journal, 749, 2, doi: 10.1088/0004-637X/749/1/2

work page doi:10.1088/0004-637x/749/1/2 2012
[15]

doi:10.1016/j.astropartphys.2015.02.007 , eprint =

Kislat, F., Clark, B., Beilicke, M., & Krawczynski, H. 2015, Astroparticle Physics, 68, 45, doi: https://doi.org/10.1016/j.astropartphys.2015.02.007

work page doi:10.1016/j.astropartphys.2015.02.007 2015
[16]

Li, H., Zhao, Q.-C., Feng, H., Tao, L., & Tsygankov, S. S. 2025, arXiv e-prints, arXiv:2504.04775, doi: 10.48550/arXiv.2504.04775

work page doi:10.48550/arxiv.2504.04775 2025
[17]

G., Jordan, C

Lyne, A. G., Jordan, C. A., Graham-Smith, F., et al. 2014, Monthly Notices of the Royal Astronomical Society, 446, 857, doi: 10.1093/mnras/stu2118

work page doi:10.1093/mnras/stu2118 2014
[18]

2018, A&A, 615, A54, doi: 10.1051/0004-6361/201731971

Mikhalev, V. 2018, A&A, 615, A54, doi: 10.1051/0004-6361/201731971

work page doi:10.1051/0004-6361/201731971 2018
[19]

2023, PASJ, 75, 1298, doi: 10.1093/pasj/psad070

Mizuno, T., Ohno, H., Watanabe, E., et al. 2023, PASJ, 75, 1298, doi: 10.1093/pasj/psad070

work page doi:10.1093/pasj/psad070 2023
[20]

Ng, C.-Y., & Romani, R. W. 2004, The Astrophysical Journal, 601, 479, doi: 10.1086/380486 P´ etri, J. 2013, MNRAS, 434, 2636, doi: 10.1093/mnras/stt1214

work page doi:10.1086/380486 2004
[21]

H., Jahoda, K., & Lyne, A

Rots, A. H., Jahoda, K., & Lyne, A. G. 2004, ApJL, 605, L129, doi: 10.1086/420842

work page doi:10.1086/420842 2004
[22]

N., & Harding , A

Timokhin, A. N., & Harding, A. K. 2019, ApJ, 871, 12, doi: 10.3847/1538-4357/aaf050

work page doi:10.3847/1538-4357/aaf050 2019
[23]

V., Chattopadhyay, T., Mithun, N

Vadawale, S. V., Chattopadhyay, T., Mithun, N. P. S., et al. 2018, Nature Astronomy, 2, 50, doi: 10.1038/s41550-017-0293-z

work page doi:10.1038/s41550-017-0293-z 2018
[24]

Wang, J., Verbiest, J. P. W., Shaifullah, G. M., et al. 2024, A&A, 687, A154, doi: 10.1051/0004-6361/202449366

work page doi:10.1051/0004-6361/202449366 2024
[25]

2025, Monthly Notices of the Royal Astronomical Society, 539, 902, doi: 10.1093/mnras/staf538

Wei, W., Xie, F., La Monaca, F., et al. 2025, Monthly Notices of the Royal Astronomical Society, 539, 902, doi: 10.1093/mnras/staf538

work page doi:10.1093/mnras/staf538 2025
[26]

C., Elsner, R

Weisskopf, M. C., Elsner, R. F., Hanna, D., et al. 2006, arXiv e-prints, astro, doi: 10.48550/arXiv.astro-ph/0611483

work page doi:10.48550/arxiv.astro-ph/0611483 2006
[27]

C., Silver, E

Weisskopf, M. C., Silver, E. H., Kestenbaum, H. L., Long, K. S., & Novick, R. 1978, ApJL, 220, L117, doi: 10.1086/182648

work page doi:10.1086/182648 1978
[28]

C., Hester, J

Weisskopf, M. C., Hester, J. J., Tennant, A. F., et al. 2000, The Astrophysical Journal Letters, 536, L81, doi: 10.1086/312733

work page doi:10.1086/312733 2000
[29]

B., & Fishman, G

Wilson, R. B., & Fishman, G. J. 1983, ApJ, 269, 273, doi: 10.1086/161039

work page doi:10.1086/161039 1983
[30]

2024, arXiv, doi: 10.48550/arxiv.2407.12779

Wong, J., Mizuno, T., Bucciantini, N., et al. 2024, arXiv, doi: 10.48550/arxiv.2407.12779

work page doi:10.48550/arxiv.2407.12779 2024
[31]

2024, ApJ, 973, 172, doi: 10.3847/1538-4357/ad6309

Wong, J., Mizuno, T., Bucciantini, N., et al. 2024, ApJ, 973, 172, doi: 10.3847/1538-4357/ad6309

work page doi:10.3847/1538-4357/ad6309 2024