Next Article in Journal
Lorentz Violation by the Preferred Frame Effects and Cosmic and Gamma Ray Propagation
Next Article in Special Issue
Promise of Persistent Multi-Messenger Astronomy with the Blazar OJ 287
Previous Article in Journal
A GMRT Narrowband vs. Wideband Analysis of the ACT−CL J0034.4+0225 Field Selected from the ACTPol Cluster Sample
Previous Article in Special Issue
The Radiative Newtonian 1 < γ ≤ 1.66 and the Paczyński–Wiita γ = 5/3 Regime of Non-Isothermal Bondi Accretion onto a Massive Black Hole with an Accretion Disc
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

γ-ray Flux and Spectral Variability of Blazar Ton 599 during Its 2021 Flare

Indian Institute of Astrophysics, Block II Koramangala, Bangalore 560034, India
*
Author to whom correspondence should be addressed.
Submission received: 15 November 2021 / Revised: 9 December 2021 / Accepted: 10 December 2021 / Published: 13 December 2021
(This article belongs to the Special Issue Global Understanding of Accretion and Ejection around Black Holes)

Abstract

:
Blazars are known to emit exceptionally variable non-thermal emission over the wide range (from radio to γ -rays) of electromagnetic spectrum. We present here the results of our γ -ray flux and spectral variability study of the blazar Ton 599, which has been recently observed in the γ -ray flaring state. Using 0.1–300 GeV γ -ray data from the Fermi Gamma-ray Space Telescope (hereinafter Fermi ), we generated one-day binned light curve of Ton 599 for a period of about one-year from MJD 59,093 to MJD 59,457. During this one year period, the maximum γ -ray flux detected was 2.24 ± 0.25 × 10 6 ph cm 2 s 1 at MJD 59,399.50. We identified three different flux states, namely, epoch A (quiescent), epoch B (pre-flare) and epoch C (main-flare). For each epoch, we calculated the γ -ray flux variability amplitude (F v a r ) and found that the source showed largest flux variations in epoch C with F v a r 35%. We modelled the γ -ray spectra for each epoch and found that the Log-parabola model adequately describes the γ -ray spectra for all the three epochs. We estimated the size of the γ -ray emitting region as 1.03 × 10 16 cm and determined that the origin of γ -ray radiation, during the main-flare, could be outside of the broad line region.

1. Introduction

Blazars are the jetted subclass of the Active Galactic Nuclei (AGNs) that are understood as AGNs with very small viewing angles to the line of sight [1,2]. Blazars emit broadband electromagnetic radiation that ranges from radio to extremely high γ -ray energies. The jets of blazars are highly Doppler boosted, resulting in flux variations over the entire accessible electromagnetic wavebands. Flat spectrum radio quasars (FSRQs) and BL Lacertae (BL Lacs) objects are the two subclass of blazars. The difference between these subclasses is determined by the equivalent width (EW) of the emission lines in their optical spectra, with FSRQs having EW > 5 Å and BL Lacs having EW < 5 Å [3]. A more physical criterion for distinguishing between FSRQs and BL Lacs was proposed by [4] which is based on the ratio of broad line region (BLR) luminosity ( L B L R ) to Eddington Luminosity ( L E d d ). FSRQs have the value L B L R / L E d d > 5 × 10 4 , whereas BL Lacs have the value L B L R / L E d d < 5 × 10 4 .
The broad-band spectral energy distributions (SEDs) of blazars comprise of two humps; the low energy hump and the high energy hump. The low energy hump peaks at optical/UV/X-ray region and the high energy hump peaks at MeV/GeV/TeV region [5,6,7]. The genesis of the low energy hump is well understood through synchrotron emission mechanism of the relativistic electrons, while the high energy hump originates through inverse Compton (IC) emission process [6]. Based on the position of the synchrotron peak frequency ( ν s y n ), blazars are further classified as low synchrotron peaked (LSP; ν s y n < 10 14 Hz), intermediate synchrotron peaked (ISP; 10 14 Hz < ν s y n < 10 15 Hz) and high synchrotron peaked (HSP; ν s y n > 10 15 Hz) blazars.
Due to the Doppler boosting in blazars’ jets, the observed emission S o b s relative to the emission in the co-moving frame S i n t is defined as e.g., [8]
S o b s = S i n t δ q
where q = 3 + α for a moving compact source and q = 2 + α for a stationary jet [9]. α is the spectral index which is defined as f ν ν α . δ is the Doppler boosting factor which is described as δ = 1/ Γ ( 1 β cos θ ), where Γ is bulk Lorentz factor ( Γ = 1/ ( 1 β 2 ) 1 / 2 ), β is the speed of jet in units of the speed of light and θ is the viewing angle between observer’s line of sight and jet’s axis. The observed time is also shortened by the effect of Doppler boosting by a factor δ 1 . These two consequences of the Doppler boosting increase the chances of observed variations in blazars over a wide range of wavelengths and make them the brightest objects in the extragalactic sky.
The study of flux variability is a valuable tool for determining the size of the emission zone in blazars. Blazars exhibit variability over a wide range of timescales, from few minutes to several years, across the full wavelength range i.e., from radio to γ -rays [10,11,12,13,14]. The variability in blazars can be explained by the shock-in-jet model [15]. In this scenario, the inhomogeneities in the jet flow produce relativistic shocks. These shocks travel through the jet plasma at relativistic speed and accelerate the particles. The other scenario is the magnetic reconnection, which is responsible for the rapid variations in blazar jets and has been investigated extensively in recent years [16,17,18,19]. Within the reconnection region, magnetic reconnection possibly enables small-scale ultra-relativistic flow, also called jet-in-jet scenario [16,19,20]. The ultra-relativistic motion of small plasmoids causes additional Doppler boosting and results in very short and bright flares.
The launch of Fermi Telescope in 2008 gave us an unprecedented opportunity to explore the γ -ray regime in blazars. Using Fermi data at γ -ray energies, flux doubling time scales have been detected as short as few minutes in both FSRQs and BL Lacs [14,21,22,23]. Causality considerations in these circumstances point to highly compact γ -ray emission zones in blazars’ jets. The detection of γ -ray photons with energy greater than 10 GeV in FSRQs leads to the conclusion that the γ -ray emission region should be located at outside the cavity formed by the broad emission line (BLR) [24]. Otherwise, pair production on UV photons emitted by (BLR) clouds should severely attenuate the γ -ray photons. Despite several studies, the physics of γ -ray variability in blazars is still a captivating subject of research.
Ton 599 (4FGL J1159.5 + 2914; [25]) is an FSRQ, located at redshift z = 0.725 [26] with R.A. = 11 h 59 m 31.8 s and Dec. = + 29 14 43.8 . Ton 599 is a strongly polarized and a highly optically violent variable quasar [27]. It was first detected in high energy γ -rays by the Energetic Gamma Ray Experiment Telescope (EGRET) in the second EGRET catalogue [28] and later by the Fermi Large Area Telescope ( Fermi -LAT) [29]. The correlation study between radio and γ -ray bands was carried out by [30] for this source to put the constraints on γ -ray emission region in parsec-scale jets. In 2017, for the first time, Ton 599 went through a protracted flaring condition spanning the full electromagnetic spectrum. During this flare, the maximum γ -ray flux observed was 1.26 × 10 6 ph cm 2 s 1 [31]. The detailed study of γ -ray flux variability during this flare was carried out by [31,32]. A multiwavelength study and the broad-band SED modelling of this flare was conducted by [33] using a two-component leptonic emission model. An EC mechanism with a dusty torus (DT) photon field producing seed photons was identified to be responsible for the GeV emission in this study. Ton 599 recently displayed a bright flare from July to September 2021, allowing us to investigate its γ -ray emission process and its consequences during this source’s flaring condition.
In this paper, we describe the γ -ray analysis performed on the source Ton 599 utilising data collected over a one-year period from September 2020 to August 2021, with the goal of constraining the γ -ray emission region in blazars. The data used in this study is presented in Section 2. Section 3 describes the γ -ray light curve, followed by results in Section 4 and discussion in Section 5. Section 6 presents a summary of the results.

