arxiv: v1 [astro-ph.he] 26 Feb 2015

Transcription

1 th Fermi Symposium : Nagoya, Japan : 20 Oct, Spectral Studies of Flaring FSRQs at GeV Energies Using Pass 8 Fermi-LAT Data R. J. Britto, S. Razzaque Department of Physics, University of Johannesburg, Auckland Park 2006, South Africa B. Lott Centre d Études Nucléaires de Bordeaux-Gradignan, Université Bordeaux 1, UMR 797/IN2P, 17 Gradignan, France on behalf of the Fermi-LAT Collaboration arxiv: v1 [astro-ph.he] 26 Feb 201 Flat spectrum radio quasars (FSRQs) are bright active galactic nuclei surrounded by gas clouds within a UVvisible intense radiation field that form the so-called broad line region (BLR). These objects emit relativistic jets from a region close to the central supermassive black hole and through the BLR. The Fermi-Large Area Telescope (Fermi-LAT) is sensitive to gamma-ray photons from 0 MeV to more than 00 GeV. We have performed spectral analysis of bright FSRQs in a. year (200801) data sample collected by Fermi-LAT, using the new Pass 8 event selection and instrument response function. Also, our study of flaring episodes in a limited time range brings interesting results while compared to the full. year data samples. 1. Modelling the BLR radiation field FSRQs constitute a class of active galactic nuclei (AGNs) with a dense BLR in which gamma rays with energy > GeV are absorbed due to electronpositron pair creation, if produced deep inside the BLR. Indeed, BLR is expected to be denser in FS- RQs compared to the BL Lac class. Operating since 2008, the Fermi satellite has amassed more than 6 years of data, continuously surveying the whole sky []. The sensitivity of Fermi- LAT is ideal for the study of the gamma-ray absorption inside FSRQs in the 0 MeV-00 GeV range. From constraints on gamma-ray absorption we may infer limits on the location of the gamma-ray emission region in the FSRQ jets. We expect > GeV photons of FSRQs to undergo absorption in the BLR, where the target photon with energy ǫ is a UV photon from the BLR radiation field. As most of these photons are expected to come from the emission lines, we use a model that includes the 6 strongest lines (NV, Lyα, OVI Lyβ, CIII NIII, NeVIII OIV, HeII Lyα) between to 1 ev [1]. We model these lines using a Breit-Wigner distribution, given by: n i ω i BW(ǫ) = 2 π[(ǫ ǫ i ) 2 +(ω i /2) 2 ], (1) where n i and ω i are the number density and width, respectively, for a given line i. Under the commonly used relations L BLR = 0.1 L disc and R BLR = L disc [, 7], where L BLR is the luminosity of the BLR, R BLR its radius, and L disc the luminosity of the accretion disc. The photon density n i of the radiation field for each line i can be written: ( )( n i [cm ] L i 17 ) cm erg s 1 ǫ i,ev RBLR 2. (2) The opacity is derived from [, 8] and is expressed as a function of E and z: [ dτ γγ dx (E,z) = r2 0 m 2 c ] 2 2 E(1+z) [ 6 ϕ ǫe(1+z) m n i ω 2 c ]dǫ i, m 2 c [(ǫ ǫ E(1+z) i ) 2 +(ω i /2) 2 ]ǫ 2 () i=1 where r 0 is the classical electron radius, m the electron mass. The τ γγ opacity in the BLR is then calculated as following: τ γγ (E,z) = a R BLR dτ γγ (E,z), () dx assuming the gamma rays are produced within a R BLR, where R BLR is the outer radius of BLR, and a < 1. Since this absorption happens at some distance from the supermassive black hole, this corrective factor that we called a represents the fraction of the BLR responsible for the absorption. In Table I are displayed the line properties of the average spectrum of quasars we used in our model, as they were given in [1]. Since the He II Lyα line has quite large uncertainties, we arbitrary fixed its EW and relative flux to be equal to the ones of N V (uncertainty represented by (*)). Very high energy gamma rays travelling from far distances undergo absorption in the extragalactic background light ( EBL, mainly composed of infrared-uv radiation). This absorption is to be considered above GeV and has been implemented in our studies, from the model presented in [6]. Evidence of absorption in the BLR for some FSRQs have been reported in [11, 12].

