First 3 mm-VLBI imaging of the two-sided jet in Cygnus A
Key Words.:galaxies: jets – galaxies: active – instrumentation: high angular resolution
Aims:We present for the first time Very-Long-Baseline Interferometry images of the radio galaxy Cygnus A at the frequency of . Thanks to the high spatial resolution of only 200 Schwarzschild radii (), such observations provide an extremely detailed view of the nuclear regions in this archetypal object and allow us to derive important constraints for theoretical models describing the launching of relativistic jets.
Methods:A pixel-based analysis of the jet outflow, which still appears two-sided on the scales probed, was performed. By fitting Gaussian functions to the transverse intensity profiles, we could determine the jet width in the nuclear region.
Results:The base of the jets appears wide. The minimum measured transverse width of is significantly larger than the radius of the Innermost Stable Circular Orbit, suggesting that the outer accretion disk is contributing to the jet launching. The existence of a faster and Doppler de-boosted inner section, powered either from the rotation of the inner regions of the accretion disk or by the spinning black hole, is suggested by the kinematic properties and by the observed limb brightening of the flow.
The launching mechanism of relativistic extragalactic jets is a much debated and crucial topic for understanding these fascinating objects. There is currently a broad consensus on the fundamental role that the strong magnetic fields anchored in the accretion disk play in the process. However, to fully reproduce the observed properties of such magnetically-driven jets, e.g. Lorentz factors as high as 50 and opening angles as small as half a degree, is still a challenge for theorists. General relativistic MHD simulations (e.g. Tchekhovskoy et al., 2011) have shown that highly collimated, relativistic jets can be efficiently powered by the black hole’s rotational energy extracted through the magnetic field, as described in the work of Blandford & Znajek (1977). This scenario is supported by a recent observational study of a large sample of blazars (Ghisellini et al., 2014), indicating that the contribution of a spinning black hole is necessary to account for the estimated jet powers. Alternatively, the rotation of the accretion disk has been proposed as a viable power source for magnetic jet launching (Blandford & Payne, 1982). Predictions for the terminal Lorentz factor of such disk winds are diverse, with values ranging from mildly relativistic (, McKinney & Gammie, 2004; De Villiers et al., 2005) to relativistic (, Komissarov et al., 2007).
A fundamental tool for discriminating between the two launching mechanisms is the direct imaging of the jet base through VLBI observations and the determination of its transverse size. This kind of study has been rarely performed because the jet base is usually not resolved. Moreover, owing to synchrotron opacity effects (Königl, 1981), the location of the VLBI core is often shifted at large distances, Schwarzschild radii (), from the black hole (e.g. Marscher et al., 2008). One exception is M 87, for which a distance of only 14-23 has been inferred for the 43 GHz core (Hada et al., 2011). In the same source, VLBI pilot studies at 1.3 (with the Event Horizon Telescope) have determined a core size of (Doeleman et al., 2012), which represents an upper limit for the true size of the jet apex. This important measurement indicates that the jet of M 87 is likely anchored close to the Innermost Stable Circular Orbit (ISCO). Performing such an analysis in other sources is clearly fundamental for testing the variety of possible scenarios. However, while future instruments with improved capabilities (e.g. phased ALMA) will certainly enrich the sample of suitable targets, this is currently quite limited.
In this letter we present, for the first time, a detailed imaging of the jet base in Cygnus A through Global VLBI observations at 86 (GMVA). Among the powerful FR II radio galaxies, Cygnus A is one of the very few two-sided sources whose jet launching region can be probed with sufficient sensitivity at high radio frequencies. At the source redshift (z=0.056)111 , assuming a cosmology with , , and assuming a black hole mass of (Tadhunter et al., 2003), GMVA observations allow scales of the order of to be probed. The large jet viewing angle enables a prominent counter-jet to be detected in this source, which represents a fundamental advantage when investigating the nuclear regions of a jet. In fact, since the transverse width can be determined continuously along the two-sided flow, the assumption-dependent back-extrapolations of the properties of the jet apex can be replaced by actual measurements.
The following analysis complements a previous VLBI study at 43 of the jet transverse profile and of the kinematic properties (Boccardi et al., 2016, hereafter B1). This revealed the presence of a pronounced transverse stratification, both in flux density (limb brightening) and bulk Lorentz factor (), typical of a flow with spine-sheath structure.
|Eb, PdB, PV|
|Eb, PdB, PV|
|Eb, PdB, PV, Mh|
2 Observations and VLBI images
The GMVA data set comprises three epochs from observations in 2009/2010 with a cadence between five and seven months (Table 1). The duration of the observations was 15-16 hours, 8 of which were spent on Cygnus A. This, together with the large number of antennas employed and with the favourable declination of the source, enabled good uv-coverage. Data were recorded in single polarization mode with a recording rate of 512 , which corresponds to an observing bandwidth of 128 . The correlation took place at the MPIfR in Bonn, Germany. The data reduction was carried out in AIPS (Greisen, 1990). Aside from the standard calibration procedures, special care was devoted during the fringe fitting. First, as described in Martí-Vidal et al. (2012), the manual phase calibration was applied. Single sub-band phases and delays were calculated using one or more scans showing strong detections from the calibrator J2015+3710. Then, the global fringe fitting was performed (Schwab & Cotton, 1983). The amplitude calibration was based on the measured system temperatures and gain curves at each telescope, and it included a correction for atmospheric opacity.
The imaging and self-calibration of amplitude and phase were performed in DIFMAP (Shepherd et al., 1994). The final clean maps obtained are presented in Fig. 1. The nuclear region is detected with a high signal-to-noise ratio S/N (between 1200-2000), with the brightest feature – which we refer to as the “core” – showing a prominent variability over the monitoring period. The total flux density varies between 0.65 and 1 and is consistent with the measurement of provided by 86 single-dish observations at IRAM (Agudo et al., 2010).
The amplitude variation during the self-calibration procedure was not greater than 30%, thus this percentage can be assumed as the uncertainty associated with the flux density measurements.
A visual inspection of the maps already shows the high morphological complexity of Cygnus A at this frequency and at this angular resolution. A prominent counter-jet is detected in all three epochs. The jet emission is relatively smooth – i.e. there are no bright features besides the core – but at the same time rich with sub-structures. It can also be noticed that, due to the low S/N in the faintest regions, the emission is not detected over the full jet cross-section in a single epoch. The jet structure can be recovered much better after stacking the three maps, as shown in Fig. 2. In this image, characterized by a reduced noise level of , the jet limb brightening already observed at 43 GHz also appears to be better defined, especially in the outer and faintest regions of the approaching jet.
3 Analysis and discussion
3.1 Transverse width profile
The analysis of the jet transverse structure in Cygnus A was first presented at the frequency of 43 in B1. In this study it was shown that the flow emanates parabolically from the position of an emission gap located at east of the brightest feature. A measurement of the jet transverse width was also provided for the nuclear region. However, the minimum size of occurring in the vicinity of the faint emission gap was about seven times smaller than the beam size in the transverse direction. In this letter we exploit the improved resolution achieved by GMVA observations for obtaining a more solid estimate of the transverse width at the base of the flow. Given the better quality and the reduced noise level of the stacked image as compared to the single epochs, we focus our analysis on the former. The stacked image was restored with a circular beam of 0.1 , which is approximately the natural resolution in the transverse direction (Table 1). Then, it was sliced pixel by pixel (1px=0.02 ) along the jet axis using the AIPS task SLICE. The transverse intensity distribution in each slice was finally fitted with a single Gaussian profile using the task SLFIT.
Figure 3 shows the dependence of the de-convolved jet width on the distance from the core in the central of the source. The distance , also expressed in , was de-projected assuming a viewing angle (as calculated in B1). The error bars are conservatively assumed equal to one fifth of the measured FWHM. The minimum jet width occurs for , where the measured transverse size is . This corresponds to about half of the nominal resolution in the transverse direction. However, as discussed by several authors (e.g. Fomalont, 1989; Lobanov, 2005; Kovalev et al., 2005), the limiting resolution of a certain brightness distribution mainly depends on the S/N of the detection. Following Lobanov (2005), the resolution limit for the case of a Gaussian profile fitted to a brightness distribution with uniform weighting is
where is the peak flux density of the transverse intensity profile, is the rms noise of the map and the equivalent circular beam (for uniform weighting) expressed as a function of the major and minor beam axes, and , as . The value of was calculated for each slice, and its profile as a function of is also reported in Fig. 3. Except for a faint area in the counter-jet ( ), the measured width is everywhere well above the resolution limit, therefore the jet is resolved in the nuclear region.
3.2 Minimum width and implications for the jet launching
The minimum jet width as measured at 86 translates into a linear size of . Unlike the case of M 87 (Doeleman et al., 2012), the jet apex appears to be much wider than the ISCO radius, which varies between 1 and 9 depending on the BH spin and on the kind of orbit (prograde or retrograde) of the BH/accretion disk system. These results suggest that at least part of the jet base in Cygnus A is anchored at large disk radii from the central BH, i.e. the outer regions of the accretion disk are involved in the jet launching process. Based on the sub-parsec scale kinematic study presented in B1, we infer that such a disk wind is characterized by mildly relativistic speeds, with the bulk Lorentz factor varying from at the jet edge to in sections of the flow closer to the jet axis.
We point out, however, that our width measurement does not exclude the existence of both a narrower and an even broader component. Concerning the former, it should be recalled that in jets seen at large viewing angles Doppler de-boosting affects not only the receding side but also the approaching one. For sufficiently high speeds, the Doppler factor can be much smaller than one, so that the flow is completely invisible. In the case of Cygnus A, the gradual fading of the accelerating inner-jet components was directly observed in our kinematic analysis. This, together with the observed limb brightening, strongly indicate that an invisible faster flow exists in the central body of the jet. This spine may be driven either from the rotation of the inner regions of the accretion disk or by the spinning black hole itself. A hypothetical broader and diffuse disk wind, instead, is likely resolved out in these high resolution images. Indeed VLBI observations at 5 (Carilli et al., 1991) support the existence of such a component. At this frequency, the jet appears as an almost cylindrical resolved flow with transverse width between 2 and 2.2 , i.e. of the order of 9000 . Although this large size may partially be the effect of scattering in the interstellar medium, there is an indication that the jet base extends to much larger disk radii than measured at 86 .
3.3 Comparison with results at 43 GHz
The discussion in the previous section has not quite addressed the problem of truly locating the central engine in Cygnus A. Given the roughly symmetric expansion of jet and counter-jet shown in Fig. 3, the most natural interpretation would be that the supermassive black hole is located in the vicinity of . However, although the jet base is assumed to be the narrowest part of the jet in the simplest models, this may not necessarily be the case. An interesting argument comes from the comparison of the 86 width profile with results at 43 . For the purposes of the comparison, the analysis was repeated at the two frequencies following the method described in Sect. 3.1 after convolving the maps with the same circular beam of 0.15 . Results are presented in Fig. 4.
We first notice that the dip in the 43 profile at , i.e. in the vicinity of the gap of emission, is not observed at 86 . Correspondingly, the 86 map (Fig. 3) does not show a pronounced emission gap.222In principle, this comparison may be affected by a misalignment of the maps due to an opacity shift of the cores. However, this is expected to be negligible at high frequencies, as confirmed by a cross-correlation of the jet parameters from which we infer a shift . The observed frequency dependence of the properties of this region supports the idea that the gap is the result of absorption, which becomes very reduced at 86 . If synchrotron opacity is responsible for the it, the gap most likely marks the true location of the central engine. A detailed spectral analysis using more sensitive data is in progress with the aim of clarifying the nature of this feature. Absorption may also explain another striking aspect of the width profiles in Fig. 4. While results at the two frequencies are quite consistent in the jet side, the 43 counter-jet appears narrower than at 86 . The puzzling asymmetry between the opening angles of the approaching and receding side observed at 43 GHz in B1, almost disappears at 86 GHz. The existence of a parsec scale free-free absorber covering the counter-jet was suggested by several studies (Krichbaum et al., 1998; Bach et al., 2004; Struve & Conway, 2010) and may well explain the observed features. The results presented in this letter hint at the presence of sub-pc scale material that is still highly opaque at 43 GHz but not at 86 .
We have presented the first VLBI images of the radio galaxy Cygnus A at the frequency of 86 . The enhanced spatial resolution of these maps allowed us to investigate the transverse width profile at the onset of the two-sided flow. The jet base appears wide with a minimum width of . This value is much higher than the ISCO radius, thus we infer that the emission is produced by a mildly relativistic, parabolically expanding disk wind. The existence of a central and faster jet, possibly driven by the black hole and invisible due to Doppler de-boosting, is however not excluded. It is, on the contrary, suggested by the kinematic properties and by the observed limb brightening of the flow. In the hypothesis that the Blandford & Znajek (1977) and Blandford & Payne (1982) mechanisms act simultaneously in relativistic jets, giving rise to a spine-sheath structure of the flow (e.g. Hardee et al., 2007; Xie et al., 2012), we indeed expect the slower disk wind to always dominate the emission in radio galaxies and the BH-driven jet to be the main observable component in blazars due to differential boosting. The comparison of the width profiles at 43 and 86 indicates that free-free opacity and/or synchrotron opacity are still strong at 43 but not at 86 . Observing at high radio frequencies is therefore fundamental for unveiling the physical properties of jets at their onset.
Acknowledgements.We thank the referee for helpful comments. We also thank Andrei Lobanov and Vassilis Karamanavis for the useful discussions, and Michael Lindqvist for his support during the observations at Onsala. 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.
- Agudo et al. (2010) Agudo, I., Thum, C., Wiesemeyer, H., & Krichbaum, T. P. 2010, ApJS, 189, 1
- Bach et al. (2004) Bach, U., Krichbaum, T. P., Middelberg, E., et al. 2004, in European VLBI Network on New Developments in VLBI Science and Technology, ed. R. Bachiller, F. Colomer, J.-F. Desmurs, & P. de Vicente, 155–156
- 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
- Boccardi et al. (2016) Boccardi, B., Krichbaum, T. P., Bach, U., et al. 2016, A&A, 585, A33
- Carilli et al. (1991) Carilli, C. L., Bartel, N., & Linfield, R. P. 1991, AJ, 102, 1691
- De Villiers et al. (2005) De Villiers, J.-P., Hawley, J. F., Krolik, J. H., & Hirose, S. 2005, ApJ, 620, 878
- Doeleman et al. (2012) Doeleman, S. S., Fish, V. L., Schenck, D. E., et al. 2012, Science, 338, 355
- Fomalont (1989) Fomalont, E. B. 1989, in Astronomical Society of the Pacific Conference Series, Vol. 6, Synthesis Imaging in Radio Astronomy, ed. R. A. Perley, F. R. Schwab, & A. H. Bridle, 213
- Ghisellini et al. (2014) Ghisellini, G., Tavecchio, F., Maraschi, L., Celotti, A., & Sbarrato, T. 2014, Nature, 515, 376
- Greisen (1990) Greisen, E. W. 1990, in Acquisition, Processing and Archiving of Astronomical Images, ed. G. Longo & G. Sedmak, 125–142
- Hada et al. (2011) Hada, K., Doi, A., Kino, M., et al. 2011, Nature, 477, 185
- Hardee et al. (2007) Hardee, P., Mizuno, Y., & Nishikawa, K.-I. 2007, Ap&SS, 311, 281
- Komissarov et al. (2007) Komissarov, S. S., Barkov, M. V., Vlahakis, N., & Königl, A. 2007, MNRAS, 380, 51
- Königl (1981) Königl, A. 1981, ApJ, 243, 700
- Kovalev et al. (2005) Kovalev, Y. Y., Kellermann, K. I., Lister, M. L., et al. 2005, AJ, 130, 2473
- Krichbaum et al. (1998) Krichbaum, T. P., Alef, W., Witzel, A., et al. 1998, A&A, 329, 873
- Lobanov (2005) Lobanov, A. P. 2005, ArXiv Astrophysics e-prints:0503225
- Marscher et al. (2008) Marscher, A. P., Jorstad, S. G., D’Arcangelo, F. D., et al. 2008, Nature, 452, 966
- Martí-Vidal et al. (2012) Martí-Vidal, I., Krichbaum, T. P., Marscher, A., et al. 2012, A&A, 542, A107
- McKinney & Gammie (2004) McKinney, J. C. & Gammie, C. F. 2004, ApJ, 611, 977
- Schwab & Cotton (1983) Schwab, F. R. & Cotton, W. D. 1983, AJ, 88, 688
- 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
- Struve & Conway (2010) Struve, C. & Conway, J. E. 2010, A&A, 513, A10
- Tadhunter et al. (2003) Tadhunter, C., Marconi, A., Axon, D., et al. 2003, MNRAS, 342, 861
- Tchekhovskoy et al. (2011) Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79
- Xie et al. (2012) Xie, W., Lei, W.-H., Zou, Y.-C., et al. 2012, Research in Astronomy and Astrophysics, 12, 817