Multiwavelength spectral and temporal analysis of VHE Blazar 1ES 1959+650: Tracing emission mechanisms across flux states

TL;DR

The paper addresses how emission mechanisms in the high-energy peak BL Lac 1ES 1959+650 operate across flux states during a period of enhanced activity, using coordinated Fermi-LAT, Swift-XRT/UVOT, and LHAASO data. It employs a one-zone leptonic SSC framework to fit state-resolved SEDs, finding systematic trends: higher flux corresponds to harder electron spectra, stronger Doppler boosting, and weaker magnetic fields, with the break energy shifting to higher values as activity increases. The γ-ray–X-ray correlation and energy-dependent variability support a leptonic, acceleration-and-cooling driven origin, while VHE data reveal that a single-zone SSC can explain most states but not the brightest VHE state (VHE–FX1), suggesting potential multi-zone or hadronic contributions. These results provide quantitative constraints on jet dynamics, particle acceleration efficiency, and energy partition in BL Lac jets, highlighting the role of shock-driven processes and Doppler boosting in producing rapid, broadband flares with implications for modeling, jet structure, and energy budgets in HBLs.

Abstract

The high-synchrotron-peaked BL Lac object 1ES\,1959+650 exhibited pronounced activity between MJD~60310 -- 60603, including a very high energy (VHE) detection reported by LHAASO. To investigate the underlying emission mechanisms, we performed a comprehensive temporal and spectral analysis using multiwavelength data from \textit{Swift}-XRT/UVOT and \textit{Fermi}-LAT, covering the optical/UV to GeV $γ$-ray bands. The source shows strong energy-dependent variability, with the largest fractional variability in $γ$-rays, followed by X-rays and UV/optical, consistent with leptonic emission scenarios. Based on the variability patterns, we identified distinct flux states (F1, F2, F3, F4, F5, VHE-FX1, and VHE-FX2). The X-ray spectra exhibit a clear ``harder-when-brighter'' trend across these states. We modeled the broadband spectral energy distributions (SEDs) using a one-zone model incorporating synchrotron and synchrotron self-Compton (SSC) emission, implemented in \textsc{xspec} using $χ^{2}$ minimization. During the VHE detection, the corresponding X-ray/optical emission likely resembled the F2 state. Modeling the VHE SED using F1-state data led to an SSC overprediction of the VHE flux, whereas all other states were well described within the one-zone framework. Systematic trends in physical parameters are observed across flux states, including spectral hardening, increasing break energy, rising bulk Lorentz factor, and decreasing magnetic field with increasing flux. These results suggest that enhanced particle acceleration efficiency and stronger Doppler boosting drive the observed flaring activity, while the decrease in magnetic field indicates conversion of magnetic energy into particle kinetic energy, consistent with shock-driven scenarios.

Figures (3)

• Figure 1: Multi-wavelength light curve of 1ES 1959+650 in different flux states. The top panel displays the 3-day binned light curve integrated over the energy range 0.1–300  GeV, including data points with $|\mathrm{TS}| > 4$, the upper middle panel displays the X-ray light curve in the energy 0.3-10 keV, and the lower panel and the bottom panel display UV and Optical light curves, respectively.
• Figure 2: Energy-dependent fractional variability in different energy bands.
• Figure 3: Broadband SEDs of 1ES 1959+650 obtained during different flux states. Top left panel: F1 state; top right panel: F2 state; upper middle left panel: VHE--FX1 state; upper middle right panel: VHE--FX2 state; lower middle left panel: F3 state; lower middle right panel: F4 state; bottom panel: F5 state. In all panels, the flux points denoted by diamond markers correspond to Swift-UVOT, circular markers to Swift-XRT, blue square markers to Fermi-LAT, and pink square markers in the VHE--FX1 and VHE--FX2 panels indicate the Crab flux point (0.5 Crab units above 1 TeV). The solid red curve represents the best-fit model spectrum.