Optical Outburst of the Blazar S4 0954+658 in Early 2015

We analyze the behavior of the BL Lac object S4 0954+658 during an unprecedented bright optical flare in early 2015. The optical flare was accompanied by a powerful γ-ray flare and the detection of very-high-energy γ-ray emission. We analyze total and polarized intensity images obtained with the VLBA at 43 GHz and discover a new bright polarized superluminal knot, which was ejected from the VLBI-core during the peak of the flare.


Introduction
The BL Lac object S4 0954+658 (z = 0.367) displays significant flux and polarization variability on both intra-day and longer timescales.Optical variability of the source has been comprehensively studied in [1,2].The multi-wavelength behavior of S4 0954+658 during a powerful outburst in March-April 2011 was analyzed in [3].Here we perform an analysis of an outburst that occurred in early 2015.We combine observations from cm-wave to γ-ray energies, and compare the light curves with the evolution of the parsec-scale jet structure of the source at 43 GHz.

Observation and Data Reduction
Optical (R-band) flux densities were obtained from photometric observations at the 0.4 m telescope of St. Petersburg State U. (LX200) and the 0.7 m telescope of the Crimean Astrophysical Observatory (AZT-8).The data analysis for these telescopes is described in [4].We also use R-band data obtained with the Perkins Telescope of Lowell Observatory.
We derived 0.1-200 GeV γ-ray flux densities by analyzing data from the Large Area Telescope (LAT) of the Fermi γ-ray Space Telescope with the standard software [5], and constructed a γ-ray light curve with seven-day binning, with a detection criterion that the maximum-likelihood test statistic (TS) should exceed 10.0.Swift/XRT 0.3-10 keV automatically processed data [6] were used to construct an X-ray light curve.S4 0954+658 is monitored approximately monthly by the Boston University (BU) group with the VLBA at 43 GHz within a sample of bright γ-ray blazars.We analyze total and polarized intensity images produced by the BU group and collaborators.Each image in Stokes I, Q, and U parameters is fit by a model consisting of a number of components with circular Gaussian brightness distributions.Identification of components in the jet across epochs is based on analysis of their flux, position angle, distance from the core, size, degree of polarization, and EVPA (Electric Vector Position Angle).We have computed kinematic parameters of knots (proper motion, velocity, and acceleration) by fitting the positions of each component over epochs by different polynomials, in the same manner as described in [7].In addition, we use data from the Owens Valley Radio Observatory (OVRO) [8] and the Submillimeter Array (SMA) [9] to construct radio-band light curves.

Results and Discussion
In early 2015, S4 0954+658 underwent a dramatic outburst at all wavelengths from gamma to radio band.Figure 1  Although the γ-ray flux falls below the detection limit during most of the period of Fermi observations, the source exhibits the highest γ-ray activity and shows substantial variability during the 2015 flare.According to the Swift data, there is a bright X-ray flare contemporaneous with the γ-ray flare.In addition to the prominent GeV γ-ray activity, VHE emission was discovered by MAGIC [10] on 13-14 February 2015 (MJD = 57067).Tanaka et al. [11] reported significant hardening of the spectrum in MeV/GeV during the flare maximum.
The flux of the 43 GHz VLBI-core, as well as the polarization degree, increased just before and after the optical flare maximum.Unfortunately, we have no VLBA data at 43 GHz during the flare maximum.Because of this, we cannot determine directly whether the core was in a high state, but the optical flare was contemporaneous with flares at 1 mm and 2 cm, which most likely originated in the core.One can see in Figure 1 that the VLBI-core light curve is in agreement with the flux variations at 1 mm and 2 cm during the entire time span of the observations.The pattern of the flux variations in radio and optical bands is quite similar, although the radio light curve is smoother.
We carefully study the optical and radio polarization behavior of S4 0954+658 near the ejection times of the components.At the quiescent state, the optical EVPA is quite stable (∼−10 • ) and tends to align with the parsec-scale jet direction at 43 GHz (∼−20 • ). Figure 1 shows that the optical EVPA is also in good agreement with the EVPA of the core.The fractional polarization of the core usually increases during optical activity (see [3]).
During 2015, several components can be seen in the jet at 43 GHz: K12,K13, A, and T (see [3] for an analysis of the jet kinematics during 2006-2012).Knot K12 was ejected at ∼ MJD 56570 (β app = 12.9 ± 0.23 c, P.A. = 13 • ± 8 • ).Knot A is a stationary component, and knot T is most likely a trailing component generated by K12 [12].Separations of knots from the VLBI-core as a function of time at 43 GHz and the trajectories of the knots are presented in Figure 2. The polarization parameters of K12 (EVPA and polarization degree) are similar to the optical polarization parameters at several epochs after its ejection from the 43 GHz VLBI-core (Figure 1).The knot maintained a high degree of polarization during later epochs (see Figure 3).
The ejection of knot K13 (MJD 57080±15) is contemporaneous with the optical flare in early 2015 (see Figure 1).The apparent velocity of the knot (β app = 23.8 ± 0.5 c) is higher than that of previous components, and the P.A. (−7 • ± 2 • ) is quite different from the mean direction of previously ejected components (∼−15 • to −25 • , see [3]).The EVPA of K13 is also close to the optical value, but, due to gaps in time coverage of the VLBA data, it is difficult to resolve the 180 • ambiguity at MJD ∼ 57200 (see Figure 1).
The observed simultaneous multi-wavelength flares and the ejection of a new superluminal component confirm the connection between the jet behavior at 43 GHz and activity at higher energies that was inferred from previous observations of S4 0954+658.The presence of such a connection implies that the multi-wavelength flare in early 2015 occurred near the 43 GHz VLBA-core.

Figure 1 .
Figure 1.(a) Light curves of S4 0954+658 from 2005 to 2016 (from top to bottom): γ-ray, X-ray, R-band optical, and radio (1 mm, 7 mm, and 2 cm).The ejections of new superluminal knots are marked with green bands.(b) Upper panel shows R-band optical light curve, fractional polarization, and EVPA.Gray dots show measurements in the R band, black dots are those of the VLBI-core at 43 GHz.Lower panel shows an enlargement of the light curve during the ejections of knots K12 and K13. Green and blue dots correspond to knots K12 and K13, respectively.

Figure 2 .
Figure 2. (a) Separations of knots from the VLBI-core as a function of time at 43 GHz; (b) 43 GHz image of S4 0954+658 with trajectories of knots superposed.

Figure 3 .
Figure 3.Total (yellow contours) and polarized (color scale) intensity images at 43 GHz; yellow line segments over the color scale show the direction of the electric vector.