Free Access
Issue
A&A
Volume 617, September 2018
Article Number A30
Number of page(s) 15
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201832624
Published online 13 September 2018

© ESO 2018

1 Introduction

Blazars are a subclass of active galactic nuclei (AGN) in which the relativistic jet presents a small viewing angle toward the observer and thus where relativistic effects on the observed emission are more extreme. Conventionally, blazars are subdivided in BL Lac objects and flat spectrum radio quasars (FSRQs) depending on the characteristic of their optical spectrum: while BL Lac objects are dominated by the featureless continuum emission from the jet, FSRQs show wide optical emission lines. The optical spectrum of a fraction of BL Lac objects presents features attributable to the host galaxy.

The blazar S4 0954+65 hosts a black hole of mass MBH ~ 3.3 × 108 M, estimated from the width of the Hα line (Fan & Cao 2004). The detection of the Hα line is not confirmed by Landoni et al. (2015; see the discussion on the redshift determination) so that the mass estimation cannot be confirmed either. This blazar presents strong variability in the optical band, already well studied by Wagner et al. (1990) and by Morozova et al. (2014). Intra night variability has been found both in optical and radio wavelengths (Wagner et al. 1993). The optical high brightness state of February 2015, presented here, is however exceptional for the object, with a brightening of more than three magnitudes in the R-band with respect to the average monitored state1. This not only spurred many alerts in the community (see Astronomer’s Telegram (ATel) #6996, #7001, #7057, #7083, #7093; Carrasco et al. 2015; Stanek et al. 2015; Spiridonova et al. 2015; Bachev 2015; Ojha et al. 2015), but also the first and only detection of the object at very high energies (VHE, E ≳ 100 GeV), thanks to observations by the MAGIC Telescopes. This detection by MAGIC and the multiwavelength (MWL) data collected alongside it are the focus of the present work.

The source GRO J0957+65, detected with the Energetic Gamma Ray Experiment Telescope (EGRET) telescope on board the Compton Gamma-Ray Observatory, has been associated through optical and radio observations with S4 0954+65 by Mukherjee et al. (1995). S4 0954+65 has been afterward always included in the released catalogs of sources detected by the Large Area Telescope (LAT) instrument on board the Fermi satellite (Abdo et al. 2010; Nolan et al. 2012; Ackermann et al. 2013, 2016; Acero et al. 2015; Ajello et al. 2017), with the exclusion of the bright source list released after the first three months of Fermi-LAT data integration.

The classification of the object, based on the available literature, is still unclear. In most of the ATels mentioned above S4 0954+65 is referenced as a FSRQ, but in most of the literature this is classified as a BL Lac object due to the small equivalent width of the emission lines in its spectrum (see, e.g., Stickel et al. 1991). Sambruna et al. (1996) classified the spectral energy distribution (SED) of S4 0954+65 as “FSRQ-like”, in a sample limited to the sources with a detection from EGRET data. It indeed presents a flatter spectral index than most BL Lac objects, in both X-ray and γ-ray bands (see Raiteri er al. 1999; and references therein). Among BL Lac objects, a further phenomenological subdivision can be made based on the frequency of the synchrotron peak, ranging from optical to X-ray frequency and identifying the classes of low-, intermediate-, or high-peaked BL Lac object (LBL, IBL, HBL, respectively). Ghisellini et al. (2011) classified this object as a LBL based on the SED. When including the kinematic features from the radio jet in the classification templates, Hervet et al. (2016) classify this as their kinematic class II, mostly composed of FSRQ. S4 0954+65 can thus be interpreted as a transitional object between FSRQ and classical BL Lac objects.

The most numerous extragalactic sources detected at VHE from Imaging Air Cherenkov Telescopes (IACTs), presently, belong to the HBL class. Therefore the VHE detection of an object such as S4 0954+65 provides a rare opportunity to study VHE emission conceivably produced in a different kind of environment. Indeed, while emission in HBL can mostly be satisfactorily modeled taking into account only processes in a compact feature in the jet, for FSRQs the inclusion of the interactions of such a feature with the surrounding ambient becomes of greater importance (see e.g., Tavecchio 2016). The structure of the broadband SED collected here will also be put in context with other common characteristics of a FSRQ classification, such as intrinsic brightness, peak of the synchrotron component and Compton dominance.

The question of redshift for S4 0954+65 is also still not settled, as claims of line detection in the optical spectrum are not always confirmed. The redshift of the source was first determined at z = 0.368 by the identification of lines by Lawrence et al. (1986, 1996). Stickel et al. (1993) obtained, from different measurements, the same redshift estimate based on line identification. None of these lines were confirmed by the observations reported in Landoni et al. (2015), who instead pose a lower limit of z ≥ 0.45. The latter results were obtained with a superior resolution spectra. At the time of the observation the magnitude in R-band of the object was 15.5, while it is known from variability studies that it could be even two magnitudes lower. In the following we will adopt the redshift z = 0.368.

The outline of this paper is as follows. In Sect. 2, we present the MAGIC telescopes and the relative data set on S4 0954+65. Section 3 reviews all the MWL data that were collected during this exceptional burst, whereas Sect. 4 discusses the implication of this burst for the source state and inner jet structure. Additional information on the MAGIC data analysis, the parameters derived from the radio data, and the full dataset for X-ray data will be found in Appendix A, B, and C respectively.

2 MAGIC observations

The Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) telescopes are an array of two IACTs located in the Island of La Palma (Spain) at an altitude of ~2200 m asl. The system is sensitive down to an energy threshold of E ~ 50 GeV (Aleksić et al. 2016) for low zenith angle observations. This is of particular relevance for the monitoring of variable sources and of those that tend to exhibit a steep spectrum at VHE. The full data have been analyzed using the standard MAGIC analysis chain and the MAGIC Standard Analysis Software (MARS; Zanin et al. 2013; Aleksić et al. 2016).