2. γ -ray Observations and Data Reduction

In this work, we used the γ -ray observations of Ton 599 taken by the Fermi -LAT for a period of ∼1 year, from 2020 September to 2021 August (MJD: 59,093–59,457). Fermi -LAT is a pair-conversion γ -ray telescope that can detect γ rays with energies ranging from 20 MeV to more than 1 TeV. It has a large field of view of about 2.4 sr and scans the full sky in every 3 h, covering 20% of the sky at any time [34]. We used the package ScienceTools v1.2.23 with the instrument response function P8R3_SOURCE_V31 for our analysis. We used the latest LAT Pass 8 data in the energy range 100 MeV to 300 GeV, where the photon-like events are classified as ‘evclass = 128, evtype = 3’. The region of interest (ROI) is specified as a circle with a radius of 10 and is centred on the γ -ray position of the source. We used a maximum zenith angle of 90 to remove γ -ray contamination from the earth’s limb. The latest isotropic model ‘iso_P8R3_SOURCE_V3_v1’ and the Galactic diffuse emission model ‘gll_iem_v07’ were used to analyze the data. The recommended condition ‘(DATA QUAL > 0)&&(LAT CONFIG= = 1)’ was then used to construct the required good time intervals. Over the time period of interest, an unbinned likelihood analysis is performed to generate the 1-day binned γ -ray light curve of Ton 599. In the light curve, the source was considered to be detected if the test statistics (TS) > 9, which corresponds to a 3 σ detection [35]. Our final one-day binned γ -ray light curve contains 256 confirmed measurements of Ton 599.

3. γ -ray Light Curve

The one-day binned γ -ray light curve of Ton 599 from 2020 September to 2021 August (MJD: 59,093–59,457) is shown in Figure 1. For the entire period of the light curve, we estimated the average γ -ray flux which was found to be 0.43 ± 0.36 × 10 6 ph cm 2 s 1 . Based on the average γ -ray flux, we visually identified three different flux states, namely, quiescent state, pre-flaring state, and the main flaring state, which were labelled as epoch A, epoch B, and epoch C, respectively. The duration of these epochs are marked by vertical lines in Figure 1 and the details of these epochs are given in Table 1. We classified epoch A as a quiescent epoch since the flux was lower than the average flux for the whole time span. During epochs B and C, the flux increased 2–3 times than the average flux, so we classified these epochs as flaring epochs. Over the whole time span of the light curve, the value of the flux was maximum during the flaring epoch C. The value of the maximum flux was determined to be 2.24 ± 0.25 × 10 6 ph cm 2 s 1 at MJD 59,399.50.

