Modelling the Time-variable Broadband Emission and Correlation Study of FSRQ S5 1044+71

TL;DR

This study analyzes 16 years of Fermi-LAT and Swift data for the FSRQ S5 1044+71 to map time-variable broadband emission and interband correlations. By applying Bayesian blocks to 3-day γ-ray light curves, the authors identify four major outbursts and measure a shortest variability timescale of 4.5 hours, constraining the γ-ray emission region to within about 0.03 pc of the central engine, near the BLR. Cross-correlation reveals strong γ-ray–X-ray coupling with X-rays lagging by ~42.5 days and a modest γ-ray–optical/UV lag, while SED modeling with a one-zone leptonic scenario requires EC scattering from both IR and BLR photon fields, with high-flux states showing harder electron spectra and reduced magnetic fields, consistent with enhanced particle acceleration and Compton dominance. The results support a compact, BLR-adjacent emission region and a leptonic origin for the γ-ray emission, with implications for jet physics and emission-site localization in powerful FSRQs. The work highlights the value of long-term, multiwavelength monitoring and time-resolved SED fitting for disentangling variability processes in blazar jets and motivates future VHE observations to further test the proposed emission geometry.

Abstract

We present a detailed temporal and spectral analysis of the blazar S5\,1044+71 using multi-wavelength data obtained from the \emph{Fermi}-LAT and Swift-XRT/UVOT telescopes. Applying the Bayesian block algorithm to the 3-day binned $γ$-ray lightcurve, we identify pronounced variability, including four major outbursts marked by significant flux enhancements. The highest flux recorded was $(1.1 \pm 0.2)\times 10^{-6}\,\text{ph}\,\text{cm}^{-2}\,\text{s}^{-1}$ on 57868.5 MJD. Each outburst comprises multiple components, and lightcurve profile analysis indicates predominantly symmetric temporal structures. The shortest variability timescale of 4.5 hours constrains the emission region to be located within 0.03 pc of the central engine, likely near the broad-line region (BLR). Additionally, two highest-energy photons were detected with energies of 46.4 GeV (on 57739.6 MJD) and 42.5 GeV (on 59161.9 MJD), observed outside the peak flaring activity. The fractional variability shows an overall increasing trend from UV/optical to $γ$-ray bands, with a noticeable dip in the X-ray range, consistent with the shape of the broadband spectral energy distribution (SED). The flux distributions during flares exhibit log-normal or double log-normal behavior, suggesting multiplicative variability processes and evolving emission zones. Cross-correlation analysis reveals a strong positive correlation between the $γ$-ray and X-ray bands, with X-rays lagging by 42.5 days. Broadband SED modeling across different flux states supports a one-zone leptonic scenario, with $γ$-ray emission produced via external Compton scattering of IR and BLR photons. High-flux states show harder electron spectra, elevated break energies, and reduced magnetic fields-features consistent with efficient particle acceleration and Compton dominance.

Figures (8)

• Figure 1: (a) Fermi-LAT lightcurve of S5 1044+71 from August 2008 to February 2024 at flux $\text{F}_{\gamma}$ (E $>$100 MeV), with 3-day binning, in units of $\text{ph}\,\text{cm}^{-2}\,\text{s}^{-1}$. The grey vertical strips indicate the time intervals selected for the flux distribution study. (b) $\gamma$-ray spectral-index ($\Gamma_{\gamma}$) as a function of time, with the horizontal dashed purple line representing the average value of index. (c) TS values ($>$ 4) for each time bin. (d) Arrival times and energies of $\gamma$-ray photons E$_{\gamma}$, in units of GeV, with significance level above $3\sigma$. All photons above yellow dashed line have E$_{\gamma}>20$ GeV.
• Figure 2: 3-day binned $\gamma$-ray lightcurves of S5 1044+71 fitted with the SOE function defined in equation (\ref{['eq:SOE']}).
• Figure 3: Flux distribution of S5 1044+71 in the $\gamma$-ray band. Epochs S1, S2, and S4 are fitted with a Gaussian PDF, whereas epoch S3 and the entire 3-day binned $\gamma$-ray lightcurve are fitted with a double Gaussian PDF, respectively. The flux are in units of $\text{ph}\,\text{cm}^{-2}\,\text{s}^{-1}$.
• Figure 4: Multi-wavelength lightcurves (MLCs) of S5 1044+71 obtained by using Fermi-LAT and Swift-XRT/UVOT observations. The observations spanned a period from 54702.5 – 60344.5 MJD. Top panel is 7-day binned $\gamma$-ray lightcurve, second, third, and fourth panels are the X-ray, UV, and optical band lightcurves, respectively. The grey vertical stripes indicate the regions where broadband spectral modelling is performed, while the region between the red vertical lines marks the time interval used for cross-correlation analysis among different bands.
• Figure 5: Fractional variability amplitude in the $\gamma$-ray, X-ray, optical, and UV bands plotted as a function of energy over the period 54702.5 – 60344.5 MJD.