Fermi Large Area Telescope and multi-wavelength observations of the flaring activity of PKS 1510-089 between 2008 September and 2009 June

We report on the multi-wavelength observations of PKS 1510-089 (a flat spectrum radio quasar at z=0.361) during its high activity period between 2008 September and 2009 June. During this 11 months period, the source was characterized by a complex variability at optical, UV and gamma-ray bands, on time scales down to 6-12 hours. The brightest gamma-ray isotropic luminosity, recorded on 2009 March 26, was ~ 2x10^48erg s^-1. The spectrum in the Fermi-LAT energy range shows a mild curvature well described by a log-parabolic law, and can be understood as due to the Klein-Nishina effect. The gamma-ray flux has a complex correlation with the other wavelengths. There is no correlation at all with the X-ray band, a weak correlation with the UV, and a significant correlation with the optical flux. The gamma-ray flux seems to lead the optical one by about 13 days. From the UV photometry we estimated a black hole mass of ~ 5.4x10^8 solar masses, and an accretion rate of ~ 0.5 solar masses/year. Although the power in the thermal and non-thermal outputs is smaller compared to the very luminous and distant flat spectrum radio quasars, PKS 1510-089 exhibits a quite large Compton dominance and a prominent big blue bump (BBB) as observed in the most powerful gamma-ray quasars. The BBB was still prominent during the historical maximum optical state in 2009 May, but the optical/UV spectral index was softer than in the quiescent state. This seems to indicate that the BBB was not completely dominated by the synchrotron emission during the highest optical state. We model the broadband spectrum assuming a leptonic scenario in which the inverse Compton emission is dominated by the scattering of soft photons produced externally to the jet. The resulting model-dependent jet energetic content is compatible with the accretion disk powering the jet, with a total efficiency within the Kerr black hole limit.

💡 Research Summary

This paper presents a comprehensive multi‑wavelength campaign on the flat‑spectrum radio quasar PKS 1510‑089 (z = 0.361) during an exceptionally active phase that lasted from September 2008 to June 2009. Using data from the Fermi Large Area Telescope (LAT) together with contemporaneous observations in the radio, optical, ultraviolet (UV), X‑ray, and γ‑ray bands, the authors investigate the source’s variability, spectral properties, and energetics with unprecedented temporal resolution.

The γ‑ray light curve, sampled on 6–12 hour intervals, shows extreme variability, with flux changes of more than an order of magnitude on timescales as short as half a day. The brightest γ‑ray episode occurred on 2009 March 26, when the isotropic γ‑ray luminosity reached ≈2 × 10⁴⁸ erg s⁻¹, placing PKS 1510‑089 among the most luminous γ‑ray flat‑spectrum radio quasars (FSRQs). The LAT spectrum in the 0.1–100 GeV range is best described by a log‑parabolic law (photon index ≈2.2, curvature parameter β≈0.08), indicating a mild curvature that the authors attribute to Klein‑Nishina (KN) effects in the inverse‑Compton (IC) scattering process.

Cross‑correlation analysis reveals a complex relationship between the γ‑ray band and lower‑energy emission. The γ‑ray flux shows no statistically significant correlation with the X‑ray (0.3–10 keV) flux, a weak positive correlation with the UV flux, and a strong positive correlation with the optical flux. Notably, the γ‑ray variations lead the optical variations by roughly 13 days, suggesting a scenario in which freshly accelerated electrons first produce high‑energy γ‑rays via external‑Compton (EC) scattering and later cool to emit synchrotron radiation in the optical band.

UV photometry, combined with a standard thin‑disk model, yields an estimate of the central black‑hole mass of M_BH ≈ 5.4 × 10⁸ M_⊙ and an accretion rate of Ṁ ≈ 0.5 M_⊙ yr⁻¹, corresponding to an Eddington ratio of ≈0.2. The spectral energy distribution (SED) displays a prominent big blue bump (BBB) that remains visible even during the optical flare of May 2009, although the optical/UV spectral index softens, indicating that the synchrotron component does not completely dominate the BBB during the highest optical state.

The authors model the broadband SED within a leptonic framework. The low‑energy component (radio to optical) is produced by synchrotron radiation of a relativistic electron population described by a log‑parabolic energy distribution (break Lorentz factor γ_break ≈ 10³, low‑energy slope s₁≈2.0, high‑energy slope s₂≈3.5). The high‑energy component is dominated by EC scattering of external photons originating from the accretion disk, the broad‑line region (BLR), and an infrared torus. The BLR, with a characteristic radius of ~0.1 pc and an energy density of ~10⁻² erg cm⁻³, provides the dominant seed photons for the γ‑ray production, consistent with the observed Compton dominance (UC/Usyn ≈ 10–20). The magnetic field is set to B ≈ 1 G, the emitting region size to R ≈ 3 × 10¹⁶ cm, and the Doppler factor to δ ≈ 20.

From the modeling, the total jet power (including kinetic, magnetic, and radiative contributions) is estimated to be ≈10⁴⁶ erg s⁻¹, roughly ten percent of the accretion‑disk luminosity. This power level is comfortably below the theoretical maximum set by the spin‑energy extraction limit of a Kerr black hole (≈0.3 Ṁc²), indicating that the jet can be powered by the accretion flow without violating energy‑budget constraints.

In summary, PKS 1510‑089, despite having a moderate accretion rate compared with the most luminous, distant FSRQs, exhibits a large Compton dominance and a conspicuous BBB, hallmarks of powerful γ‑ray quasars. The observed 13‑day lag between γ‑ray and optical emission, the log‑parabolic curvature of the LAT spectrum, and the successful EC‑dominated SED modeling together provide strong evidence that the high‑energy output is governed by external‑photon fields and that Klein‑Nishina effects shape the γ‑ray spectrum. The study underscores the importance of coordinated, high‑cadence multi‑wavelength monitoring for disentangling the complex interplay between jet dynamics, particle acceleration, and the surrounding radiation fields in blazars.