4. Results

4.1. γ -ray Flux Variability

We used the fractional variability amplitude ( F v a r ) to characterize the flux variability properties of Ton 599. The F v a r is commonly used to quantify the intrinsic variations in blazar light curves and is defined as, e.g., [12,36,37]
F v a r = S 2 σ e r r 2 ¯ x ¯ 2
In the above equation S 2 is the sample variance and σ e r r 2 ¯ is mean square error. These are defined as
S 2 = 1 N 1 i = 1 N ( x i x ¯ ) 2
and
σ e r r 2 ¯ = 1 N i = 1 N σ e r r , i 2
Here N is the total number of data points in the light curve and x ¯ is the arithmetic mean of the light curve.
The uncertainty in F v a r is given by the following equation
e r r ( F v a r ) = 1 2 N σ e r r 2 ¯ x ¯ 2 F v a r 2 + σ e r r 2 ¯ N 1 x ¯ 2
For each epoch, we calculated the value of F v a r separately by considering the average flux ( x ¯ ) and the data points (N) within the epoch only. The values of N, x ¯ , and F v a r for each epoch are listed in Table 1. The errors calculated using Equation (5) are given in Table 1. For epoch A, the error in the value of F v a r is relatively large which also indicates that it is a quiescent state. We considered the source to be variable in epochs if F v a r > 3 × e r r ( F v a r ). According to our criteria for variations, the source was not variable in epoch A (quiescent state), while significant γ -ray flux variations were noticed during epochs B and C with F v a r values of 22% and 35%, respectively.

4.2. γ -ray Spectral Fitting

We modeled the γ -ray spectrum of Ton 599 for each epoch to investigate the inherent distribution of electrons that causes the γ -ray emission during the epoch. We used the power-law (PL) and the log parabola (LP) models to fit the γ -ray spectra of Ton 599 using the maximum likelihood analysis. The PL model is defined as follows [38]:
d N ( E ) d E = N E E Γ
where N is the prefactor (normalization of the energy spectrum), E is the pivot energy (∼523.51 MeV) given in the 4FGL catalog, and Γ is the photon index, while, the LP model is given by
d N ( E ) d E = N E E b α β l n E E b
In this equation N is the normalization, E b is the pivot energy (same as in PL), α is photon index at E b and β is the curvature index that defines the curvature of the spectrum around the peak.
For each epoch, the model fitted γ -ray spectra of Ton 599 is shown in Figure 2, and the values of best fitted model parameters are given in Table 2. In Figure 2, the uncertainties are large at higher energies because of relatively low photon counts at these energies.
To determine whether the γ -ray spectrum has a curvature or not, and which model (LP or PL) best describes the γ -ray spectrum of Ton 599, we calculated T S c u r v e (curvature of the test statistics) [38], which is defined as:
T S c u r v e = 2 ( l o g L L P l o g L P L )
In the above equation, L represents the likelihood function. The value of T S c u r v e for each epoch is given in Table 2. We employed the T S c u r v e > 16 threshold (i.e., 4 σ level; [35]) to determine the presence of statistically significant curvature in the γ -ray spectrum. We found that the LP model best describes the γ -ray spectra of Ton 599 for all three epochs.

4.3. Location of the γ -ray Emission Region

To determine the size of the γ -ray emitting region, the flux doubling time scale is usually estimated, e.g., [31,39]. We calculated the flux doubling time scale for the one-day binned γ -ray light curve during the main-flaring epoch (epoch C) of Ton 599 as follows:
F ( t 2 ) = F ( t 1 ) × 2 Δ t / τ d
In the above expression, F( t 1 ) and F( t 2 ) are the flux values at times t 1 and t 2 respectively, Δ t = t 2 t 1 and τ d represents the flux doubling time.
We found a flux doubling timescale of ∼13.2 h during the epoch C for the blazar Ton 599. Using the flux doubling time scale, we estimated the size of the γ -ray emitting region for Ton 599 using the following expression:
r c τ d δ / ( 1 + z )
A gamma-ray Doppler factor of delta = 12.5 was calculated for Ton 599 by [40] using the multiwavelength data and the known radio Doppler factors with the assumption that the observed boosted emission is from the SSC model. This value is consistent with the lower limits (10.75 and 13.45) of delta obtained for the two bright flares of Ton 599 by [32]. Using δ = 12.5, the size of the γ -ray emission region is estimated to be 1.03 × 10 16 cm.
According to [41], the location of the γ -ray emission region with respect to the central super massive black hole (SMBH) can be approximated by using the relation R = r/ ϕ , where r is the size of the γ -ray emitting blob and ϕ is the jet opening angle.
An intrinsic opening angle of 0.58 (∼0.01 radian) was determined for Ton 599 by [42]. They estimated the jet opening angles for Fermi detected AGNs using 15.4 GHz Very Long Baseline Array (VLBA) observations following two different methods: (a) by analyzing transverse jet profiles in the image plane and (b) by model fitting the data in the (u, v) plane. Using these values of r and ϕ , the location of the γ -ray emitting region for Ton 599 was found to be at 1.03 × 10 18 cm from the SMBH.

5. Discussion

