| Issue |
A&A
Volume 700, August 2025
|
|
|---|---|---|
| Article Number | L22 | |
| Number of page(s) | 6 | |
| Section | Letters to the Editor | |
| DOI | https://doi.org/10.1051/0004-6361/202555706 | |
| Published online | 20 August 2025 | |
Letter to the Editor
Evidence of a non-equipartition energy regime in 1803+784
Core-shift and Faraday rotation measurements from simultaneous multi-frequency polarimetric VGOS observations
1
Observatorio Astronómico Nacional (OAN-IGN), Calle Alfonso XII 3, 28014 Madrid, Spain
2
Dpt. Astronomia i Astrofísica, Universitat de València, C/ Dr. Moliner 50, 46120 Burjassot, Spain
3
Observatori Astronòmic, Universitat de València, C/ Cat. José Beltrán 2, 46980 Paterna, Spain
4
Centro de Desarrollos Tecnológicos, Observatorio de Yebes (IGN), 19141 Yebes, Guadalajara, Spain
⋆ Corresponding author.
Received:
28
May
2025
Accepted:
30
July
2025
Context. Compact jets from active galactic nuclei (AGN) are commonly assumed to be in equipartition between particle and magnetic-field energy densities at the regions where the radio emission dominates at centimetre wavelengths. This assumption has significant implications for both jet physics and the accuracy of VLBI-based astrometry and geodesy.
Aims. We tested the validity of the energy equipartition hypothesis in AGN cores at centimetre wavelengths by analysing the blazar 1803+784 using simultaneous broadband full-polarization observations with the VLBI Global Observing System (VGOS).
Methods. We present VGOS observations of the blazar 1803+784 covering the 3–11 GHz frequency range. The data were processed using a dedicated calibration pipeline, followed by model fitting and multi-frequency imaging analysis. We measured the frequency-dependent core shift and mapped the spectral index and rotation measure (RM) across the source.
Results. We find a core-shift power-law index of kr = 0.73−0.19+0.12, significantly deviating from the expected equipartition value of kr = 1. This indicates that either the equipartition condition or the conical jet geometry, or both, are not fulfilled in the centimetre-wavelength core region. The wide frequency coverage of VGOS also allows us to decouple the Faraday rotation of the core into an internal component (RMI = 121 ± 8 rad m−2, produced in the core region) and an external component (RME = −44 ± 9 rad m−2, associated with a distant, extended medium that may also affect the polarization in downstream regions of the jet at larger scales).
Conlcusions. These results demonstrate the power of VGOS for high-fidelity simultaneous multi-frequency polarimetric studies of compact AGN jets, and underline the need to account for non-equipartition effects in both jet astrophysics and geodetic VLBI.
Key words: techniques: interferometric / techniques: polarimetric / galaxies: jets / quasars: general
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Compact jets from active galactic nuclei (AGN), particularly those associated with blazars and radio-loud quasars, are widely used in geodetic very long baseline interferometry (VLBI). However, the physical conditions and structure of these jets remain under active study. Jet models typically assume equipartition between the energy densities of relativistic particles and magnetic fields (e.g. Plavin et al. 2019; O’Sullivan & Gabuzda 2008; O’Sullivan & Gabuzda 2009). This assumption has been tested in some sources and found to be consistent with observations, as shown in Hada et al. (2011) and Sokolovsky et al. (2011). The equipartition hypothesis implies a specific frequency dependence for the apparent core position, characterized by a power-law shift with exponent kr ≈ 1 (Blandford & Königl 1979; Lobanov 1998):
When kr = 1, the frequency-dependent shift in the core position does not introduce a group delay. Therefore, group-delay VLBI astrometry is not affected by frequency-dependent core shifts (Porcas 2009). However, when kr ≠ 1, the group-delay position becomes frequency-dependent and deviates from the physical jet origin. This can introduce systematic errors in astrometric VLBI solutions. Kovalev et al. (2008) showed that such core shifts can bias the alignment between radio and optical reference frames, compromising the precision of celestial reference frame realizations such as the International Celestial Reference Frame (ICRF; Ma et al. 1998; Charlot et al. 2020) and the Gaia catalogue (Gaia Collaboration 2016, 2023; Lindegren et al. 2018).
The advent of the VLBI Global Observing System (VGOS), the next generation of geodetic VLBI observations, developed under the coordination of the International VLBI Service for Geodesy and Astrometry (IVS)1 (Nothnagel et al. 2017), allows a robust test of the equipartition hypothesis on a sample of radio-loud AGN routinely used for geodetic studies. Designed to achieve station position accuracies at the millimetre level, and with the capability for simultaneous multi-frequency observations, VGOS employs ultra-wideband receivers capable of covering approximately 2–15 GHz. The broad frequency coverage enables the measurement of group delays with picosecond-level precision, a critical requirement for high-accuracy geodesy and astrometry (Petrachenko et al. 2009; Niell et al. 2018).
Verifying whether AGN cores adhere to the expected kr = 1 regime is not only important for understanding jet physics, but is also crucial for ensuring the accuracy of geodetic measurements (e.g. Porcas 2009; Xu et al. 2022) and for aligning radio and optical reference frames (Kovalev et al. 2008; Xu et al. 2021a).
In this Letter we present an analysis of high-sensitivity VGOS observations of the prototypical blazar 1803+784. These multi-frequency data allowed us to derive a precise measurement of the core shift as a function of frequency. Additionally, we exploited the broad frequency coverage of VGOS to measure both internal and external Faraday rotation in the jet of 1803+784.
Finally, we produced high-resolution multi-frequency images of this compact jet, and created detailed maps of the spectral index and rotation measure. Furthermore, accurately characterizing the source structure in each band and its alignment through the core shift makes it possible to correct the structure effects described by Xu et al. (2021b), as demonstrated by the direct application of this pipeline to VGOS data in Jaron et al. (2024).
2. Observations
For this work we used observations from the IVS experiment VO2187, observed on 6–7 July 2022, involving eight antennas: Goddard (GS), Ishioka (IS), Kokee (K2), McDonald (MG), Onsala East (OE), Onsala West (OW), Westford (WF), and Yebes (YJ). Full details of the frequency set-up and observation strategy can be found in Pérez-Díez et al. (2024).
The data is divided into four frequency bands labelled A, B, C, and D, following the naming convention adopted by the VGOS community. Each band consists of eight spectral windows (spws) of 32 MHz and spans 0.5 GHz. The bands are centred at 3.25 GHz (A), 5.5 GHz (B), 6.75 GHz (C), and 10.5 GHz (D). In this Letter we focus on the source that was most extensively observed, the prototypical blazar 1803+784.
3. Methods
3.1. Calibration
The calibration strategy used in this work is described in Pérez-Díez et al. (2024). The first step consisted in converting the linear-polarization visibilities into circular basis using PolConvert (Martí-Vidal et al. 2016). After this, we applied a wide-band global fringe-fitting algorithm that solves simultaneously for dispersive and non-dispersive delays across the full VGOS frequency range. This method is optimized for the VGOS multi-band design and uses ionospheric models from GNSS IONEX maps produced by Jet Propulsion Laboratory (JPL)2 (Martire et al. 2024) to improve convergence.
The amplitude calibration was refined through a hybrid mapping procedure using Difmap (Shepherd 1997), combining iterative CLEAN deconvolution and self-calibration. We first performed several iterations of phase self-calibration, with decreasing solution intervals (ranging from 30 minutes to 1 minute), using total intensity (Stokes I) CLEAN models. Then, a few extra iterations of amplitude and phase self-calibration were performed, also based on Stokes I, using longer solution intervals (30–20 minutes). The instrumental polarization (D-terms) was refined using the PolSolve software (Martí-Vidal et al. 2021), which performs a global fit of the radio interferometer measurement equation (RIME), formulated by Hamaker et al. (1996), to minimize instrumental polarization leakage. This step ensures reliable full-Stokes calibration across the entire VGOS bandwidth, correcting for residual leakage effects that may remain after the initial polarization conversion.
3.2. Model fitting
To characterize the source structure, we performed model fitting to the visibilities following a procedure based on Martí-Vidal et al. (2014). The source brightness model is a combination of several Gaussian components: a central elliptical Gaussian for the optically thick core, and two Gaussians to represent the extended optically thin jet emission. The fitting was initially performed in total intensity (Stokes I) at C band (6.75 GHz), which provides the best compromise between resolution and signal-to-noise ratio.
To model the frequency-dependent core shift, the positions and sizes (but not the flux densities) of every optically thin jet component were kept fixed across bands and the position of the core component was left free. In this fit, the Gaussian shape parameters of the core were also held fixed to minimize systematics related to jet blending at the lower frequencies.
Lastly, the polarization structure (Stokes Q, U, and V) was modelled using the same Gaussian components fitted to Stokes I. In this case, only the flux densities of the individual components were allowed to vary to ensure consistency between the total intensity and polarimetric models across all frequencies.
3.3. Multi-frequency imaging
To reconstruct the source structure across the wide frequency coverage provided by VGOS, we employed a regularized maximum likelihood (RML) approach using the eht-imaging software (Chael et al. 2016, 2018). This technique uses closure phases, log closure amplitudes, and visibility amplitudes as observables. This combination is robust against station-based calibration errors, and thus simplifies the imaging process.
The imaging was carried out in multi-frequency mode, recently implemented in ehtim (Chael et al. 2023). This algorithm performs a log-log Taylor expansion of the image intensity as a function of frequency around a reference frequency ν0, which enables the simultaneous reconstruction of a total intensity map and a spectral index map, thus producing aligned images across frequencies. Maps for the other Stokes parameters were deconvolved by imaging each band separately with ehtim.
4. Results
4.1. Core shift and non-equipartition regime
We have measured a core shift of approximately 0.20 milliarcseconds (mas) between 3 and 11 GHz. In Fig. 1 we show the core positions relative to the one at the highest frequency (10.66 GHz). The shift is decomposed into its right ascension (RA) and declination (Dec) components, both exhibiting consistent trends with frequency. The total core-shift vector, Δr, follows a clear power-law behaviour with frequency. The wide bandwidth of VGOS, with eight spws per band, allows for relatively small uncertainties in the measurements. The uncertainty in band A is larger due to greater scatter among spws, possibly caused by its bigger beam.
![]() |
Fig. 1. Measured relative core shift as a function of frequency for 1803+784. The total shift vector, along with its components in right ascension and declination, are shown. The error bars represent 3σ uncertainties. The best-fit model is shown as a dashed green line. |
By fitting the frequency dependence of the core shift described in Eq. (1), we obtain a best-fit power-law index of
. This result represents a significant deviation from the standard equipartition condition (kr ∼ 1), which assumes a balance between the energy densities of relativistic particles and magnetic fields in a conical, freely expanding jet (Lobanov 1998). A value of kr < 1 indicates that the equipartition condition, the conical jet geometry, or both may not hold at the region we are observing at these frequencies. These findings underline the necessity of carefully evaluating the core-shift behaviour in blazars when deriving physical jet parameters, such as magnetic field strength, and when conducting geodetic studies.
4.2. Rotation measure
The wide simultaneous frequency coverage provided by VGOS allows us, for the first time, to separate the internal and external contributions of the Faraday rotation observed in 1803+784. While a distant plasma screen introduces an external rotation measure, RME, that rotates the “electric vector position angle (EVPA), χ, as Δχ = RME Δλ2, an internal (and/or surrounding) screen around the optically thick core produces a rotation measure RMI that depends on frequency as RMI ∝ νa, as derived in Jorstad et al. (2007). The exponent a describes the radial dependence of the electron density, assuming a power-law decrease with distance from the black hole: ne ∝ r−a (Hovatta et al. 2019). Reported values of the exponent a in the literature typically range from 0.4 to 4, with average values of around 2 (Jorstad et al. 2007; O’Sullivan & Gabuzda 2009; Kravchenko et al. 2017; Hovatta et al. 2019). This average value of a = 2 suggests that the Faraday rotation takes place in a sheath surrounding a conically expanding jet. Lower values could be attributed to a more collimated jet structure (Hovatta et al. 2019), different geometry of the jet, flaring activity, or a possible filamentary structure of the Faraday screen (Kravchenko et al. 2017).
By fitting the polarization angle, χ, of the VLBI core as a function of RME and RMI (see Fig. 2), we measure an external Faraday rotation of RMEcore = −44 ± 9 rad m−2 and an internal rotation of RMIcore = 121 ± 8 rad m−2 at the reference frequency of 5.5 GHz. When fitting for the exponent of RMIcore with ν, the result is a = 2.18 ± 0.95. Due to large uncertainties, we cannot confirm this value, but it seems to be consistent with the theoretical expectation for a common expanding jet.
![]() |
Fig. 2. Fit of the polarization angle as a function of λ2 for the core and jet regions of 1803+784. The shaded regions represent the 3σ confidence intervals of the model fits. Band D data have been averaged because of the low S/N in this band. |
On the other hand, the plasma in the jet components is optically thin, so there is no core shift and the effects of RMEjet and RMIjet both follow a λ2 relation. In this case, the fitted rotation measure is
rad m−2, and χ follows a clean λ2 dependence, consistent with the presence of a total optically thin Faraday screen (Fig. 2).
If we consider that the effect of the distant plasma screen does not change between the core and the inner jet, we can further infer the RMIjet contribution by taking RMEjet ∼ RMEcore. With this assumption, the internal rotation measure in the optically thin jet is approximately RMIjet ∼ 93 ± 10 rad m−2, which is lower than RMIcore. This is in agreement with a jet magnetic field that decreases with distance to the black hole, as is implied in the model, RM ∝ νp, used to fit the χ at the core. Hence, the Faraday rotation at the core and the jet are both self-consistent with the picture of having two Faraday screens, one distant and extended and a second one located in the immediate jet surroundings (and/or internal to the jet itself).
4.3. Multi-frequency image
Figure 3 presents the multi-frequency imaging results for 1803+784. Panels (a)–(c) show the total intensity (Stokes I) maps at three frequency bands (A, C, and D), with contour lines superimposed over colour raster images. The contours are drawn at 0.08, 0.2, 0.5, 1, 5, 10, 20, 40, 60, and 80% of the peak intensity in each band. The colour raster highlights the detailed jet structure at each frequency.
![]() |
Fig. 3. Multi-frequency polarimetric image of 1803+784. Panels (a)–(c) show the total intensity maps at bands A, C, and D in logarithmic scale, normalized to their peak. The dashed vertical line indicates the core position at band D, marking the core shift. Panel (e) displays all contours from the four bands superposed, along with the core and jet Gaussian models. Panel (f) shows the spectral index map (α). Panel (g) presents the rotation measure (RM) map. Only regions above 0.8% of the peak intensity are included. |
A vertical dashed grey line marks the RA position of the core in band D. A shift of the core to the west with increasing frequency is clearly noticeable and compatible with the one obtained in the model fitting.
Panel (e) overlays the total intensity contours from all bands, allowing a direct visualization of the core shift. The Gaussian components derived from the model fitting described in Sect. 3.2 are also plotted, with colours corresponding to each band: red (A), orange (B), green (C), and purple (D).
Panel (f) displays the spectral index (α) map, computed from the multi-frequency imaging described in Sect. 3.3. The core region exhibits an inverted spectrum (α > 0) consistent with synchrotron self-absorption, while the jet shows a steep spectrum (α < 0), as expected for optically thin synchrotron emission.
Finally, panel (g) presents the rotation measure (RM) map derived from the multi-band polarimetric data. The RM structure reveals a complex distribution, with higher RM values near the core and lower values along the jet.
5. Discussion and conclusions
Our results provide the first direct evidence for a significant departure from equipartition conditions in the compact jet of 1803+784, based on simultaneous multi-frequency and full-polarization VLBI observations. The measured core-shift power-law index (
) strongly deviates from the standard equipartition prediction (kr ∼ 1; Lobanov 1998), indicating that either the energy density balance between particles and magnetic fields, the conical jet geometry, or both assumptions, are not fully satisfied in the innermost jet regions during this epoch.
Kovalev et al. (2008) obtained a core-shift value of 0.72 mas between 2.3 and 8.6 GHz for 1803+784. Extrapolating our fit to a frequency of 2.3 GHz, a core-shift value of 0.61 mas is obtained, although this difference may be due to the fact that the observations are separated by 17 years.
More recent studies of another source, NGC 315, obtain kr values similar to those obtained in this work. Boccardi et al. (2021) fits core-shift values to six different frequencies (in a range similar to ours), obtaining a value of kr = 0.84 ± 0.06. However, a reanalysis of the data by Ricci et al. (2022), where they add an additional point at 43 GHz, fits two kr values for different frequency ranges: kr = 0.93 ± 0.01 for frequencies lower than 10 GHz and kr = 0.57 ± 0.17 for higher frequencies. According to their interpretation, these two distinct regimes may test the parabolic and conical zones of the jet, that is, the transition between a non-equipartition zone and an equipartition zone.
Following this hypothesis, we performed the same 10 GHz separation. Fitting without data at D band (10.5 GHz) we obtain a value of
for the three frequencies lower than 10 GHz. This value is closer to unity, without reaching equipartition, suggesting that the scenario proposed by Ricci et al. (2022) may also apply to 1803+784, but the transition would occur at a lower frequency. Observations at frequencies below our band A would be necessary to confirm this. We also note that the observations of Boccardi et al. (2021), Ricci et al. (2022) are separated by more than 20 years, while the VGOS observations presented here are truly simultaneous.
The ability to separate internal and external Faraday rotation contributions, thanks to VGOS’s wide instantaneous bandwidth, represents another major advance. Our measurement of a positive internal RM and a negative external RM near the core supports the presence of two distinct Faraday screens, likely associated with a stratified plasma layer surrounding the base of the jet. Our inferred external Faraday rotation (RME = −44 ± 9 rad m−2 after correcting for internal contributions) is broadly consistent with previous measurements at larger scales (i.e. farther downstream from the core) by Gabuzda & Chernetskii (2003).
The spectral index and RM maps derived from our multi-band imaging show an excellent agreement with previous studies such as the MOJAVE programme (Hovatta et al. 2012) This confirms the reliability of VGOS observations and our calibration pipeline in obtaining astrophysical results.
These results highlight the potential of VGOS for high-fidelity multi-frequency full-polarization studies of compact AGN jets. In particular, they demonstrate the importance of accounting for non-equipartition effects when deriving jet parameters or performing high-precision VLBI astrometry.
Acknowledgments
This work has been partially supported by the Generalitat Valenciana GenT Project CIDEGENT/2018/021 and by the MICINN Research Projects PID2019-108995GB-C22 and PID2022-140888NB-C22. This work has been supported by the grant PRE2020-092200 funded by MCIN/AEI/ 10.13039/501100011033 and by ESF invest in your future.
References
- Blandford, R. D., & Königl, A. 1979, ApJ, 232, 34 [Google Scholar]
- Boccardi, B., Perucho, M., Casadio, C., et al. 2021, A&A, 647, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chael, A. A., Johnson, M. D., Narayan, R., et al. 2016, ApJ, 829, 11 [Google Scholar]
- Chael, A. A., Johnson, M. D., Bouman, K. L., et al. 2018, ApJ, 857, 23 [Google Scholar]
- Chael, A., Issaoun, S., Pesce, D. W., et al. 2023, ApJ, 945, 40 [NASA ADS] [CrossRef] [Google Scholar]
- Charlot, P., Jacobs, C. S., Gordon, D., et al. 2020, A&A, 644, A159 [EDP Sciences] [Google Scholar]
- Gabuzda, D. C., & Chernetskii, V. A. 2003, MNRAS, 339, 669 [Google Scholar]
- Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hada, K., Doi, A., Kino, M., et al. 2011, Nature, 477, 185 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaker, J. P., Bregman, J. D., & Sault, R. J. 1996, A&AS, 117, 137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hovatta, T., Lister, M. L., Aller, M. F., et al. 2012, AJ, 144, 105 [Google Scholar]
- Hovatta, T., O’Sullivan, S., Martí-Vidal, I., Savolainen, T., & Tchekhovskoy, A. 2019, A&A, 623, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jaron, F., Baldreich, L., Johm, J., et al. 2024, in Proceedings of the 16th EVN Symposium, eds. E. Ros, P. Benke, S. A. Dzib, & I. Rottmann, 145 [Google Scholar]
- Jorstad, S. G., Marscher, A. P., Stevens, J. A., et al. 2007, AJ, 134, 799 [NASA ADS] [CrossRef] [Google Scholar]
- Kovalev, Y. Y., Lobanov, A. P., Pushkarev, A. B., & Zensus, J. A. 2008, A&A, 483, 759 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kravchenko, E. V., Kovalev, Y. Y., & Sokolovsky, K. V. 2017, MNRAS, 467, 83 [NASA ADS] [Google Scholar]
- Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lobanov, A. P. 1998, A&A, 330, 79 [NASA ADS] [Google Scholar]
- Ma, C., Arias, E. F., Eubanks, T. M., et al. 1998, AJ, 116, 516 [Google Scholar]
- Martire, L., Runge, T. F., Meng, X., et al. 2024, J. Geod., 98, 44 [Google Scholar]
- Martí-Vidal, I., Vlemmings, W. H. T., Muller, S., & Casey, S. 2014, A&A, 563, A136 [Google Scholar]
- Martí-Vidal, I., Roy, A., Conway, J., & Zensus, A. J. 2016, A&A, 587, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martí-Vidal, I., Mus, A., Janssen, M., de Vicente, P., & González, J. 2021, A&A, 646, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Niell, A., Barrett, J., Burns, A., et al. 2018, Radio Sci., 53, 1269 [NASA ADS] [CrossRef] [Google Scholar]
- Nothnagel, A., Artz, T., Behrend, D., & Malkin, Z. 2017, J. Geod., 91, 711 [Google Scholar]
- O’Sullivan, S. P., & Gabuzda, D. C. 2008, The Role of VLBI in the Golden Age for Radio Astronomy, 9, 16 [Google Scholar]
- O’Sullivan, S. P., & Gabuzda, D. C. 2009, MNRAS, 400, 26 [Google Scholar]
- Pérez-Díez, V., Martí-Vidal, I., Albentosa-Ruiz, E., et al. 2024, A&A, 688, A151 [Google Scholar]
- Petrachenko, B., Niell, A., Behrend, D., et al. 2009, Design Aspects of the VLBI2010 System., NASA/TM-2009-214180, 2009 [Google Scholar]
- Plavin, A. V., Kovalev, Y. Y., Pushkarev, A. B., & Lobanov, A. P. 2019, MNRAS, 485, 1822 [Google Scholar]
- Porcas, R. W. 2009, A&A, 505, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ricci, L., Boccardi, B., Nokhrina, E., et al. 2022, A&A, 664, A166 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shepherd, M. C. 1997, Astron. Data Anal. Software Syst. VI, 125, 77 [Google Scholar]
- Sokolovsky, K. V., Kovalev, Y. Y., Pushkarev, A. B., & Lobanov, A. P. 2011, A&A, 532, A38 [CrossRef] [EDP Sciences] [Google Scholar]
- Xu, M. H., Lunz, S., Anderson, J. M., et al. 2021a, A&A, 647, A189 [EDP Sciences] [Google Scholar]
- Xu, M. H., Savolainen, T., Zubko, N., et al. 2021b, J. Geophys. Res.: Solid Earth, 126, e2020JB021238 [Google Scholar]
- Xu, M. H., Savolainen, T., Anderson, J. M., et al. 2022, A&A, 663, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Core-shift fitting
To determine the power-law index kr describing the frequency dependence of the core shift, we fitted the averaged total shift values Δr as a function of frequency using the model Δr(ν) = A ⋅ ν−1/kr. The fit was performed using a weighted least-squares minimization.
To estimate the uncertainty of kr, we generated Monte Carlo simulations of the fit parameters, sampling from the multi-variate normal distribution defined by the best-fit covariance matrix. The asymmetric uncertainties were derived from the 16th and 84th percentiles of the resulting distribution.
The data used for the fit is presented in Table A.1. For each band, the spw measurements were averaged, and the fitting was carried out using these averaged values, which are emphasized in boldface.
Core-shift measurements with associated uncertainties.
Appendix B: Rotation measure data
To estimate the Faraday rotation measures (RM) for the core and jet components, we fitted the observed electric vector position angles (EVPA) as a function of wavelength squared (λ2) using a model that accounts for both external and internal Faraday rotation. For the core component, the fit includes an intrinsic rotation term with a power-law index describing the frequency dependence of the internal Faraday rotation.
The fits were performed by minimizing the residuals using a non-linear least-squares method. To quantify the uncertainties of the fitted parameters, we generated Monte Carlo simulations by randomly perturbing the EVPA data within their measurement uncertainties. The final uncertainties were estimated from the standard deviation of the resulting parameter distributions.
The individual EVPA measurements used for the fits are summarized in Table B.1 (for the core) and Table B.2 (for the jet).
EVPA data for the Core component.
EVPA data for the Jet component.
All Tables
All Figures
![]() |
Fig. 1. Measured relative core shift as a function of frequency for 1803+784. The total shift vector, along with its components in right ascension and declination, are shown. The error bars represent 3σ uncertainties. The best-fit model is shown as a dashed green line. |
| In the text | |
![]() |
Fig. 2. Fit of the polarization angle as a function of λ2 for the core and jet regions of 1803+784. The shaded regions represent the 3σ confidence intervals of the model fits. Band D data have been averaged because of the low S/N in this band. |
| In the text | |
![]() |
Fig. 3. Multi-frequency polarimetric image of 1803+784. Panels (a)–(c) show the total intensity maps at bands A, C, and D in logarithmic scale, normalized to their peak. The dashed vertical line indicates the core position at band D, marking the core shift. Panel (e) displays all contours from the four bands superposed, along with the core and jet Gaussian models. Panel (f) shows the spectral index map (α). Panel (g) presents the rotation measure (RM) map. Only regions above 0.8% of the peak intensity are included. |
| In the text | |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.



