2025ApJ...993L..26Y

We present the long-term spectral and timing evolution of the ultraluminous X-ray pulsar NGC 7793 P13 from 2011 to 2024 based on archival data from XMM-Newton, Chandra, NICER, and NuSTAR, including unpublished data after 2020. This data set enables us to investigate the observational properties across a long modulation of ∼10 yr. Although previous studies suggested an increasing trend in flux in 2020, we find that the pulsar decreased its flux to $2.{6}_{-0.7}^{+0.7}\times 1{0}^{-14}\,{\rm{erg}}\,{{\rm{cm}}}^{-2}\,{{\rm{s}}}^{-1}$ in 2021, and rebrightened to $3.{9}_{-0.1}^{+0.1}\times 1{0}^{-12}\,{\rm{erg}}\,{{\rm{cm}}}^{-2}\,{{\rm{s}}}^{-1}$ by 2024. Moreover, in the last two years, the spin-up rate was twice as large as that before 2020. However, the pulsed fraction was roughly as expected from an anticorrelation with the flux confirmed before 2020. Furthermore, we perform systematic phase-resolved spectroscopy to investigate the spectral evolution. The results show that, during the flux-decaying epoch, spectral hardness in the off-pulse phase softened while that in its on-pulse phase remained almost unchanged. This softening was not observable in the rebrightening epoch. Our results could provide new insights into the accretion geometry of neutron stars in supercritical accretion systems.