Multiwavelength picture of the blazar S5 0716+714 during its brightest outburst

TL;DR

This paper presents a comprehensive multiwavelength analysis of the brightest outburst of the blazar S5 0716+714 in 2015, combining data from radio to VHE gamma rays to study its spectral energy distribution and variability.

Contribution

It provides the first detailed multiwavelength spectral and temporal analysis of S5 0716+714 during its brightest optical and gamma-ray outburst, integrating data from numerous observatories.

Findings

Detection of VHE gamma-ray emission with 13σ significance during the outburst

Construction of the broad-band spectral energy distribution during the flare

Analysis of time variability across multiple wavelengths

Abstract

S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13{\sigma} significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Mets\"ahovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst.

Figures

Full text

\articlenumber

x \doinum10.3390/—— \pubvolumexx

\externaleditorAcademic Editors: \history

\TitleMultiwavelength Picture of the Blazar S5 0716+714 during Its Brightest Outburst \AuthorMarina Manganaro 1,∗, Giovanna Pedaletti 2, Marlene Doert 3, Denis Bastieri 4, Vandad Fallah Ramazani 5, Dario Gasparrini 6, Elina Lindfors 5, Benoit Lott 7, Mireia Nievas 8, Bindu Rani 9, David J. Thompson 9, for the MAGIC and Fermi-LAT Collaborations; Emmanouil Angelakis 10, George Borman 11, Mark Gurwell 12, Talvikki Hovatta 13, Ryosuke Itoh 14, Svetlana Jorstad 15, Alex Kraus 10, Thomas P. Krichbaum 10, Paul Kuin 16, Anne Lähteenmäki 13,17, Valeri Larionov 18, Amy Yarleen Lien 9, Ioannis Myserlis 10, Merja Tornikoski 13, Ivan Troitsky 18, J. Anton Zensus 10 \AuthorNamesMarina Manganaro, Giovanna Pedaletti, Marlene Doert, Denis Bastieri, Vandad Fallah Ramazani, Dario Gasparrini, Elina Lindfors, Benoit Lott, Mireia Nievas, Bindu Rani, David J. Thompson, for the MAGIC and Fermi-LAT Collaborations; Emmanouil Angelakis, George Borman, Mark Gurwell, Talvikki Hovatta, Ryosuke Itoh, Svetlana Jorstad, Alex Kraus, Thomas P. Krichbaum, Paul Kuin, Anne Lähteenmäki, Valeri Larionov, Amy Yarleen Lien, Ioannis Myserlis, Merja Tornikoski, Ivan Troitsky, J. Anton Zensus

\corresCorrespondence: manganaro@iac.es; Tel.: +34-922-605-551

\conferencetitle"Blazars through sharp MW eyes"

1 Introduction

The BL Lac object S5 0716+714 is among the brightest and most active blazars Rani2013 . Because of the featureless optical continuum, it is hard to estimate its redshift. A value of $z=0.31\pm 0.08$ is derived from the photometric detection of the host galaxy nilsson2008 . It was detected in VHE for the first time by MAGIC in 2008 Anderhub2009 . Because of its extreme variability, S5 0716+714 has been the subject of several optical monitoring campaigns (e.g., wagner1996 ; rani2011 ; montagni2006 ; Rani2013 ; rani2015 and references therein). S5 0716+714 is also among the bright blazars in the Fermi-LAT (Large Area Telescope) Bright AGN Sample (LBAS) abdo2010LBAS , whose GeV spectra are governed by a broken power law. The combined GeV–TeV spectrum of the source displays absorption-like features in 10–100 GeV energy range senturk2013 . The blazar S5 0716+714 was observed in an unprecedented outburst phase during January 2015. The source was detected at its historic high brightness at optical and IR bands. In January 2015, MAGIC, triggered by the high optical state and by high energy photons detected by Fermi-LAT, detected the source with a significance of $\sim 13\sigma$ Atel6999 .

2 Instruments Involved

