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.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii. The data, showing predominantly negative RM sign without significant difference of the RMs on the northern and southern jet edges, suggest that the spatial extent of the Faraday screen is much larger than the jet. We apply models of hot accretion flows, thought to be prevalent in active galactic nuclei having relatively low luminosity such as M87, and find that the decrease of RM is described well by a gas density profile $\rho \propto r^{-1}$. This behavior matches the theoretically expected signature of substantial winds, nonrelativistic un-collimated gas outflows from hot accretion flows, which is consistent with the results of various numerical simulations. The pressure profile inferred from the density profile is flat enough to collimate the jet, which can result in gradual acceleration of the jet in a magneto-hydrodynamical process. This picture is in good agreement with the observed gradual collimation and acceleration of the M87 jet inside the Bondi radius. The dominance of negative RMs suggests that jet and wind axis are misaligned such that the jet emission exposes only one side of the toroidal magnetic fields permeating the winds.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world.