pith. machine review for the scientific record. sign in

arxiv: 2604.20095 · v2 · submitted 2026-04-22 · 🌌 astro-ph.HE · astro-ph.GA

Recognition: unknown

Detections of nearly bias-free core shifts with 5-30 μas precisions at 8-43 GHz in BL Lacertae

Niu Liu , Jun Yang , Xiaopeng Cheng , Ai-Ling Zeng , Wen Chen , Xiao-Long Yang , Xiaoyu Hong , Xia-Xuan Zhang , Jia-Cheng Liu , Zi Zhu

Authors on Pith no claims yet

Pith reviewed 2026-05-09 23:59 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords BL Lacertaecore shiftVLBI astrometryrelativistic jetactive galactic nucleusradio frequencyenergy equipartitionblazar
0
0 comments X

The pith

VLBI astrometry detects a 250 microarcsecond core shift in BL Lacertae between 8.4 and 43.2 GHz, with the shift index consistent with equipartition.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper measures the frequency-dependent positions of the compact core in the jet of BL Lacertae using high-precision VLBI observations at five frequencies from 8.4 to 43.2 GHz. By using inverse phase-referencing to a nearby calibrator, they achieve position precisions of 5 to 30 microarcseconds and find a total core shift of 250 microarcseconds between the lowest and highest frequencies. These shifts indicate an optically thick region in the upstream jet, and the derived power-law index for the frequency dependence is consistent with the standard model where magnetic and particle energies are in equipartition. Such measurements help locate the jet launching point and constrain the physical conditions near the central black hole in active galaxies.

Core claim

The observations reveal that the apparent core position shifts by up to 250 μas from 8.4 to 43.2 GHz, scaling as frequency to the power of minus one over k_r, where k_r is measured as 1.18 with uncertainties allowing consistency with 1. This supports the presence of an optically thick region upstream in the jet and energy equipartition between particles and magnetic fields, even during a flaring state.

What carries the argument

The inverse phase-referencing technique with a nearby steep-spectrum calibrator, which enables nearly bias-free astrometric positions, combined with fitting the core shifts to a power-law model in frequency.

Load-bearing premise

The nearby calibrator introduces no residual position biases and the flaring jet structure follows the assumed power-law dependence for core shifts.

What would settle it

Repeated VLBI observations at the same frequencies during a quiescent state of BL Lacertae yielding a core-shift index outside the current uncertainty range would challenge the result.

Figures

Figures reproduced from arXiv: 2604.20095 by Ai-Ling Zeng, Jia-Cheng Liu, Jun Yang, Niu Liu, Wen Chen, Xiao-Long Yang, Xiaopeng Cheng, Xiaoyu Hong, Xia-Xuan Zhang, Zi Zhu.

