A phase-resolved analysis of the mode-changing pulsar PSR J2021+4026 shows hints of a multipolar magnetosphere

A. Fiori,M. Razzano,A. K. Harding,M. Kerr,R. P. Mignani,P. M. Saz Parkinson

DOI: https://doi.org/10.1051/0004-6361/202348924

2024-02-26

Abstract:The radio-quiet pulsar PSR J2021+4026 is a peculiar pulsar showing repeated and quasi-periodic mode changes. Its flux shows repeated variations between two states at intervals of 3.5$ years. These events occur over timescales $<$100 days and are correlated with sudden changes in the spin-down rate. Multiwavelength observations also revealed an X-ray phase shift relative to the profile for one of the events. PSR J2021+4026 is currently the only known isolated pulsar showing significant variability, and thus it has been the object of thorough investigations. The goal of our work is to study the mode changes of PSR J2021+4026 with improved detail. By accurately characterizing variations in the spectrum and pulse profile, we aim to relate the observations to theoretical models. We also aim to interpret the mode changes in terms of variations in the structure of a multipolar dissipative magnetosphere. We continually monitored the rotational evolution and the flux of PSR J2021+4026 using more than 13 years of data with a binned likelihood approach. We investigated the features of the phase-resolved spectrum and pulse profile, and from these we inferred the macroscopic conductivity, the electric field parallel to the magnetic field and the curvature radiation cutoff energy. These physical quantities are related to the spin-down rate and the flux and therefore are relevant to the theoretical interpretation of the mode changes. We introduced a simple magnetosphere model that combines a dipole field with a strong quadrupole component. We simulated magnetic field configurations to determine the positions of the polar caps for different sets of parameters. We clearly detect the previous mode changes and confirm a more recent mode change that occurred around June 2020. We provide a full set of best-fit parameters for the phase-resolved spectrum and the pulse profile obtained in five distinct time intervals. We computed the relative variations in the best-fit parameters, finding typical flux changes between $13$ and $20$ . Correlations appear between the flux and the spectral parameters as the peak of the spectrum shifts by sim 10 toward lower energies when the flux decreases. The analysis of the pulse profile reveals that the pulsed fraction of the light curve is larger when the flux is low. Finally, the magnetosphere simulations show that some configurations could explain the observed multiwavelength variability. However, self-consistent models are required to reproduce the observed magnitudes of the mode changes.

astronomy & astrophysics