-
Modeling blazar broadband emission with convolutional neural networks -- II. External Compton model
Authors:
N. Sahakyan,
D. Bégué,
A. Casotto,
H. Dereli-Bégué,
P. Giommi,
S. Gasparyan,
V. Vardanyan,
M. Khachatryan,
A. Pe'er
Abstract:
In the context of modeling spectral energy distributions (SEDs) for blazars, we extend the method that uses a convolutional neural network (CNN) to include external inverse Compton processes. The model assumes that relativistic electrons within the emitting region can interact and up-scatter external photon originating from the accretion disk, the broad-line region, and the torus, to produce the o…
▽ More
In the context of modeling spectral energy distributions (SEDs) for blazars, we extend the method that uses a convolutional neural network (CNN) to include external inverse Compton processes. The model assumes that relativistic electrons within the emitting region can interact and up-scatter external photon originating from the accretion disk, the broad-line region, and the torus, to produce the observed high-energy emission. We trained the CNN on a numerical model that accounts for the injection of electrons, their self-consistent cooling, and pair creation-annihilation processes, considering both internal and all external photon fields. Despite the larger number of parameters compared to the synchrotron self-Compton model and the greater diversity in spectral shapes, the CNN enables an accurate computation of the SED for a specified set of parameters. The performance of the CNN is demonstrated by fitting the SED of two flat-spectrum radio quasars, namely 3C 454.3 and CTA 102, and obtaining their parameter posterior distributions. For the first source, the available data in the low-energy band allowed us to constrain the minimum Lorentz factor of the electrons, $γ_{\rm min}$, while for the second source, due to the lack of these data, $γ_{\rm min} = 10^2$ was set. We used the obtained parameters to investigate the energetics of the system. The model developed here, along with one from Bégué et al. (2023), enables self-consistent, in-depth modeling of blazar broadband emissions within leptonic scenario.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Modeling blazar broadband emission with convolutional neural networks -- I. Synchrotron self-Compton model
Authors:
Damien Bégué,
Narek Sahakyan,
Hüsne Dereli Bégué,
Paolo Giommi,
Sargis Gasparyan,
Mher Khachatryan,
Andrea Casotto,
Asaf Pe'er
Abstract:
Modeling the multiwavelength spectral energy distributions (SEDs) of blazars provides key insights into the underlying physical processes responsible for the emission. While SED modeling with self-consistent models is computationally demanding, it is essential for a comprehensive understanding of these astrophysical objects. We introduce a novel, efficient method for modeling the SEDs of blazars b…
▽ More
Modeling the multiwavelength spectral energy distributions (SEDs) of blazars provides key insights into the underlying physical processes responsible for the emission. While SED modeling with self-consistent models is computationally demanding, it is essential for a comprehensive understanding of these astrophysical objects. We introduce a novel, efficient method for modeling the SEDs of blazars by the mean of a convolutional neural network (CNN). In this paper, we trained the CNN on a leptonic model that incorporates synchrotron and inverse Compton emissions, as well as self-consistent electron cooling and pair creation-annihilation processes. The CNN is capable of reproducing the radiative signatures of blazars with high accuracy. This approach significantly reduces computational time, thereby enabling real-time fitting to multi-wavelength datasets. As a demonstration, we used the trained CNN with MultiNest to fit the broadband SEDs of Mrk 421 and 1ES 1959+650, successfully obtaining their parameter posterior distributions. This novel framework for fitting the SEDs of blazars will be further extended to incorporate more sophisticated models based on external Compton and hadronic scenarios, allowing for multi-messenger constraints in the analysis. The models will be made publicly available via a web interface, the Markarian Multiwavelength Datacenter, to facilitate self-consistent modeling of multi-messenger data from blazar observations.
△ Less
Submitted 6 November, 2023;
originally announced November 2023.
-
Gradient boosting decision trees classification of blazars of uncertain type in the fourth Fermi-LAT catalog
Authors:
N. Sahakyan,
V. Vardanyan,
M. Khachatryan
Abstract:
The deepest all-sky survey available in the $γ$-ray band - the last release of the Fermi-LAT catalogue (4FGL-DR3) based on the data accumulated in 12 years, contains more than 6600 sources. The largest population among the sources is blazar subclass - 3743, $60.1\%$ of which are classified as BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs), while the rest are listed as blazar…
▽ More
The deepest all-sky survey available in the $γ$-ray band - the last release of the Fermi-LAT catalogue (4FGL-DR3) based on the data accumulated in 12 years, contains more than 6600 sources. The largest population among the sources is blazar subclass - 3743, $60.1\%$ of which are classified as BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs), while the rest are listed as blazar candidates of uncertain type (BCU) as their firm optical classification is lacking. The goal of this study is to classify BCUs using different machine learning algorithms which are trained on the spectral and temporal properties of already classified BL Lacs and FSRQs. Artificial Neural Networks, \textit{XGBoost} and LightGBM algorithms are employed to construct predictive models for BCU classification. Using 18 input parameters of 2219 BL Lacs and FSRQs, we train (80\% of the sample) and test (20\%) these algorithms and find that LightGBM model, state-of-the-art classification algorithm based on gradient boosting decision trees, provides the highest performance. Based on our best model, we classify 825 BCUs as BL Lac candidates and 405 as FSRQ candidates, however, 190 remain without a clear prediction but the percentage of BCUs in 4FGL is reduced to 5.1\%. The $γ$-ray photon index, synchrotron peak frequency, and high energy peak frequency of a large sample are used to investigate the relationship between FSRQs and BL Lacs (LBLs, IBLs, and HBLs).
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
Modeling the time variable spectral energy distribution of the blazar CTA 102 from 2008 to 2022
Authors:
N. Sahakyan,
D. Israyelyan,
G. Harutyunyan,
S. Gasparyan,
V. Vardanyan,
M. Khachatryan
Abstract:
We present long-term multiwavelength observations of blazar CTA 102 ($z=1.037$). Detailed temporal and spectral analyses of $γ$-ray, X-ray and UV/optical data observed by {\it Fermi}-LAT, Swift XRT, NuSTAR and Swift-UVOT over a period of 14 years, between August 2008 and March 2022, was performed. We found strong variability of source emission in all the considered bands, especially in the $γ$-ray…
▽ More
We present long-term multiwavelength observations of blazar CTA 102 ($z=1.037$). Detailed temporal and spectral analyses of $γ$-ray, X-ray and UV/optical data observed by {\it Fermi}-LAT, Swift XRT, NuSTAR and Swift-UVOT over a period of 14 years, between August 2008 and March 2022, was performed. We found strong variability of source emission in all the considered bands, especially in the $γ$-ray band it exhibited extreme outbursts when the flux crossed the level of $10^{-5}\:{\rm photon\:cm^{-2}\:s^{-1}}$. Using the Bayesian Blocks algorithm, we split the adaptively binned $γ$-ray light curve into 347 intervals of quiescent and flaring episodes and for each period built corresponding multiwavelength spectral energy distributions (SEDs), using the available data. Among the considered SEDs, 117 high-quality (quasi) contemporaneous SEDs which have sufficient multiwavelength data, were modeled using JetSeT framework within a one-zone leptonic synchrotron and inverse Compton emission scenario assuming the emitting region is within the broad-line-region and considering internal and external seed photons for the inverse Compton up-scattering. As a result of modeling, the characteristics of the relativistic electron distribution in the jet as well as jet properties are retrieved and their variation in time is investigated. The applied model can adequately explain the assembled SEDs and the modelling shows that the data in the bright flaring periods can be reproduced for high Doppler boosting and magnetic field. The obtained results are discussed in the context of particle cooling in the emitting region.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Multiwavelength Study of High-Redshift Blazars
Authors:
N. Sahakyan,
D. Israyelyan,
G. Harutyunyan,
M. Khachatryan,
S. Gasparyan
Abstract:
High-redshift blazars are among the most powerful objects in the Universe. The spectral and temporal properties of thirty-three distant blazars ($z>2.5$) detected in the high energy $γ$-ray band are investigated by analyzing the {\it Fermi}-LAT and {\it Swift} UVOT/ XRT data. The considered sources have soft time averaged $γ$-ray spectra ($Γ_{\rm γ}\geq2.2$) whereas those that have been observed i…
▽ More
High-redshift blazars are among the most powerful objects in the Universe. The spectral and temporal properties of thirty-three distant blazars ($z>2.5$) detected in the high energy $γ$-ray band are investigated by analyzing the {\it Fermi}-LAT and {\it Swift} UVOT/ XRT data. The considered sources have soft time averaged $γ$-ray spectra ($Γ_{\rm γ}\geq2.2$) whereas those that have been observed in the X-ray band have hard X-ray spectra ($Γ_{\rm X}=1.01-1.86$). The $γ$-ray flux of high-redshift blazars ranges from $4.84\times10^{-10}$ to $1.50\times10^{-7}\:{\rm photon\:cm^{-2}\:s^{-1}}$ and the luminosity is within $(0.10-5.54)\times10^{48}\:{\rm erg\:s^{-1}}$ which during the $γ$-ray flares increases up to $(0.1-1)\times10^{50}\:{\rm erg\:s^{-1}}$. In the X-ray band, only the emission of PKS 0438-43, B2 0743+25 and TXS 0222+185 is found to vary in different Swift XRT observations whereas in the $γ$-ray band, the emission is variable for fourteen sources: the flux of B3 1343+451 and PKS 0537-286 changes in sub-day scales, that of PKS 0347-211 and PKS 0451-28 in day scales, while the$γ$-ray variability of the others is in week or month scales. The properties of distant blazar jets are derived by modeling the multiwavelength spectral energy distributions within a one-zone leptonic scenario assuming that the X-ray and $γ$-ray emissions are produced from inverse Compton scattering of synchrotron and dusty torus photons. From the fitting, the emission region size is found to be $\leq0.05$ pc and the magnetic field and the Doppler factor are correspondingly within $0.10-1.74$ G and $10.0-27.4$. By modeling the optical-UV excess, we found that the central black hole masses and accretion disk luminosities are within $L_{\rm d}\simeq(1.09-10.94)\times10^{46}\:{\rm erg \: s^{-1}}$ and $(1.69-5.35)\times10^{9}\:M_\odot$, respectively.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.