2 γ 2 th Fermi Symposium : Nagoya, Japan : 20 Oct, 201 τ γ Opacity in the BLR using the Breit-Wigner line model C. / z = 0.89 / L BLR =.00 N V (.00 ev) Ly α O VI + Ly β C III + N III (.20 ev) (12.0 ev) (12.6 ev) Ne VIII + O IV (1.90 ev) He II Ly α (0.81 ev) erg.s -1 F (erg cm s -1 ) C. /0 Aug Jan 201 / Pass 7 vs Pass 8 PASS 7 PASS E (MeV) Figure 1: Opacity τ γγ(e,z) versus Energy, for C.. The opacity sum on the 6 lines is drawn in plain black. Table I Properties of main lines of the average spectrum of quasars as compiled in [1], and He II Lyα as we defined it for this study. Line ǫ (ev) EW (ev) Relative flux NV Lyα OVI Lyβ CIII NIII NeVIII OIV HeII Lyα * 0.22* 2. Data processing and model fitting We have analysed data of 7 bright gamma-ray FSRQs. Plots of the spectral energy distributions (SEDs) under the label. years have been processed from August 2008 till 0 April 201. The sources we present in this paper are listed in Table II, LII and BII being respectively the Galactic latitude and longitude in decimal degrees. Data were processed using the Pass 8 data representation (P8 SOURCE V), and the Science Tools version v Signal is reconstructed from each source using the unbinned likelihood tool 1, applied to LAT data in the GeV energy range, within a region of interest (ROI) of radius. A sourceregion extended to an additional annulus accounted for all the point sources of the Second Fermi-LAT source catalog [9], and for the Galactic diffuse emission (template years P8 V2 scaled) and the isotropic diffuse emission (isotropic source years P8V). We computed the SEDs for all the sources of the selected sample with the Pass 8 data representation. Additionally, for the two brightest objects of our 1 LogParabola Power Law with exp cutoff Figure 2: Comparison of the SED of C., produced with the Pass 7 (magenta) and Pass 8 (green) data sets. FSRQ sample, i.e. C. and PKS 1-089, we also computed the SEDs with the PASS 7 reprocessed dataset (P7REP SOURCE V1) and verified the consistencyoftheresultswith respecttothepass8ones. Although some bin-to-bin fluctuations appear due to energy wise event migrations, the two SEDs (Pass 7 and Pass 8) for both the sources are compatible (as shown in Figure 2 for C.). A first set of fits was performed from 0 MeV till the highest energy data point (excluding upper limits), using a log-parabola (LP: dn(e)/de = N 0 (E/E 0 ) α β log(e/e0)), with E 0 kept fixed at MeV, and where log is the natural logarithm), a broken power law (BPL: dn(e)/de = N 0 (E/E b ) Γi, with i = 1 if E < E b and i = 2 if E > E b ), and a power law with an exponential cutoff (PLEC: dn(e)/de = N 0 (E/E p ) ΓPLEC exp( E/E c ), with E p kept fixed at 12.7 MeV). Other sets of fits were performed by adding exponential factors to model the EBL absorption [6] and the opacity of the BLR. As the fits presentedin the Sections and arebinning dependant, the values of the fit parameters vary from one choice of binning to another. Narrow data binning could held spurious fluctuations, while wide data-binning could hide features. In order to estimate this systematic effect, we do the following: a first LP fit is performed on the SED, while keeping all parameters fixed to the values obtained by the unbinned likelihood analysis, and a χ 2 /ndf is returned. A second LP fit is performed with N 0, α and β kept free,

