Next Article in Journal
Cosmological Effects of Quantum Vacuum Condensates
Next Article in Special Issue
Optical/Infrared Polarised Emission in X-ray Binaries
Previous Article in Journal
Exploring String Axions with Gravitational Waves
Previous Article in Special Issue
Linear Polarization Properties of Parsec-Scale AGN Jets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetic Field Studies in BL Lacertae through Faraday Rotation and a Novel Astrometric Technique

1
Instituto de Astrofísica de Andalucía-CSIC, Glorieta de la Astronomía s/n, E-18008 Granada, Spain
2
International Centre for Radio Astronomy Research, M468, The University of Western Australia, 35 Stirling Hwy, Crawley, Western Australia 6009, Australia
3
CSIRO Astronomy and Space Science, 26 Dick Perry Avenue, Kensington Western Australia 6151, Australia
4
Observatorio Astronómico Nacional (IGN), Alfonso XII, 3 y 5, E-28014 Madrid, Spain
*
Author to whom correspondence should be addressed.
Galaxies 2017, 5(4), 97; https://doi.org/10.3390/galaxies5040097
Submission received: 25 September 2017 / Revised: 5 December 2017 / Accepted: 5 December 2017 / Published: 12 December 2017
(This article belongs to the Special Issue Polarised Emission from Astrophysical Jets)

Abstract

:
It is thought that dynamically important helical magnetic fields twisted by the differential rotation of the black hole’s accretion disk or ergosphere play an important role in the launching, acceleration, and collimation of active galactic nuclei (AGN) jets. We present multi-frequency astrometric and polarimetric Very Long Baseline Array (VLBA) images at 15, 22, and 43 GHz, as well as Faraday rotation analyses of the jet in BL Lacertae as part of a sample of AGN jets aimed to probe the magnetic field structure at the innermost scales to test jet formation models. The novel astrometric technique applied allows us to obtain the absolute position at mm wavelengths without any external calibrator.

1. Introduction

One of the main questions to answer in active galactic nuclei (AGN) jets is the nature of the core, defined as the upstream end of the jet. There are two complementary models. In one of these, the core is the place where the jets becomes optically thin, so the position of the core shifts with observing frequency. This was observed at centimeter wavelength [1]. On the other hand, there is evidence that most gamma-ray flares are simultaneous with a major outburst in the core of the jet observed at millimeter wavelength. These flares can be explained as a superluminal component going through the mm-core, identified as a recollimation shock [2,3,4,5]. In this case, the position of the core is fixed at a jet location.
It is broadly accepted that magnetic fields play a relevant role in the dynamics of the innermost regions of relativistic jets. Most of the models propose a helical geometry, although there are still uncertainties regarding the actual configuration.
With the aim of obtaining a better understanding of the structure of the magnetic field and analyzing the nature of the core jet, we have developed a polarized multi-frequency Very Long Baseline Array (VLBA) program in which we have observed a sample of AGN (BL Lacertae, 3C120, 3C273, CTA102, 0716+714, 3C111). In this work, we present the results for BL Lacertae.

2. Observation

Very Long Baseline Interferometry at radio and millimeter frequencies has provided ultra-high-resolution observations which allow us to study the innermost part of AGN jets. The Very Long Baseline Array (VLBA) have receivers in a wide range of frequencies (from 1.2 GHz to 96 GHz), allowing us to perform multifrequency studies in polarization.
We present polarized intensity VLBA images of BL Lacertae at 15, 22, and 43 GHz, obtained in July 2013. These data were taken as a part of an astrometric multifrequency program which includes BL Lacs, flat spectrum radio quasars (FSRQs), and radio galaxies observations with the VLBA. Calibration of the data was performed within the AIPS software following the standard procedure for polarimetric observations. The data were edited, self-calibrated, and imaged both in total and polarized intensity with a combination of the AIPS and Difmap software packages [6]. Calibration of the electric vector position angle (EVPA) at 15 and 22 GHz was performed by comparison of the integrated polarization measured from the VLBA images and that obtained from MOJAVE team observations [7] and VLBA-BU-BLAZAR program [8], respectively, at contemporaneous epochs. With these values, we have calculated the rotation measure (RM) on the southern part of the source (which is an optically thin region and has very low values of RM) and obtained the absolute orientation of the EVPA at 22 GHz.

3. Astrometric Technique

