On the Cooling Trend of SGR 0526-66

We present a systematic analysis of all archival Chandra observations of the soft-gamma repeater SGR 0526-66. Our results show that the X-ray flux of SGR 0526-66 decayed by about 20% between 2000 and 2009. We employ physically motivated X-ray spectral models and determine the effective temperature and the strength of the magnetic field at the surface as kT = 0.354_{-0.024}^{+0.031} keV and B = (3.73^{+0.16}_{-0.08})x10^{14} G, respectively. We find that the effective temperature remains constant within the statistical uncertainties and attribute the decrease in the source flux to a decrease in the emitting radius. We also perform timing analysis to measure the evolution of the spin period and the period derivative over the nine year interval. We find a period derivative of .P = (4.0 +/- 0.5)x10^{-11} ss^{-1}, which allows us to infer the dipole magnetic field strength and compare it with the one determined spectroscopically. Finally, we compare the effective temperature of SGR 0526-66 with the expected cooling trends from magnetized neutron stars and suggest an initial magnetic field strength of 10^{15-16} G for the source.

💡 Research Summary

This paper presents a comprehensive study of the long‑term X‑ray variability and cooling behavior of the soft‑gamma repeater SGR 0526‑66, located in the Large Magellanic Cloud. Using all available Chandra ACIS‑S observations from 2000, 2001, and a series of four pointings in 2009, together with a deep XMM‑Newton EPIC‑pn exposure from 2007, the authors construct a uniform data set spanning nine years. Data reduction was performed with CIAO 4.2 and CALDB 4.3.0, applying barycentric corrections and extracting source events from a 2″ radius while using annular regions for background and for the surrounding supernova remnant N49.

A crucial step is the independent determination of the interstellar hydrogen column density (N_H) from the SNR spectrum. By fitting the SNR with two plane‑parallel shocked plasma components (vpshock), they obtain N_H = (0.15 ± 0.03) × 10²² cm⁻², which is then fixed in all subsequent fits of the magnetar spectrum. The magnetar’s X‑ray emission is modeled with the physically motivated STEMS (Surface Thermal Emission and Magnetospheric Scattering) model. This model incorporates a strongly magnetized hydrogen atmosphere, resonant cyclotron scattering in the magnetosphere, and four free parameters: effective surface temperature (kT), surface magnetic field strength (B), scattering optical depth (τ), and the average particle velocity (β).

Initial fits allowing all parameters to vary independently for each observation yielded large uncertainties due to limited counts. The authors therefore linked B, τ, and β across all epochs, allowing only kT and the model normalization (which translates to the emitting radius) to vary. The simultaneous fit of the six spectra gives a consistent magnetic field B = (3.73 +0.08 − 0.16) × 10¹⁴ G, τ ≈ 5.5, β ≈ 0.52, and an effective temperature kT = 0.355 +0.031 − 0.024 keV. The temperature remains statistically unchanged over the nine‑year interval.

The unabsorbed 0.5–6.5 keV flux declines from 1.33 × 10⁻¹² erg s⁻¹ cm⁻² in 2000 to ≈1.05 × 10⁻¹² erg s⁻¹ cm⁻² in 2009, a ~20 % decrease. Because the temperature is constant, the flux reduction is attributed to a shrinking emitting area: the inferred blackbody radius drops from ~13.3 km to ~11.8 km, implying a ~10 % contraction of the hot spot on the neutron‑star surface.

Timing analysis was performed using the Z²₂ test on the same Chandra data sets plus the XMM‑Newton observation. Pulsations near 8 s were detected with high significance in the early Chandra observations and the XMM‑Newton data; later Chandra pointings yielded marginal detections. Four reliable period measurements were obtained, and a linear fit yields a spin‑down rate ˙P = (4.02 ± 0.49) × 10⁻¹¹ s s⁻¹. Assuming magnetic dipole braking, this corresponds to a dipolar field B_dipole ≈ 3.9 × 10¹⁴ G, in excellent agreement with the spectroscopic B.

Finally, the authors compare the measured temperature and magnetic field with theoretical magnetar cooling curves that include field decay (e.g., Viganò et al. 2013). The observed point lies on the cooling trajectory expected for a neutron star born with an initial field of 10¹⁵–10¹⁶ G and an age of roughly 5 kyr, consistent with the Sedov age of the associated SNR N49 (~4.8 kyr). This suggests that SGR 0526‑66 has undergone significant magnetic field decay while maintaining a relatively stable surface temperature, a behavior characteristic of magnetars transitioning from an active bursting phase to a quiescent state.

In summary, the paper demonstrates that (1) the ~20 % X‑ray flux decline over nine years is driven by a reduction in the emitting radius rather than cooling, (2) the surface temperature and magnetic field have remained essentially constant, matching dipole‑braking estimates, and (3) the source’s thermal state aligns with magneto‑thermal evolution models that require an initially ultra‑strong field. The authors recommend continued high‑resolution X‑ray monitoring and more frequent timing observations to refine the spin‑down measurement and to track possible future changes in the emitting area or temperature.