The MAGIC collaboration supports a program of targets of opportunity (ToO), triggered by MWL monitoring. The ToO program was activated for observations of S4 0954+65 at the end of January 2015 after the first hints of enhanced optical state (triggered by the Tuorla monitoring in R-band, see Sect. 3.3). We observed the source with the MAGIC telescopes for 2 nights (MJD 57049-57050, 2015 January 27 and 28), for a total of 1 h high-quality dark time data, but obtained no detection. We resumed the ToO observations in February after the Tuorla monitoring revealed a very exceptional flux state, later confirmed by other monitoring programs (see Sect. 3.3). We obtained a detection at a significance of ~ 7.4σ from observations during 2015 February 14 (MJD 57067, ATel #7080 Mirzoyan et al. 2015). We continued observing S4 0954+65, barring adverse atmospheric conditions, until full moon days when standard MAGIC observations are not possible due to the elevated level of background light (last day of observation, with already large moonlight contamination, on 2015 March 1, MJD 57082). A detailed breakdown of the observation conditions and relative results can be found in Appendix A.

The total excess from the dark-time data is consistent with a point source emission (see Fig. 1). No other significant emission is found in the field of view apart from the one coincident with S4 0954+65 at the center.

The SED points presented in Sect. 4 below are derived for the day of the flare (MJD 57067, 2015 February 14), using only data taken in dark conditions (that allow for the lowest threshold and lowest systematic uncertainty, Appendix A). We followed the standard MAGIC unfolding procedure (Albert et al. 2007) to obtain the intrinsic spectrum.

The γ-ray emission from sources at high redshift is absorbed via photon–photon pair production on photons from the extragalactic background light (EBL; see e.g., Finke et al. 2010; Domínguez et al. 2011). S4 0954+65 redshift is assumed to be z = 0.368. The spectral shape of the intrinsic emission, that is, after the correction for the EBL absorption, can be fitted with a simple power law (1)

with normalization at E0 = 0.15 TeV and spectral index Γ = 3.98 ± 0.67stat ± 0.15sys. The quoted systematic uncertainties are derived from the standard evaluation in MAGIC data presented by Aleksić et al. (2016). We note that the calculated systematic uncertainty on N0 does not contain the uncertainty on the energy scale, that is about 15%. The unfolded MAGIC spectrum is shown in Fig. 2. The unfolded observed spectrum, that is, without correcting for the EBL absorption, can be described also by a simple power law with at E0 = 0.15 TeV and spectral index Γ = 4.58 ± 0.66stat ± 0.15sys.

thumbnail Fig. 1

Distribution of the squared angular distance (θ2) between the reconstructed event direction and the nominal source direction. The filled histogram is the background estimation, obtained from sky regions within the field of view with similar detector acceptance. We show only data taken in dark condition (condition 1, see Appendix A). The standard MAGIC low energy (LE) cuts are applied to the data (see Appendix A and Table A.1). The vertical line corresponds to the optimal cut (θ2 = 0.02 deg2) for point source analysis in LE cuts, used to derive significance values.

Open with DEXTER
thumbnail Fig. 2

Spectrum for the VHE MAGIC detection. MAGIC data are for flare night only (2015 February 14, MJD 57067.14). Violet filled circles are for the unfolded observed points, while open circles are deabsorbed for EBL absorption (EBL model by Domínguez et al. 2011). The solid line is the fit for the observed points and the dashed line is the fit for the de-absorbed ones, with details in the text.

Open with DEXTER

3 The multiwavelength coverage

All the data presented in this section are collected to produce the light curves and SED, whose interpretation is later presented in Sect. 4.

3.1 Fermi-LAT

The LAT on board the Fermi satellite scans the entire sky every 3 h. From the data of the first four years of operation, S4 0954+65 was detected with an average significance of 27.2σ in the energy range from 100 MeV to 300 GeV as reported in the Fermi-LAT Third Source Catalog (Acero et al. 2015; 3FGL). A dedicated analysis from MJD 56952 (2014 October 22) to MJD 57208 (2015 July 05) is presented in this work. We selected Pass 8 source class events within a 10° circular region centered on the position of S4 0954+65, in the energy range 0.1–500 GeV. The spectral analysis was performed through an unbinned likelihood fit, using the ScienceTools software package version v11-05-00 along with the instrument response functions P8R2_SOURCE_V6. The model of the likelihood fit includes a Galactic diffuse emission model and an isotropic component2. In addition, we included the sources in the 3FGL catalog within a 20° circular region centered on S4 0954+65. The spectral indexes and fluxes of the 3FGL sources located within a region of 10° from S4 0954+65 were left free to vary, while the sources in the region from 10° to 20° were fixed to their catalog values. The results were obtained from two iterations of maximum-likelihood analysis, after the sources with a test statistics (Mattox et al. 1996) TS < 10 were removed. The strongest source located beyond 10° from S4 0954+65 is at an angular distance of 10.8°. This source has a variability index of 42.4 in the 3FGL catalog, that allows us to treat it as a nonvariable source and thus to fix its spectral index and flux to the values reported in the 3FGL catalog.

The light curve was calculated in day timescale bins, modeling the source with a single power-law spectrum (as it is also described in the 3FGL). Both the flux and spectral index of S4 0954+65 were left free during the likelihood fits, while the rest of the point sources were fixed and only the diffuse Galactic and isotropic models were allowed to vary. In case of TS < 4, an upper limit on the flux was calculated fixing the spectral index to 2.38 as given in the 3FGL catalog. The results are shown in Fig. 3. The figure also shows the light curve calculated in a 15-day bin as comparison. The light curve was obtained with the same procedure described above for the 1-day binning. During the HE flare in November 2014 (MJD 56976, ATel #6709; Krauss 2014) the LAT spectral index is compatible with its 3FGL value of 2.38 ± 0.04, averaged from four years of data. Moreover, the visibility of the source by MAGIC was at an unfavorable zenith angle of 60° (implying a high energy threshold). Therefore, no ToO observation was activated with MAGIC for this flare. MAGIC observations were activated later on during the strong flare on February 2015 when the LAT detected a hardening of the spectrum as shown previously by Tanaka et al. (2016) where the LAT analysis using Pass 7 reprocessed data is presented.

The spectral analysis for the MWL SED corresponds to 1-day integration centered in the MAGIC observation (MJD 57067.14, 2015 February 14). From a first likelihood fit we found the best spectral fit was a power-law spectral index of 1.87 ± 0.09 (significantly harder than its average 3FGL value) and was fixed in the model for the spectral points calculation. Moreover, all the sources included in the model except the diffuse Galactic and isotropic models were also fixed. The source was detected during this period with a TS of 379.7. A curved spectral model is not significantly favored in this day (TS for a log parabola fit is TS LP = 380.10 to be compared with a simple power law fit with TS PWL = 379.74).

3.2 Swift dataset

The 22 multi epochs event-list obtained by the X-ray Telescope (XRT; Burrows et al. 2004) on board the Neil Gehrels Swift Observatory satellite in the period of 2014 November 17 (MJD 56978.96395) to 2015 March 11 (MJD 57092.26632) with a total exposure time of ~11.12 h were processed using the procedure described by Fallah Ramazani et al. (2017). All these observations had been performed in photon counting (PC) mode, with an average integration time of 1.8 ks each. The equivalent Galactic hydrogen column density is fixed to the value of nH = 5.17 × 1020[cm−2] (Kalberla et al. 2005).

The average integral photon X-ray flux (0.3–10 keV) in this period is 1.64 × 10−11erg cm−2 s−1. The X-ray flux is peaking at MJD 57070.76523 with F(0.3−10keV) = 3.18 × 10−11erg cm−2 s−1 which is a factor of about two higher than the average flux of the analyzed period. The average flux outside the flare period (2006–2015) is F(0.3−10keV) = 4.3 × 10−12erg cm−2 s−1, that we derived from a sample of XRT data comprising 25 X-ray exposures in the XRT database, not including the 22 multi epochs event-list described above. This indicates that the source was clearly in its X-ray high state during the VHE γ-ray detection. The X-ray spectral index during the analyzed period varies between 1.15 ± 0.06 ≤ ΓX ≤ 1.82 ± 0.1. It is notable that the softest spectral index was obtained a night prior to the VHE γ-ray flare while the spectra starts to harden after 2015 February 14 and reach its historical hardest spectra ten days after the VHE γ-ray flare. The X-ray spectra on the night before and after the VHE γ-ray flare can be well described with a power-law with spectral index of ΓX, Feb. 13 = 1.82  ± 0.05 (χ2 ∕d.o.f. = 1.024/41) and ΓX, Feb. 15 = 1.49 ± 0.07 (1.025/24 χ2 ∕d.o.f.) respectively.The full dataset analysis is given in Appendix C.

The Swift satellite hosts an additional instrument, the Ultraviolet/Optical Telescope (UVOT; Poole et al. 2008). The data taken during the period of interest for this work have already been presented by Tanaka et al. (2016). They follow the behavior of the optical light curve that we present below. Therefore, they are not reproduced again nor shown in Fig. 3. The UVOT bands are however important for the SED modeling presented in Sect. 4 and will therefore be included there for MJD 57067 (2015 February 14, day of the VHE detection). The dataset presented by Tanaka et al. (2016) suffers from an incorrect exposure calculation by a factor of two, related to the deadtime correction, and thus a lower reconstructed flux. We, therefore have performed a reanalysis here for the two exposures taken with UVOT on MJD 57066.76. Data reduction has been done on all the available filters (v,   b,   u,  w1,  m2,  w2), following the standard UVOT data analysis prescriptions3. We present both exposures separately, due to the high variability in this night (e.g., for the V -band there is a variation of ~0.3 magnitudes in ~1.5 h).

3.3 The optical domain

Optical data were collected with: 35 cm KVA telescope (La Palma Island, Spain) used in the Tuorla monitoring program; 1.8 m Perkins telescope of Lowell Observatory (Flagstaff, Arizona); 70 cm telescope AZT-8 at the Crimean Astrophysical Observatory (Nauchny, Russia); 40 cm telescope LX-200 of St. Petersburg State University (St. Petersburg, Russia); IAC80/Camelot at the Teide Observatory (Tenerife, Spain). The data analysis from KVA was performed with the semi-automatic pipeline using the standard analysis procedures (Nilsson et al. in prep). The differential photometry was performed using the comparison star magnitudes from Villata et al. (1997). For the Perkins telescope see Jorstad et al. (2010) and references therein. The details of observationsand data reductions with AZT-8 and LX-200 are given by Larionov et al. (2008). IAC80/Camelot data were automatically processed by the pipeline Redcam and calibrated astrometrically using XParallax, both available at the telescope. Instrumental magnitudes for IAC80/Camelot data were extracted using Sextractor (Bertin & Arnouts 1996) and calibration of the source magnitude was obtained with respect to the reference stars provided by Raiteri er al. (1999).

All the telescopes mentioned above provide R-band photometry. We have applied the calibration of Mead et al. (1990) for all optical measurements to transform magnitudes into flux densities, and dereddened the flux according to the absorption by Schlafly & Finkbeiner (2011). The host galaxy is not detected for this object.

From the Perkins, AZT-8+ST7 and LX-200 telescopes we collect also polarization information. In Fig. 3 we show the optical photometry data and time evolution of the fractional linear polarization and the electrical vector position angle (EVPA) in R-band. The EVPA measurements have been arranged such to minimize the impact of the ± 180° ambiguity, that is, adding or subtracting 180° whenever two subsequent measurements differ by more than 90°.

In the same timeframe of the VHE detection and the optical flare, a substantial change in the optical EVPA can be identified (see Fig. 3). The EVPA rotation starts just before the optical and VHE flare and reaches a total change of roughly 100°. The optical flare in February 2015 is a factor of about 3 larger in flux than the 2011 flare (see Morozova et al. 2014), that was already exceptional and concurrent with a series of γ-ray flares evident in Fermi-LAT data. During the most extreme flare in 2011, the EVPA rotated by about 300°.

thumbnail Fig. 3

MWLlight curves and polarization evolution of S4 0954+65 ranging from MJD 56970 (2014 November 9) to MJD 57200 (2015June 27). The energy range of each panel and the corresponding instrument can be found in the legend. Please refer to the text for details on the data taking and reduction for each instrument.

Open with DEXTER
thumbnail Fig. 4

Sequence of total (contours) and polarized (segments) intensity images of S4 0954+658 at 43 GHz, convolved with a beam of 0.24 × 0.15 mas2 at PA = –10°. The global total intensity peak is 1606 mJy beam−1 and the global polarized intensity peak is 104 mJy beam−1; black line segments within each image show the direction of polarization and their length is proportional to the polarized intensity. The black horizontal line indicates the position of the core, A0, and gray, blue, and red circles show the locations of knots K14a,b and K15, respectively. The size of the circles is proportional to the estimated average size in each epoch. The detailed characteristics of the knots can be found in Table 1 and B.1.

Open with DEXTER
Table 1

Characteristics of the new radio knots observed from the jet of S4 0954+65.

3.4 Radio and millimeter ranges

The source S4 0954+65 was monitored at 3.5 mm (86 GHz) and 1.3 mm (229 GHz) wavelengths from the IRAM 30 m Millimeter Radiotelescope under the Polarimetric Monitoring of AGN at Millimeter Wavelengths4 (POLAMI) program. The program monitors the four Stokes parameters of a sample of the brightest 40 northern blazars with a cadence better than a month (see Agudo et al. 2018a,b; Thum et al. 2018). Results from the observations are presented in Fig. 3. The data reduction, calibration, and flagging procedures were described in detail by Agudo et al. (2018a), see also Agudo et al. (2010, 2014). Figure 3 includes also the 1.3 mm flux density data that were obtained at the Submillimeter Array (SMA) located in Hawaii. S4 0954+65 is included in an ongoing monitoring program at the SMA to determine the fluxes of compact extragalactic radio sources that can be used as calibrators at millimeter wavelengths (Gurwell et al. 2007). Observations of available potential calibrators are from time to time observed for 3–5 min, and the measured source signal strength calibrated against known standards, typically solar system objects (Titan, Uranus, Neptune, or Callisto). Data from this program are updated regularly and are available at the SMA website5. The largest flux in the considered period is at MJD 57072-57076, showing an increase of the flux between 1 and 3 mm wavelengths. We note, however, the lack of exactly simultaneous data to the MAGIC peak detection (MJD 57067).

S4 0954+65 is monitored monthly by the Boston University (BU) group with the Very Long Baseline Array (VLBA) at 43 GHz within a sample of bright γ-ray blazars through the VLBA-BU-BLAZAR program6. The VLBA data are calibrated and imaged in the same manner as discussed by Jorstad et al. (2005, 2017). The VLBA imaging monitoring program allows us to study the kinematics of the inner jet at pc scale. The inner jet has been monitored also for months after the VHE flare (see Fig. 4). In addition to the stable core at mm wavelengths (dubbed A0, see Fig. 4) it was possible to identify the emergence of three new knots whose characteristics are tabulated in Table 1. The nomenclature of the knots follows in sequential order from the beginning of the VLBA monitoring program. Previous knots characteristics can be found in Morozova et al. (2014).

Of particular interest is knot K15, which is very compact, with a FWHM average size of 0.05 ± 0.01 mas and presents the largest apparent speed of (25.27 ± 1.20)c, cf. Fig. 5. The zero-epoch separation of this knot is consistent with the VHE flare considering its 18-day uncertainty. The intensity of the core is increasing in the epoch of MJD 57067 observation, but no significant change in the core polarization can be appreciated. The detailed information on the time evolution of the radio knot can be found in Table B.1, while the polarization evolution details are shown in Table B.2. No new knot appears in November 2014, when the source was high in the HE band as observed by Fermi-LAT but without optical enhancement. The zero epoch-separation from the core of knots K14a,b are not coincident with the high state in Fermi-LAT data of November 2014, but happen months before. We analyzed Fermi-LAT data for the period included within the error band for K14a,b zero epoch-separation and found no particular enhancements.

We also note the position angle of K15 with respect to the core, (PA = − 5.9° ± 1.9°). This is different than the values reconstructed from previous knots, ranging from roughly PA = − 15° to PA = − 25° in Morozova et al. (2014), that are in turn consistent with the values for K14a,b. The mean jet direction is at PA ≃−20°. A difference in PA and in apparent speed could be simply related to a small difference in the angle to the observer. However, the highest apparent speed can be used to estimate the Doppler factor, considering the upper limit to largest possible viewing angle θobs < arcsin(1∕βapp) and ultimately leading to δapp ~ βapp. Applying this to the above mentioned knots (averaging the apparent speed to βapp ~ 13c for K14a/b): θobs,K15 < 2.3° and δapp,K15 ~ 25; θobs,K14 < 4.4° and δapp,K14 ~ 13.

The 37 GHz observations were made with the 13.7 m diameter telescope at Aalto University Metsähovi Radio Observatory. A detailed description of the data reduction and analysis is given by Teraesranta et al. (1998). The error estimate in the flux density includes contributions from the measurement RMS and the uncertainty of the absolute calibration. The S4 0954+65 observations were done as part of the regular monitoring program and the GASP-WEBT campaign. There are no strictly simultaneous 37 GHz data to the MAGIC detection, however an increase in flux can be seen when comparing observation taken 1 day before (2015 February 13, MJD 57066.15, Fν = 1.27 ± 0.07 Jy) and one day after the MAGIC detection (2015 February 15, MJD 57068.15, Fν = 1.65 ± 0.09 Jy).

The Owens Valley Radio Observatory (OVRO) 40 m uses off-axis dual-beam optics and a cryogenic pseudo-correlation receiver with a 15.0 GHz center frequency and 3 GHz bandwidth. Calibration is achieved using a temperature-stable diode noise source to remove receiver gain drifts and the flux density scale is derived from observations of 3C 286 assuming the Baars et al. (1977) value of 3.44 Jy at 15.0 GHz. The systematic uncertainty of about 5% in the flux density scale is not included in the error bars. Complete details of the reduction and calibration procedure are found in Richards et al. (2011). The long-term monitoring program at OVRO monitors the variability of this source at 15 GHz over a longer time than what shown here. While it is obvious that the source was variable also during February 2015, it is not an exceptionally bright flux state of the source in the radio band. From a decade long monitoring, the source shows brighter levels (highest at F15 GHz = 2.53 Jy) and fainter levels (lowest at F15 GHz = 0.85 Jy).

Both 15 and 37 GHz data seem to be in agreement with the behavior seen from mm wavelength data. Again we note the lack of strictly simultaneous data to the MAGIC peak detection (MJD 57067).

thumbnail Fig. 5

Apparent distance from the radio core A0 of the new emerging knots, K14a,b and K15, as a function of time. The images fromwhich the apparent distances are calculated can be found in Fig. 4. We note that the K15 knot presents the highest apparent speed.

Open with DEXTER

4 Discussion

The coverage of flaring states at VHE is helpful in understanding jet dynamics. We present a discussion of the SED for the day of the flare (2015 February 14). We did not attempt SED modeling for other days, for which the MAGIC data would provide only nonconstraining upper limits to emission at VHE. The day of the VHE detection is instead put in context with a longer time span behavior in the MWL dataset. However, the VHE sampling of the state is too scarce to attempt a numerical correlation study of the light curves.

4.1 Light-curve phenomenology

The MWL light curves of the source for all the instruments involved in the present work are reported in Fig. 3, and cover a period of seven months, from MJD 56970 (2014 November 19) to MJD 57200 (2015 June 27). The red region in Fig. 3 indicates the time window where the knot K15 was ejected in the VLBA analysis, as reported in Table 1: a time range of 36 days centered in MJD 57081 (2015 February 28). The VHE detection and the enhanced activity in the other bands are found inside the K15 ejection time window, making this event important for the understanding of the whole scenario. The spectral index at HE as inferred from the Fermi-LAT data is harder than the average spectral index of Γ = 2.38 ± 0.04 from the 3FGL catalog dataspan. In the presented timeframe, the X-ray emission peaks around the observation on MJD 57070.76434(2015 February 17), with a delay with respect to the detection in VHE. The approximately 3 hours of observationsin VHE in the same night did not lead to a detection (see Table A.1). However, during the period of enhanced MWL activity, there is a clear hardening of the X-ray spectrum. Hardening at both X-ray and γ-ray energies points toward the emergence of a new component in the nonthermal spectrum.

The optical band is very bright during the VHE detection, reaching peaks of more than 20 mJy of flux density when the average behavior of the source isfound around a few mJy (see the optical monitoring from Tuorla observatory). The optical emission is polarized by a fraction of ≳ 10% and the polarization angle rotates by ~100° during the flare: Blinov et al. (2015) have shown that from a systematic monitoring (Robopol monitoring) of both γ-ray loud and γ-ray quiet sources, only the former class of object displays polarization angle rotation similar to the one seen here for S4 0954+65. Blinov et al. (2015) studied the change of EVPA as a function of time for smooth changes of > 90°. Requesting the same smoothness requirements, no smooth rotation of >90° can be identified in the dataset presented here, see Fig. 6. A variation of ΔEVPA ≃ 105° that is not smooth can, however, be identified between MJD 57060 and MJD 57075. This variation would imply a change of the EVPA curve slope of ΔEVPAΔt = 7 deg day−1, compatible with the bulk of the variations studied by Blinov et al. (2015). The rotations of the polarization angle are often physically linked to high flaring states of the objects in the γ-ray band. While individual occurrences of γ-ray flares and rotations cannot be firmly linked to each other, there is a low probability that all the occurrences are due to chance coincidence (from MonteCarlo simulations in Blinov et al. 2015). This hypothesis is still confirmed from 3 years of Robopol monitoring data in Blinov et al. (2018). Kiehlmann et al. (2017) also study whether a simple stochastic variation can account for the observed rotations in the Robopol monitoring. While their model is failing to recover all the observational characteristics in the monitoring, it also highlights a larger discrepancy from the expectations of stochastic model with respect to the occurrence of large variations of EVPA (>90°), however not significant. Smooth variations seem also to be more firmly linked to deterministic processes and not to a random walk effect (Kiehlmann et al. 2016). Robopol monitoring data are also used in Angelakis et al. (2016), to study the difference in the amount of polarization seen on average in γ-ray loud and γ-ray quiet sources. The median fraction variability of the S4 0954+65 dataset presented here is 16.4%. This value can be compared with the average 10% for the γ-ray loud subset of the Robopol monitoring and a value of 17.1% for S4 0954+65 computed for the observations on year 2013 and 2014. According to the interpretation by Angelakis et al. (2016), a higher fractional polarization is also expected in LSP/ISP blazars, due to the fact that in such sources the optical synchrotron emission relates to the peak synchrotron emission. Therefore, the particles associated with this emission are the most energetic, with faster cooling and thus probing a small volume of the emission region near the acceleration region, where it is expected to have a stronger ordered (helical) magnetic field, leading to higher polarization fraction.

Images at 43 GHz show the emergence of new knots. In Morozova et al. (2014), a series of optical flares of S4 0954+65 in 2011 are studied, and the emission of knots is found correlated to the simultaneous flaring of the optical and HE bands. The maximum flux in the 2011 state is a factor of 3 lower in optical than the state presented here. The polarization fraction in this 2011 flare was similar to that seen in the present work. In Morozova et al. (2014) the chance coincidence of high optical state and knot emission has very low probability.

The phenomenology of the 2015 flare described here agrees very well with the model put forward by Marscher et al. (2008) and applied to the S4 0954+65 dataset of Morozova et al. (2014). In that model, the flare is due to a newly appearing knot accelerating at the base of the jet and propagating through an helical flow streamline. The helical streamline can be expected due to the anchoring of the accelerating flow to the rotating base of the accretion disk or black hole magnetosphere, depending on modeling. The magnetic field topology in the jet is also helical and ordered. Geometrical effects and the propagation through the helical magnetic field account for the rotation of the EVPA.

In Zhang et al. (2014), a model is proposed where the EVPA rotation is also related to the propagation through an helical magnetic field, but the streamline of propagation is not necessarily helical itself. In this model the magnitude of the swing can depend on the assumptions on the settings for the flare, specifically the magnetic field strength and orientation, the acceleration efficiency and the continuous injection of freshly accelerated particles.

The model described in Marscher et al. (2008) allows the emission at radio wavelengths in a flaring state which is not simultaneous with the VHE flare. In this scenario the radio activity could be delayed several days, even months, with respect to the VHE detection. This is expected if synchrotron self absorption is involved, and hence the emission region is located closer to the central engine than the radio core (A0 in Fig. 4). The peak of radio emission is expected to be lagging behind and appear when the disturbance has propagated further down the jet, where the absorption is not an issue. The X-ray emission peak, then, could also be delayed with respect to the optical outburst. As the X-ray emission is probably due to IC of an external soft photon field by electrons in the jet (see above), the X-ray variability traces both the accelerated particle distribution and a change in the soft photon field. This retraces similar interpretation drawn for flares of other sources where the dataset was however richer and more detailed (Marscher et al. 2008, 2010; Aleksić et al. 2014; Ahnen et al. 2017b).

thumbnail Fig. 6

Light curves for R-band polarimetry of S4 0954+65. Please refer to the text for details on the data taking and reduction for each instrument.

Open with DEXTER

4.2 Emission model for the flare SED

The SED of blazars are dominated by their nonthermal emission and can usually be described by two broad components. The low energy nonthermal emission is explained as synchrotron emission, while the high energy emission is most commonly modeled through inverse Compton (IC) emission, where soft photons are upscattered to γ-ray energies by electrons within the jet emitting region. The origin of the soft photon field itself can vary for different blazar subclasses. In particular, for most of the classical BL Lac objects, the VHE emission can be reasonably modeled through synchrotron self-Compton emission (SSC, see e.g., Rees 1967; Maraschi et al. 1992). Instead, for the case of FSRQs, the modeling of the emission usually requires the inclusion of external soft photon fields from, for example, the infrared dusty torus or the optical-ultraviolet emission from the broad line region (BLR) for the IC process (see e.g., Tavecchio 2016).

A broadband SED is compiled for 2015 February 14 (MJD 57067). We collect, from the MWL sample described in Sect. 3, the data closest in time to the MAGIC observation. Fermi-LAT data points are obtained from a 1-day integration centered on the MAGIC observation. The specific dates of other wavelength observations are given in the caption of Fig. 7.

Tanaka et al. (2016) model the SED of S4 0954+65 during a similar integration time as the 2015 flare studied in this work. The data shown in Fig. 7 include, in addition to what is shown by Tanaka et al. (2016), the VHE data from the MAGIC observation, the AZT-8+ST7 and POLAMI data. Moreover, the Fermi-LAT data are reanalyzed as described in Sect. 3 to be centered at the MAGIC observation time and benefit from the latest Fermi-LAT Pass 8. The Swift-XRT and Swift-UVOT data are also reanalyzed for this work.

Tanaka et al. (2016) report that a SSC modeling of the data is challenging, requiring very low magnetic field (B ~ 1 μG in contrast to the B ~ 1 G expected in blazar jet components). Alternatively, an External Compton (EC) modeling was able to reproduce the data. In their model, the soft photon field for the EC model was the dusty torus from the source. In Fig. 7, we plot the model from Tanaka et al. (2016). This model reproduces the Fermi-LAT and MAGIC data, although their paper did not include any MAGIC data. However, the model fails to reproduce properly the optical observations. Such underestimation at optical frequencies in the model of Tanaka et al. (2016) is driven by a misreconstruction of the UVOT fluxes, explained in Sect. 3. With the reanalyzed UVOT dataset presented here, we use a new model, using the same code and most of the same assumptions as in Tanaka et al. (2016), including a redshift of z = 0.368. The code is explained in detail in Finke et al. (2008) and Dermer et al. (2009). We note that the presented SED model curves already include the effect of EBL absorption, meaning that the intrinsic emission is absorbed according to the EBL model by Finke et al. (2010). The new EC model provides a good description of the MWL data and is shown in Fig. 7. The parameters of both models are reported in Table 2. The break in the underlying electron population is similar to what expected by classical cooling, with the slope of the electron distribution before of the break (s1) and after the break (s2) differing by s2s1 = 1.2. We also tested a model with a larger break (s2s1 = 1.6) and longer variability timescale (105 s) which was better able to reproduce the UV band data. However, the larger s2s1 disagrees with expectations from synchro-Compton losses. Further, this model underestimates the MAGIC data and overestimates the emission at 1 mm. We therefore do not present this model here. The use of VHE spectral information is crucial to model the falling part of the high energy peak of the blazars SEDs, which is crucial to constrain the most energetic electrons within the leptonic framework scenario (SSC and EC models).

As mentioned in the introduction, the classification of a blazar can be aided by the study of its SED characteristics. According to the SED model presented above, the peak of the synchrotron emission is at νsyn ~ 8 × 1014 Hz, makingit an intermediate synchrotron peaked BL Lac object (Ackermann et al. 2015)7. The Compton dominance, calculated comparing the luminosity at the peak of the synchrotron emission to that of the IC peak, is LICLsyn ~ 7. Such Compton dominance value is at least 3.5 times the values obtained by Finke (2013) for long-term blazar studies.

thumbnail Fig. 7

Spectral energy distribution for the VHE MAGIC detection. Red symbols are strictly simultaneous to the VHE detection, blue symbols are for data taken during the same day and black symbols are for the closest observations. MAGIC spectral data (red circles) are for flare night only (2015 February 14, MJD 57067.14). Red filled circles are for the unfolded observed data points. The red shaded band shows the region of additional systematic uncertainty. Fermi-LAT data are the Pass 8 data for 2015 February 14 (1-day integration centered on the MAGIC observation, blue squares). Swift-XRT data are for 2015 February 13 (MJD 57066.70992, blue squares). Swift-UVOT data are given for the twoseparate exposure taken on 2015 February 13 (MJD 57066.76, blue triangles and dark blue triangles). R-band data are for 14th Feb (Tuorla, MJD 57067.16375 and AZT-8+ST7 MJD 57067.1, red diamonds). POLAMI data are for the 18th February (MJD 57071.5, black stars at 100 and 300 GHz). OVRO data are for 2015 February 10 and 19 (black circles at 15 GHz). Metsähovi for 2015 February 13 and 15 (MJD 57066.15, MJD 57068.68, black squares at 37 GHz). The gray data are for NED (light) and SSDC (dark) SED historical data points. The model from Tanaka et al. (2016; gray dashed curve) as well as the model presented here (black solide curve) include an emission component from synchrotron plus inverse Compton on a dusty torus (see text for details and Table 2 for the values of the physical parameters). The effect of the EBL attenuation is included in the modeling using the model by Finke et al. (2010) and a redshift of z = 0.368.

Open with DEXTER
Table 2

SED model parameters.

5 Conclusions

The census of extragalactic objects that present VHE emission is still limited. We present here the first detection at VHE of the blazar S4 0954+65 obtained through observations with the MAGIC Telescopes. The observations were conducted during an exceptional flare of the source in February 2015, originally identified in the optical band. We collected MWL simultaneous data to better characterize the state of the source.

The HE emission is also found in elevated state from the analysis of Fermi-LAT data, which reveal the hardest state of the HE emission to be concurrent with the detection at VHE. The X-ray emission peak is delayed by few days with respect to the VHE detection and shows a trend of spectral hardening during the period presented here. The radio and mm wavelength emission reveal a moderate elevation of the flux, that is however not exceptional in the long term behavior of the source.

The source is classified in the literature as a BL Lac, but we have shown here that it presents similarities with the FSRQ class. Results from the monitoring of optical polarization and 43 GHz jet component analysis were compared to archival observation of S4 0954+65 and of statistical behaviour of other sources. Three main measurements were considered: the day of the VHE detection of S4 0954+65 is included in the error box for the zero epoch separation of knot K15; the optical polarization fraction is increasing in the same period; a rotation of optical EVPA of ~ 100° can be identified, also in the same period, possibly related to the helical structure of the magnetic field in the acceleration region. We discuss how these measurements point to a common behavior with ISP/LSP sources. Both, the best emission model (EC on dust torus) and the MWL light curve behavior show points of contact with other sources that are either clear FSRQ (like PKS 1510-089) or are transitional objects (like BL Lac itself). This is also supported from the moderate Compton dominance in the SED model presented here and the fact that the synchtrotron peak show that the source can be classified as ISP source. The work presented here reiterates the importance of VHE γ-ray and detailed MWL studies of blazars during different flux states to test their intrinsic characteristics and shed light on the physical processes taking place within their jets.

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 ERDF under the Spanish MINECO (FPA2015-69818-P, FPA2012-36668, FPA2015-68378-P, FPA2015-69210-C6-2-R, FPA2015-69210-C6-4-R, FPA2015-69210-C6-6-R, AYA2015-71042-P, AYA2016-76012-C3-1-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 the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, the Polish National Research Centre grant UMO-2016/22/M/ST9/00382 and by the Brazilian MCTIC, CNPq, and FAPERJ. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK), and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. This work performed in part under DOE Contract DE-AC02-76SF00515. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Part of this work is based on archival data, software, or online services provided by the Space Science Data Center - ASI. 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 The research at Boston University was supported by NASA Fermi Guest Investigator program grant 80NSSC17K0694 and US National Science Foundation grant AST-1615796. The VLBA is an instrument of the Long Baseline Observatory. The Long Baseline Observatory is a facility of the National Science Foundation operated by Associated Universities, Inc. This paper is partly based on observations carried out with the IRAM 30 m Telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). IA acknowledges support by a Ramón y Cajal grant of the Ministerio de Economía, Industria y Competitividad (MINECO) of Spain. The research at the IAA–CSIC was supported in part by the MINECO through grants AYA2016–80889–P, AYA2013–40825–P, and AYA2010–14844, and by the regional government of Andalucía through grant P09–FQM–4784. St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. 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.

Appendix A Additional information on MAGIC data reduction

Table A.1

MAGIC data summary for the observation of S4 0954+65 from 2015 January 27 to 2015 March 1.

The MAGIC telescopes are supported by an extensive weather monitor program. Atmospheric transmission at different heights within the MAGIC field of view is obtained with the use of a LIDAR (for details on this see Fruck & Gaug 2015). For data quality selection we consider the transmission measured at a height of 9 km, with T9 km = 1 representing a perfectly clear sky and T9 km = 0 a complete opacity. MAGIC can carry out observations also during partial moonlight, with the drawback of having a higher energy threshold and larger systematic errors due to a higher contamination from the elevated night sky background (NSB), see Ahnen et al. (2017a). The brightness of the NSB can be monitored from the average current in the camera (DC). S4 0954+65 was observed in a zenith range ranging from 35° to 50°, for a total of 12.5 h of data, of which approximtely one hour was lost due to bad weather. In the following we refer to the different observation conditions of our data set as follows:

  • good dark data: data taken with dark sky (DC < 1.5 μA) and good atmospheric condition (T9 km > 0.85), used for detection and spectral reconstruction;

  • dark data needing atmospheric correction: data taken with dark sky (DC < 1.5 μA) but under nonoptimal weather conditions (0.55 < T9 km < 0.85), used for detection and spectral reconstruction after atmospheric correction;

  • good low moon data: data taken with elevated NSB due to moonlight (1.5 μ < DC < 4 μA) and good atmospheric condition (T9 km > 0.85), used only for detection in this particular dataset;

  • good moon data: data taken during high NSB due to moonlight (DC > 4 μA) and good atmospheric condition (T9 km > 0.85), used only for detection.

The subsample of dataset selected with condition (1) (9.48 h of good quality data) has been analyzed with the standard MAGIC analysis chain (Zanin et al. 2013). The subsample of dataset selected with condition (2) (0.32 h of data) follows the same analysis chain until the estimation of the energy for the events and evaluation of the flux. For this last step, the estimated energy and the effective area are corrected taking into account the enhanced atmospheric absorption (for validation of the procedure see Fruck & Gaug 2015).

The subsample of dataset selected with condition (3) is applicable only at the day of 14th February, with the first VHE detection. The detection can be claimed from dark data alone (i.e., selected with condition 1), but an extra 0.81 h of data were taken under low moonlight. The data are presented here for completeness, but are not used for spectral reconstruction so as not to increase the systematic error and energy threshold.

The subsample of dataset selected with condition (4) (0.94 h of data) requires a special analysis that takes care of the effect of moonlight on data taking, reconstruction and analysis. Details of the procedure can be found in Ahnen et al. (2017a).

The detailed breakdown of significances and estimated VHE fluxes is given in Table A.1. Numbers are presented for the so-called low energy (LE) and full range (FR) cuts. The LE cuts are optimized for an energy range of E ≳ 100 GeV and are particularly appropriate for steep spectrum sources, while FR cuts are optimized for an energy range of E ≳ 250 GeV. The cuts are applied on 2 parameters: the “size” parameter, integrated charge (in photoelectrons) in the cleaned shower image; the “hadronness” parameter, computed from the gamma-hadron separation random forest (RF), with a value ranging from zero for the most gamma-like images to one for the most hadron-like images. Indeed the standard MAGIC analysis chain relies on RF techniques to discriminate among gamma and hadronic shower and to better reconstruct the event directions. Lookup tables are used for energy estimation. This is achieved starting from a parametrization of the shower images in the detector. The significance of signal is then calculated with Eq. (17) from Li & Ma (1983) and using five regions of equal size and distance to the center of camera as the signal region for background estimation. Fluxes are calculated above an energy threshold of 150 GeV, which corresponds to the peak of the differential energy distribution of the excess events as a function of estimated event energy. The high energy threshold is due to the high zenith angle of the observation. We note that for data of condition 4, strong moon, we apply an additional minimum cut in the “size” parameter (“size” > 150 phe) of the reconstructed Cherenkov image as prescribed by the moonlight-adapted analysis. This increases the energy threshold to a value of ~250 GeV. In case of nondetection, we provide 95% confidence level upper limits to the flux, calculated following Rolke et al. (2005), considering a systematic error on flux estimation of 30% (Aleksić et al. 2016).

Appendix B Additional VLBA derived parameters

Table B.1

Time evolution of characteristics of the new radio knots observed from the jet of S4 0954+65.

The detailed information on the time evolution of the radio knot can be found in Table B.1, while the polatization evolution details are shown in Table B.2.

Table B.2

Time evolution of polarization parameters (percentage and angle) for the core A0 observed from the jet of S4 0954+65.

Appendix C Swift-XRT full dataset

Table C.1

S4 0954+65 Swift-XRT exposures.

Table C.1 collects all the analyzed exposures for the Swift-XRT dataset described in Sect. 3. Fluxeshave been extracted from a 20 pixel circular aperture. A different aperture was used on 2015 February 17 (MJD 57070.76), due to pile-up effects.

References


7

Intermediate-synchrotron-peaked blazar (ISP) are defined with rest-frame synchrotron peak frequencies of 1014Hz < νsyn < 1015Hz.

All Tables

Table 1

Characteristics of the new radio knots observed from the jet of S4 0954+65.

Table 2

SED model parameters.

Table A.1

MAGIC data summary for the observation of S4 0954+65 from 2015 January 27 to 2015 March 1.

Table B.1

Time evolution of characteristics of the new radio knots observed from the jet of S4 0954+65.

Table B.2

Time evolution of polarization parameters (percentage and angle) for the core A0 observed from the jet of S4 0954+65.

Table C.1

S4 0954+65 Swift-XRT exposures.

All Figures

thumbnail Fig. 1

Distribution of the squared angular distance (θ2) between the reconstructed event direction and the nominal source direction. The filled histogram is the background estimation, obtained from sky regions within the field of view with similar detector acceptance. We show only data taken in dark condition (condition 1, see Appendix A). The standard MAGIC low energy (LE) cuts are applied to the data (see Appendix A and Table A.1). The vertical line corresponds to the optimal cut (θ2 = 0.02 deg2) for point source analysis in LE cuts, used to derive significance values.

Open with DEXTER
In the text
thumbnail Fig. 2

Spectrum for the VHE MAGIC detection. MAGIC data are for flare night only (2015 February 14, MJD 57067.14). Violet filled circles are for the unfolded observed points, while open circles are deabsorbed for EBL absorption (EBL model by Domínguez et al. 2011). The solid line is the fit for the observed points and the dashed line is the fit for the de-absorbed ones, with details in the text.

Open with DEXTER
In the text
thumbnail Fig. 3

MWLlight curves and polarization evolution of S4 0954+65 ranging from MJD 56970 (2014 November 9) to MJD 57200 (2015June 27). The energy range of each panel and the corresponding instrument can be found in the legend. Please refer to the text for details on the data taking and reduction for each instrument.

Open with DEXTER
In the text
thumbnail Fig. 4

Sequence of total (contours) and polarized (segments) intensity images of S4 0954+658 at 43 GHz, convolved with a beam of 0.24 × 0.15 mas2 at PA = –10°. The global total intensity peak is 1606 mJy beam−1 and the global polarized intensity peak is 104 mJy beam−1; black line segments within each image show the direction of polarization and their length is proportional to the polarized intensity. The black horizontal line indicates the position of the core, A0, and gray, blue, and red circles show the locations of knots K14a,b and K15, respectively. The size of the circles is proportional to the estimated average size in each epoch. The detailed characteristics of the knots can be found in Table 1 and B.1.

Open with DEXTER
In the text
thumbnail Fig. 5

Apparent distance from the radio core A0 of the new emerging knots, K14a,b and K15, as a function of time. The images fromwhich the apparent distances are calculated can be found in Fig. 4. We note that the K15 knot presents the highest apparent speed.

Open with DEXTER
In the text
thumbnail Fig. 6

Light curves for R-band polarimetry of S4 0954+65. Please refer to the text for details on the data taking and reduction for each instrument.

Open with DEXTER
In the text
thumbnail Fig. 7

Spectral energy distribution for the VHE MAGIC detection. Red symbols are strictly simultaneous to the VHE detection, blue symbols are for data taken during the same day and black symbols are for the closest observations. MAGIC spectral data (red circles) are for flare night only (2015 February 14, MJD 57067.14). Red filled circles are for the unfolded observed data points. The red shaded band shows the region of additional systematic uncertainty. Fermi-LAT data are the Pass 8 data for 2015 February 14 (1-day integration centered on the MAGIC observation, blue squares). Swift-XRT data are for 2015 February 13 (MJD 57066.70992, blue squares). Swift-UVOT data are given for the twoseparate exposure taken on 2015 February 13 (MJD 57066.76, blue triangles and dark blue triangles). R-band data are for 14th Feb (Tuorla, MJD 57067.16375 and AZT-8+ST7 MJD 57067.1, red diamonds). POLAMI data are for the 18th February (MJD 57071.5, black stars at 100 and 300 GHz). OVRO data are for 2015 February 10 and 19 (black circles at 15 GHz). Metsähovi for 2015 February 13 and 15 (MJD 57066.15, MJD 57068.68, black squares at 37 GHz). The gray data are for NED (light) and SSDC (dark) SED historical data points. The model from Tanaka et al. (2016; gray dashed curve) as well as the model presented here (black solide curve) include an emission component from synchrotron plus inverse Compton on a dusty torus (see text for details and Table 2 for the values of the physical parameters). The effect of the EBL attenuation is included in the modeling using the model by Finke et al. (2010) and a redshift of z = 0.368.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.