This project was designed to obtain astrometric information from centimeter to millimeter wavelength. Obtaining astrometric measurements at millimeter wavelength is not a simple issue. The technique applied in this work is explained in more detail in [9]. One method that has worked well for millimeter wavelength is the source/frequency phase referencing (SFPR) [10], but this technique requires a second calibrator within about 10 from the target. To overcome this situation, we have developed a new approach to the SFPR, called multi-frequency phase referencing (MFPR), in which we do not need a second calibrator at millimetric wavelength [9].
We have followed a specific observation strategy consisting of four blocks, in order to calculate the dispersive (i.e., instrumental and ionospheric propagation effects) and non-dispersive residual delay errors (e.g., tropospheric propagation effects and geometric errors) without losing astrometric information [9]: (1) A usual phase-reference block, observing the target interleaved with a calibrator at 5, 8, 15, and 22 GHz. This block is used to perform the usual phase reference astrometry at centimeter wavelength. (2) An ionospheric block, in which we observe the target (in this case BL Lacertae) at 1.3, 5, and 22 GHz. These observations are used to calculate the ionospheric propagation effects. (3) A frequency-phase-transfer block, observing the target at low frequency (22 GHz) interleave with high frequencies (43 and 86 GHz). This block is necessary to scale the phase corrections at high frequency based on the low frequency ones. (4) A repeat of the ionospheric block. These four steps are cycled around for the duration of the experiment.
In order to eliminate the non-dispersive residual errors, the data at higher frequencies (43 and 86 GHz) are calibrated scaling the phase corrections from near-simultaneous observations at a lower frequency. This is possible as long as the frequency pair (e.g., 22 and 43 GHz) have an integer frequency ratio. After that, it is necessary to calibrate the dispersive residual errors. We have used the corrections of the prime calibrator to remove the instrumental effects in BL Lacertae. The novel part of this technique consists of the way to calculate the ionospheric propagation effects.
This novel technique calculates the total electron content (TEC) in the target direction using a linear fit in ν 2 to the semi-simultaneous delay measurements. We have calculated this fit for each antenna, intermediate frequency channel (IF), and interval time. This allows us to eliminate the ionospheric contribution from the data.

4. Results

Figure 1, Figure 2 and Figure 3 show polarimetric VLBA images of BL Lacertae at 15, 22, and 43 GHz, respectively. The structure of BL Lacertae is dominated by the nuclear region plus an extended and bright structure towards the south. The electric vector position angle (EVPA) along the southern jet, and also along the nuclear region, are mainly aligned with the jet direction. This behavior is consistent with the presence of a helical magnetic field geometry, where the toroidal component is dominant over the poloidal one all along the pc-scale jet [11]. This behavior is consistent with previous results presented in [12].
The rotation measure (RM) depends on the electron density and the component of the magnetic field in the direction of the line-of-sight. For this reason, RM maps give relevant information about the magnetic structure in the radio source or in the immediate vicinity. Figure 4 shows the rotation measure map of BL Lacertae obtained by combining the 15, 22, and 43 GHz data, superposed with contours representing the total intensity at 15 GHz. BL Lacertae are also affected by a small Faraday rotation contribution from our own galaxy, which has been accounted for using the values in [13]. We have estimated the errors in the RM as the uncertainty in the chi-square fit to the λ 2 dependence of the EVPAs. The EVPA errors for each observing frequency are calculated from error propagation of the uncertainties in Q and U, added in quadrature to the error in the calibration of the absolute orientation.
The values of RM increase towards the nuclear regions, as we can see in slice 1 along the jet core. This result is in agreement with previous studies of BL Lacertae (e.g., [12]), and could be produced by a higher electron density in the internal regions.
By inspecting slices 2 and 3 in Figure 4, we note a gradient of RM across the jet width, as it was observed in a number of sources [12,14,15,16].
According to [15], a minimum resolution of two observing beams is required for a reliable determination of RM gradients in jet profiles. Our achieved resolution is just above two observing beams for the slices in the core area (slices 1 and 2), and larger for slice 3 in the jet, confirming the reliability of the measured RM gradients.
The gradient across the jet width could be a signature of a helical magnetic field in which the toroidal component dominates over the poloidal one. The results presented in this work are in agreement with the structure of the RM presented by [17], in which the RM shows a gradient with position angle from the centroid of the core with positive RM values in the area upstream of the centroid and negative downstream. This change in the sign of the RM supports the model of the helical magnetic field and has also been observed in a few other sources [15,17].
The astrometry results are presented in more detail in [9]. The absolute position of the core between cm-mm wavelengths departs from the expected in case the core is located at the optically thin–thick transition. That is, the core shift between 22 and 43 GHz is smaller than the prediction from measurements of the core shift at cm wavelengths. This supports the model on which the mm-VLBI core corresponds to a recollimation shock.