Figure 1
Figure 1. Figure 1: Core shift measurements of BL Lacertae. (a) Two-dimensional core shifts at 8.4, 12.4, 15.2, and 23.6 GHz relative to 43.2 GHz. The dashed line denotes the best-fit core-shift model, and the gray region shows the scatter from 1000 randomly sampled MCMC realizations. (b) 43 GHz VLBA image. (c) Projected core displacements along the jet direction, adopting a position angle of 197◦ . The red dashed line shows … view at source ↗
Figure 2
Figure 2. Figure 2: Left: Final phase solutions at 43.2 GHz over the first observing hour during the first hour of observations. Red symbols denote the phase solutions of BL Lacertae, and blue symbols show the interpolated phases transferred from BL Lacertae to J2203+4208. Right: Residual phases for BL Lacertae after application of the final phase solutions at 43.2 GHz during the same time interval. B. MODEL-FITTING RESULTS … view at source ↗
Figure 3
Figure 3. Figure 3: VLBA maps of J2203+4208 at 8.4 GHz, 12.4 GHz, 15.2 GHz, 23.6 GHz, and 43.2 GHz, respectively. The contours start from 3σ off-source noise level. Gaussian model was adopted to fit the calibrated visibility data of J2203+4208 in order to determine its integrated flux densities and relative positions with respect to the reference position using the modelfit task. The results of model￾fitting for J2203+4208 ar… view at source ↗
Figure 4
Figure 4. Figure 4: Radio spectrum of the reference source J2203+4208. The blue open circles denote the integrated flux densities obtained in this work, while the gray zone indicates the 1-σ confidence interval [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: VLBA maps of BL Lacertae at 8.4 GHz, 12.4 GHz, 15.2 GHz, 23.6 GHz, and 43.2 GHz, respectively. The contours start from 3σ off-source noise level. Greisen, E. W. 2003, in Astrophysics and Space Science Library, Vol. 285, Information Handling in Astronomy - Historical Vistas, ed. A. Heck, 109, doi: 10.1007/0-306-48080-8_7 Hada, K., Doi, A., Kino, M., et al. 2011, Nature, 477, 185, doi: 10.1038/nature10387 Hu… view at source ↗
Figure 6
Figure 6. Figure 6: Posterior probability distributions of the core shift parameters derived from the MCMC analysis. Pearson correlation coefficients between parameters are indicated, with red and blue denoting positive and negative correlations, respectively. Marscher, A. 2010, in The Jet Paradigm: From Microquasars to Quasars, ed. T. Belloni (Berlin, Heidelberg: Springer Berlin Heidelberg), 173–201, doi: 10.1007/978-3-540-7… view at source ↗
read the original abstract

When a radio jet is partially optically thick in the launching region, its apparent compact core may display frequency-dependent positional shifts. High-precision astrometric measurements of core shifts enable astronomers to pinpoint the jet's origin and place tight constraints on the magnetic field. BL Lacertae, the archetypal BL Lac object, hosts a highly variable and well-collimated jet. To independently constrain its innermost core shifts, we conducted very long baseline interferometric (VLBI) observations at 8.4, 12.4, 15.2, 23.6, and 43.2 GHz. By exploiting a nearby (13.3 arcmin) steep-spectrum calibrator (NVSS J220340+420839) through inverse phase-referencing VLBI astrometry, we detect nearly unbiased two-dimensional core shift measurements with state-of-the-art precisions of 5-30 $\mu$as, which are significant at $>3\sigma$ confidence. The core shift between 8.4 and 43.2 GHz reaches 250 $\mu$as. The apparent core shifts scale with frequency as $\nu^{-1/k_r}$, implying the existence of an optically thick region in the upstream of jet. The derived core-shift index, $k_r\!=\!1.18^{+0.59}_{-0.34}$, is consistent, within uncertainties, with the canonical $k_r\!=\!1$ expected under energy equipartition between the jet particle and magnetic field energy densities, while allowing for modest deviations given that BL Lacertae was captured in a flaring state.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper claims to have obtained high-precision, nearly bias-free two-dimensional core shift measurements in BL Lacertae at five radio frequencies (8.4–43.2 GHz) using inverse phase-referencing VLBI astrometry with a nearby calibrator. The largest measured shift is 250 μas, and the frequency dependence is fitted to yield a core-shift index k_r = 1.18^{+0.59}_{-0.34}, consistent with the canonical value of 1 under equipartition assumptions, despite the source being in a flaring state.

Significance. Should the core positions prove robust against potential contamination from flaring-induced jet components, these results would offer state-of-the-art astrometric constraints on the jet launching region in a prototypical BL Lac object. The reported precisions of 5–30 μas and >3σ significances highlight the capability of the inverse phase-referencing technique for such measurements.

major comments (1)
  1. [Core-shift analysis and power-law fitting] The derivation of k_r relies on the assumption that the measured core positions follow r ∝ ν^{-1/k_r} as expected for a steady, self-absorbed conical jet. Given that the observations occurred during a flaring state (as stated in the abstract), it is necessary to demonstrate that transient components do not introduce frequency-dependent or independent offsets that could alter the best-fit k_r within its reported uncertainties.
minor comments (2)
  1. Providing the individual core position measurements and their error budgets in a table would allow readers to reproduce the k_r fit and assess the impact of any potential biases.
  2. [Abstract] The term 'nearly bias-free' is used; a short clause explaining the role of the steep-spectrum calibrator in mitigating position biases would enhance accessibility.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive feedback and positive assessment of the potential significance of our core-shift measurements. We address the major comment below and agree that further discussion is needed to strengthen the manuscript.

read point-by-point responses
  1. Referee: [Core-shift analysis and power-law fitting] The derivation of k_r relies on the assumption that the measured core positions follow r ∝ ν^{-1/k_r} as expected for a steady, self-absorbed conical jet. Given that the observations occurred during a flaring state (as stated in the abstract), it is necessary to demonstrate that transient components do not introduce frequency-dependent or independent offsets that could alter the best-fit k_r within its reported uncertainties.

    Authors: We agree that the flaring state warrants explicit discussion to confirm that transient components have not biased the fitted k_r. The inverse phase-referencing technique with the nearby calibrator NVSS J220340+420839 is designed to deliver nearly unbiased astrometric positions by calibrating out common-mode errors, yielding the reported 5–30 μas precisions and >3σ detections in two dimensions. The measured 250 μas shift between 8.4 and 43.2 GHz and the resulting k_r = 1.18^{+0.59}_{-0.34} already incorporate substantial uncertainties that comfortably include the canonical value of 1. In the revised manuscript we will add a dedicated subsection discussing possible contamination from flaring-induced jet components, arguing that the observed monotonic frequency dependence, the consistency of the two-dimensional shifts across all five frequencies, and the external-calibrator reference together limit any such offsets to levels smaller than the quoted uncertainties. We will also note that any residual frequency-dependent bias would likely increase the scatter in the power-law fit rather than systematically shift k_r outside the reported range. revision: yes

Circularity Check

0 steps flagged

No significant circularity in core-shift measurements and k_r derivation

full rationale

The paper reports direct VLBI astrometric observations at five frequencies (8.4–43.2 GHz) using inverse phase-referencing to a nearby calibrator, yielding measured two-dimensional core shifts with 5–30 μas precision. The core-shift index k_r is obtained by fitting these measured positions to the standard functional form r ∝ ν^{-1/k_r}. This is ordinary parameter estimation from independent data rather than a self-definitional loop, a fitted input relabeled as prediction, or any load-bearing self-citation. The subsequent statement that the fitted value is consistent with the canonical k_r=1 is a post-hoc comparison, not an input assumption that forces the result. No uniqueness theorems, ansatzes smuggled via citation, or renaming of known results appear in the derivation chain. The central claims rest on the raw positional measurements and their statistical uncertainties, which remain falsifiable against external VLBI data.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The result is a direct measurement; the only fitted quantity is k_r from the power-law model. No new entities are postulated. Standard VLBI assumptions (phase stability, calibrator position accuracy) are invoked but not derived here.

free parameters (1)
  • k_r
    Core-shift index fitted to the observed frequency-dependent shifts.
axioms (1)
  • domain assumption The jet core position follows a power-law shift with frequency of the form ν^{-1/k_r}.
    Invoked to interpret the measured shifts and derive k_r.

pith-pipeline@v0.9.0 · 5638 in / 1305 out tokens · 28342 ms · 2026-05-09T23:59:23.432648+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

57 extracted references · 56 canonical work pages · 1 internal anchor

  1. [1]

    , year = 2008, month = apr, volume =

    The inner jet of an active galactic nucleus as revealed by a radio-to- -ray outburst. , year = 2008, month = apr, volume =. doi:10.1038/nature06895 , adsurl =

  2. [2]

    Planck 2018 results. VI. Cosmological parameters. , keywords =. doi:10.1051/0004-6361/201833910 , archivePrefix =. 1807.06209 , primaryClass =

  3. [3]

    E., et al

    SciPy 1.0: fundamental algorithms for scientific computing in Python. Nature Medicine , keywords =. doi:10.1038/s41592-019-0686-2 , archivePrefix =. 1907.10121 , primaryClass =

  4. [4]

    Computing in Science and Engineering , keywords =

    Matplotlib: A 2D Graphics Environment. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2007.55 , adsurl =

  5. [5]

    P., Tollerud, E

    Astropy: A community Python package for astronomy. , keywords =. doi:10.1051/0004-6361/201322068 , archivePrefix =. 1307.6212 , primaryClass =

  6. [6]

    W., Lang, D., & Goodman, J

    emcee: The MCMC Hammer. , keywords =. doi:10.1086/670067 , archivePrefix =. 1202.3665 , primaryClass =

  7. [7]

    The Astronomical Journal , author =

    The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package. , keywords =. doi:10.3847/1538-3881/aabc4f , archivePrefix =. 1801.02634 , primaryClass =

  8. [8]

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package. , keywords =. doi:10.3847/1538-4357/ac7c74 , archivePrefix =. 2206.14220 , primaryClass =

  9. [9]

    , keywords =

    The NRAO VLA Sky Survey. , keywords =. doi:10.1086/300337 , adsurl =

  10. [10]

    C., & Varoquaux, G

    The NumPy Array: A Structure for Efficient Numerical Computation. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2011.37 , archivePrefix =. 1102.1523 , primaryClass =

  11. [11]

    2016, in ELPUB, 87–90, doi: 10.3233/978-1-61499-649-1-87

    Jupyter Notebooks a publishing format for reproducible computational workflows. IOS Press , year = 2016, pages =. doi:10.3233/978-1-61499-649-1-87 , adsurl =

  12. [12]

    , keywords =

    Relativistic jets as compact radio sources. , keywords =. doi:10.1086/157262 , adsurl =

  13. [13]

    A., Chandler, C

    The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science Case and Survey Design. , keywords =. doi:10.1088/1538-3873/ab63eb , archivePrefix =. 1907.01981 , primaryClass =

  14. [14]

    AJ , eprint =

    Polarimetric Observations of 15 Active Galactic Nuclei at High Frequencies: Jet Kinematics from Bimonthly Monitoring with the Very Long Baseline Array. , keywords =. doi:10.1086/444593 , archivePrefix =. astro-ph/0502501 , primaryClass =

  15. [15]

    BASS. XXII. The BASS DR2 AGN Catalog and Data. , keywords =. doi:10.3847/1538-4365/ac6c05 , archivePrefix =. 2207.12432 , primaryClass =

  16. [16]

    Astronomical Data Analysis Software and Systems VI , year = 1997, editor =

    Difmap: an Interactive Program for Synthesis Imaging. Astronomical Data Analysis Software and Systems VI , year = 1997, editor =

  17. [17]

    , keywords =

    Ultracompact jets in active galactic nuclei. , keywords =. doi:10.48550/arXiv.astro-ph/9712132 , archivePrefix =. astro-ph/9712132 , primaryClass =

  18. [18]

    The Radio Fundamental Catalog. I. Astrometry. , keywords =. doi:10.3847/1538-4365/ad8c36 , archivePrefix =. 2410.11794 , primaryClass =

  19. [19]

    Information Handling in Astronomy - Historical Vistas , year = 2003, series =

    AIPS, the VLA, and the VLBA. Information Handling in Astronomy - Historical Vistas , year = 2003, editor =. doi:10.1007/0-306-48080-8_7 , adsurl =

  20. [20]

    , year = 2014, month = jun, volume =

    Dynamically important magnetic fields near accreting supermassive black holes. , year = 2014, month = jun, volume =. doi:10.1038/nature13399 , adsurl =

  21. [21]

    MOJAVE. XV. VLBA 15 GHz Total Intensity and Polarization Maps of 437 Parsec-scale AGN Jets from 1996 to 2017. , keywords =. doi:10.3847/1538-4365/aa9c44 , archivePrefix =. 1711.07802 , primaryClass =

  22. [22]

    , keywords =

    Significant core shift variability in parsec-scale jets of active galactic nuclei. , keywords =. doi:10.1093/mnras/stz504 , archivePrefix =. 1811.02544 , primaryClass =

  23. [23]

    , keywords =

    Tracing Milky Way scattering by compact extragalactic radio sources. , keywords =. doi:10.1093/mnras/stac1898 , archivePrefix =. 2201.04359 , primaryClass =

  24. [24]

    MOJAVE: Monitoring of Jets in Active galactic nuclei with VLBA Experiments. IX. Nuclear opacity. , keywords =. doi:10.1051/0004-6361/201219173 , archivePrefix =. 1207.5457 , primaryClass =

  25. [25]

    , keywords =

    Time variability of the core-shift effect in the blazar 3C 454.3. , keywords =. doi:10.1051/0004-6361/202243435 , archivePrefix =. 2209.13301 , primaryClass =

  26. [26]

    , year = 2011, month = sep, volume = 477, pages =

    An origin of the radio jet in M87 at the location of the central black hole. , year = 2011, month = sep, volume =. doi:10.1038/nature10387 , adsurl =

  27. [27]

    A VLBA survey of the core shift effect in AGN jets. I. Evidence of dominating synchrotron opacity. , keywords =. doi:10.1051/0004-6361/201016072 , archivePrefix =. 1103.6032 , primaryClass =

  28. [28]

    , keywords =

    Testing particle acceleration in blazar jets with continuous high-cadence optical polarization observations. , keywords =. doi:10.1051/0004-6361/202451037 , archivePrefix =. 2406.15554 , primaryClass =

  29. [29]

    , keywords =

    Precise radio astrometry and new developments for the next-generation of instruments. , keywords =. doi:10.1007/s00159-020-00126-z , archivePrefix =. 2010.02156 , primaryClass =

  30. [30]

    , keywords =

    Magnetic field strength and spectral distribution of six parsec-scale active galactic nuclei jets. , keywords =. doi:10.1111/j.1365-2966.2009.15428.x , archivePrefix =. 0907.5211 , primaryClass =

  31. [31]

    , keywords =

    On the possible core shift break in relativistic jets. , keywords =. doi:10.1093/mnras/stae2515 , archivePrefix =. 2411.02925 , primaryClass =

  32. [32]

    , keywords =

    Dissecting the AGN Disk-Jet System with Joint VLBI-Gaia Analysis. , keywords =. doi:10.3847/1538-4357/aaf650 , archivePrefix =. 1808.05115 , primaryClass =

  33. [33]

    A&A , author =

    Single-epoch VLBI imaging study of bright active galactic nuclei at 2 GHz and 8 GHz. , keywords =. doi:10.1051/0004-6361/201219352 , archivePrefix =. 1205.5559 , primaryClass =

  34. [34]

    , keywords =

    A connection between -ray and parsec-scale radio flares in the blazar 3C 273. , keywords =. doi:10.1093/mnras/stx710 , archivePrefix =. 1703.07976 , primaryClass =

  35. [35]

    Studies of the Jet in Bl Lacertae. I. Recollimation Shock and Moving Emission Features. , keywords =. doi:10.1088/0004-637X/787/2/151 , archivePrefix =. 1404.0976 , primaryClass =

  36. [36]

    Probing the Innermost Regions of AGN Jets and Their Magnetic Fields with RadioAstron. I. Imaging BL Lacertae at 21 Microarcsecond Resolution. , keywords =. doi:10.3847/0004-637X/817/2/96 , archivePrefix =. 1512.04690 , primaryClass =

  37. [37]

    , year = 1968, month = may, volume =

    BL Lac identified as a Radio Source. , year = 1968, month = may, volume =. doi:10.1038/218663a0 , adsurl =

  38. [38]

    , keywords =

    VLBI study of 1038+528A and B : discovery of wavelength dependence of peak brightness location. , keywords =. doi:10.1086/161592 , adsurl =

  39. [39]

    Relativistic Jets in the Radio Reference Frame Image Database. II. Blazar Jet Accelerations from the First 10 Years of Data (1994-2003). , keywords =. doi:10.1088/0004-637X/758/2/84 , archivePrefix =. 1208.4399 , primaryClass =

  40. [40]

    , keywords =

    The core shift effect in the blazar 3C 454.3. , keywords =. doi:10.1093/mnras/stt2133 , archivePrefix =. 1307.4100 , primaryClass =

  41. [41]

    Galaxies , keywords =

    The VLBA-BU-BLAZAR Multi-Wavelength Monitoring Program. Galaxies , keywords =. doi:10.3390/galaxies4040047 , adsurl =

  42. [42]

    , keywords =

    Aligning VLBI images of active galactic nuclei at different frequencies. , keywords =. doi:10.1111/j.1365-2966.2008.13087.x , archivePrefix =. 0809.3313 , primaryClass =

  43. [43]

    Spectral distributions in compact radio sources. I. Imaging with VLBI data. , keywords =. doi:10.1051/aas:1998446 , archivePrefix =. astro-ph/9804112 , primaryClass =

  44. [44]

    , keywords =

    Radio astrometry with chromatic AGN core positions. , keywords =. doi:10.1051/0004-6361/200912846 , archivePrefix =. 0909.0933 , primaryClass =

  45. [45]

    , keywords =

    Erratic Flaring of BL Lac in 2012-2013: Multiwavelength Observations. , keywords =. doi:10.3847/0004-637X/816/2/53 , adsurl =

  46. [46]

    , keywords =

    Core shift effect in blazars. , keywords =. doi:10.1093/mnras/stx847 , archivePrefix =. 1704.03229 , primaryClass =

  47. [47]

    , keywords =

    Opacity, variability, and kinematics of AGN jets. , keywords =. doi:10.1093/mnras/stz885 , archivePrefix =. 1809.05536 , primaryClass =

  48. [48]

    , year = 1965, month = jan, volume =

    A 610.5 MHz survey of the sky between declinations +40 and +44 degrees. , year = 1965, month = jan, volume =. doi:10.1086/109815 , adsurl =

  49. [49]

    Astronomy Reports , keywords =

    Frequency-Dependent Core Shifts in Ultracompact Quasars. Astronomy Reports , keywords =. doi:10.1134/S1063772918110094 , archivePrefix =. 1809.10011 , primaryClass =

  50. [50]

    , keywords =

    High-precision Astrometric Millimeter Very Long Baseline Interferometry Using a New Method for Multi-frequency Calibration. , keywords =. doi:10.3847/1538-4357/834/2/177 , archivePrefix =. 1612.02958 , primaryClass =

  51. [51]

    Extragalactic Jets: Theory and Observation from Radio to Gamma Ray , year = 2008, editor =

    Phase-Referencing on BL Lac. Extragalactic Jets: Theory and Observation from Radio to Gamma Ray , year = 2008, editor =. doi:10.48550/arXiv.0711.3564 , archivePrefix =. 0711.3564 , primaryClass =

  52. [52]

    , keywords =

    The jet collimation profile at high resolution in BL Lacertae. , keywords =. doi:10.1051/0004-6361/202039616 , archivePrefix =. 2102.08952 , primaryClass =

  53. [53]

    , keywords =

    Relativistic jets as X-ray and gamma-ray sources. , keywords =. doi:10.1086/158638 , adsurl =

  54. [54]

    Jets in Active Galactic Nuclei

    Marscher, A.P. Jets in Active Galactic Nuclei. The Jet Paradigm: From Microquasars to Quasars. 2010. doi:10.1007/978-3-540-76937-8_7

  55. [55]

    , keywords =

    Toward Microarcsecond Astrometry for the Innermost Wobbling Jet of the BL Lacertae Object OJ 287. , keywords =. doi:10.3847/2041-8213/acf71c , archivePrefix =. 2309.03635 , primaryClass =

  56. [56]

    , year = 2022, month = sep, volume =

    Rapid quasi-periodic oscillations in the relativistic jet of BL Lacertae. , year = 2022, month = sep, volume =. doi:10.1038/s41586-022-05038-9 , adsurl =

  57. [57]

    The Radio Fundamental Catalog , publisher =

    Petrov, Leonid , keywords =. The Radio Fundamental Catalog , publisher =. doi:10.25966/DHRK-ZH08 , url =