Magnetar outburst models with cooling simulations

Magnetar outbursts are among the most noteworthy manifestations of magnetism in neutron stars. They are episodes in which the X-ray luminosity of a strongly magnetised neutron star swiftly rises by several orders of magnitude to then decay over the course of several months. In this work, we present simulations of outbursts as a consequence of localised heat deposition in a magnetised neutron star crust, and the subsequent surface cooling. In particular, we employed a magnetothermal evolution code adapted to the study of short-term phenomena; that is, one including in its integration domain the outer layers of the star, where heat diffusion is faster. This choice entailed the development and use of heat blanketing envelope models that are thinner than those found in the literature as the surface boundary condition. We find that such envelopes can support a higher surface temperature than the thicker ones (albeit for less time), which can account for the typical luminosities observed in outbursts even when coming from small hotspots (few km in radius). We study several parameters related to the energetics and geometry of the heating region, concluding that the cooling of a crustal hotspot found in the outer part of the crust can account for the luminosity evolution observed in outbursts both in terms of peak luminosity and timescales. Finally, we discuss the key observables that must be studied in future observations to better constrain the nature of the underlying mechanism.