In this section we give our interpretation of the findings of above-mentioned analyses, as well as a discussion of them.

5.1. γ -ray Flux Variability

In the γ -ray band, blazars exhibit remarkable flux variability. For all the three epochs of the source Ton 599, we estimated the flux variability amplitude using one-day binned light curves. We found that during the main-flaring epoch (epoch C), the source showed largest variations with F v a r = 0.35 ± 0.03. During the pre-flaring epoch (epoch B), the source was also variable with F v a r = 0.22 ± 0.04 and during the quiescent epoch (epoch A), the source was not significantly variable ( F v a r = 0.36 ± 0.12) within the error bar. During the main-flare (epoch C), the source showed a maximum flux of 2.24 ± 0.25 × 10 6 ph cm 2 s 1 at MJD 59,399.50, which is larger than the maximum flux observed during the 2017 flare [31]. We found that the photon index value is 1.99 ± 0.10 during the main-flaring epoch (epoch C), when flux is highest, harder than the 4FGL value of 2.19 ± 0.01.
The observed γ -ray variability in the blazars could be attributed to both intrinsic and extrinsic factors. The distribution of relativistic electrons responsible for the emission is one of the intrinsic effects. These relativistic electrons, which can be accelerated to a Lorentz factor of up to 10 6 , are responsible for non-thermal emission from blazar jets via synchrotron or inverse Compton emission processes. The distribution of these relativistic electrons and seed photon responsible for the inverse Compton emission process (synchrotron self Compton (SSC); [15,43] and external Compton (EC); [44,45]) are the intrinsic factors that cause short-term variability in jets. Extrinsic factors, such as the moving blob’s high bulk Lorentz factor ( Γ ) ∼ 50, are in addition to intrinsic factors. Though the large bulk Lorentz factor is a favourable explanation in the case of BL Lacs, because the seed photons come from inside the jets and produce γ -rays through inverse Compton. On the other hand, the seed photons in FSRQ, emanate from outside the jets and produce γ -rays. However, pair-production through dense broad line region (BLR) can attenuate γ -ray emission. As a result, in the case of FSRQs, the large bulk Lorentz factor that causes γ -ray variability is not a plausible explanation [46]. The γ -ray flux variability in blazars jets can be explained by the magnetic reconnection (jet-in-jet) scenario, in which the mini jets generated in the jets can have a large bulk Lorentz factor and they can produce γ -ray flux variability on short time scale [16,19].

5.2. γ -ray Spectra

For each epoch, we modelled the γ -ray spectra of Ton 599 using the power law and the log parabola models. We used the T S c u r v e value to determine the best fitted model. Based on the T S c u r v e value, we found that the LP model best describes the γ -ray spectra of Ton 599 during all the three epochs indicating that the GeV γ -ray spectrum of Ton 599 is curved.
FSRQs usually have a curved GeV γ -ray spectrum, e.g., [47,48,49]. Several theories have been proposed, in the literature, to explain the curved γ -ray spectra of FSRQs. The attenuation of γ -ray photons by pair-production (inside the BLR) could explain the curved γ -ray spectrum [50]. A viable scenario for defining the curve in the γ -ray spectrum is the Klein-Nishina effect on the inverse Compton scattering of BLR photons through relativistic electrons present in the jet [51]. However, several investigations conducted to locate the γ -ray emission region discovered that the observed γ -ray spectra in FSRQs are not caused by IC scattering of BLR photons and that the γ -ray emission site is located outside the BLR [19,52]. The curved γ -ray spectra of FSRQs could also be explained by the curved energy distribution of the electrons emitting the radiation [53,54].

5.3. Location of the γ -ray Emitting Region

We calculated the minimum γ -ray flux doubling timescale of 13.2 h during epoch C for Ton 599. Using this minimum doubling timescale, we estimated the size of the γ -ray emitting region as 1.03 × 10 16 cm. We have also identified the location of the γ -ray emitting blob to be at a distance of 1.03 × 10 18 cm from the SMBH. The size of the BLR for TON 599 was estimated as ∼ (2.11–2.45) × 10 17 cm by [32]. So, the location of the γ -ray emitting blob, found in this work, is outside of the BLR. During the 2017 γ -ray flare of this source, using two component leptonic model, Ref. [33] found that the seed photons for the GeV emission produced by the dusty torus (DT) and the GeV emitting region is located outside of the BLR. Our results are consistent with the findings of [33].

6. Summary

We investigated the γ -ray flux and spectral variability of the blazar Ton 599 during MJD 59,093 to MJD 59,457. For our study, we have chosen three epochs of different flux states: quiescent, pre-flaring, and main-flaring. The outcomes of the analysis of these epochs are summarised below.
  • We estimated the flux variability amplitude for the specified epochs. The largest γ -ray flux variations were found for the main-flaring epoch (C) with F v a r = 0.35 ± 0.03. The source also showed variations in epoch B with F v a r = 0.22 ± 0.04. However, no significant flux variation was observed in the quiescent epoch.
  • The γ -ray spectra were well fit by the LP model for all the three epochs.
  • We estimated the size of the γ -ray emitting region as 1.03 × 10 16 cm and the location of the γ -ray emitting blob could be outside of the BLR.

Author Contributions

