Connection between inner jet kinematics and broadband flux variability in the BL Lac object S5 0716+714
Key Words.:galaxies: active – BL Lacertae objects: individual: S5 0716+714 – radio continuum: galaxies – jets: galaxies – gamma-rays
We present a high-frequency very long baseline interferometry (VLBI) kinematical study of the BL Lac object S5 0716+714 over the time period of September 2008 to October 2010. The aim of the study is to investigate the relation of the jet kinematics to the observed broadband flux variability. We find significant non-radial motions in the jet outflow of the source. In the radial direction, the highest measured apparent speed is 37 c, which is exceptionally high, especially for a BL Lac object. Patterns in the jet flow reveal a roughly stationary feature 0.15 mas downstream of the core. The long-term fits to the component trajectories reveal acceleration in the sub-mas region of the jet. The measured brightness temperature, T, follows a continuous trend of decline with distance, T, which suggests a gradient in Doppler factor along the jet axis. Our analysis suggest that a moving disturbance (or a shock wave) from the base of the jet produces the high-energy (optical to -ray) variations upstream of the 7 mm core, and then later causes an outburst in the core. Repetitive optical/-ray flares and the curved trajectories of the associated components suggest that the shock front propagates along a bent trajectory or helical path. Sharper -ray flares could be related to the passage of moving disturbances through the stationary feature. Our analysis suggests that the -ray and radio emission regions have different Doppler factors.
The combination of very long baseline interferometry (VLBI) images and broadband flux density variability is a unique way to probe the emission mechanisms near the base of jets in the blazar class of active galactic nuclei (AGN). The broadband flares are often found to be connected to the ejection of new moving emission components into the jet (e.g., Marscher et al., 2008; Jorstad et al., 2013; Schinzel et al., 2012; Krichbaum et al., 2001, and references therein). Moreover, mm-VLBI observations offer a unique possibility to study the structural evolution in the parsec-scale region, which has been proposed to be the site of much of the high-energy emission (e.g., Rani et al., 2013c, 2014; Fuhrmann et al., 2014; Marscher et al., 2008; Schinzel et al., 2012). Therefore, these observations have provided new constraints on the physical parameters of the emission regions, i.e., size, brightness temperature, magnetic field, and motion.
The blazar S5 0716+714 (z 0.3, Nilsson et al., 2008; Danforth et al., 2012) is a BL Lac object with a featureless optical spectrum. It is one of the most intensively studied blazars because of its extreme variability properties across the entire electromagnetic spectrum (e.g., Villata et al., 2008; Fuhrmann et al., 2008; Rani et al., 2010b, a, 2013a, 2013b; Larionov et al., 2013). The broadband flux variability of the source is quite complex, with rapid flaring activity (on a timescale of a few hours to days) superimposed on top of a broad and slow variability trend on a timescale of year (Rani et al., 2013a; Raiteri et al., 2003). VLBI studies of the source show a core-dominated jet pointing towards the north (Bach et al., 2005; Britzen et al., 2009), while Very Large Array observations show a halo-like jet misaligned by 90 on kiloparsec scales. Britzen et al. (2009) suggested an apparent stationarity of jet components relative to the core; more recent studies, however, have reported motion as fast as 40 c (Rastorgueva et al., 2011; Larionov et al., 2013; Lister et al., 2013). Non-radial motion and wiggling component trajectories have often been observed in the inner mas jet region of the source (Britzen et al., 2009; Rastorgueva et al., 2011; Rani et al., 2014).
In Rani et al. (2013a) (hereafter Paper I), we presented the densely sampled multi-frequency observations of the source between April 2007 and January 2011. These observations allowed us to study the broadband flaring behavior of the source and to probe the physical processes, location, and size of the emission regions. The intense optical and -ray monitoring revealed fast repetitive variations (60–70 days) superimposed on a long-term variability trend with a time scale of 350 days, which propagated down to radio wavelengths with an observed time lag of 65 days. A detailed investigation of the optical flares found a variability amplitude proportional to the flux level, which can be explained by a variable Doppler factor. We also found that the shock-in-jet model for the evolution of radio flares requires geometrical variations in addition to intrinsic variations of the source. As a possible scenario to explain the observations, we suggested that the geometry significantly affects the long-term flux variations, which could be caused by a relativistic shock tracing a spiral path through the jet. To explain the multi-frequency behavior of an optical--ray outburst in 2011, Larionov et al. (2013) have also suggested a shock wave propagating along a helical path in the blazar’s jet.
In this paper, we use high-resolution multi-frequency VLBI observations to investigate the inner jet kinematics of S5 0716+714, with a focus on the major radio/optical/-ray flares over the time period of September 2008 to October 2010. We focus on the morphological evolution of the source to investigate its relation to the broadband flux variations reported in Paper I, in particular the high-energy emission. The paper is structured as follows. Section 2 provides a brief description of observations and data reduction. In Section 3, we report and discuss our results. Summary and conclusions are given in Section 4.
2 Multi-frequency VLBI data: observations and data reduction
To explore the inner jet kinematics of the source, we used the 7 mm (43 GHz) and 3 mm (86 GHz) VLBI data obtained between September 2008 and October 2010, with observations at 26 epochs during this period. The 7 mm data are from the Boston University monthly monitoring program of bright -ray blazars with the Very Long Baseline Array (VLBA)111VLBA-BU-BLAZARS, http://www.bu.edu/blazars. The data reduction was performed using standard tasks of the Astronomical Image Processing System (AIPS) and Difmap (Shepherd et al., 1994) software. The imaging of the source (including amplitude and phase self-calibration) was done using the CLEAN algorithm (Högbom, 1974) and SELFCAL procedures in Difmap (Shepherd et al., 1994). Several iterations of phase corrections followed by amplitude adjustments were adopted for the self-calibration process. More details of the data reduction can be found in Jorstad et al. (2005).
The 3 mm observations were performed using the Global mm VLBI Array (GMVA). The GMVA is currently an array consisting of 14 antennas in Europe and the United States, including the 8 VLBA stations equipped with 3 mm receivers, plus the Effelsberg, Onsala, Metsähovi, Pico Valeta, Plateau de Bure, and Yebes antennas. Between October 2008 and May 2010, observations were taken approximately every six months. Amplitude calibration is difficult at 3 mm, with atmospheric fluctuations often having a major effect. We find that amplitudes are frequently too low after the nominal calibration in AIPS. To correct for this, amplitudes are scaled by comparing the flux density on the shortest baselines with single dish measurements. Details of the 3 mm data reduction can be found in Hodgson et al. (2014).
We modeled the observed brightness distribution (the visibility amplitude and phases of the observed radio brightness of the source) by multiple circular Gaussian components, thus deriving positions, flux densities, and sizes of the distinct bright features in the jet. The model fits were carried out for all epochs within the DIFMAP package, starting with a point-like model and fixing the position of the brightest component to (0,0). Jet components were added until the addition of an extra component did not lead to a significant improvement in the value of the fit to the uv data. The uncertainties of the model component parameters were determined by comparing the parameter ranges obtained after performing model fits with a different number of model components. We used at least four different model fits for each epoch to obtain the parameter uncertainties. However, the uncertainties depend also on the self-calibration and data editing, and on the brightness and size of the components. In addition, the uncertainties increase with increasing distance from the core. This is accounted for by following an independent approach for the error estimation proposed by Krichbaum et al. (1998). The fitted model parameters for all of the epochs are listed in Table A, where we only list the formal errors obtained from the first method. Figures 1 and 2 show the fitted circular Gaussian components superposed on the clean images.
To investigate the kinematics in the jet of S5 0716+714, the individual model components were identified following the assumption that the changes of the flux density, distance from the VLBI core, position angle, and size should be small for the time period between adjacent epochs. In order to prevent a potentially large systematic error arising from the incorrect cross-identification of moving features from epoch to epoch, the simplest scheme was adopted while identifying the jet-features. A self-consistent cross-identification is proposed using all available model-fit parameters. These cross-identifications are not necessarily unique, especially if the source evolves in a more complex manner than assumed here.
3 Results and discussion
In this section we present the flux and spatial evolution of the bright radio emission in the jet of S5 0716+714. Special attention is given to the kinematics of the inner jet ( 2 mas from the core). We also investigate possible correlations between the jet kinematics and the broadband flux variations reported in Paper I.
3.1 Jet kinematics
Figures 1 and 2 display a sequence of VLBI maps convolved with a natural beam. The model fit parameters are given in Table A. In the following sub-sections, we will discuss in detail the apparent motions of components and brightness temperature gradient in the parsec-scale jet.
3.1.1 Component motion
For the kinematic study of individual components, we choose the VLBI core as a reference point,
fixed to coordinates (0,0). However, it is important to note that the absolute position of the
reference point i.e. core may change e.g. due to changes in opacity or instabilities in the jet; but,
the results presented here are not sensitive to that as we are interested in the relative motion of components
w.r.t. the core. The VLBI data were fitted by using different numbers of components
during different activity stages over the time span considered here. During
this period, we identified a total of 12 components - C1, K1, A1, A2, A3, and NC1 to NC7 - in addition to
the core, C0 (see Fig. 1).
Figure 3 plots the evolution of the distance of different knots from the core, and
Fig. 4 shows their trajectories in the XY-plane projected on the sky.
Stationary Features :
The sequence of VLBI images allow us to investigate the component motion along the jet as a function of time. Most of the components exhibit significant motion down the jet, except A1, A2, and A3. The latter are comparatively stable in position, although they exhibit significant scatter in position angle exceeding the uncertainties (see Fig. 4). Nevertheless, their radial distance from the core remains at 0.15 mas. Stationary features are a common characteristic in AGN jets (e.g., Fromm et al., 2013a; Jorstad et al., 2001; Britzen et al., 2010). In straight jets, these can be produced by recollimation shocks, instabilities [magnetohydrodynamic (MHD), and/or Kelvin-Helmholtz (KH)], or magnetic pinches (see Marscher, 2009; Hardee, 2006; Meier et al., 2001, for details). Bends in the jet can also cause quasi-stationary features, either because the jet turns more into the line of sight, thus increasing the Doppler beaming factor, or due to the formation of a shock that deflects the flow (Alberdi et al., 1993). Numerical simulations indicate that when a moving knot passes through a standing re-collimation shock, the components blend into a single feature, then split up after the collision with no lasting changes in the proper motion of the moving knot (Gomez et al., 1997; Fromm et al., 2012). During the interaction, the feature may move a short distance downstream before returning to its previous position. Therefore, it is quite possible that A1, A2, and A3 are the same standing feature, which becomes disturbed by moving features. We find that the kinematics of three new components NC2, NC3, and NC6, are significantly different before and after the interaction. The components move more slowly () as they approach the stationary feature, then faster () after the interaction. Such behavior is consistent with the stationary feature representing a bend in the jet (Alberdi et al., 1993). Indeed, a bend at 0.15 mas is evident in the 3 mm maps of the source (see Fig. 2). This suggests that the observed stationary features are standing oblique shocks that cause the flow of the jet to bend at 0.15 mas owing to a transverse pressure gradient or an impact with an interstellar cloud.
Moving Features :
The other components, C1, K1, and NC1 to NC7, exhibited significant motion in the radial direction. Component C1 represents very faint emission at a core separation of 1 mas. It seems to be the remnants of a component ejected in earlier epochs. The bright feature K1 could be a trailing component that forms behind a strong shock (Agudo et al., 2001) as it is identified at four epochs after the ejection of a new component (NC1).
As seen in Fig. 4, all components follow curved trajectories in the XY-plane projected on the sky. In addition, the trajectories differ from one component to the next. The wiggling trajectories in the jet might be a signature of helical motion. Helical jet models (e.g., Gomez et al., 1994; Hardee, 2006) represent the jet as an inhomogeneous flow of plasma interacting with ambient matter. Magnetohydrodynamic or Kelvin-Helmholtz instabilities (Hardee, 2006; Perucho et al., 2006) can explain the bending and helical structures, and possibly also less regular wiggles observed in parsec-scale radio jets.
In order to investigate component motions in the jet of S5 0716+714, we use the XY positions of the components relative to the core to fit the trajectories with polynomials of different order. The resulting fitted radial separation of each component, , is shown in Fig. 5. We find that a linear function is sufficient to fit the trajectories of components NC1, NC4, NC5, and NC7. However, a second-order polynomial is required to fit the trajectories of components NC2, NC3, and NC6, which corresponds to an apparent acceleration. Two multiple straight lines needed to fit different portions of the trajectories of these components in Fig. 5 provides a clear demonstration of the inadequacy of simple linear fits. The dashed curves represent the quadratic fits.
For each component, the fits yield an average proper motion and a mean speed. Using the estimated angular speed , we have computed the kinematic parameters of the jet, e.g., the apparent speed, , and Doppler factor, . We derive the apparent speed of the components, , from the angular speed, , using
where is the luminosity distance and the redshift of the source. The luminosity distance corresponding to z=0.3 is = 1600 Mpc for a CDM cosmology with = 0.27, = 0.73, and = 71 km s Mpc (Spergel et al., 2003). The angular and apparent speed of the individual components are given in Table 1. The calculated values range from to , the latter of which is unusually high for a blazar, especially for a BL Lac object.
For the non-ballistic components, NC2, NC3, and NC6, the quadratic fits represent significant acceleration along the jet axis.
Formally, we obtain = 4.901.55 mas yr for NC2,
= 2.450.80 mas yr for NC3, and = 0.940.36 mas yr for NC6.
Such apparent acceleration is a commonly
observed feature of many objects in the MOJAVE sample (Lister et al., 2013), for which a statistically significant
tendency for acceleration in the base of jets has been found.
A standing shock slows the flow while deflecting it, after which the flow can accelerate back to
the velocity it had upstream of the shock
(Gomez et al., 1997). However, since the apparent speed was relatively slow in the upstream region, this seems not to fit the case of S5 0716+714.
The data are more consistent with the scenario that, after a moving component passes through the stationary feature, it
bends towards or away from our line-of-sight. It then appears to move faster because its velocity vector subtends an angle closer to that which maximizes apparent motion, , where
is the velocity of the component divided by the speed of light.
An alternative possibility is that the underlying flow of the jet systematically accelerates outward.
Theoretical models involve strong magnetic fields associated
with the putative supermassive black hole/accretion disk system that play a key role in the initial acceleration
and collimation of the jet (Blandford & Payne, 1982; Blandford & Znajek, 1977; Meier et al., 2001). The conversion
of Poynting flux to flow energy is gradual and may persist out to parsec scales (e.g., Sikora et al., 2005).
Under this explanation, however, it is difficult to explain the straight-line kinematic evolution of some of the components.
Variability Doppler factor:
A variability Doppler factor using VLBI jet kinematics was defined by Jorstad et al. (2005) as
where is the angular size of the component, defined as 1.6 for a Gaussian with FWHM diameter measured at the epoch of maximum brightness, and is the luminosity distance. The timescale of variability is defined as = (Burbidge et al., 1974), where is the time separation in years between maximum () and minimum () flux densities. The estimated values of for the individual components, listed in Table 1, range from 6 to 21. Using the broadband flux density and spectral variability study of the source over the same time period (Paper I, Sect. 3.4), we obtain a range of self-consistent values for the Doppler factor. The different independent approaches suggest that 20, which is consistent with upper range of the values of obtained here.
3.1.2 Ejection Epochs
Back-extrapolation of the components’ motion allows us to estimate the time of zero separation from the core (i.e., the “ejection” time, T). As shown in the previous section, the components follow curved trajectories; therefore, we back-extrapolate the component trajectories in the XY-plane to determine the epoch of (0,0) separation. The error estimates are not very straightforward to calculate for components following non-ballistic trajectories. We use the estimated error on the fitted parameters for a given function to calculate its uncertainty. Back-extrapolation of the two envelopes of the function provide us with the error in T. The calculated ejection times for the individual components are listed in Table 1.
and denote, respectively, lower and upper values of the component speeds.
3.1.3 Brightness temperature gradient in the jet
The redshift-corrected brightness temperature (T) of the bright emission features can be approximated using the following relation (Jorstad et al., 2005):
where is the component flux density in Jy, is the FWHM of the circular Gaussian component in mas, and is the observing frequency in GHz; for the calculations here, we have only used 43 GHz data. The calculated brightness temperatures () of the components generally decline with radial separation from the core (see Fig. 8). This decaying trend can be approximated by a power law, with = 2.430.36, where is the radial distance along the jet axis. The solid line in Fig. 8 represents the fitted power law.
We place the observed brightness temperature gradient in the context of the shock-in-jet model (Marscher & Gear, 1985). In the common picture of this model, a relativistic shock propagates down a conical jet, slowly expanding adiabatically while maintaining shock conditions during propagation. While the shock propagates down the jet, it undergoes three major evolutionary stages dominated by inverse Compton, synchrotron, and adiabatic energy losses. As a result, the observed brightness temperature declines as a power law, T . The value of can be derived from the spectral evolution of the radio emission (Lobanov & Zensus, 1999; Fromm et al., 2013b). Following and S() (Lobanov & Zensus, 1999), we obtain
where , , and parametrize the variations along the jet axis of the Doppler factor, , magnetic field, , and power law distribution of energy of the emitting electrons, . Assuming a constant Doppler factor, i.e., = 0, and using a typical value of = 2 (corresponding to optically thin synchrotron spectral index = 0.5), and = 1 to 2 (for toroidal and poloidal magnetic fields, respectively), we obtain = 1.3 to 1.6, = 4.2 to 5.7, and = 3.2 to 4.7.
The observed brightness temperature gradient in S5 0716+714 (Fig. 8) has a slope equal to 2.430.36. The slope of the observed intensity gradient rules out the simple assumptions of a constant Doppler factor. We therefore consider a variable Doppler factor along the jet axis. We find that the observed value of is consistent with the adiabatic loss phase if , with 0.12 to 0.40, which corresponds to a moderate variation in the Doppler factor. However, with relatively larger variations in the Doppler factor 0.47 to 0.74, the observed value of would agree with synchrotron loss phase as well.
The evolution of the flux density and frequency ( and ) at the spectral turnover from synchrotron self-absorption of the radio flares, discussed in Paper I (Section 3.3.2), also suggests a variation of along the jet axis. In that study, we found that changes as with during the rise and during the decay of the first radio flare. The evolution of the second flare was governed by during the rising phase and during the decay (see Paper I for details). Therefore, it is evident that the two flares require a different dependence of on distance from the core, which is also inconsistent with what we obtained for the brightness temperature gradient ( to 0.74). This can be explained in two ways. First, the single-dish radio flux correlates well with the VLBI core flux (see Section 3.2 for details). This implies that the radio flares originate within the core. The core is more compact and farther upstream than the jet components, with a comparatively higher brightness temperature, and can follow a different dependence of on . The individual components trace the intensity gradient further downstream of the jet, and therefore can individually have a different – dependence. A second possibility is that the brightness temperature gradient plot corresponds to an averaged behavior of the two radio flares.
3.2 Jet kinematics and broadband emission
For the comparison of the observed broadband flaring behavior with the jet kinematics, we have used the flux density light curves from Paper I, where the details of observations and data reduction are given. Broadband flux light curves are shown in Fig. 7. Panel (a) presents the weekly averaged -ray photon flux light curve. Panel (b) shows the X-ray light curve, while panel (c) displays the optical V-band light curve. Panel (d) presents the radio (single-dish) and 43 GHz core222We used a mean value of the measured optically thin spectral index, = 0.4 (Rani et al., 2013a) to scale the 86 GHz flux density measurements to those at 43 GHz. flux density curves plotted in the same panel. Panel (e) exhibit flux density light curves of the individual components shown in different colors, with the same color used to shade their respective ejection times. The major outburst between JD’ 840 and 1150 is marked with a horizontal arrow in panel (b). The repetitive optical flares with a typical duration of 120–140 days are label as “0” to “9”, and the broadband flares are denoted “G” for -rays, “X” for X-rays, “O” for optical and “R” for radio followed by the number adjacent to them; e.g., we have flares G0 to G9 at -rays, X6 at X-rays, O0 to O9 at optical, and R6 and R9 at radio frequencies. The sharp -ray flares are labeled as “A” to “J”, and have a duration of 10–20 days.
Although VLBI observations of S5 0716+714 are missing at the peak of the two major radio flares (JD’ 1150–1200 and 1300), there is still an apparent one-to-one correlation between the core flux and single-dish radio flux measurements. This suggests that the single-dish radio flares are produced in the core, which agrees with the high mean core-to-jet flux ratio, which ranges from 1.3 to 21.4 with a mean of 4.94.0. The broadband flux cross-correlations (Paper I) also indicate that the flux variations at -ray (and optical) and single-dish radio frequencies are correlated such that the peak of the major optical/-ray outburst leads the peak of the radio outburst by a two-month time period.
3.2.1 Component ejections and their connection with broadband flares
Because of the superposition of different modes of flaring activity, the broadband flaring behavior is quite complex in the source. The major outburst (marked with a horizontal arrow in Fig. 7) is accompanied by faster repetitive optical/-ray flares (labeled as “0” to “9”) and sharp -rays flares (labeled as “A” to “J”). The jet kinematic study during this period suggests significant changes in the jet morphology of the source; in total, seven new components were ejected. It is therefore plausible that the ejection of new components could be connected to the broadband flares. We have identified a stationary feature at a distance of 0.15 mas from the core (Section 3.1.1). Interaction of a moving component with the stationary feature could amplify the magnetic field and accelerate particles to produce a multi-waveband flare (Marscher, 2014; Fromm et al., 2012).
|Component||Flares observed at the time of|
|zero separation||component crosses|
|from the core||stationary feature|
In Fig. 7, the shaded areas represent the ejection times of new jet components and the green dashed lines mark the average epochs when a moving component crosses the stationary feature at 0.15 mas. In Table LABEL:tab3, we list the new components (col. 1), along with the flares that coincide with their respective ejection times (col. 2). Flares observed closer in time to the passage of a moving component across the stationary feature are listed in col. 3. We note that most of the sharp -ray flares (B, C, D, G, I, and J) are observed when a moving component crosses the stationary feature, while most of the repetitive faster flares (0, 2, 3, 5, and 6) coincide with the ejections of new components from the core.
In Paper I, we found that the high-energy (optical and -ray) emission is significantly correlated with the radio flares such that the major optical/-ray outburst propagates down to radio frequencies (R6 flare) with a time delay of 65 days, following a power-law dependence on frequency with a slope 0.3. We also noticed that the evolution of radio flare R6 is in accordance with a generalized shock model (Valtaoja et al., 1992). This suggests that the origin of the major outburst could be related to a disturbance/shock that moves down the jet.
We interpret the broadband flaring events in the following manner (see Fig. 8). A disturbance at the base of the jet could cause a shock wave to form and propagate down the jet (Valtaoja et al., 1992; Marscher & Gear, 1985), which could be related to the onset of the major optical/-ray outburst (JD’ 840) and disturbs the jet outflow, resulting in the ejection of two new components (NC2 and NC3, Fig. 8a). The core is optically thick at that point, as the shock forms upstream of the core; for this reason, we do not see much activity in the core. Repetitive optical/-ray flares and the curved trajectories of the associated components favor the propagation of the shock along a bent trajectory, possibly a helical path. The moving shock reaches its maximum strength at JD’ 1150 (Fig. 8 b), close to the peak of the major optical/-ray outburst (see Section 3.3, Paper I for details), which also coincides with the onset of radio flare R6. At this point, the shock reaches the core, which starts to become optically thin. After this, further disturbances passing through the core correspond to the ejections of three new components (NC3, NC4, and NC5). Major radio outburst R6 is observed 65 days later, accompanied by the ejection of component NC7 (Fig. 8 c). As the moving features cross the the stationary feature at 0.15 mas, at first, the stationary feature gets destructed after the interaction, and later it appears again. The times of passage of moving components through the stationary feature appear to coincide with the sharp -ray flares in most cases. Therefore, the sharp -ray flares seem to be produced via the interaction of moving features with the stationary one.
X-ray flare X6 coincides with a similar increase in the flux of the 43 GHz core. Optical flare O7 peaks about a month later, which roughly coincides with the expected passage through the stationary feature of the same disturbance that caused flare X6. However, no bright moving knot is apparent in the next VLBA image days after O7. Also, the dramatic decrease in the millimeter-wavelength flux during this interval implies that the disturbance is strongly quenched following flare O7. There are two anomalies in this sequence of events: there is no optical flare as the disturbance crosses the core, and there is no -ray activity during the entire time span covering flares X6 and O7. The first indicates that, for reasons that are unclear (e.g., non-favorable magnetic field orientation relative to a shock; Summerlin & Baring, 2012), electrons with energies high enough to produce optical and -ray emission were not produced in the core region at this time. The second implies that there was a dearth of seed photons for inverse-Compton scattering 0.15 mas from the core during flare O7.
3.2.2 Doppler beaming of radio and gamma-ray emission regions
The observed time lag between -ray and 43 GHz core flares can be used to estimate the separation between the two emission regions by using the following equation (Fuhrmann et al., 2014):
where is the viewing angle of the source, is the apparent jet speed, and is the observed time lag. Using the observed value = 8232 days (Rani et al., 2014), we have = (0.0690.026)() pc.
The measured apparent speeds in the source cover a wide range. The current observations suggest slower motion (8 c) on sub-mas scales (0.3 mas), which later get accelerated to speeds as fast as 40 c on mas scales. The source also has a prominent bent jet structure at sub-mas scales (see Fig. 2); as a result, the sub-mas and mas regions could have different viewing angles. For = 6-8 c, the viewing angle333, , is 6-9; this gives = 2.9-4.4 pc. Using = 40 c, we have = 1.4 and 110 pc. Even for a larger viewing angle, =10, we obtain 16 pc. Pushkarev et al. (2012) estimated that the 15 GHz core is located at a distance of 6.68 pc from the jet base/apex. Using the technique to align the optically-thin emission based on model-fitted components and morphological similarities, we derive a shift of 0.0420.013 mas (0.190.06 parsec) between the 15 GHz and 43 GHz core, which puts the location of the 43 GHz core at a distance of 6.5 pc from the jet-apex. Therefore, the distance between -ray emission region and 43 GHz core cannot be larger than 6.5 pc. For a typical viewing angle , 6.5 pc gives 8 c. However, we note that for a larger viewing angle ( 10), could be as high as 10-12 but not unless the viewing angle is larger than 30, which is unlikely to be the case for S5 0716+714. The results therefore suggest that the high-energy and radio emission regions have different apparent speeds and probably different viewing angles as well. This implies that the -ray and radio emission regions have different Doppler factors.
4 Summary and conclusion
The monthly sampled mm-VLBI monitoring of S5 0716+714 allows us to study the inner jet kinematics and to investigate its relation to the observed broadband flux variability. Our study reveals significant non-radial motions in the jet outflow, including variations in the orientation of the sub-parsec-scale jet (see Rani et al. 2014 for details) and possibly helical trajectories of moving components.
In radial directions, the individual components exhibit extreme apparent speeds as high as 37 c. An apparent speed of 20 c seems to be more typical in S5 0716+714 on parsec scales. Recently, Lister et al. (2013) reported an apparent speed of 43.61.3 c. These values would even be extreme for a quasar. Therefore, moving components in the jets of BL Lac objects can be just as fast as in quasars, despite flux-limited surveys at lower resolution (e.g. Lister et al., 2009) finding that the mean apparent velocity of BL Lac objects is lower than for quasars. Further mm-VLBI studies are needed to determine the extent to which angular resolution and the variety of spectral classes of BL Lac objects (low- to high-frequency peaks of the synchrotron spectral energy distributions) skew the apparent velocity distribution. It is possible that the fast BL Lacs are actually quasars in disguise, with their broad emission lines swamped by the non-thermal optical continuum. On the other hand, emission lines have never been observed in S5 0716+714.
The bright components in 0716+714 follow a continuous decay trend in the observed brightness temperature (T), which implies a variation in Doppler factor along the jet axis. The spectral evolution of radio flares (see Paper I for details) also suggests variations in Doppler factor. A significant correlation between -ray flux and variations in orientation of the sub-parsec scale jet (Rani et al., 2014) supports the possibility that the observed flux variations in S5 0716+714 are strongly influenced by Doppler factor variations, and therefore may have a geometrical origin.
Our study suggests a connection between jet kinematics and the observed broadband flaring activity. We find that the origin of major optical/-ray/radio outbursts in the source seems to be related to a disturbance, such as shock wave, propagating down the jet. The repetitive, faster optical/-ray flares, superposed on the major outbursts, and non-ballistic motions of the associated jet components can be interpreted via helical jet models (as discussed by, e.g., Blandford & Payne, 1982; Blandford & Znajek, 1977; Camenzind & Krockenberger, 1992) and/or magnetohydrodynamic instabilities (Hardee, 2006). Our analysis favors that the major optical/-ray flares in S5 0716+714 are produced upstream of the 7 mm VLBI core. The sharp -ray flares, however, seem to be produced when the moving components cross a stationary feature 0.15 mas from the core. Finally, our analysis suggests that the upstream -ray emission region has a different Doppler factor than observed in the parsec-scale radio jet.
Acknowledgements.The Fermi/LAT Collaboration acknowledges the generous support of a number of agencies and institutes that have supported the Fermi/LAT Collaboration. 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 research has made use of data obtained with the Global Millimeter VLBI Array (GMVA), which consists of telescopes operated by the MPIfR, IRAM, Onsala, Metsahovi, Yebes and the VLBA. The data were correlated at the correlator of the MPIfR in Bonn, Germany. The VLBA is an instrument of the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This study makes use of 43 GHz VLBA data from the VLBA-BU Blazar Monitoring Program (VLBA-BU-BLAZAR; http://www.bu.edu/blazars/VLBAproject.html), funded by NASA through the Fermi Guest Investigator Program. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated by Associated Universities, Inc. 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. BR is thankful to Svetlana Jorstad, Benoit Lott, Biagina Boccardi, and Shoko Koyama for useful discussion and comments on the paper. We thank the referee for constructive comments.
- Agudo et al. (2001) Agudo, I., Gómez, J.-L., Martí, J.-M., et al. 2001, ApJ, 549, L183
- Alberdi et al. (1993) Alberdi, A., Marcaide, J. M., Marscher, A. P., et al. 1993, ApJ, 402, 160
- Bach et al. (2005) Bach, U., Krichbaum, T. P., Ros, E., et al. 2005, A&A, 433, 815
- Blandford & Payne (1982) Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883
- Blandford & Znajek (1977) Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433
- Britzen et al. (2009) Britzen, S., Kam, V. A., Witzel, A., et al. 2009, A&A, 508, 1205
- Britzen et al. (2010) Britzen, S., Kudryavtseva, N. A., Witzel, A., et al. 2010, A&A, 511, A57
- Burbidge et al. (1974) Burbidge, G. R., Jones, T. W., & Odell, S. L. 1974, ApJ, 193, 43
- Camenzind & Krockenberger (1992) Camenzind, M. & Krockenberger, M. 1992, A&A, 255, 59
- Danforth et al. (2012) Danforth, C. W., Nalewajko, K., France, K., & Keeney, B. A. 2012, ArXiv e-prints
- Fromm et al. (2012) Fromm, C. M., Perucho, M., Ros, E., et al. 2012, International Journal of Modern Physics Conference Series, 8, 323
- Fromm et al. (2013a) Fromm, C. M., Ros, E., Perucho, M., et al. 2013a, A&A, 551, A32
- Fromm et al. (2013b) Fromm, C. M., Ros, E., Perucho, M., et al. 2013b, A&A, 551, A32
- Fuhrmann et al. (2008) Fuhrmann, L., Krichbaum, T. P., Witzel, A., et al. 2008, A&A, 490, 1019
- Fuhrmann et al. (2014) Fuhrmann, L., Larsson, S., Chiang, J., et al. 2014, MNRAS, 441, 1899
- Gomez et al. (1994) Gomez, J. L., Alberdi, A., & Marcaide, J. M. 1994, A&A, 284, 51
- Gomez et al. (1997) Gomez, J. L., Marti, J. M. A., Marscher, A. P., Ibanez, J. M. A., & Alberdi, A. 1997, ApJ, 482, L33
- Hardee (2006) Hardee, P. E. 2006, in American Institute of Physics Conference Series, Vol. 856, Relativistic Jets: The Common Physics of AGN, Microquasars, and Gamma-Ray Bursts, ed. P. A. Hughes & J. N. Bregman, 57–77
- Högbom (1974) Högbom, J. A. 1974, A&AS, 15, 417
- Jorstad et al. (2005) Jorstad, S. G., Marscher, A. P., Lister, M. L., et al. 2005, AJ, 130, 1418
- Jorstad et al. (2001) Jorstad, S. G., Marscher, A. P., Mattox, J. R., et al. 2001, ApJ, 556, 738
- Jorstad et al. (2013) Jorstad, S. G., Marscher, A. P., Smith, P. S., et al. 2013, ApJ, 773, 147
- Krichbaum et al. (1998) Krichbaum, T. P., Alef, W., Witzel, A., et al. 1998, A&A, 329, 873
- Krichbaum et al. (2001) Krichbaum, T. P., Graham, D. A., Witzel, A., et al. 2001, in Astronomical Society of the Pacific Conference Series, Vol. 250, Particles and Fields in Radio Galaxies Conference, ed. R. A. Laing & K. M. Blundell, 184
- Larionov et al. (2013) Larionov, V. M., Jorstad, S. G., Marscher, A. P., et al. 2013, ApJ, 768, 40
- Lister et al. (2013) Lister, M. L., Aller, M. F., Aller, H. D., et al. 2013, ArXiv e-prints
- Lister et al. (2009) Lister, M. L., Cohen, M. H., Homan, D. C., et al. 2009, AJ, 138, 1874
- Lobanov & Zensus (1999) Lobanov, A. P. & Zensus, J. A. 1999, ApJ, 521, 509
- Marscher (2009) Marscher, A. P. 2009, ArXiv:0909.2576
- Marscher (2014) Marscher, A. P. 2014, ApJ, 780, 87
- Marscher & Gear (1985) Marscher, A. P. & Gear, W. K. 1985, ApJ, 298, 114
- Marscher et al. (2008) Marscher, A. P., Jorstad, S. G., D’Arcangelo, F. D., et al. 2008, Nature, 452, 966
- Meier et al. (2001) Meier, D. L., Koide, S., & Uchida, Y. 2001, Science, 291, 84
- Nilsson et al. (2008) Nilsson, K., Pursimo, T., Sillanpää, A., Takalo, L. O., & Lindfors, E. 2008, A&A, 487, L29
- Perucho et al. (2006) Perucho, M., Lobanov, A. P., Martí, J.-M., & Hardee, P. E. 2006, A&A, 456, 493
- Pushkarev et al. (2012) Pushkarev, A. B., Hovatta, T., Kovalev, Y. Y., et al. 2012, A&A, 545, A113
- Raiteri et al. (2003) Raiteri, C. M., Villata, M., Tosti, G., et al. 2003, A&A, 402, 151
- Rani et al. (2010a) Rani, B., Gupta, A. C., Joshi, U. C., Ganesh, S., & Wiita, P. J. 2010a, ApJ, 719, L153
- Rani et al. (2010b) Rani, B., Gupta, A. C., Strigachev, A., et al. 2010b, MNRAS, 404, 1992
- Rani et al. (2013a) Rani, B., Krichbaum, T. P., Fuhrmann, L., et al. 2013a, A&A, 552, A11
- Rani et al. (2013b) Rani, B., Krichbaum, T. P., Lott, B., Fuhrmann, L., & Zensus, J. A. 2013b, Advances in Space Research, 51, 2358
- Rani et al. (2014) Rani, B., Krichbaum, T. P., Marscher, A. P., et al. 2014, A&A, 571, L2
- Rani et al. (2013c) Rani, B., Lott, B., Krichbaum, T. P., Fuhrmann, L., & Zensus, J. A. 2013c, A&A, 557, A71
- Rastorgueva et al. (2011) Rastorgueva, E. A., Wiik, K. J., Bajkova, A. T., et al. 2011, A&A, 529, A2+
- Schinzel et al. (2012) Schinzel, F. K., Lobanov, A. P., Taylor, G. B., et al. 2012, A&A, 537, A70
- Shepherd et al. (1994) Shepherd, M. C., Pearson, T. J., & Taylor, G. B. 1994, in Bulletin of the American Astronomical Society, Vol. 26, Bulletin of the American Astronomical Society, 987–989
- Sikora et al. (2005) Sikora, M., Begelman, M. C., Madejski, G. M., & Lasota, J.-P. 2005, ApJ, 625, 72
- Spergel et al. (2003) Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175
- Summerlin & Baring (2012) Summerlin, E. J. & Baring, M. G. 2012, ApJ, 745, 63
- Valtaoja et al. (1992) Valtaoja, E., Terasranta, H., Urpo, S., et al. 1992, A&A, 254, 71
- Villata et al. (2008) Villata, M., Raiteri, C. M., Larionov, V. M., et al. 2008, A&A, 481, L79
Appendix A Model fit component results
|822||2.0440.102||00||00||0.0340.002||core 7 mm|
|856||2.1670.108||00||00||0.0340.002||core 7 mm|
JD = JD - 2454000,
S : The integrated flux in the component,
r : The radial distance of the component center from the center of the map,
: The position angle of the center of the component, and
: The FWHM of the component.
|Table 1 continued.|
: Identification of the individual components. If a component appeared only in a single epoch,
we labeled it with X.