MAGIC is a stereoscopic system consisting of two 17 m diameter Imaging Atmospheric Cherenkov Telescopes located at the Observatorio del Roque de los Muchachos, on the Canary Island of La Palma (Spain). The current integral sensitivity for the zenith range $35^{\circ}<zd<50{{}^{\circ}}$ above 100 GeV is $1.76\%\pm 0.03\%$ of the Crab Nebula’s flux in 50 h. The Field of View (FoV) of MAGIC is $3.5^{\circ}$ Aleksic2015_Upgrade2 . The Large Area Telescope (LAT) on board of Fermi is an imaging high-energy $\gamma$-ray telescope in the energy range from about 20 MeV to more than 300 GeV Atwood2009 . The LAT’s Field of View covers about 20% of the sky at any time, and it scans continuously, covering the whole sky every three hours. Multiwavelength observations were performed collecting data from a large number of instruments. Our dataset includes data from Swift-XRT Burrows2005 in the X-ray and Swift-UVOT Roming2005 in optical-UV; for the optical band we show data from the Tuorla Blazar monitoring program (http://users.utu.fi/kani/1m), the 1.8 meter Perkins telescope (https://lowell.edu/research/research-facilities/1-8-meter-perkins/), the Steward Observatory (http://www.bu.edu/blazars/), the 70-cm AZT-8 reflector of the Crimean Astrophysical Observatory (http://craocrimea.ru/ru), the 40-cm LX-200 telescope in St. Petersburg, Russia, HagenTorn2006 and Kanata telescope at the Hiroshima Observatory (http://hasc.hiroshima-u.ac.jp/telescope/kanatatel-e.html); radio band was covered by Metsähovi Terasranta , the Owens Valley Radio Observatory (OVRO) 40 m telescope Richards2011 and the Effelsberg 100 m radio telescope (http://www.mpifr-bonn.mpg.de/en/effelsberg). Radio data at 230 GHz band are obtained using the Submillimeter Array (http://sma1.sma.hawaii.edu/) (SMA).

3 Observations and Results

3.1 MAGIC Analysis

MAGIC started to observe S5 0716+714 on 19 January 2015. Unfortunately, the following two nights of strong wind during the S5 0716+714 visibility window prevented observation of the source, which was resumed on 23 January. From that day, the flux constantly increased up to the maximum flux ever observed in the VHE range for this source, on 27 January. The next activity of S5 0716+714 in the VHE was detected by MAGIC on the 13 February, and this time lasting three days only, up to the 16 February. In the present work, we will refer to January data (from MJD 57040 to MJD 57050) as phase A and to February data (from MJD 57066 to MJD 57069) as phase B. The observations were performed in wobble mode with a $0.4\deg$ offset and four symmetric positions with respect to the camera center. The observations were performed in 12 nights in January and February 2015 for a total duration of 17.57 h at the zenith range of $40^{\circ}<zd<50^{\circ}$). The corresponding distribution of the squared angular distance $\theta^{2}$ between the reconstructed source position and the nominal source position (points) or the background estimation position (shaded area) for the whole period of observation is shown in Figure 1a. The vertical dashed line shows the value of $\theta^{2}$ up to which the number of excess events and significance are integrated. The cuts used are optimized for steep spectra and low energy thresholds. The analysis energy threshold is $\sim$ 65 GeV, measured as the peak of the energy distribution of a sample of Monte Carlo simulated events for a source with the spectral shape of S5 0716+714. The analysis of the data collected was performed within the standard MAGIC analysis framework MARS Zanin2013_MARS ; Aleksic2015_Upgrade2 .

3.2 Fermi-LAT Analysis

The HE $\gamma$-ray (0.1–300 GeV) observations were obtained in a survey mode by the Fermi-LAT Large Area Telescope. The LAT data were analyzed using the standard ScienceTools (software version v10.01.01) and instrument response function CALDB. Photons in the source event class were selected for the analysis. We analyzed a region of interest of 10*∘* in radius centred at the position of S5 0716+714 using a maximum-likelihood algorithm mattox1996 . We included all the 54 sources of the 3FGL catalogue Acero2015 within 10*∘* in the unbinned likelihood analysis. Model parameters for sources within 5 degrees of the ROI are kept free, while we kept the model parameters for the rest fixed to their catalogue values.

3.3 Multiwavelength Light Curves

The unprecedented outburst of S5 0716+714 was covered in all the wavelengths: such a rich collection of simultaneous data allows us to study with great precision the broadband spectral energy distribution of the source, deeply investigating its nature. Multiwavelength flux light curves of S5 0716+714 during the period from MJD 57010 to MJD 57075 are shown in Figure 2. The top panel shows the MAGIC light curves: in the VHE band, the no-variability hypothesis has been discarded because it resulted in a $\chi^{2}/ndf=67/7$ for January and $\chi^{2}/ndf=22/3$ for the February data for a constant fit. Intranight variability in VHE was not detected. The shadowed areas indicate Phase A (from MJD 57040 to MJD 57050) and Phase B (from MJD 57066 to MJD 57069), the two high states in VHE, and the corresponding activity in the other bands.

When MAGIC discovered S5 0716+714 in 2008 (Anderhub2009 ) the radio band was in a quiescent level, while here, even if delayed with respect to the optical peak, it is possible to discern a response in the radio frequency. The VHE integral flux measured in 2008 by MAGIC was 10 and six times lower with respect to the one detected during the phase A and phase B outbursts, respectively. The two phases show quite different behavior: phase A is characterized by a high state of the source clearly visible in every energy band, especially in optical where the highest flux ever detected for S5 0716+714 is reported. HE and optical curves show a double peaked structure and the polarization angle undergoes a fast change between the two optical peaks. After an intermediate time period when activity is limited to the radio band, phase B is characterized by VHE and X-ray high flux, while the other bands reached a lower state after the outburst of phase A. The fractional variability of both phases shows an increasing trend with the energy. The phase B flare is clearly visible in the VHE and X-ray band only. Such differences in behaviour between the two activity states is very interesting to study, and together with the jet analysis of the source and the modeling of the broadband spectral energy distribution, which will be presented in preparation , it will improve our understanding of the BL Lac emission mechanism Rani2016 .

4 Discussion

The unprecedented outburst of the BL Lac object S5 0716+714 in January 2015 was observed in all the wavebands. The MAGIC detection of the source, together with simultaneous data in the HE range from Fermi-LAT, will allow us to constrain the Inverse Compton part of the broadband spectrum. The rich multiwavelength coverage presented here is part of a more comprehensive study of the source in preparation preparation . The impressive activity in phase A was followed by some days of decreasing flux, and then a new high-state was observed in the VHE and X-ray bands (phase B). The structure of the light curve in phase A showed a double-peaked shape, and the electric vector position angle (EVPA) of the linearly polarized optical emission rotated by $~{}\sim 360^{\circ}$ in between the two peaks. Such swings in EVPA have been frequently observed in connection with $\gamma$-ray flaring epochs (e.g., Marscher2010 ; Aleksic2014 ; Blinov2015 ). A detailed correlation study and modelling of the spectral energy distribution are in progress.

Acknowledgements.

We would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The financial support of the German BMBF and MPG, the Italian INFN and INAF, the Swiss National Fund SNF, the he ERDF under the Spanish MINECO (FPA2015-69818-P, FPA2012-36668, FPA2015-68278-P,FPA2015-69210-C6-2-R, FPA2015-69210-C6-4-R, FPA2015-69210-C6-6-R, AYA2013-47447-C3-1-P, AYA2015-71042-P, ESP2015-71662-C2-2-P, CSD2009-00064), and the Japanese JSPS and MEXT is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia “Severo Ochoa” SEV-2012-0234 and SEV-2015-0548, and Unidad de Excelencia “María de Maeztu” MDM-2014-0369, by grant 268740 of the Academy of Finland, by the Croatian Science Foundation (HrZZ) Project 09/176 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, and by the Polish MNiSzW grant 745/N-HESS-MAGIC/2010/0. The Fermi-LAT Collaboration acknowledges support for LAT development, operation and data analysis from NASA and DOE (United States), CEA/Irfu and IN2P3/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. We thank the Swift team duty scientists and science planners. The Metsähovi team acknowledges the support from the Academy of Finland to our observing projects (numbers 212656, 210338, 121148, and others). The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. The OVRO 40-m monitoring program is supported in part by NASA grants NNX08AW31G, NNX11A043G and NNX14AQ89G, and NSF grants AST-0808050 and AST-1109911. St.-Petersburg University team was supported by research grant 6.38.335.2015 and RFBR grant 15-02-00949. \authorcontributionsMarina Manganaro, Giovanna Pedaletti and Mireia Nievas performed the MAGIC analysis of the source. Bindu Rani, Denis Bastieri, Dario Gasparrini took care of the Fermi-LAT analysis. This proceeding is part of a joint MAGIC-Fermi-LAT work involving also other multiwavelength authors. The preliminary results here included will be part of a paper in preparation which is led by Marina Manganaro, Marlene Doert and Bindu Rani. All the other authors listed provided the analysis of their data and participated in the discussion of the multiwavelength curve. Elina Lindfors is the PI of the MAGIC program of ToO observations of flaring AGNs based on optical, X-ray and $\gamma$-ray triggers. Marina Manganaro wrote this proceeding. \conflictofinterestsThe authors declare no conflict of interest.