Conceptualization, B.R.; Data curation, B.R. and A.P.; Formal analysis, B.R.; Investigation, B.R.; Project administration, A.P.; Resources, B.R.; Visualization, B.R.; Writing—original draft, B.R.; Writing—review & editing, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in this work is publicly available from the Fermi-LAT data archive https://fermi.gsfc.nasa.gov/ssc/data/access/, accessed on 1 September 2021.

Acknowledgments

We thank the anonymous referees for their constructive comments, which helped us to improve the manuscript. This publication made use of data from the Fermi Gamma-ray Space Telescope, which accessed from the Fermi Science Support Center. We also acknowledge the use of High Performance Computing Facility (Nova cluster) at Indian Institute of Astrophysics.

Conflicts of Interest

The authors declare no conflict of interest.

Note

1

References

  1. Urry, C.M.; Padovani, P. Unified Schemes for Radio-Loud Active Galactic Nuclei. Publ. Astron. Soc. Pac. 1995, 107, 803. [Google Scholar] [CrossRef] [Green Version]
  2. Padovani, P.; Alexander, D.M.; Assef, R.J.; De Marco, B.; Giommi, P.; Hickox, R.C.; Richards, G.T.; Smolčić, V.; Hatziminaoglou, E.; Mainieri, V.; et al. Active galactic nuclei: What’s in a name? Astrophys. Space Phys. Rev. 2017, 25, 2. [Google Scholar] [CrossRef] [Green Version]
  3. Stocke, J.T.; Morris, S.L.; Gioia, I.M.; Maccacaro, T.; Schild, R.; Wolter, A.; Fleming, T.A.; Henry, J.P. The Einstein Observatory Extended Medium-Sensitivity Survey. II. The Optical Identifications. Astrophys. J. Suppl. 1991, 76, 813. [Google Scholar] [CrossRef]
  4. Ghisellini, G.; Tavecchio, F.; Foschini, L.; Ghirland, A.G. The transition between BL Lac objects and flat spectrum radio quasars. Mon. Not. R. Astron. Soc. 2011, 414, 2674–2689. [Google Scholar] [CrossRef]
  5. Fossati, G.; Maraschi, L.; Celotti, A.; Comastri, A.; Ghisellini, G. A unifying view of the spectral energy distributions of blazars. Mon. Not. R. Astron. Soc. 1998, 299, 433–448. [Google Scholar] [CrossRef]
  6. Abdo, A.A.; Ackermann, M.; Agudo, I.; Ajello, M.; Aller, H.D.; Aller, M.F.; Angelakis, E.; Arkharov, A.A.; Axelsson, M.; Bach, U.; et al. The Spectral Energy Distribution of Fermi Bright Blazars. Astrophys. J. 2010, 716, 30–70. [Google Scholar] [CrossRef] [Green Version]
  7. Mao, P.; Urry, C.M.; Massaro, F.; Paggi, A.; Cauteruccio, J.; Künzel, S.R. A Comprehensive Statistical Description of Radio-through-Gamma-Ray Spectral Energy Distributions of All Known Blazars. Astrophys. J. Suppl. 2016, 224, 26. [Google Scholar] [CrossRef]
  8. Lin, C.; Fan, J.H.; Xiao, H.B. The intrinsic γ-ray emissions of Fermi blazars. Res. Astron. Astrophys. 2017, 17, 066. [Google Scholar] [CrossRef] [Green Version]
  9. Lind, K.R.; Blandford, R.D. Semidynamical models of radio jets: Relativistic beaming and source counts. Astrophys. J. 1985, 295, 358–367. [Google Scholar] [CrossRef]
  10. Gaur, H.; Gupta, A.C.; Strigachev, A.; Bachev, R.; Semkov, E.; Wiita, P.J.; Peneva, S.; Boeva, S.; Slavcheva-Mihova, L.; Mihov, B.; et al. Optical flux and spectral variability of blazars. Mon. Not. R. Astron. Soc. 2012, 425, 3002–3023. [Google Scholar] [CrossRef] [Green Version]
  11. Rani, B.; Krichbaum, T.P.; Fuhrmann, L.; Böttcher, M.; Lott, B.; Aller, H.D.; Aller, M.F.; Angelakis, E.; Bach, U.; Bastieri, D.; et al. Radio to gamma-ray variability study of blazar S5 0716+714. Astron. Astrophys. 2013, 552, A11. [Google Scholar] [CrossRef] [Green Version]
  12. Pandey, A.; Gupta, A.C.; Wiita, P.J. X-ray Intraday Variability of Five TeV Blazars with NuSTAR. Astrophys. J. 2017, 841, 123. [Google Scholar] [CrossRef] [Green Version]
  13. Prince, R.; Majumdar, P.; Gupta, N. Long-term Study of the Light Curve of PKS 1510-089 in GeV Energies. Astrophys. J. 2017, 844, 62. [Google Scholar] [CrossRef] [Green Version]
  14. Shukla, A.; Mannheim, K.; Patel, S.R.; Roy, J.; Chitnis, V.R.; Dorner, D.; Rao, A.R.; Anupama, G.C.; Wendel, C. Short-timescale γ-ray Variability in CTA 102. Astrophys. J. Lett. 2018, 854, L26. [Google Scholar] [CrossRef]
  15. Marscher, A.P.; Gear, W.K. Models for high-frequency radio outbursts in extragalactic sources, with application to the early 1983 millimeter-to-infrared flare of 3C 273. Astrophys. J. 1985, 298, 114–127. [Google Scholar] [CrossRef]
  16. Giannios, D. Reconnection-driven plasmoids in blazars: Fast flares on a slow envelope. Mon. Not. R. Astron. Soc. 2013, 431, 355–363. [Google Scholar] [CrossRef] [Green Version]
  17. Sironi, L.; Petropoulou, M.; Giannios, D. Relativistic jets shine through shocks or magnetic reconnection? Mon. Not. R. Astron. Soc. 2015, 450, 183–191. [Google Scholar] [CrossRef]
  18. Werner, G.R.; Uzdensky, D.A.; Cerutti, B.; Nalewajko, K.; Begelman, M.C. The Extent of Power-law Energy Spectra in Collisionless Relativistic Magnetic Reconnection in Pair Plasmas. Astrophys. J. Lett. 2016, 816, L8. [Google Scholar] [CrossRef] [Green Version]
  19. Shukla, A.; Mannheim, K. Gamma-ray flares from relativistic magnetic reconnection in the jet of the quasar 3C 279. Nat. Commun. 2020, 11, 4176. [Google Scholar] [CrossRef]
  20. Giannios, D.; Uzdensky, D.A.; Begelman, M.C. Fast TeV variability from misaligned minijets in the jet of M87. Mon. Not. R. Astron. Soc. 2010, 402, 1649–1656. [Google Scholar] [CrossRef] [Green Version]
  21. Aharonian, F.; Akhperjanian, A.G.; Bazer-Bachi, A.R.; Behera, B.; Beilicke, M.; Benbow, W.; Berge, D.; Bernlöhr, K.; Boisson, C.; Bolz, O.; et al. An Exceptional Very High Energy Gamma-Ray Flare of PKS 2155-304. Astrophys. J. Lett. 2007, 664, L71–L74. [Google Scholar] [CrossRef]
  22. Arlen, T.; Aune, T.; Beilicke, M.; Benbow, W.; Bouvier, A.; Buckley, J.H.; Bugaev, V.; Cesarini, A.; Ciupik, L.; Connolly, M.P.; et al. Rapid TeV Gamma-Ray Flaring of BL Lacertae. Astrophys. J. 2013, 762, 92. [Google Scholar] [CrossRef] [Green Version]
  23. Ackermann, M.; Anantua, R.; Asano, K.; Baldini, L.; Barbiellini, G.; Bastieri, D.; Becerra Gonzalez, J.; Bellazzini, R.; Bissaldi, E.; Blandford, R.D.; et al. Minute-timescale > 100 MeV γ-Ray Variability during the Giant Outburst of Quasar 3C 279 Observed by Fermi-LAT in 2015 June. Astrophys. J. Lett. 2016, 824, L20. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, H.T.; Bai, J.M. Absorption of 10-200 GeV Gamma Rays by Radiation from Broad-Line Regions in Blazars. Astrophys. J. 2006, 653, 1089–1097. [Google Scholar] [CrossRef] [Green Version]
  25. Abdollahi, S.; Acero, F.; Ackermann, M.; Ajello, M.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; et al. Fermi Large Area Telescope Fourth Source Catalog. Astrophys. J. Suppl. 2020, 247, 33. [Google Scholar] [CrossRef] [Green Version]
  26. Hewett, P.C.; Wild, V. Improved redshifts for SDSS quasar spectra. Mon. Not. R. Astron. Soc. 2010, 405, 2302–2316. [Google Scholar] [CrossRef] [Green Version]
  27. Fan, J.H.; Tao, J.; Qian, B.C.; Gupta, A.C.; Liu, Y.; Yuan, Y.H.; Yang, J.H.; Wang, H.G.; Huang, Y. Optical Photometrical Observations and Variability for Quasar 4C 29.45. Publ. Astron. Soc. Jpn. 2006, 58, 797–808. [Google Scholar] [CrossRef] [Green Version]
  28. Thompson, D.J.; Bertsch, D.L.; Dingus, B.L.; Esposito, J.A.; Etienne, A.; Fichtel, C.E.; Friedlander, D.P.; Hartman, R.C.; Hunter, S.D.; Kendig, D.J.; et al. The Second EGRET Catalog of High-Energy Gamma-Ray Sources. Astrophys. J. Suppl. 1995, 101, 259. [Google Scholar] [CrossRef]
  29. Abdo, A.A.; Ackermann, M.; Ajello, M.; Allafort, A.; Antolini, E.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; et al. The First Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope. Astrophys. J. 2010, 715, 429–457. [Google Scholar] [CrossRef] [Green Version]
  30. Ramakrishnan, V.; León-Tavares, J.; Rastorgueva-Foi, E.A.; Wiik, K.; Jorstad, S.G.; Marscher, A.P.; Tornikoski, M.; Agudo, I.; Lähteenmäki, A.; Valtaoja, E.; et al. The connection between the parsec-scale radio jet and γ-ray flares in the blazar 1156+295. Mon. Not. R. Astron. Soc. 2014, 445, 1636–1646. [Google Scholar] [CrossRef] [Green Version]
  31. Prince, R. Multi-frequency Variability Study of Ton 599 during the High Activity of 2017. Astrophys. J. 2019, 871, 101. [Google Scholar] [CrossRef] [Green Version]
  32. Patel, S.R.; Chitnis, V.R.; Shukla, A.; Rao, A.R.; Nagare, B.J. Temporal Variability and Estimation of Jet Parameters for Ton 599. Astrophys. J. 2018, 866, 102. [Google Scholar] [CrossRef]
  33. Patel, S.R.; Chitnis, V.R. Leptonic modelling of Ton 599 in flare and quiescent states. Mon. Not. R. Astron. Soc. 2020, 492, 72–78. [Google Scholar] [CrossRef]
  34. Atwood, W.B.; Abdo, A.A.; Ackermann, M.; Althouse, W.; Anderson, B.; Axelsson, M.; Baldini, L.; Ballet, J.; Band, D.L.; Barbiellini, G.; et al. The Large Area Telescope on the Fermi Gamma-Ray Space Telescope Mission. Astrophys. J. 2009, 697, 1071–1102. [Google Scholar] [CrossRef] [Green Version]
  35. Mattox, J.R.; Bertsch, D.L.; Chiang, J.; Dingus, B.L.; Digel, S.W.; Esposito, J.A.; Fierro, J.M.; Hartman, R.C.; Hunter, S.D.; Kanbach, G.; et al. The Likelihood Analysis of EGRET Data. Astrophys. J. 1996, 461, 396. [Google Scholar] [CrossRef]
  36. Vaughan, S.; Edelson, R.; Warwick, R.S.; Uttley, P. On characterizing the variability properties of X-ray light curves from active galaxies. Mon. Not. R. Astron. Soc. 2003, 345, 1271–1284. [Google Scholar] [CrossRef]
  37. Pandey, A.; Gupta, A.C.; Wiita, P.J. X-ray Flux and Spectral Variability of Six TeV Blazars with NuSTAR. Astrophys. J. 2018, 859, 49. [Google Scholar] [CrossRef] [Green Version]
  38. Nolan, P.L.; Abdo, A.A.; Ackermann, M.; Ajello, M.; Allafort, A.; Antolini, E.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; et al. Fermi Large Area Telescope Second Source Catalog. Astrophys. J. Suppl. 2012, 199, 31. [Google Scholar] [CrossRef] [Green Version]
  39. Paliya, V.S. Fermi-Large Area Telescope Observations of the Exceptional Gamma-Ray Flare from 3C 279 in 2015 June. Astrophys. J. Lett. 2015, 808, L48. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, L.Z.; Fan, J.H.; Cheng, K.S. The Multiwavelength Doppler Factors for a Sample of Gamma-Ray Loud Blazars. Publ. Astron. Soc. Jpn. 2002, 54, 159–169. [Google Scholar] [CrossRef] [Green Version]
  41. Foschini, L.; Ghisellini, G.; Tavecchio, F.; Bonnoli, G.; Stamerra, A. Search for the shortest variability at gamma rays in flat-spectrum radio quasars. Astron. Astrophys. 2011, 530, A77. [Google Scholar] [CrossRef] [Green Version]
  42. Pushkarev, A.B.; Kovalev, Y.Y.; Lister, M.L.; Savolainen, T. Jet opening angles and gamma-ray brightness of AGN. Astron. Astrophys. 2009, 507, L33–L36. [Google Scholar] [CrossRef]
  43. Ghisellini, G.; Maraschi, L. Bulk acceleration in relativistic jets and the spectral properties of blazars. Astrophys. J. 1989, 340, 181–189. [Google Scholar] [CrossRef]
  44. Begelman, M.C.; Sikora, M.; Giommi, P.; Barr, P.; Garilli, B.; Gioia, I.M.; Maccacaro, T.; Maccagni, D.; Schild, R.E. Inverse Compton scattering of ambient radiation by a cold relativistic jet - A source of beamed, polarized continuum in blazars? Astrophys. J. 1987, 322, 650–661. [Google Scholar] [CrossRef]
  45. Sikora, M.; Begelman, M.C.; Rees, M.J. Comptonization of diffuse ambient radiation by a relativistic jet: The source of gamma rays from blazars? Astrophys. J. 1994, 421, 153–162. [Google Scholar] [CrossRef]
  46. Sbarrato, T.; Foschini, L.; Ghisellini, G.; Tavecchio, F. Study of the variability of blazars gamma-ray emission. Adv. Space Res. 2011, 48, 998–1003. [Google Scholar] [CrossRef] [Green Version]
  47. Abdo, A.A.; Ackermann, M.; Ajello, M.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Bechtol, K.; et al. Spectral Properties of Bright Fermi-Detected Blazars in the Gamma-Ray Band. Astrophys. J. 2010, 710, 1271–1285. [Google Scholar] [CrossRef]
  48. Paliya, V.S.; Sahayanathan, S.; Stalin, C.S. Multi-Wavelength Observations of 3C 279 During the Extremely Bright Gamma-Ray Flare in 2014 March-April. Astrophys. J. 2015, 803, 15. [Google Scholar] [CrossRef]
  49. Rajput, B.; Stalin, C.S.; Sahayanathan, S.; Rakshit, S.; Mandal, A.K. Temporal correlation between the optical and γ-ray flux variations in the blazar 3C 454.3. Mon. Not. R. Astron. Soc. 2019, 486, 1781–1795. [Google Scholar] [CrossRef]
  50. Coogan, R.T.; Brown, A.M.; Chadwick, P.M. Localizing the γ-ray emission region during the 2014 June outburst of 3C 454.3. Mon. Not. R. Astron. Soc. 2016, 458, 354–365. [Google Scholar] [CrossRef] [Green Version]
  51. Cerruti, M.; Dermer, C.D.; Lott, B.; Boisson, C.; Zech, A. Gamma-Ray Blazars near Equipartition and the Origin of the GeV Spectral Break in 3C 454.3. Astrophys. J. Lett. 2013, 771, L4. [Google Scholar] [CrossRef] [Green Version]
  52. Costamante, L.; Cutini, S.; Tosti, G.; Antolini, E.; Tramacere, A. On the origin of gamma-rays in Fermi blazars: Beyondthe broad-line region. Mon. Not. R. Astron. Soc. 2018, 477, 4749–4767. [Google Scholar] [CrossRef]
  53. Tramacere, A.; Giommi, P.; Perri, M.; Verrecchia, F.; Tosti, G. Swift observations of the very intense flaring activity of Mrk 421 during 2006. I. Phenomenological picture of electron acceleration and predictions for MeV/GeV emission. Astron. Astrophys. 2009, 501, 879–898. [Google Scholar] [CrossRef] [Green Version]
  54. Dermer, C.D.; Yan, D.; Zhang, L.; Finke, J.D.; Lott, B. Near-equipartition Jets with Log-parabola Electron Energy Distribution and the Blazar Spectral-index Diagrams. Astrophys. J. 2015, 809, 174. [Google Scholar] [CrossRef] [Green Version]