Acknowledgments

This research has been supported by the Spanish Ministry of Economy and Competitiveness grant AYA2016-80889-P. 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.

Author Contributions

J.L.G., R.D. and M.J.R. designed the experiment and the observations. S.N.M. analyzed the data, made the plots, and wrote the paper. R.D. and M.J.R. contributed to the analysis of the data and J.L.G. contributed in the writing and interpretation of the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hada, K.; Doi, A.; Kino, M.; Nagai, H.; Hagiwara, Y.; Kawaguchi, N. An origin of the radio jet in M87 at the location of the central black hole. Nature 2011, 477, 185–187. [Google Scholar] [CrossRef] [PubMed]
  2. Marscher, A.P.; Jorstad, S.G.; D’Arcangelo, F.D.; Smith, P.S.; Williams, G.G.; Larionov, V.M.; Oh, H.; Olmstead, A.R.; Aller, M.F.; Aller, H.D.; et al. The inner jet of an active galactic nucleus as revealed by a radio-to-gamma-ray outburst. Nature 2008, 452, 966–969. [Google Scholar] [CrossRef] [PubMed]
  3. Marscher, A.P.; Jorstad, S.G.; Larionov, V.M.; Aller, M.F.; Aller, H.D.; Lähteenmäki, A.; Agudo, I.; Smith, P.S.; Gurwell, M.; Hagen-Thorn, V.A. Probing the inner jet of the quasar PKS 1510–089 with multi-waveband monitoring during strong gamma-ray activity. Astrophys. J. Lett. 2010, 710, L126. [Google Scholar] [CrossRef]
  4. Jorstad, S.G.; Marscher, A.P.; Smith, P.S.; Larionov, V.M.; Agudo, I.; Gurwell, M.; Wehrle, A.E.; Lähteenmäki, A.; Nikolashvili, M.G.; Schmidt, G.D.; et al. A tight connection between gamma-ray outbursts and parsec-scale jet activity in the quasar 3c 454.3. Astrophys. J. 2013, 773, 147. [Google Scholar] [CrossRef]
  5. Casadio, C.; Gómez, J.L.; Jorstad, S.G.; Marscher, A.P.; Larionov, V.M.; Smith, P.S.; Gurwell, M.A.; Lähteenmäki, A.; Agudo, I.; Molina, S.N.; et al. A multi-wavelength polarimetric study of the blazar CTA 102 during a gamma-ray flare in 2012. Astrophys. J. 2015, 813, 51. [Google Scholar] [CrossRef]
  6. Shepherd, M.C. Difmap: An Interactive Program for Synthesis Imaging. ASP Conf. Ser. 1997, 125, 77. [Google Scholar]
  7. Lister, M.L.; Cohen, M.H.; Homan, D.C.; Kadler, M.; Kellermann, K.I.; Kovalev, Y.Y.; Ros, E.; Savolainen, T.; Zensus, J.A. MOJAVE: Monitoring of jets in active galactic nuclei with VLBA experiments. VI. Kinematics analysis of a complete sample of blazar jets. Astron. J. 2009, 138, 1874–1892. [Google Scholar] [CrossRef]
  8. Jorstad, S.; Marscher, A. The VLBA-BU-BLAZAR Multi-Wavelength Monitoring Program. Galaxies 2016, 4, 47. [Google Scholar] [CrossRef]
  9. Dodson, R.; Rioja1, M.J.; Molina, S.N.; Gómez, J.L. High-precision astrometric millimeter very long baseline interferometry using a new method for multi-frequency calibration. Astrophys. J. 2017, 834, 10. [Google Scholar] [CrossRef]
  10. Rioja, M.J.; Dodson, R. High-precision astrometric millimeter very long baseline interferometry using a new method for atmospheric calibration. Astron. J. 2011, 141, 15. [Google Scholar] [CrossRef]
  11. Lyutikov, M.; Pariev, V.I.; Gabuzda, D.C. Polarization and structure of relativistic parsec-scale AGN jets. Mon. Not. Roy. Astron. Soc. 2005, 360, 869–891. [Google Scholar] [CrossRef]
  12. O’Sullivan, S.P.; Gabuzda, D.C. Three-dimensional magnetic field structure of six parsec-scale active galactic nuclei jets. Mon. Not. Roy. Astron. Soc. 2009, 393, 429–456. [Google Scholar] [CrossRef]
  13. Taylor, A.R.; Stil, J.M.; Sunstrum, C. A rotation measure image of the sky. Astrophys. J. 2009, 702, 1230–1236. [Google Scholar] [CrossRef]
  14. Asada, K.; Inoue, M.; Uchida, Y.; Kameno, S.; Fujisawa, K.; Iguchi, S.; Mutoh, M. A Helical Magnetic Field in the Jet of 3C 273. Publ. Astron. Soc. Jpn. 2002, 54, L39. [Google Scholar] [CrossRef]
  15. Hovatta, T.; Lister, M.L.; Aller, M.F.; Aller, H.D.; Homan, D.C.; Kovalev, Y.Y.; Pushkarev, A.B.; Savolainen, T. MOJAVE: Monitoring of Jets in Active Galactic Nuclei with VLBA Experiments. VIII. Faraday Rotation in Parsec-scale AGN Jets. Astron. J. 2012, 144, 105. [Google Scholar] [CrossRef]
  16. Gómez, J.L.; Marscher, A.P.; Jorstad, S.G.; Agudo, I.; Roca-Sogorb, M. Faraday rotation and polarization gradients in the jet of 3C 120: Interaction with the external medium and a helical magnetic field? Astrophys. J. 2008, 681, L69. [Google Scholar] [CrossRef]
  17. Gómez, J.L.; Lobanov, A.P.; Bruni, G.; Kovalev, Y.Y.; Marscher, A.P.; Jorstad, S.G.; Mizuno, Y.; Bach, U.; Sokolovsky, K.V.; Anderson, J.M. Probing the innermost regions of agn jets and their magnetic fields with Radioastron. I. Imaging BL Lacertae at 21 μ as resolution. Astrophys. J. 2016, 817, 14. [Google Scholar] [CrossRef]