3 th Fermi Symposium : Nagoya, Japan : 20 Oct, 201 Table II Characteristics of the 7 FSRQs used in this paper. Sources are ordered by decreasing flux in the 1-0 GeV energy range. Values of the luminosity in the broad line region (L BLR) are taken or derived from [, 1]. Name RA Dec LII BII Redshift Flux 2FGL Flux 1FHL Photon index L BLR (J2000.0) (J2000.0) (deg.) (deg.) z 1-0 GeV -0 GeV 1FGL ( erg.s 1 ) C e-8 1.e PKS e-8 7.e C e-8 7.e C e-8.7e PKS e e B e-8.27e PKS (B) e-8 2.9e E b kept fixed to the value returned by the unbinned likelihood analysis. We compare the fits to validate that the binned fit is compatible with the unbinned fit, though results differ. These systematics could be overcome by implementing an unbinned analysis for all fitted models in future. Our modelling study is then done by fitting the SEDs using the LP, BPL and PLEC function that all now include EBL. These fits will be reported in the Sections and, and compared to the fits that include both EBL and BLR absorption. The latter fits are written as LPτ, BPLτ, and PLECτ. The observed spectrum F obs (E) will then be ultimately written: F obs (E) = e τebl(e,z) e a τγγ(e,z) F int (E), () where F int (E) is the LP, BPL orplec fitting function. Parametera is kept free in the [, 1] range to account for the fraction of radius of the BLR in which gamma rays may be absorbed. The fitting procedure using the absorption models is implemented by interpolation of both the τ EBL and τ γγ graphs. While comparing each EBL + BLR absorption fit (LPτ, BPLτ or PLECτ) with the EBL + no BLR absorption fit (LP, BPL or PLEC), and if both fits have a χ 2 /ndf < 1, we obtain a p-value which indicates the discrepancy between the fit with model and the fit without model, for a given function.. Results on the. years of data We present now the results on the 7 bright FSRQs during the. years period previously defined in Section 2. Data under the label. years are processed from August 2008 till 0 April 201. In Figure are presented the SEDs of C., PKS 1-089, TXS and PKS Some of the fits are not visible, mainly the BPL and PLEC. The reason is that they are overwritten by the BPLτ and PLECτ functions, for which the parameter a is very small. We would consider having evidence for absorption in the BLR if we get all of the following: atleastonegoodqualityfitamongoneofthe fits with EBL+BLR absorption (LPτ, BPLτ and PLECτ); parameter a with a relatively small error bar; a small p-value or a bad fit of the corresponding function with only EBL absorption (LP, BPL or PLEC). In Table III are displayed the fit parameters of the 7 sources, along with the p-values used to compare the models (with versus without absorption). In blue bold face are the parameters that suggest a possible BLR absorption, as some of the above conditions are partiallymet forc.(with ap-valueof6.2 /.9 σ C.L.), and for TXS 120+1(with a p-value of 9. /2.6σ C.L.). Wehavenohintofabsorption for the other sources we studied.. Results on high state/flaring episodes Data were analysed during flaring/high state periods for the following sources: C. (high state and giant flare) during 02 Nov-0 Dec 20 (MJD 02.-.)[1]; PKS 1-089during 19 Feb-0 Apr 2012 (MJD ), along with MAGIC data from the same period (MAGIC data points taken from [2]); PKS 12-1 during 0 Sep Jul 201 (MJD ), as a combined series of successive radio flares). During these outburst episodes, the gamma-ray emission region can have a different location compared

4 th Fermi Symposium : Nagoya, Japan : 20 Oct, 201 C. (Pass 8 -. years) PKS (Pass 8 -. years) z = 0.89 L BLR =.00 z = 0.60 L BLR = Log-parabola x exp( -aτ γ γ - τ EBL ) Log-parabola x exp( -aτ γ γ - τ EBL ) BPL x exp(-a τ ) γ γ BPL x exp(-a τ ) γ γ PLEC x exp(-a τ γ ) γ PLEC x exp(-a τ γ ) γ TXS (Pass 8 -. years) PKS 12-1 (Pass 8 -. years) z = 1.8 L BLR = 8.00 z = 1.22 L BLR = Log-parabola x exp( -aτ γ γ - τ EBL ) Log-parabola x exp( -aτ γ γ - τ EBL ) BPL x exp(-a τ γ ) γ BPL x exp(-a τ γ ) γ PLEC x exp(-a τ γ ) γ PLEC x exp(-a τ γ ) γ Figure : SEDs of C., PKS 1-089, TXS and PKS BPL and PLEC fits are often hidden beneath BPLτ and PLECτ. to the quiescent state. In this section we present the results obtained on these three high states (Figure and Table IV). Under the same criteria than the ones used in Section, we still find no evidence of BLR absorption, though we still have a hint of it for C. with a p-value of around 2.2 for the discrepancy between the BPL and BPLτ fits (significance of about σ). Due to the unusual shape of the SED of PKS 12-1 during this series of flares, all fits have a large χ 2 value. Though we dispose of less photons in the data analysis of flaring episodes, during strong and long flares, it could be possible to constrain the location of the gamma-ray production region if it is deep enough inside the BLR.