Figure 1. One-day binned γ -ray light curve of Ton 599 (upper panel) with the TS values corresponding to each flux value (lower panel). The vertical black lines refer to the quiescent epoch (Epoch A), whereas vertical blue lines refer to the γ -ray flaring epochs (Epoch B and Epoch C). In the upper panel the horizontal cyan line represents the average flux from September 2020 to August 2021.
Figure 1. One-day binned γ -ray light curve of Ton 599 (upper panel) with the TS values corresponding to each flux value (lower panel). The vertical black lines refer to the quiescent epoch (Epoch A), whereas vertical blue lines refer to the γ -ray flaring epochs (Epoch B and Epoch C). In the upper panel the horizontal cyan line represents the average flux from September 2020 to August 2021.
Galaxies 09 00118 g001
Figure 2. γ -ray spectrum for the epochs (AC). The name of the epoch is mentioned in each plot.
Figure 2. γ -ray spectrum for the epochs (AC). The name of the epoch is mentioned in each plot.
Galaxies 09 00118 g002
Table 1. The epochs considered for the γ -ray light curve study are listed in the table below. Here, N, x ¯ , and F v a r denote the total number of data points, the average γ -ray flux, and the fractional variability amplitude, respectively, for the epoch. The γ -ray flux is in the units of 10 6 ph cm 2 s 1 .
Table 1. The epochs considered for the γ -ray light curve study are listed in the table below. Here, N, x ¯ , and F v a r denote the total number of data points, the average γ -ray flux, and the fractional variability amplitude, respectively, for the epoch. The γ -ray flux is in the units of 10 6 ph cm 2 s 1 .
EpochsMJDCalendar Date (dd-mm-yyyy)Peak FluxN x ¯ F var Remarks
StartEndStartEnd
Epoch A59,14459,22622-10-202012-01-20210.37 ± 0.16420.15 ± 0.010.36 ± 0.12Quiescent
Epoch B59,31559,36011-04-202124-05-20211.16 ± 0.18440.71 ± 0.020.22 ± 0.04Flaring
Epoch C59,37459,41109-06-202116-07-20212.24 ± 0.25371.02 ± 0.030.35 ± 0.03Flaring
Table 2. Details of the PL and LP model fits for the three epochs. Here, γ -ray flux is in the units of 10 6 ph cm 2 s 1 .
Table 2. Details of the PL and LP model fits for the three epochs. Here, γ -ray flux is in the units of 10 6 ph cm 2 s 1 .
EpochPLLP TS curve
Γ Flux−Log L α β Flux−Log L
A−2.381 ± 0.0040.135 ± 0.002146,206.1602.045 ± 0.0851.075 ± 0.4730.111 ± 0.012146,191.32629.67
B−2.142 ± 0.0010.800 ± 0.00292,649.8202.051 ± 0.0210.973 ± 0.1170.735 ± 0.02092,634.26631.108
C−2.044 ± 0.0021.084 ± 0.00584,267.7001.923 ± 0.0020.958 ± 0.0110.980 ± 0.00384,245.25844.884
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rajput, B.; Pandey, A. γ-ray Flux and Spectral Variability of Blazar Ton 599 during Its 2021 Flare. Galaxies 2021, 9, 118. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9040118

AMA Style

Rajput B, Pandey A. γ-ray Flux and Spectral Variability of Blazar Ton 599 during Its 2021 Flare. Galaxies. 2021; 9(4):118. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9040118

Chicago/Turabian Style

Rajput, Bhoomika, and Ashwani Pandey. 2021. "γ-ray Flux and Spectral Variability of Blazar Ton 599 during Its 2021 Flare" Galaxies 9, no. 4: 118. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9040118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop