VLBA monitoring of Mrk 421 at 15 GHz and 24 GHz during 2011
Key Words.:Galaxies: active – BL Lacertae objects: Mrk 421 – Galaxies: jets
Context:High-resolution radio observations are ideal for constraining the value of physical parameters in the inner regions of active-galactic-nucleus jets and complement results on multiwavelength (MWL) observations. This study is part of a wider multifrequency campaign targeting the nearby TeV blazar Markarian 421 (z=0.031), with observations in the sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and rays (Fermi-LAT, MAGIC, VERITAS).
Aims:We investigate the jet’s morphology and any proper motions, and the time evolution of physical parameters such as flux densities and spectral index. The aim of our wider multifrequency campaign is to try to shed light on questions such as the nature of the radiating particles, the connection between the radio and -ray emission, the location of the emitting regions and the origin of the flux variability.
Methods:We consider data obtained with the Very Long Baseline Array (VLBA) over twelve epochs (one observation per month from January to December 2011) at 15 GHz and 24 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch.
Results:The structure of Mrk 421 is dominated by a compact (0.13 mas) and bright component, with a one-sided jet detected out to 10 mas. We identify 5-6 components in the jet that are consistent with being stationary during the 12-month period studied here. Measurements of the spectral index agree with those of other works: they are fairly flat in the core region and steepen along the jet length. Significant flux-density variations are detected for the core component.
Conclusions:From our results, we draw an overall scenario in which we estimate a viewing angle and a different jet velocity for the radio and the high-energy emission regions, such that the respective Doppler factors are and .
|Observation||MJD||Map peak||Beam||1 rms||Notes|
|date||(mJy/beam)||(mas mas, )||(mJy/beam)|
|15 GHz||24 GHz||15 GHz||24 GHz||15 GHz||24 GHz|
|2011/01/14||55575||348||319||1.05 0.65, 15.3||0.79 0.47, 8.84||0.19||0.18||No MK, no NL|
|2011/02/25||55617||391||338||1.16 0.74, 14.4||0.64 0.39, 6.52||0.35||0.18||NL snowing|
|2011/03/29||55649||386||359||1.06 0.66, 5.37||0.65 0.39, 4.26||0.17||0.24||No HK|
|2011/04/25||55675||367||308||0.92 0.50, 3.79||0.61 0.34, 3.88||0.28||0.33||-|
|2011/05/31||55712||355||297||0.93 0.51, 5.56||0.64 0.35, 8.41||0.25||0.29||-|
|2011/06/29||55741||262||208||0.89 0.50, 7.17||0.56 0.32, 12.3||0.17||0.37||No LA|
|2011/07/28||55770||220||197||0.91 0.56, 0.89||0.60 0.37, 1.14||0.20||0.30||-|
|2011/08/29||55802||275||200||0.97 0.55, 0.26||0.63 0.35, 2.70||0.18||0.27||No HK|
|2011/09/28||55832||264||238||1.06 0.67, 16.2||0.72 0.47, 18.4||0.26||0.24||No MK|
|2011/10/29||55863||261||167||1.06 0.69, 0.09||0.73 0.44, 3.73||0.26||0.14||HK snowing|
|2011/11/28||55893||283||201||1.02 0.59, 18.0||0.70 0.42, 14.3||0.28||0.17||No PT, FD, MK|
|2011/12/23||55918||295||287||0.89 0.49, 1.82||0.61 0.35, 5.74||0.15||0.21||No HK|
Markarian 421 (R.A.= , Dec.= 12′ 31.79906″, J2000) is one of the nearest () and brightest BL Lac objects in the sky. It was the first extragalactic source detected at TeV energies by the Cherenkov telescope at Whipple Observatory (Punch et al. 1992). The spectral energy distribution (SED) of this object, dominated by non-thermal emission, has two smooth broad components: one at lower energies, from radio band to the soft X-ray domain, and another at higher energies peaking at -ray energies (Abdo et al. 2011). The low-frequency peak is certainly due to synchrotron emission from relativistic electrons in the jet interacting with the magnetic field, whereas the high-frequency peak is probably due to the inverse Compton scattering of the same population of relativistic electrons with synchrotron low-energy photons (synchrotron-self-compton model, SSC, see Abdo et al. 2011; Tavecchio et al. 2001). In this framework, multiwavelength (MWL) coordinated campaigns are a fundamental tool for understanding the physical properties of the source, e.g. by studying variability, which is present at all frequencies, but particularly TeV energies where Gaidos et al. (1996) measured a doubling time of minutes. The accurate MWL study and SED modeling performed by Abdo et al. (2011) revealed some interesting results, such as the size of the emitting region and the magnetic field , which in the context of the leptonic scenario they constrained to be and 0.05 G. However, the details of the physical processes responsible for the observed emission are still poorly constrained. Because of the considerable variability and the broadband spectrum, multiwavelength long-term observations are required for a good comprehension of the emission mechanisms.
This study is part of a new multi-epoch and multi-instrument campaign, which also involves observations in the sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and rays (Fermi-LAT, MAGIC, VERITAS), as well as at the cm wavelengths with low resolution observations (e.g. F-GAMMA, Medicina). The aim of this observational effort is to shed light on fundamental questions such as the nature of the radiating particles, the connection between the radio and -ray emission, the location of the emitting regions, and the origin of the flux variability. Very long baseline interferometry (VLBI) plays an important role in addressing these scientific questions because it is the only technique that can resolve (at least partially) the inner structure of the jet. Therefore, cross-correlation studies of Very Long Baseline Array (VLBA) data with data from other energy ranges (in particular rays) can provide us with important information about the structure of the jet and the location of the blazar emission.
At radio frequencies, Mrk 421 clearly shows a one-sided jet structure aligned at a small angle with respect to the line of sight (Giroletti et al. 2006). In this work, we present new VLBA observations to study in detail the inner jet structure on parsec scales. We are able to investigate the evolution of shocks that arise in the jet, by means of the model-fitting technique. In earlier works (Piner et al. 1999; Piner & Edwards 2004), the jet components show only subluminal apparent motion, which seems to be a common characteristic of TeV blazars. Thanks to accurate measurements of changes on parsec scales, by the VLBA, we can find valid constraints on the geometry and kinematics of the jet.
This paper is structured as follows: in Section 2 we introduce the dataset, in Section 3 we report the results of this work (model fits, flux density variations, apparent speeds, jet sidedness, spectral index), and in Section 4 we discuss results giving our own interpretation in the astrophysical context. We used the following conventions for cosmological parameters: km sec Mpc, and , in a flat Universe. We defined the spectral index such that .
We observed Mrk 421 throughout 2011 with the VLBA. The source was observed once per month, for a total of 12 epochs, at three frequencies: 15, 24, and 43 GHz. In this paper, we present the complete analysis of the whole 15 and 24 GHz datasets. We also observed, at regular intervals, three other sources (J0854+2006, J1310+3220, and J0927+3902) used as fringe finders and calibrators for the band pass, the instrumental (feed) polarization, and the electric-vector position angle. At each epoch, Mrk 421 was observed for nearly 40 minutes at each frequency, spread into several scans of about 3 minutes each, interspersed with calibrator sources in order to improve the (u,v)-coverage. Calibrators were observed for about 10 minutes each, generally spread into three scans of 3 minutes.
For calibration and fringe-fitting, we used the AIPS software package (Greisen 2003), and for image production the standard self-calibration procedures included in the DIFMAP software package (Shepherd 1997), which uses the CLEAN algorithm proposed by Högbom (1974). In some epochs, one or more antennas did not work properly because of technical problems; for a complete report (see Table 1).
We show in Fig. 1 two images of Mrk 421 at 15 and 24 GHz, which were produced by stacking all the images respectively at 15 and 24 GHz created with DIFMAP. The alignment of the images was checked by comparing the pixel position of the peak. In both images, we set the lowest contour equal to about three times the off-source residual rms noise level.
All the 12 images at each frequency show a similar structure, consisting of a well-defined and well-collimated one-sided jet structure emerging from a compact nuclear region (core-dominated source). This is the typical structure of a BL Lac object (Giroletti et al. 2004a). The jet extends for roughly 4.5 mas (2.67 pc)1111 mas corresponds to 0.59 pc., with a position angle (PA) of (measured from north through east). This morphology agrees with the results of other studies of similar angular resolution (Marscher 1999).
Since the morphology was very stable from epoch to epoch, the stacking did not smear any details of the structure. A sample couple of 15 and 24 GHz images is shown in Fig. 2.
3.2 Model fits
For each epoch, we used the model-fitting routine in DIFMAP to fit the visibility data of the source in the plane with either elliptical or circular Gaussian components. In this way, we were able to investigate in detail the inner jet structure and its evolution. For all epochs at 15 GHz, a good fit was obtained with five Gaussian components, while at 24 GHz we needed six components. At both frequencies, we identified the core with the brightest, innermost, and most compact feature. We label the other components C1, C2, C3, and C4, starting from the outermost (C1) to the innermost (C4). The higher angular resolution achieved at 24 GHz resolves the second innermost 15 GHz component (C4 located at 0.45 mas from the core) into two features (C4b at 0.3 mas and C4a 0.7 mas from the core) (see Fig. 2).
Thanks to the extremely fine time-sampling, we were able to make an attempt to identify the same component in each epoch. Overall, the components extend out to a region of about 5 mas. In this way, with a limited number of components, it was possible to analyze the proper motions and flux density levels at various times. From Fig. 3, we can clearly see that the data occupy well-defined regions in the radius vs. time plot, and that this behavior helps us to identify the individual components across epochs.
All details of the model fit analysis are shown in Table 4. We calculated the uncertainties in the position (error bars in Fig. 3) using the ratio of the size of each component to the signal-to-noise ratio (S/N). In the case of very bright, compact components, the nominal error value is too small and replaced by a conservative value equal to 10% of the beam size (Orienti et al. 2011). On the other hand, when the calculated error is very large (i.e. comparable to the component radius, as in the case of a very extended component with a very low flux density), we replace it with a threshold value equal to the maximum value of the calculated errors for the same component at different epochs.
By comparing with Piner et al. (2010) and Piner & Edwards (2005), we can argue that the components C2, C3, C4a, and C4b of the present analysis quite likely correspond respectively to the components C5, C6, C7, and C8 of the aforementioned works. We also suggest that our component C1 corresponds to their component C4a. We will consider a more accurate identification of our components with those of previous analyses in a forthcoming paper, where the 43 GHz data analysis will also be presented.
3.3 Flux density variability
Using the results of the model-fit technique, we analyzed the temporal evolution of the flux density for each component of the source. The brightest component represents the core; at 15 GHz, it has a mean value of around 350 mJy, which decreases along the jet until values of about 10 mJy. By comparing the flux density of each component at the various epochs, it emerges that there are no significant variations in the flux densities of the C1-C4 components. The flux density of each component remains roughly constant at various times within the uncertainties calculated; in any case, there is no indication of flaring activity. Any small variations may be artifacts brought about by our fitting procedures: for instance, the flux density of the inner components may be underestimated in some cases, because part of it was incorporated into the core component flux, or the flux density of the most extended features may be underestimated at epochs missing some short baseline data because of telescope failures. This is a remarkable result because it also allows us to identify components based on their flux density values and confirms our choice of components based on their positions.
In Fig. 4, we show the light curve for Mrk 421 during 2011 at 15 GHz (upper panel) and at 24 GHz (lower panel), considering the total flux density (squares) and the core flux density (triangles). The light curve reveals an interesting feature: in the second part of the year (starting at MJD 55700), we clearly note a decrease in the total flux density. From the complete flux-density analysis, we found that the core is the component responsible for the decrease, while the extended region does not display any significant variations. To further exclude calibration effects, we performed the same analysis on the three calibrators. In Fig. 4, we present the light curves of the calibrator J1310+3220 (diamonds) at 15 GHz (upper panel) and 24 GHz (lower panel). From this comparison, we clearly see that the trend of the light curves for the two sources is very different. We can assert that the flux density decrease observed for the core of Mrk 421 is a real feature. Error bars were calculated by considering a calibration error of about of the flux density and a statistical error equal to three times the map rms noise.
3.4 Spectral index analysis
Thanks to our multi-frequency data set, we were able to carry out a detailed analysis of the spectral index distribution on the parsec scale. We conducted a quantitative assessment of the spectral index of the components by comparing the flux density at the two frequencies at each epoch and averaging the results to reduce the statistical fluctuations. For the core, we obtained a value of , and for the outermost components (C1, C2, and C3, considered altogether), we obtained . The uncertainties were calculated from the theory of the propagation of errors using the formula
where and represent the flux density and the respective uncertainty (see the last paragraph of section 3.3). These values agree with the results obtained from the spectral index image, which we produced with the following procedure: we first produced new images for each epoch and frequency using the same range and the same restoring beam ( mas mas). We then produced the spectral index image by combining the two images, resulting respectively from the stacking of all images at 15 GHz and at 24 GHz, using the AIPS task COMB, clipping the pixels with a S/N 3 using the input images. By summing images we were able to increase the S/N, and results are more reliable; the stability of the individual features guarantees that we do not lose structural information in the averaging process. The alignment of the images at the two frequencies was checked by comparing the pixel position of the peak. When a shift was present, we used the LGEOM task in AIPS to align the images. The resulting image presents the typical flat spectrum in the core region, with a steepening along the jet radius. However, despite the higher S/N obtained with image stacking, a significant spread in the spectral index values is present in the diffuse jet region, hence we do not show the image here.
3.5 Apparent speeds
From our model-fit, we infer a small or no displacement for the jet components. To verify and quantify this statement, we determined the speeds of each component by means of linear fits to the separation of the individual features from the core at different epochs. For the three outer components (C1, C2, and C3), we used combined data at 15 and 24 GHz, since the positions of each component at the two frequencies are generally consistent within the error bars. For the two inner components (C4a and C4b), we only used data at 24 GHz. Results are shown in Table 3.
We found low values for the apparent speeds, in agreement with previous studies (e.g. Piner & Edwards 2005). The two innermost components (C4a and C4b) are essentially stationary, with an upper limit to their separation velocity of c. In addition C2 and C3 are consistent with being stationary, while the outermost component C1 has a low-significance (1.5) subluminal motionc. If this trend of increasing velocity at larger radii were real and if the apparent speeds shown in Table 3 represented the bulk apparent speed of the plasma in the jet, we could speculate that some mechanism involving an acceleration acts in the outer region of the jet.
3.6 Jet/counter-jet ratio
We estimated the ranges of viewing angles and of from the jet/counter-jet brightness ratio. Assuming that the source has two symmetrical jets of the same intrinsic power, we used the equation
where and are, respectively, the jet and counter-jet brightnesses and represents the spectral index defined in the introduction; we adopted the exponent, since the jet is smooth and does not contain well-defined compact blobs. For the jet brightness, we used mJy/beam, measured at 24 GHz, in the image resulting from the stacking of all the 12 epoch images, in the jet region located at mas from the core. For the counter-jet, which is not visible, we used an upper limit provided by the rms noise level measured in the image, which resulted in mJy/beam; this consequently yields a lower limit to both and . With a value of , in agreement with our spectral index images, we obtained and then . Therefore, the minimum allowed jet bulk velocity is (corresponding to a bulk Lorentz factor ) and the maximum viewing angle is .
4 Discussion and conclusions
The Doppler factor defined as
is a key element in the study of blazars, since it affects various parameters such as the observed brightness, the SED peak frequency, the variability timescale, and more. Modeling of the SED and study of the variability in different wavebands generally require large values of the blazar Doppler factors; in the case of Mrk421, Gaidos et al. (1996) estimated , from the observed TeV variability time of about 30 min and Abdo et al. (2011) required a Doppler factor between 20 and 50 to reproduce the broadband SED. In turn, VLBI observations can also constrain by posing limits on and , as provided by the various arguments discussed in this work.
When closely spaced repeated observations are available, the study of the proper motion is a useful tool in determining the ranges for and . Surprisingly, several works (e.g. Giroletti et al. 2004b; Piner et al. 2008, 2010) have reported subluminal motions, sometimes consistent with the component being stationary, in the jets of all the 10 TeV blazars for which proper motion studies have been performed in the literature. Thanks to the large number of dual-frequency observations, the fine time-sampling, and the high quality of the data provided by the good -coverage, we performed a robust identification of the Gaussian components and constrained their motion to be consistent with no displacement at all. At the same time, the high sensitivity and in particular the stacked image place significant constraints on the jet/counter-jet ratio.
The first immediate consequence is that we can reject the hypothesis that the small is solely due to a projection effect, since it would require an unrealistically narrow viewing angle: in the case of component C4b, the upper limit to the observed motion implies a viewing angle to reproduce the observed jet/counter-jet ratio (and even smaller to agree with the high energy limits). If the jets’ distribution is isotropic on the sky, the real number of misaligned sources (parent population) is incompatible with these very small values of . For example, in the Bologna Complete Sample selected by Giovannini et al. (2005) at low frequency (and thus free from Doppler favoritism bias), one would expect fewer than 0.03 sources with (see also Tavecchio 2005; Henri & Saugé 2006).
On the other hand, since larger values of the viewing angle capable of reproducing the observed lack of proper motion are incompatible with the jet/counter-jet ratio, we conclude that the pattern velocity cannot be representative of the bulk flow velocity. The Gaussian components obtained in our model fit provide a good description of the visibility data but do not represent well-defined, high-contrast jet features (see Lyutikov & Lister 2010). In our interpretation, the low apparent speeds found imply that the proper motion of Mrk421 does not provide any information about the jet bulk velocity; even on the basis of the sole jet brightness ratio, untenable viewing angles would be necessary to match the pattern and bulk velocities.
What then are the real values of the viewing angle and the jet bulk velocity in Mrk421? To reproduce the observed jet asymmetry, we needed to consider a range of velocities and angles . We could exclude the upper range of the values, since this would not reproduce the high Doppler factors required by high energy observations (Gaidos et al. 1996; Abdo et al. 2011); in particular, we were unable to achieve when , respectively.
Smaller angles thus seem to be favored since they were the only ones consistent with high values of ; however, such small angles would still represent a challenge to the observed radio properties. We were able to estimate the intrinsic power of the radio core , by debeaming the observed monochromatic luminosity of the core with the equation
With a value of W Hz at 15 GHz, , and , we obtained for the intrinsic power of the core W Hz; this value was at the very low end of the typical range of intrinsic power found for different samples of radio galaxies (e.g. Liuzzo et al. 2011), suggesting that lower values of provide a more typical core power.
Moreover, the 12 monthly observations have not revealed any dramatic flux-density variability in the core of the source, which further points to a lower Doppler factor for the radio jet. We estimated the variability brightness temperature of the core () with the formula proposed by Hovatta et al. (2009)
where is the observed frequency in GHz, is the redshift, is the luminosity distance in meters, is the difference between the maximum value for the core flux density and the minimum value, and is the variability time. With the values provided by our observations and 90 days, we obtaind a value of K, which does not require any significant beaming.
where is the flux density of the component measured in Jy, and are the full widths at half maximum of the major and minor axes respectively of the component measured in mas, is the redshift, and is the observation frequency in GHz. The resulting ’s are on the order of a few K, only slightly exceeding the limit derived by Readhead (1994) from equipartition arguments.
Taken together, the lack of superluminal features, the low core dominance, and the weak variability suggest a scenario in which no strong beaming is required in the radio jet. This is not uncommon in TeV blazars (Piner & Edwards 2004; Piner et al. 2008), but unprecedentedly firm observational support for it has been provided by our intensive campaign. Low values of the Doppler factor, e.g. , can reproduce the observational radio properties, including the jet brightness asymmetry.
We conclude that the Doppler factor must be different in the radio band than the -ray band. Since we do not expect that the viewing angle changes significantly, this leads us to the necessity of a velocity structure in the jet, as previously discussed by e.g. Chiaberge et al. (2000), Georganopoulos & Kazanas (2003), and Ghisellini et al. (2005). Our images do not provide strong evidence in favour of either a radial or transverse velocity structure, although previous works have revealed a limb brightening in Mrk421, on both a milliarcsecond scale at 43 GHz (Piner et al. 2010) and at mas at 5 GHz (Giroletti et al. 2006). This would favor the presence of a transverse velocity structure across the jet axis. This structure consists of two components: a fast inner spine and a slower outer layer. Different Doppler factors were obtained depending on whether we measured the speed of the spine or the layer.
A viable scenario for Mrk421 is that the viewing angle is between and , which is consistent with the statistical counts of low-power radio sources and the possibility of reaching the high Doppler factors required by SED modeling and high-energy variability. The jet velocity is structured, with a typical Lorentz factor of in the radio region (yielding ), and in the high-energy emission region. For example, assuming , , and , we obtained a value of and and successfully reproduced all the observational properties of the source.
In summary, the detailed analysis presented in this paper has largely confirmed with improved quality expectations based on the knowledge so far achieved for TeV blazars. We have also estimated with a good level of significance some important and fundamental parameters (, , , ) that characterize the physical processes in blazars. However, there is still much to be understood and we expect to obtain other significant results from the analysis extended to other wavelengths, particularly in the -ray domain. Additional works on the dataset presented in this paper are planned, and will deal with, e.g., the 43 GHz images and the polarization properties. Moreover, we intend to combine our dataset with those of other works (e.g. Piner & Edwards 2005 or the MOJAVE survey, Lister et al. 2009) to increase the temporal coverage of the observations and obtain even tighter constraints over a longer time frame.
Acknowledgements.This work is based on observations obtained through the BG207 VLBA project, which makes use of the Swinburne University of Technology software correlator, developed as part of the Australian Major National Research Facilities Programme and operated under licence (Deller et al. 2011). The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. For this paper we made use of the NASA/IPAC Extragalactic Database NED which is operated by the JPL, Californian Institute of Technology, under contract with the National Aeronautics and Space Administration. We acknowledge financial contribution from grant PRIN-INAF-2011. This research is partially supported by KAKENHI (24540240). KVS and YYK are partly supported by the Russian Foundation for Basic Research (project 11-02-00368), and the basic research program “Active processes in galactic and extragalactic objects” of the Physical Sciences Division of the Russian Academy of Sciences. YYK is also supported by the Dynasty Foundation. The research at Boston University was supported in part by NASA through Fermi grants NNX08AV65G, NNX08AV61G, NNX09AT99G, NNX09AU10G, and NNX11AQ03G, and by US National Science Foundation grant AST-0907893. We thank Dr. Claire Halliday for the language editing work which improved the text of the present manuscript.
- Abdo et al. (2011) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011, ApJ, 736, 131
- Chiaberge et al. (2000) Chiaberge, M., Celotti, A., Capetti, A., & Ghisellini, G. 2000, A&A, 358, 104
- Cornwell (2009) Cornwell, T. J. 2009, A&A, 500, 65
- Deller et al. (2011) Deller, A. T., Brisken, W. F., Phillips, C. J., et al. 2011, PASP, 123, 275
- Gaidos et al. (1996) Gaidos, J. A., Akerlof, C. W., Biller, S., et al. 1996, Natur, 383, 319
- Georganopoulos & Kazanas (2003) Georganopoulos, M., & Kazanas, D. 2003, ApJ, 594, L27
- Ghisellini et al. (2005) Ghisellini, G., Tavecchio, F., & Chiaberge, M. 2005, A&A, 432, 401
- Giovannini et al. (2005) Giovannini, G., Taylor, G. B., Feretti, L., et al. 2005, ApJ, 618, 635
- Giroletti et al. (2004a) Giroletti, M., Giovannini, G., Taylor, G. B., & Falomo, R. 2004a, ApJ, 613, 752
- Giroletti et al. (2004b) Giroletti, M., Giovannini, G., Feretti, L., et al. 2004b, ApJ, 600, 127
- Giroletti et al. (2006) Giroletti, M., Giovannini, G., Taylor, G. B., & Falomo, R. 2006, ApJ, 646, 801
- Greisen (2003) Greisen, E. W. 2003, ASSL, 285, 109
- Henri & Saugé (2006) Henri, G., & Saugé, L. 2006, ApJ, 640, 185
- Hovatta et al. (2009) Hovatta, T., Valtaoja, E., Tornikoski, M., Lähteenmäki, A. 2009, A&A, 494, 527
- Högbom (1974) Högbom, J. A. 1974, A&AS, 15, 417
- Lister et al. (2009) Lister, M. L., Aller, H. D., Aller, M. F., et al. 2009, AJ, 137, 3718
- Liuzzo et al. (2011) Liuzzo, E., Giovannini, G., & Giroletti, M. 2011, arXiv:1110.6360
- Lyutikov & Lister (2010) Lyutikov, M., & Lister, M. 2010, ApJ, 722, 197
- Marscher (1999) Marscher, A. P. 1999, APh, 11, 19
- Orienti et al. (2011) Orienti, M., Venturi, T., Dallacasa, D., et al. 2011, MNRAS, 417, 359
- Piner et al. (1999) Piner, B. G., Unwin, S. C., Wehrle, A. E., Edwards, P. G., Fey, A. L., & Kingham, K. A. 1999, ApJ, 525, 176
- Piner & Edwards (2004) Piner, B. G., & Edwards, P. G. 2004, ApJ, 600, 115
- Piner & Edwards (2005) Piner, B. G., & Edwards, P. G. 2005, ApJ, 622, 168
- Piner et al. (2008) Piner, B. G., Pant, N., & Edwards, P. G. 2008, ApJ, 678, 64
- Piner et al. (2010) Piner, B. G., Pant, N., & Edwards, P. G. 2010, ApJ, 723, 1150
- Punch et al. (1992) Punch, M., Akerlof, C. W., Cawley, M. F., et al. 1992, Natur, 358, 477
- Readhead (1994) Readhead, A. C. S. 1994, ApJ, 426, 51
- Shepherd (1997) Shepherd, M. C. 1997, ASPC, 125, 77
- Tavecchio et al. (2001) Tavecchio, F., Maraschi, L., Pian, E., et al. 2001, ApJ, 554, 725
- Tavecchio (2005) Tavecchio, F. 2005, Proceedings of the MG10 Meeting held at CBPF, Rio de Janeiro, Brazil, 20-26 July 2003, p. 512
- Tingay et al. (1998) Tingay, S. J., Murphy, D. W., Lovell, J. E. J., et al. 1998, ApJ, 497, 594
|EpochFor more details on various epochs see Table 1.||Frequency||Component ID||SFlux density in mJy.||Estimated errors for the component flux density.||rr and PA are the polar coordinates of the component’s center with respect to the core. The position angle (PA) is measured from North through East.||Estimated errors in the component position.||PAr and PA are the polar coordinates of the component’s center with respect to the core. The position angle (PA) is measured from North through East.||aa and b are the FWHM of the major and minor axes of the Gaussian component.||Position angle of the major axis measured from North through East.|