5 th Fermi Symposium : Nagoya, Japan : 20 Oct, 201 C. (Pass 8 - Nov-Dec20) PKS (Pass 8 - Feb-Apr 2012) PKS 12-1 (Pass 8 - Sep-Jul 201) -8 z = 0.89 L BLR =.00-8 z = 0.60 L BLR =.62 z = 1.22 L BLR = MAGIC Log-parabola x exp( -aτ γ γ - τ EBL ) Log-parabola x exp( -aτ γ γ - τ EBL ) BPL x exp(-a τ γ ) γ BPL x exp(-a τ γ ) γ PLEC x exp(-a τ ) γ γ Log-parabola x exp( -aτ γ γ - τ EBL ) PLEC x exp(-a τ ) γ γ Figure : SEDs of C., PKS and PKS 12-1 during the outburst periods. PLEC fits are hidden beneath PLECτ for C. and PKS The flare of PKS was studied along with MAGIC data above 90 GeV and the combined LAT-MAGIC SED was fitted only with LP and LPτ function (LP hidden beneath LPτ).. Conclusions and perspectives We find that the gamma-gamma absorption in the BLR is not significant enough to claim discovery for the models of BLR and spectral functions we have investigated. There are hints of absorption in case of C. and TXS with significance of the order of σ. An implication of our results could be that the gamma-ray emission zone in FSRQs might be located outside or at the outer edge of the BLR. However, further investigation on binning effects on the SED fits are required. Future work is also expected to improve the modelling of the BLR. Acknowledgments The Fermi-LAT Collaboration acknowledges support for LAT development, operation and data analysis from NASA and DOE (United States), CEA/Irfu and IN2P/CNRS (France), ASI and INFN (Italy), MEXT, KEK, and JAXA (Japan), and the K.A. Wallenberg Foundation, the Swedish Research Council and the National Space Board (Sweden). Science analysis support in the operations phase from INAF (Italy) and CNES (France) is also gratefully acknowledged. Thanks to Pfesesani Van Zyl for providing the PKS 12-1 dates during radio flares and to Julian Sitarek and the MAGIC Collaboration for providing the MAGIC data points. References [1] A. A. Abdo et al, 2011, ApJL, 7:L26 [2] J. Aleksic et al, 201, A&A, 69, A6 [] W. B. Atwood et al, 2009, ApJ, 697, 71 [] J. A. Baldwin & H. Netzer, 1978, ApJ, 226, [] R. W. Brown, K. O. Mikaelian and R. J. Gould, 197, ApJ L, 1, L20-L20 [6] J. D. Finke, S. Razzaque & C. D. Dermer, 20, ApJ, 712, 28 [7] G. Ghisellini & F. Tavecchio, 2009, MNRAS, 97, 98 [8] R. G. Gould & G. P. Schréder, 1967, Physical Review, 1, [9] P. L. Nolan, et al, 2012, ApJ S, 199, 1) [] L. Pacciani et al, 201, ApJ, 790, [11] J. Poutanen & B. Stern, 20, ApJ L, 717:L118 [12] B. Stern & J. Poutanen, 201, ApJ, 79, 8 [1] R. C. Telfer, W. Zheng, G. A. Kriss and A. F. Davidsen, 2002, ApJ, 6:77-78 [1] D. R. Xiong & X. Zhang, 201, MNRAS, 1,, p.7-9

