Search

We study the collimation and acceleration of jets in the nearby giant radio galaxy NGC 315, using multifrequency Very Long Baseline Array observations and archival High Sensitivity Array and Very Large Array data. We find that the jet geometry transitions from a semi-parabolic shape into a conical/hyperbolic shape at a distance of ≈105 gravitational radii. We constrain the frequency-dependent position of the core, from which we locate the jet base. The jet collimation profile in the parabolic region is in good agreement with the steady axisymmetric force-free electrodynamic solution for the outermost poloidal magnetic field line anchored to the black hole event horizon on the equatorial plane, similar to the nearby radio galaxies M87 and NGC 6251. The velocity field derived from the asymmetry in brightness between the jet and counterjet shows gradual acceleration up to the bulk Lorentz factor of Γ ∼ 3 in the region where the jet collimation occurs, confirming the presence of the jet acceleration and collimation zone. These results suggest that the jets are collimated by the pressure of the surrounding medium and accelerated by converting Poynting flux to kinetic energy flux. We discover limb brightening of the jet in a limited distance range where the angular resolution of our data is sufficient to resolve the jet transverse structure. This indicates that either the jet has a stratified velocity field of fast-inner and slow-outer layers or the particle acceleration process is more efficient in the outer layer owing to the interaction with the surroundings on parsec scales. © 2021. The American Astronomical Society. All rights reserved., J.P. acknowledges financial support from the Korean National Research Foundation (NRF) via Global PhD Fellowship grant 2014H1A2A1018695 and support through the EACOA Fellowship awarded by the East Asia Core Observatories Association, which consists of the Academia Sinica Institute of Astronomy and Astrophysics, the National Astronomical Observatory of Japan, Center for Astronomical Mega-Science, Chinese Academy of Sciences, and the Korea Astronomy and Space Science Institute. This work is supported by the Ministry of Science and Technology of Taiwan grant MOST 109-2112-M-001-025 and 108-2112-M-001-051 (K.A). G.-Y.Z. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709). M.K. acknowledges financial support from JSPS KAKENHI grant JP18K03656 and JP18H03721. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated by Associated Universities, Inc. This work is based in part on observations made with the KaVA, which is operated by the the Korea Astronomy and Space Science Institute and the National Astronomical Observatory of Japan.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification., Comment: Accepted for publication in ApJ; This is a new version after updating Figure 11 in the manuscript

PKS 1510--089 is a bright and active $\gamma$-ray source that showed strong and complex $\gamma$-ray flares in mid-2015 during which the Major Atmospheric Gamma Imaging Cherenkov telescopes detected variable very high energy (VHE; photon energies $>$100 GeV) emission. We present long-term multi-frequency radio, optical, and $\gamma$-ray light curves of PKS 1510--089 from 2013 to 2018, and results of an analysis of the jet kinematics and linear polarization using 43 GHz Very Long Baseline Array data observed between late 2015 and mid-2017. We find that a strong radio flare trails the $\gamma$-ray flares in 2015, showing an optically thick spectrum at the beginning and becoming optically thin over time. Two laterally separated knots of emission are observed to emerge from the radio core nearly simultaneously during the $\gamma$-ray flares. We detect an edge-brightened linear polarization near the core in the active jet state in 2016, similar to the quiescent jet state in 2008--2013. These observations indicate that the $\gamma$-ray flares may originate from compression of the knots by a standing shock in the core and the jet might consist of multiple complex layers showing time-dependent behavior, rather than of a simple structure of a fast jet spine and a slow jet sheath., Comment: 11 pages, 7 figures, To appear in ApJ

