Abstract:Context. We have investigated the multiwavelength emission of PSR J2021+4026, the only isolated gamma-ray pulsar known to be variable, which in October 2011 underwent a simultaneous change in gamma-ray flux and spin-down rate, followed by a second mode change in February 2018. Multiwavelength monitoring is crucial to understand the physics behind these events and how they may have affected the structure of the magnetosphere. Aims.The monitoring of pulse profile alignment is a powerful diagnostic tool for constraining magnetospheric reconfiguration. We aim to investigate timing or flux changes related to the variability of PSR J2021+4026 via multiwavelength observations, including gamma-ray observations from Fermi-LAT, X-ray observations from XMM-Newton, and a deep optical observation with the Gran Telescopio Canarias.Methods. We performed a detailed comparison of the timing features of the pulsar in gamma and X-rays and searched for any change in phase lag between the phaseogram peaks in these two energy bands. Although previous observations did not detect a counterpart in visible light, we also searched for optical emission that might have increased due to the mode change, making this pulsar detectable in the optical. Results.We have found a change in the gamma-to X-ray pulse profile alignment by 0.21$\pm$0.02 in phase, which indicates that the first mode change affected different regions of the pulsar magnetosphere. No optical counterpart was detected down to g'=26.1 and r'=25.3. Conclusions.We suggest that the observed phase shift could be related to a reconfiguration of the connection between the quadrupole magnetic field near the stellar surface and the dipole field that dominates at larger distances. This is consistent with the picture of X-ray emission coming from the heated polar cap and with the simultaneous flux and frequency derivative change observed during the mode changes.

Changes in the nature of the spectral continuum and the stability of the cyclotron line in the X-ray pulsar GRO J2058+42

Study of the X-ray Pulsar XTE J1946+274 with NuSTAR

Long-term spin-down and low luminosity regime in the Be/X-ray binary pulsar GX 304-1

RX J0440.9+4431: another supercritical X-ray pulsar

Drop in the hard pulsed fraction and a candidate cyclotron line in IGR J16320-4751 seen by NuSTAR

Temporal Variations And Spectral Properties of Be/X-ray Pulsar GRO J1008-57 Studied by INTEGRAL

Study of the X-ray Pulsar IGR J21343+4738 based on NuSTAR, Swift, and SRG data

Cyclotron Line Variability

NuSTAR detection of a broad absorption line in IGR J06074+2205

Detection of a cyclotron line in the Be X-ray pulsar IGR J06074$+$2205

Spectrum formation in X-ray pulsars at very low mass accretion rate: Monte-Carlo approach

On the cyclotron absorption line and evidence of the spectral transition in SMC X-2 during 2022 giant outburst

Evidence for a cyclotron absorption line and spectral transition in EXO 2030+375 during 2021 giant outburst

Investigation of the Timing and Spectral Properties of an Ultraluminous X-Ray Pulsar NGC 7793 P13

Temporal and spectral variations of the X-ray pulsar Cen X-3 observed by NuSTAR

X-ray emission from magnetized neutron star atmospheres at low mass accretion rates. I. Phase-averaged spectrum

Peculiarities of super-Eddington flares from the X-ray pulsar LMC X-4 with NuSTAR

The quiescent state of the accreting X-ray pulsar SAX J2103.5+4545

Multiwavelength observations of PSR J2021+4026 across a mode change reveal a phase shift in its X-ray emission