6 Parameter Model C. PKS C +21. C 279 PKS 0 TXS PKS12 LP unbinned α 2.28 ± ± ± ± ± ± ± LP unbinned β ± ± ± ± ± ± ± PL E b 298 ±. 28 ± ±.8 28 ±.2 8 ± ± ± 2.8 function LP α 2.20 ± ± ± ± ± ± ± 0.01 (without BLR LP β ± ± ± ± ± ± ± absorption BPL Γ ± ± ± ± ± ± ± 0.0 model) BPL Γ ± ± ± ± ± ± ± BPL E b 107 ± ± ± ± ± ± ± 7.6 PLEC Γ PLEC 2.21 ± ± ± ± ± ± ± 0.01 PLEC E c 820 ± ± ± ± ± ± ± LPτ α 2.2 ± ± ± ± ± ± ± 0.01 LPτ β ± ± ± ± ± ± ± function BPLτ Γ ± ± ± ± ± ± ± (with BLR BPLτ Γ ± ± ± ± ± ± ± absorption BPLτ E b 122 ± ± ± ± 7.9 ± ± ± 1.2 model) PLECτ Γ PLEC 2.21 ± ± ± ± ± ± ± 0.01 PLECτ E c 821 ± ± ± ± ± ± ± LPτ ± ± ± ± ± ± ± a BPLτ ± ± ± ± ± ± ± PLECτ ± ± ± ± ± ± ± LP unbinned 6.19 (0) 7. (17) (16).288 (16) 1.78 (16) (1) 6.17 (16) χ 2 (ndf) LP (27) (1) (1).609 (1) (1) 6.8 (12) (1) without abs. BPL.67 (26) 9.96 (1) (12) (12) (12) (11) (12) model PLEC (26) 1.11 (1) (12) 7.22 (12) (12).18 (11) (12) χ 2 (ndf) LPτ (26) 7.2 (1) (12).82 (12) 6.26 (12).181 (11) 20.9 (12) with abs. BPLτ (2) 7.22 (12) (11).69 (11).282 (11) () 17.7 (11) model PLECτ.017 (26) (1) (12) 7.2 (12) (12).19 (11) (12) χ 2 LP/PLτ 1.17e+01.60e e e e e e+00 (ndf) BPL/BPLτ 2.08e e e e e e e+00 PLEC/PLECτ 1.812e e e-0 1.6e e e-0.02e-0 LP/LPτ 6.166e-0 8.6e e e e e-02 - p-value BPL/BPLτ e e e e-0 - PLEC/PLECτ Table III Fitting parameters and derived significances for the 7 bright FSRQs of our sample. We mention (*) when a reached the lower edge of the fitting interval. 6 th Fermi Symposium : Nagoya, Japan : 20 Oct, 201

7 th Fermi Symposium : Nagoya, Japan : 20 Oct, Table IV Fitting parameters and derived significances for the flares we studied in our FSRQ sample. We mention (*) when a reached the lower edge of the fitting interval. Parameter Model C. PKS PKS12 (Nov-Dec 20) (Feb-Apr 2012) (Sep 2012-Jul 201) LP unbinned α 2.12 ± ± 0.0 LP unbinned β ± ± PL E b 286 ± ±. function LP α 2.1 ± ± 0.01 (without BLR LP β ± ± absorption BPL Γ ± ± 0.00 model) BPL Γ ± ± 0.0 BPL E b 921 ± ± 1.6 PLEC Γ PLEC 2.18 ± ± PLEC E c 9981 ± ± LPτ α 2.17 ± ± 0.01 LPτ β ± ± function BPLτ Γ ± ± (with BLR BPLτ Γ ± ± 0.06 absorption BPLτ E b 921 ± ± 0.2 model) PLECτ Γ PLEC 2.18 ± ± PLECτ E c 998 ± ± 11.7 LPτ ± ± a BPLτ 0.09 ± ± PLECτ ± ± LP unbinned. (1) (17) 79.6 (1) χ 2 (ndf) LP.01 (11) (1) 7.81 (12) without abs. BPL 1.18 () 1.17 (1) 1.8 (11) model PLEC () 7.61 (1) 77.6 (11) χ 2 (ndf) LPτ () (1) 7.89 (11) with abs. BPLτ.780 ( 9) () model PLECτ 6.91 () (11) χ 2 LP/PLτ.11e e 7.8e-0 (ndf) BPL/BPLτ 9.78e e+01 PLEC/PLECτ 1.86e e-0 LP/LPτ 7.769e p-value BPL/BPLτ 2.19e PLEC/PLECτ 0.000e

8 Format Problems, source edited Bad or no PostScript file Bad Fonts Figure Problems Unreadable disks/missing files Paper too long Equation problems Author requested changes Pages numbered/empty etc