We study the linear polarization of the radio cores of eight blazars simultaneously at 22, 43, and 86 GHz with observations obtained by the Korean VLBI Network (KVN) in three epochs between late 2016 and early 2017 in the frame of the Plasma-physics of Active Galactic Nuclei (PAGaN) project. We investigate the Faraday rotation measure (RM) of the cores; the RM is expected to increase with observing frequency if core positions depend on frequency due to synchrotron self-absorption. We find a systematic increase of RMs at higher observing frequencies in our targets. The RM--$\nu$ relations follow power-laws with indices distributed around 2, indicating conically expanding outflows serving as Faraday rotating media. Comparing our KVN data with contemporaneous optical polarization data from the Steward Observatory for a few sources, we find indication that the increase of RM with frequency saturates at frequencies of a few hundreds GHz. This suggests that blazar cores are physical structures rather than simple $\tau=1$ surfaces. A single region, e.g. a recollimation shock, might dominate the jet emission downstream of the jet launching region. We detect a sign change in the observed RMs of CTA 102 on a time scale of $\approx$1 month, which might be related to new superluminal components emerging from its core undergoing acceleration/deceleration and/or bending. We see indication for quasars having higher core RMs than BL Lac objects, which could be due to denser inflows/outflows in quasars., Comment: 22 pages, 10 figures, 3 tables; to appear in ApJ (submitted 2018 March 7; accepted 2018 May 11)

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources., 36 pages, 9 figures, published by AJ

We present the observational results of the Gamma-ray blazar, 3C 66A, at 2.3, 8.4, and 22 GHz at 4 epochs during 2004-05 with the VLBA. The resulting images show an overall core-jet structure extending roughly to the south with two intermediate breaks occurring in the region near the core. By model-fitting to the visibility data, the northmost component, which is also the brightest, is identified as the core according to its relatively flat spectrum and its compactness. As combined with some previous results to investigate the proper motions of the jet components, it is found the kinematics of 3C 66A is quite complicated with components of inward and outward, subluminal and superluminal motions all detected in the radio structure. The superluminal motions indicate strong Doppler boosting exists in the jet. The apparent inward motions of the innermost components last for at least 10 years and could not be caused by new-born components. The possible reason could be non-stationarity of the core due to opacity change., 24 pages, 4 figures

We study a warping instability of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on a circular orbit. Such a circumbinary disk is subject to not only tidal torques due to the binary gravitational potential but also radiative torques due to radiation emitted from an accretion disk around each black hole. We find that a circumbinary disk initially aligned with the binary orbital plane is unstable to radiation-driven warping beyond the marginally stable warping radius, which is sensitive to both the ratio of vertical to horizontal shear viscosities and the mass-to-energy conversion efficiency. As expected, the tidal torques give no contribution to the growth of warping modes but tend to align the circumbinary disk with the orbital plane. Since the tidal torques can suppress the warping modes in the inner part of circumbinary disk, the circumbinary disk starts to be warped at radii larger than the marginally stable warping radius. If the warping radius is of the order of 0.1 pc, a resultant semi-major axis is estimated to be of the order of 10^-2 pc to 10^-4 pc for 10^7 Msun black hole. We also discuss the possibility that the central objects of observed warped maser disks in active galactic nuclei are binary supermassive black holes with a triple disk: two accretion disks around the individual black holes and one circumbinary disk surrounding them., 23 pages, 4 figures, 2 tables, accepted for publication in ApJ

To explain the multi-wavelength light curves (from radio to X-ray) of HST-1 in the M87 jet, we propose an hourglass model that is a modified two-zone system of Tavecchio & Ghisellini (hereafter TG08): a slow hourglass-shaped or Laval nozzle-shaped layer connected by two revolving exponential surfaces surrounding a fast spine, through which plasma blobs flow. Based on the conservation of magnetic flux, the magnetic field changes along the axis of the hourglass. We adopt the result of TG08---the high-energy emission from GeV to TeV can be produced through inverse Compton by the two-zone system, and the photons from radio to X-ray are mainly radiated by the fast inner zone system. Here, we only discuss the light curves of the fast inner blob from radio to X-ray. When a compressible blob travels down the axis of the first bulb in the hourglass, because of magnetic flux conservation, its cross section experiences an adiabatic compression process, which results in particle acceleration and the brightening of HST-1. When the blob moves into the second bulb of the hourglass, because of magnetic flux conservation, the dimming of the knot occurs along with an adiabatic expansion of its cross section. A similar broken exponential function could fit the TeV peaks in M87, which may imply a correlation between the TeV flares of M87 and the light curves from radio to X-ray in HST-1. The Very Large Array (VLA) 22 GHz radio light curve of HST-1 verifies our prediction based on the model fit to the main peak of the VLA 15 GHz radio light curve., 14 pages, 2 figures, accepted for publication in AJ