Figure 1. Very Long Baseline Array (VLBA) images of BL Lacertae at 15 GHz. Contours are the total intensity, grey scale represents the linearly polarized intensity, and the short sticks indicate the electric vector position angle.
Figure 1. Very Long Baseline Array (VLBA) images of BL Lacertae at 15 GHz. Contours are the total intensity, grey scale represents the linearly polarized intensity, and the short sticks indicate the electric vector position angle.
Galaxies 05 00097 g001
Figure 2. Same as Figure 1 at 22 GHz.
Figure 2. Same as Figure 1 at 22 GHz.
Galaxies 05 00097 g002
Figure 3. Same as Figure 1 at 43 GHz.
Figure 3. Same as Figure 1 at 43 GHz.
Galaxies 05 00097 g003
Figure 4. Rotation measure (RM) map, obtained by combining images at 15-22-43 GHz, and transverse RM slices of BL Lacertae. Contours show the 15 GHz image, grey scale indicates the rotation measure, and bars plot the electric vector position angles (EVPAs).
Figure 4. Rotation measure (RM) map, obtained by combining images at 15-22-43 GHz, and transverse RM slices of BL Lacertae. Contours show the 15 GHz image, grey scale indicates the rotation measure, and bars plot the electric vector position angles (EVPAs).
Galaxies 05 00097 g004

Share and Cite

MDPI and ACS Style

Molina, S.N.; Gómez, J.L.; Dodson, R.; Rioja, M.J. Magnetic Field Studies in BL Lacertae through Faraday Rotation and a Novel Astrometric Technique. Galaxies 2017, 5, 97. https://doi.org/10.3390/galaxies5040097

AMA Style

Molina SN, Gómez JL, Dodson R, Rioja MJ. Magnetic Field Studies in BL Lacertae through Faraday Rotation and a Novel Astrometric Technique. Galaxies. 2017; 5(4):97. https://doi.org/10.3390/galaxies5040097

Chicago/Turabian Style

Molina, Sol N., José L. Gómez, Richard Dodson, and María J. Rioja. 2017. "Magnetic Field Studies in BL Lacertae through Faraday Rotation and a Novel Astrometric Technique" Galaxies 5, no. 4: 97. https://doi.org/10.3390/galaxies5040097

APA Style

Molina, S. N., Gómez, J. L., Dodson, R., & Rioja, M. J. (2017). Magnetic Field Studies in BL Lacertae through Faraday Rotation and a Novel Astrometric Technique. Galaxies, 5(4), 97. https://doi.org/10.3390/galaxies5040097

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop