SVOM discovery of a strong X-ray outburst of the blazar 1ES~1959+650 and multi-wavelength follow-up with the Neil Gehrels Swift observatory

TL;DR

SVOM detected the first X-ray outburst of a blazar, a 1ES 1959+650 event, and, together with Swift, conducted a dense multi-wavelength campaign to track the flare's temporal and spectral evolution. Spectral modeling with log-parabolic fits reveals a harder-when-brighter behavior and a synchrotron peak that shifts to higher energies as flux increases, while evidence for a single acceleration mechanism is inconclusive. The data favor a mixed acceleration scenario combining stochastic processes and possibly shock-related acceleration, as indicated by the positive S_p–E_p correlation and lack of a clear E_p–beta signature or HR–flux hysteresis. The study demonstrates SVOM–Swift synergy for time-domain blazar science and demonstrates the value of broad, continuous coverage from 0.3 to ~50 keV to constrain particle acceleration in jets.

Abstract

On December 6, 2024, 1ES 1959+650, one of the X-ray brightest blazars known, underwent a high-amplitude X-ray outburst detected by SVOM, the first such discovery with this mission. The source was subsequently monitored with SVOM and Swift from December 2024 to March 2025. We report the detection and multi-wavelength follow-up of this event, and describe the temporal and spectral evolution observed during the campaign. Data from SVOM/MXT, SVOM/ECLAIRs, and Swift/XRT were analyzed with log-parabola models to track flux and spectral variability. The source was detected in a bright state over the 0.3-50 keV range. During the three months of monitoring, the X-ray flux varied significantly, showing episodes of spectral hardening at high flux levels. The spectral curvature evolved more irregularly and did not show a clear trend with flux. A shift of the Spectral Energy Distribution (SED) synchrotron peak to higher energies is seen when the flux increases. This constitutes the first blazar outburst discovered in X-rays by SVOM. The coordinated follow-up with Swift provided continuous coverage of the flare and highlights the strong complementarity of the two missions for time-domain studies of blazars. The flare shows no clear signatures of either Fermi I or Fermi II acceleration, suggesting a mixed Fermi I/II scenario.

Figures (6)

• Figure 1: Evolution of the spectral parameters and the optical to X-ray fluxes. Are presented from top to bottom : intrinsic flux in the 0.3--10.0 keV band, the photon index, the curvature parameter, the energy of the SED synchrotron peak and the associated flux, the UV (W1, M2, W2 filters) and optical (V, B, U filters, as well as the VT R and B band) fluxes. The dashed blue line in the $\beta$ panel represent $\beta=0.4$, considered as a threshold value between low and high curvature. The red stars represent points that are unconstrained whereas red triangle represent $1\sigma$ upper limits. The blue and magenta backgrounds are respectively highlighting Period 1 and 2.
• Figure 2: Example of joint-fit of the XRT, MXT and ECLAIRs for the observations taken on December 16, 2024. The bottom panel shows the residuals, computed as the difference between the observed data ($D$) and the model ($M$), normalized by the error ($\sigma$).
• Figure 3: $E_p-\beta$ plan for the entire period (1+2). The purple and blue lines represent what we expect for Fermi II and Fermi I acceleration processes with an additional multiplicative constant fitted only using data from period 1. In both cases, the colored regions represent the error at the $1\sigma$ level. The red stars represent points that are unconstrained whereas red triangles represent $1\sigma$ upper limits. Errors on $E_{\rm p}$ at the $1\sigma$ level were obtained through error propagation using the $1\sigma$ uncertainties of the parameters $\alpha$ and $\beta$. When the lower bound of the confidence interval included zero, the value was treated as a $1\sigma$ upper limit. In addition, values were considered unconstrained when their uncertainties were at least one order of magnitude larger than the measured value itself, and when both limits of the confidence interval were poorly constrained.
• Figure 4: $S_p-E_p$ plan. The dashed line represents the case $S_p \propto E_p^{0.6}$. The orange and green lines represent the model $S_p = C \times E_p^q$ fitted to the data from period 1 and period 2, respectively. In both cases, the colored region represent the error at the $1\sigma$ level. The red stars represent points that are unconstrained whereas red triangles represent $1\sigma$ upper limits.
• Figure 5: $HR-F_{0.3-10 \ \mathrm{keV}}$ plan for the entire follow-up campaign (period 1 and 2). The number near the data points represent the order to follow, ranging from 0 to 30. The areas shaded in orange represent the three loops initially identified and finally reject by the bootstrap algorithm. The three loops are : 0-6 CW, 0-14 CW; 0-15 CW.