We report on high angular resolution polarimetric observations of the nearby radio galaxy M87 using the Very Long Baseline Array at 24 GHz ($\lambda=$1.3 cm) and 43 GHz ($\lambda=$7 mm) in 2017-2018. New images of the linear polarization substructure in the nuclear region are presented, characterized by a two-component pattern of polarized intensity and smooth rotation of the polarization plane around the 43 GHz core. From a comparison with an analogous dataset from 2007, we find that this global polarization pattern remains stable on a time interval of 11 yr, while showing smaller month-scale variability. We discuss the possible Faraday rotation toward the M87 nucleus at centimeter to millimeter wavelengths. These results can be interpreted in a scenario where the observed polarimetric pattern is associated with the magnetic structure in the confining magnetohydrodynamic wind, which also serves as the source of the observed Faraday rotation.

The supermassive black holes (SMBHs) of narrow-line Seyfert 1 galaxies (NLS1s) are at the lowest end of mass function of active galactic nuclei (AGNs) and preferentially reside in late-type host galaxies with pseudobulges, which are thought to be formed by internal secular evolution. On the other hand, the population of radio-loud NLS1s presents a challenge for the relativistic jet paradigm that powerful radio jets are exclusively associated with very high mass SMBHs in elliptical hosts, which are built-up through galaxy mergers. We investigated distorted radio structures associated with the nearest gamma-ray emitting, radio-loud NLS1 1H 0323+342. This provides supporting evidence for the merger hypothesis based on the past optical/near-infrared observations of its host galaxy. The anomalous radio morphology consists of two different structures, the inner curved structure of currently active jet and the outer linear structure of low-brightness relics. Such a coexistence might be indicative of the stage of an established black hole binary with precession before the black holes coalesce in the galaxy merger process. 1H 0323+342 and other radio-loud NLS1s under galaxy interactions may be extreme objects on the evolutionary path from radio-quiet NLS1s to normal Seyfert galaxies with larger SMBHs in classical bulges through mergers and merger-induced jet phases.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3\sigma$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region.

Observations of gamma-ray bursts (GRBs) at Very High Energy (VHE) offer a unique opportunity to investigate particle acceleration processes, magnetic fields and radiation fields in these events. Very Long Baseline Interferometry (VLBI) observations have been proven to be a powerful tool providing unique information on the source size of the GRBs at mas scales, as well as their accurate positions and possible expansion speeds. This paper reports on the follow-up observations of GRB 190114C, the first ever GRB detected with high significance at TeV photon energies by the MAGIC telescope, conducted with the East Asia VLBI Network (EAVN) at 22 GHz on three epochs, corresponding to 6, 15 and 32 days after the burst. The derived maps do not show any significant source above 5 sigma. The inferred upper limits on the GRB 190114C flux density at 22 GHz are used here to constrain the allowable two-dimensional parameter space for the afterglow emission. We find that our limits are consistent with most afterglow parameter combinations proposed so far in the literature. This is the first effort for the EAVN to search and monitor a radio transient in the Target of Opportunity mode. In addition to the useful constraints on GRB 190114C radio emission, experience gained from these observations is very helpful for future routine operation of EAVN transient program.

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 investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between $\sim300$ to $4 \times 10^7$ Schwarzschild Radii($R_{\rm S}$) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the VLBI Space Observatory Programme (VSOP), the Very Long Baseline Array (VLBA), and the Very Large Array (VLA). The jet-width profile of receding jets are apparently consistent with that of approaching jets throughout the measuring distance ranges, indicating symmetry at least up to the sphere of gravitational influence of the central black hole. The power-law index $a$ of the jet-width profile ($w_{\rm{jet}} \propto r^{a}$, where $w_{\rm jet}$ is the jet width, $r$ is the distance from the central engine in the unit of $R_{\rm S}$) apparently shows a transition from $a \sim 0$ to $a \sim 1$, i.e., the cylindrical-to-conical jet structures, at a distance of $\sim1\times10^{4} \ R_{\mathrm{S}}$. The cylindrical jet shape at the small distances is reminiscent of the innermost jets in 3C 84. Both the central engines of NGC 1052 and 3C 84 are surrounded by dense material, part of which is ionized and causes heavy free-free absorption.

We propose a new method of estimating a mass of a super massive black hole residing in the center of an active galaxy. The active galaxy M87 offers a convenient test case for the method due to the existence of a large amount of observational data on the jet and ambient environment properties in the central area of the object. We suggest that the observed transition of a jet boundary shape from a parabolic to a conical form is associated with the flow transiting from the magnetically dominated regime to the energy equipartition between plasma bulk motion and magnetic field. By coupling the unique set of observations available for the jet kinematics, environment and boundary profile with our MHD modelling under assumption on the presence of a dynamically important magnetic field in the M87 jet, we estimate the central black hole mass and spin. The method leads us to believe that the M87 super massive black hole has a mass somewhat larger than typically accepted so far.

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.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

The mechanisms producing fast variability of the $\gamma$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $\gamma$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $\gamma$-ray bands. In this paper we present the MAGIC VHE $\gamma$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $\gamma$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $\gamma$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $\gamma$-rays is very similar to two other observed VHE $\gamma$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $\gamma$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $\gamma$-ray flare, but all have some tension with the multiwavelength observations.

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.

The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet.

This paper theoretically investigates the relations between the structure of relativistic jets and produced synchrotron images, by using a steady axisymmetric force-free jet model. We especially focus on the limb-brightened jets that are largely symmetric to the jet axes and observed in some active galactic nuclei such as M87, Mrk~501, Cyg~A, and 3C84. We find that symmetrically limb-brightened images can be produced when magnetic field lines of the jet penetrate a fast-spinning black hole as motivated by the Blandford-Znajek mechanism. On the other hand, jets with magnetic field lines that pass through a slowly spinning black hole or the Keplerian accretion disk produce highly asymmetric radio images. In addition, the edge of a counter jet tends to be luminous in the accretion-disk model even for rather small viewing angles, which may be problematic for some observed jets. We also suggest that the site of particle accelerations in relativistic jets can be constrained by fitting the radio images to observations. This kind of study focusing on the jet images far away from the central engine is complimentary to those concentrating directly on the innermost region with upcoming data of the Event Horizon Telescope.

Radio jets in active galaxies have been expected to interact with circumnuclear environments in their early phase evolutions. By performing the multi-epoch monitoring observation with the KVN and VERA Array (KaVA) at 43~GHz, we investigate the kinematics of the notable newborn bright component C3 located at the tip of the recurrent jet of 3C~84. During 2015 August-September, we discover the flip of C3 and the amount of the flip is about 0.4~milli-arcsecond in angular scale, which corresponds to 0.14 parsec in physical scale. After the flip of C3, it wobbled at the same location for a few months and then it restarted to propagate towards the southern direction. The flux density of C3 coherently showed the monotonic increase during the observation period. The flip is in good agreement with hydrodynamical simulations of jets in clumpy ambient medium. We estimate the number density of the putative clump based on the momentum balance between the jet thrust and the ram pressure from the clump and it is about $10^{3-5}~{\rm cm^{-3}}$. We briefly discuss possible origins of the clump.

We present a kinematic study of the subparsec-scale radio jet of the radio galaxy 3C 84/NGC 1275 with the VLBI Exploration of Radio Astrometry (VERA) array at 22 GHz for 80 epochs from 2007 October to 2013 December. The averaged radial velocity of the bright component "C3" with reference to the radio core is found to be $0.27 pm 0.02c$ between 2007 October and 2013 December. This constant velocity of C3 is naturally explained by the advancing motion of the head of the mini-radio lobe. We also find a non-linear component in the motion of C3 with respect to the radio core. We briefly discuss possible origins of this non-linear motion.

We investigated the detailed radio structure of the jet of 1H 0323+342 using high-resolution multi-frequency Very Long Baseline Array observations. This source is known as the nearest $\gamma$-ray emitting radio-loud narrow-line Seyfert 1 (NLS1) galaxy. We discovered that the morphology of the inner jet is well characterized by a parabolic shape, indicating the jet being continuously collimated near the jet base. On the other hand, we found that the jet expands more rapidly at larger scales, resulting in a conical-like shape. The location of the "collimation break" is coincident with a bright quasi-stationary feature at 7 mas from core (corresponding to a deprojected distance of the order of $\sim$100pc), where the jet width locally contracts together with highly polarized signals, suggesting a recollimation shock. We found that the collimation region is coincident with the region where the jet speed gradually accelerates, suggesting the coexistence of the jet acceleration and collimation zone, ending up with the recollimation shock, which could be a potential site of high-energy $\gamma$-ray flares detected by the Fermi-LAT. Remarkably, these observational features of the 1H 0323+342 jet are overall very similar to those of the nearby radio galaxy M87 and HST-1 as well as some blazars, suggesting that a common jet formation mechanism might be at work. Based on the similarity of the jet profile between the two sources, we also briefly discuss the mass of the central black hole of 1H 0323+342, which is also still highly controversial on this source and NLS1s in general.

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.

Very Long Baseline Interferometry observations at 86 GHz reveal an almost hollow jet in M87 with a forked morphology. The detailed analysis presented here indicates that the spectral luminosity of the central spine of the jet in M87 is a few percent of that of the surrounding hollow jet $200 -400 \mu\rm{as}$ from the central black hole. Furthermore, recent jet models in indicate that a hollow "tubular" jet can explain a wide range of plausible broadband spectra originating from jetted plasma located within $\sim 30\mu\rm{as}$ of the central black hole, including the 230 GHz correlated flux detected by the Event Horizon Telescope. Most importantly, these hollow jets from the inner accretion flow have an intrinsic power capable of energizing the global jet out to kiloparsec scales. Thus motivated, this paper considers new models of the event horizon magnetosphere (EHM) in low luminosity accretion systems. Contrary to some models, the spine is not an invisible powerful jet. It is an intrinsically weak jet. In the new EHM solution, the accreted poloidal magnetic flux is weak and the background photon field is weak. It is shown how this accretion scenario naturally results in the dissipation of the accreted poloidal magnetic flux in the EHM not the accumulation of poloidal flux required for a powerful jet. The new solution indicates less large scale poloidal magnetic flux (and jet power) in the EHM than in the surrounding accretion flow and cannot support significant EHM driven jets.

We report the discovery of a local convergence of a jet cross section in the quasi-stationary jet feature in the gamma-ray-emitting narrow-line Seyfert 1 galaxy (NLS1) 1H 0323+342. The convergence site is located at approximately 7 mas (corresponding to the order of 100 pc in deprojection) from the central engine. We also found limb-brightened jet structures at both the upstream and downstream of the convergence site. We propose that the quasi-stationary feature showing the jet convergence and limb-brightening occurs as a consequence of recollimation shock in the relativistic jets. The quasi-stationary feature is one of the possible gamma-ray-emitting sites in this NLS1, in analogy with the HST-1 complex in the M87 jet. Monitoring observations have revealed that superluminal components passed through the convergence site and the peak intensity of the quasi-stationary feature, which showed apparent coincidences with the timing of observed gamma-ray activities.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

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.

We investigate the nature of the compact, and possibly variable nuclear radio source in the centre of WISE J0716-19, the proposed host galaxy of fast radio burst, FRB 150418. We observed WISE J0716-19 at 5.0 GHz with the European VLBI Network four times between 2016 March 16 and June 2. At three epochs, we simultaneously observed the source with e-MERLIN at the same frequency. We detected a compact source in the EVN data in each epoch with a significance up to ~8sigma. The four epochs yielded consistent results within their uncertainties, for both peak surface intensity and positions. The mean values for these quantities are I_peak = (115+-9) {\mu}Jy/beam and r.a. = 07:16:34.55496(7), dec. = -19:00:39.4754(8), respectively. The e-MERLIN data provided ~3-5sigma detections, at a position consistent with those of the EVN data. The presence of emission on angular scales intermediate between the EVN and e-MERLIN is consistent with being null. The brightness temperature of the EVN core is Tb~10^8.5K, close to the value required by Akiyama & Johnson (2016) to explain the radio properties of WISE J0716-19 in terms of interstellar induced short-term variability. Our observations provide direct, independent evidence of the existence of a nuclear compact source in WISE J0716-19, a physical scenario with no evident connection with FRB 150418. However, the EVN data do not show indication of the variability observed with the VLA.

We present full-polarization observations of the compact, steep-spectrum radio quasar 3C~286 made with the ALMA at 1.3~mm. These are the first full-polarization ALMA observations, which were obtained in the framework of Science Verification. A bright core and a south-west component are detected in the total intensity image, similar to previous centimeter images. Polarized emission is also detected toward both components. The fractional polarization of the core is about 17\%, this is higher than the fractional polarization at centimeter wavelengths, suggesting that the magnetic field is even more ordered in the millimeter radio core than it is further downstream in the jet. The observed polarization position angle (or EVPA) in the core is $\sim$\,$39^{\circ}$, which confirms the trend that the EVPA slowly increases from centimeter to millimeter wavelengths. With the aid of multi-frequency VLBI observations, we argue that this EVPA change is associated with the frequency-dependent core position. We also report a serendipitous detection of a sub-mJy source in the field of view, which is likely to be a submillimeter galaxy.

For short-wavelength VLBI observations, it is difficult to measure the phase of the visibility function accurately. The closure phases are reliable measurements under this situation, though it is not sufficient to retrieve all of the phase information. We propose a new method, Phase Retrieval from Closure Phase (PRECL). PRECL estimates all the visibility phases only from the closure phases. Combining PRECL with a sparse modeling method we have already proposed, imaging process of VLBI does not rely on dirty image nor self-calibration. The proposed method is tested numerically and the results are promising.

High-resolution millimeter wave very-long-baseline interferometry (mm-VLBI) is an ideal tool for probing the structure at the base of extragalactic jets in detail. The TeV blazar Mrk 501 is one of the best targets among BL Lac objects for studying the nature of off-axis jet structures because it shows different jet position angles at different scales. The aim of this study is to investigate the properties of the off-axis jet structure through high-resolution mm-VLBI images at the jet base and physical parameters such as kinematics, flux densities, and spectral indices. We performed Very Long Baseline Array (VLBA) observations over six epochs from 2012 February to 2013 February at 43 GHz. Quasi-simultaneous Global Millimeter VLBI Array (GMVA) observations at 86 GHz were performed in May 2012. We discover a new jet component at the northeast direction from the core in all the images at 43 and 86 GHz. The new component shows the off-axis location from the persistent jet extending to the southeast. The 43 GHz images reveal that the scattering of the positions of the NE component is within ~0.2 mas. The 86 GHz data reveals a jet component located 0.75 mas southeast of the radio core. We also discuss the spectral indices between 43 and 86 GHz, where the northeast component has steeper spectral index and the southeast component has comparable or flatter index than the radio core does.

We report on results from new high-sensitivity, high-resolution 86GHz (3.5 millimeter) observations of the jet base in the nearby radio galaxy M87, obtained by the Very Long Baseline Array in conjunction with the Green Bank Telescope. The resulting image has a dynamic range exceeding 1500 to 1, the highest ever achieved for this jet at this frequency, resolving and imaging a detailed jet formation/collimation structure down to ~10 Schwarzschild radii (Rs). The obtained 86GHz image clearly confirms some important jet features known at lower frequencies, i.e., a wide-opening angle jet base, a limb-brightened intensity profile, a parabola-shape collimation profile and a counter jet. The limb-brightened structure is already well developed at < 0.2mas (< 28Rs, projected) from the core, where the corresponding apparent opening angle becomes as wide as ~100 degrees. The subsequent jet collimation near the black hole evolves in a complicated manner; there is a "constricted" structure at tens Rs from the core, where the jet cross section is locally shrinking. We suggest that an external pressure support from the inner part of radiatively-inefficient accretion flow may be dynamically important in shaping/confining the footprint of the magnetized jet. We also present the first VLBI 86GHz polarimetric experiment for this source, where a highly polarized (~20%) feature is detected near the jet base, indicating the presence of a well-ordered magnetic field. As a by-product, we additionally report a 43/86 GHz polarimetric result for our calibrator 3C 273 suggesting an extreme rotation measure near the core.

We report multifrequency phase-referenced observations of the nearby radio galaxy NGC 4261, which has prominent two-sided jets, using the Very Long Baseline Array at 1.4-43 GHz. We measured radio core positions showing observing frequency dependences (known as "core shift") in both approaching jets and counter jets. The limit of the core position as the frequency approaches infinity, which suggests a jet base, is separated by 82$\pm$16 ${\mu}$as upstream in projection, corresponding to (310$\pm$60)Rs (Rs: Schwarzschild radius) as a deprojected distance, from the 43 GHz core in the approaching jet. In addition, the innermost component at the counter jet side appeared to approach the same position at infinity of the frequency, indicating that cores on both sides are approaching the same position, suggesting a spatial coincidence with the central engine. Applying a phase referencing technique, we also obtained spectral index maps, which indicate that emission from the counter jet is affected by free-free absorption (FFA). The result of the core shift profile on the counter jet also requires FFA because the core positions at 5-15GHz cannot be explained by a simple core shift model based on synchrotron self-absorption (SSA). Our result is apparently consistent with the SSA core shift with an additional disk-like absorber over the counterjet side. Core shift and opacity profiles at the counter jet side suggest a two-component accretion: a radiatively inefficient accretion flow at the inner region and a truncated thin disk in the outer region. We proposed a possible solution about density and temperature profiles in the outer disk on the basis of the radio observation.

We investigate the location of the radio jet bases ("radio cores") of blazars in radio images, and their stationarity by means of dense very long baseline interferometry (VLBI) observations. In order to measure the position of a radio core, we conducted 12 epoch astrometric observation of the blazar Markarian 421 with the VLBI Exploration of Radio Astrometry at 22 GHz immediately after a large X-ray flare, which occurred in the middle of 2011 September. For the first time, we find that the radio core is not stationary but rather changes its location toward 0.5 mas downstream. This angular scale corresponds to the de-projected length of a scale of $10^5$ Schwarzschild radii (Rs) at the distance of Markarian~421. This radio-core wandering may be a new type of manifestation associated with the phenomena of large X-ray flares.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature located 0.75 mas southeast of the radio core, which is consistent with the previous result. The location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show flat-to-steep, steep, and flat-to-invert, respectively.

The relativistic jet in M87 offers a unique opportunity for understanding the detailed jet structure and emission processes due to its proximity. In particular, the peculiar jet region HST-1 at ~1 arcsecond (or 80 pc, projected) from the nucleus has attracted a great deal of interest in the last decade because of its superluminal motion and broadband radio-to-X-ray outbursts, which may be further connected to the gamma-ray productions up to TeV energies. Over the last five years, we have been doing an intensive monitoring of HST-1 with EVN at 5GHz in order to examine the detailed structural evolution and its possible connection to high-energy activities. While this program already yielded interesting results in terms of the detailed mas-scale structure, proper motion measurements and structural variations, the recent HST-1 brightness is continuously decreasing at this frequency. To counter this, we have shifted our monitoring frequency to 1.7GHz from October 2013. This strategy successfully recovered the fainter emission that was missed in the last 5GHz session. Moreover, we again discovered the sudden emergence of a new component at the upstream edge of HST-1, demonstrating that the use of EVN 1.7GHz is indeed powerful to probe the current weak nature of HST-1. Here we report early results from the 1.7GHz monitoring as well as further progress on the long-term kinematic study.

The radio galaxy 3C84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGN) and one of the best laboratories to study the radio properties of subparsec scale jets. We discuss here the past and present activity of the nuclear region within the central 1pc and the properties of subparsec-sized components C1, C2 and C3. We compare these results with the high resolution space-VLBI image at 5GHz obtained with the RadioAstron satellite and we shortly discuss the possible correlation of radio emission with the gamma-ray emission.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at 120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi Large Area Telescope. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure in the radio band. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central black hole.

We explore the degree of magnetization at the jet base of M87 by using the observational data of the event horizon telescope (EHT) at 230~GHz obtained by Doeleman et al. By utilizing the method in Kino et al., we derive the energy densities of magnetic fields ($U_{B}$) and electrons and positrons ($U_{\pm}$) in the compact region detected by EHT (EHT-region) with its full-width-half-maximum size $40~{\rm \mu as}$. First, we assume that an optically-thick region for synchrotron self absorption (SSA) exists in the EHT-region. Then, we find that the SSA-thick region should not be too large not to overproduce the Poynting power at the EHT-region. The allowed ranges of the angular size and the magnetic field strength of the SSA-thick region are $21~{\rm \mu as} \le \theta_{\rm thick}\le 26.3~{\rm \mu as}$ and $50~{\rm G}\le B_{\rm tot}\le 124~{\rm G}$, respectively. Correspondingly $U_{B}\gg U_{\pm} $ is realized in this case. We further examine the composition of plasma and energy density of protons by utilizing the Faraday rotation measurement ($RM$) at 230~GHz obtained by Kuo et al. Then, we find that $U_{B}\gg U_{\pm}+U_{p} $ still holds in the SSA-thick region. Second, we examine the case when EHT-region is fully SSA-thin. Then we find that $U_{B}\gg U_{\pm}$ still holds unless protons are relativistic. Thus, we conclude that magnetically driven jet scenario in M87 is viable in terms of energetics close to ISCO scale unless the EHT-region is fully SSA-thin and relativistic protons dominated.

We investigate the position of the radio core in a blazar by multi-epoch astrometric observations at 43 GHz. Using the VLBI Exploration of Radio Astrometry (VERA), we have conducted four adjacent observations in February 2011 and another four in October 2011, and succeeded in measuring the position of the radio core in the TeV blazar Mrk 501 relative to a distant compact quasar NRAO 512. During our observations, we find that (1) there is no positional change within ~0.2 mas or ~2.0 pc de-projected with $\pm1\sigma$ error for the weighted-mean phase-referenced positions of Mrk 501 core relative to NRAO 512 over four adjacent days, and (2) there is an indication of position change for 3C 345 core relative to NRAO 512. By applying our results to the standard internal shock model for blazars, we constrain the bulk Lorenz factors of the ejecta.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure at radio. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central supermassive black hole.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We report our intensive radio monitoring observations of the jet in M87 with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi-LAT. During this period, an elevated level of the M87 flux is reported at VHE gamma rays. We detected a remarkable increase of the radio flux density from the unresolved jet base (radio core) with VERA at 22 and 43GHz coincident with the VHE activity. Meanwhile, we confirmed with EVN at 5GHz that HST-1 (an alternative gamma-ray production candidate site) remained quiescent in terms of its flux density and structure. These results in the radio bands strongly suggest that the VHE gamma-ray activity in 2012 originates in the jet base within 0.03pc or 56 Schwarzschild radii from the central supermassive black hole. We further conducted VERA astrometry for the M87 core during the flaring period, and detected core shifts between 22 and 43GHz. We also discovered a clear frequency-dependent evolution of the radio core flare at 43, 22 and 5GHz; the radio flux density increased more rapidly at higher frequencies with a larger amplitude, and the light curves clearly showed a time-lag between the peaks at 22 and 43GHz. This indicates that a new radio-emitting component was created near the black hole in the period of the VHE event, and then propagated outward with progressively decreasing synchrotron opacity. By combining these results, we estimated an apparent speed of the newborn component, and derived a sub-luminal speed of less than ~0.2c. This value is significantly slower than the super-luminal (~1.1c) features that appeared from the core during the prominent VHE flaring event in 2008, suggesting that the stronger VHE activity can be associated with the production of the higher Lorentz factor jet.

We explore energy densities of magnetic field and relativistic electrons in the M87 jet. Since the radio core at the jet base is identical to the optically thick surface against synchrotron self absorption (SSA), the observing frequency is identical to the SSA turnover frequency. As a first step, we assume the radio core as a simple uniform sphere geometry. Using the observed angular size of the radio core measured by the Very Long Baseline Array at 43 GHz, we estimate the energy densities of magnetic field ($U_{B}$) and relativistic electrons ($U_{e}$) based on the standard SSA formula. Imposing the condition that the Poynting power and relativistic electron one should be smaller than the total power of the jet, we find that (i) the allowed range of the magnetic field strength ($B_{tot}$) is from 1 G to 15 G, and that (ii) $1 times 10^{-5} < U_{e}/U_{B} < 6 times 10^{2}$ holds. The uncertainty of $U_{e}/U_{B}$ comes from the strong dependence on the angular size of the radio core and the minimum Lorentz factor of non-thermal electrons ($gamma_{e,min}$) in the core. It is still open that the resultant energetics is consistent with either the magnetohydrodynamic jet or with kinetic power dominated jet even on ~10 Schwarzschild radii scale.

We present a study of sub-pc scale radio structure of the radio galaxy 3C 84/NGC 1275 based on the Very Long Baseline Array (VLBA) data at 43 GHz. We discover a limb-brightening in the "restarted" jet associated with the 2005 radio outburst. In the 1990s, the jet structure was ridge-brightening rather than limb-brightening, despite the observations being done with similar angular resolution. This indicates that the transverse jet structure has changed recently. This change in the morphology shows an interesting agreement with the $\gamma$-ray flux increase, i.e., the $\gamma$-ray flux in 1990s was at least seven times lower than the current one. One plausible explanation for the limb-brightening is the velocity structure of the jet in the context of the stratified jet, which is a successful scenario to explain the $\gamma$-ray emission in some active galactic nuclei (AGNs). If this is the case, the change in apparent transverse structure might be caused by the change in the transverse velocity structure. We argue the possibility that the transition from ridge-brightening to limb-brightening is related to the $\gamma$-ray time variability on the timescale of decades. We also discuss the collimation profile of the jet.

To obtain a better understanding of the location and mechanisms for the production of the gamma-ray emission in jets of AGN we present a detailed study of the HST-1 structure, 0.8 arcsec downstream the jet of M87, previously identified as a possible candidate for TeV emission. HST-1 shows a very peculiar structure with superluminal as well as possible stationary sub-components, and appears to be located in the transition from a parabolic to a conical jet shape, presumably leading to the formation of a recollimation shock. This scenario is supported by our new RHD simulations in which the interaction of a moving component with a recollimation shock leads to the appearance of a new superluminal component. To discern whether HST-1 is produced by a recollimation shock or some other MHD instability, we present new polarimetric 2.2 and 5 GHz VLBA, as well as 15, 22 and 43 GHz JVLA observations obtained between November 2012 and March 2013.

Abridged: We investigated the detailed radio structure of the nucleus of the Sombrero galaxy using high-resolution, quasi-simultaneous, multi-frequency, phase-referencing VLBA observations. We obtained the VLBI images toward this nucleus, with unprecedented sensitivities and resolutions, at the seven frequencies between 1.4 and 43 GHz, where those at 15, 24 and 43 GHz are the first clear VLBI detections. At 43 GHz, the nuclear structure was imaged on a linear scale under 100 Schwarzschild radii. For the first time, we have discovered the presence of the extended structure in this nucleus, which is directing from the radio core in two sides toward northwest/southeast directions. The nuclear structure shows a clear spatial gradient on the radio spectra, which is similar to that commonly seen in more luminous AGN with powerful relativistic jets. Moreover, the radio core shows a frequency-dependent size with an elongated shape, and the position of the core also tends to be frequency dependent. A set of these new findings provide evidence that the central engine of the Sombrero is powering radio jets. Based on the observed brightness ratio of jet-to-counter jet, core position shift and its comparison with a theoretical model, we constrained the following fundamental physical parameters for the M 104 jets: (1) the northern side is the approaching jet, whereas the southern side is receding: (2) the inclination angle of the jet is relatively close to our line-of-sight, probably less than ~25degrees: (3) the jet intrinsic velocity is highly sub-relativistic at a speed less than ~0.2c. The derived pole-on nature of the M 104 jet is in accordance with the previous argument that M 104 contains a true type II AGN, i.e., the broad line region of this nucleus is actually absent or intrinsically weak, if the plane of the presumed circumnuclear torus is perpendicular to the axis of the radio jets.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We investigated the detailed inner jet structure of M87 using the Very Long Baseline Array data at 2, 5, 8.4, 15, 23.8, 43, and 86 GHz, especially focusing on the multi-frequency properties of the radio core at the jet base. First, we measured a size of the core region transverse to the jet axis, defined as $W_{\rm c}$, at each frequency \nu, and found a relation between $W_{\rm c}$ and $\nu$ as $W_{\rm c}(\nu) \propto \nu^{-0.71\pm0.05}$. Then, by combining $W_{\rm c}(\nu)$ and the frequency dependence of the core position $r_{\rm c}(\nu)$, which was obtained by our previous study, we have constructed a collimation profile of the innermost jet $W_{\rm c}(r)$ down to ~10 Schwarzschild radii ($R_{\rm s}$) from the central black hole. We found that $W_{\rm c}(r)$ smoothly connects with the width profile of the outer edge-brightened, parabolic jet, and then follows a similar radial dependence down to several tens of $R_{\rm s}$. Closer to the black hole, the measured radial profile suggests a possible change of the jet collimation shape from the outer parabolic one, in which the jet shape tends to become more radially-oriented. This could be related to a magnetic collimation process or/and interaction with surrounding materials at the jet base. The present results shed light on the importance of higher-sensitivity/resolution imaging studies for M87 at 86, 43 and also 22 GHz, and should be examined more rigorously.

The nearby low-luminosity active galactic nucleus (LLAGN) NGC 4258 has a weak radio continuum emission at the galactic center. Quasi-simultaneous multi-frequency observations using the Very Large Array (VLA) from 5 GHz (6 cm) to 22 GHz (1.3 cm) showed inverted spectra in all epochs, which were intra-month variable, as well as complicated spectral features that cannot be represented by a simple power law, indicating multiple blobs in nuclear jets. Using the Nobeyama Millimeter Array (NMA), we discovered a large amplitude variable emission at 100 GHz (3 mm), which had higher flux densities at most epochs than those of the VLA observations. A James Clerk Maxwell Telescope (JCMT) observation at 347 GHz (850 micron) served an upper limit of dust contamination. The inverted radio spectrum of the nucleus NGC 4258 is suggestive of an analogy to our Galactic center Sgr A*, but with three orders of magnitude higher radio luminosity. In addition to the LLAGN M 81, we discuss the nucleus of NGC 4258 as another up-scaled version of Sgr A*.

The nearby low-luminosity active galactic nucleus (LLAGN) NGC 4258 has a weak radio continuum component at the galactic center. We investigate its radio spectral properties on the basis of our new observations using the Nobeyama Millimeter Array at 100 GHz and archival data from the Very Large Array (VLA) at 1.7-43 GHz and the James Clerk Maxwell telescope at 347 GHz. The NGC 4258 nuclear component exhibits (1) an intra-month variable and complicated spectral feature at 5-22 GHz and (2) a slightly inverted spectrum at 5-100 GHz (a spectral index of ~0.3) in time-averaged flux densities, which are also apparent in the closest LLAGN M81. These similarities between NGC 4258 and M81 in radio spectral natures in addition to previously known core shift in their AU-scale jet structures produce evidence that the same mechanism drives their nuclei. We interpret the observed spectral property as the superposition of emission spectra originating at different locations with frequency-dependent opacity along the nuclear jet. Quantitative differences between NGC 4258 and M81 in terms of jet/counter jet ratio, radio loudness, and degree of core shift can be consistently understood by fairly relativistic speeds (bulk Lorentz factors of >~ 3) of jet and their quite different inclinations. The picture established from the two closest LLAGNs is useful for understanding the physical origin of unresolved and flat/inverted spectrum radio cores that are prevalently found in LLAGNs, including Sgr A*, with starved supermassive black holes in the present-day universe.

We investigate observational properties of a kilo-parsec scale knot in radio-loud quasar 3C 380 by using two epoch archival data obtained by Very Long Baseline Interferometry (VLBI) at 5 GHz on 1998 July and 2001 April. We succeed in obtaining the highest spatial resolution image of the bright knot K1 located at 732 milliarcseconds, or more than 20 kpc de-projected, downstream from the nucleus three times better than previously obtained highest resolution image by Papageorgiou et al. (2006). Our images reveal, with new clarity, "inverted bow-shock" structure in K1 facing the nucleus and its morphology resembles a conical shock wave. By comparing the two epoch images directly, we explore the kinematics of K1 and obtain the upper limit of apparent velocity, 0.25 mas/yr or 9.8 c of K1 for the first time. The upper limit of apparent velocity is marginally smaller than superluminal motions seen in the core region. Further new epoch VLBI observations are necessary to measure the proper motion at K1.

We report on the detailed radio status of the M87 jet during the Very-High-Energy (VHE) gamma-ray flaring event in April 2010, obtained from high-resolution, multi-frequency, phase-referencing VLBA observations. We especially focus on the properties for the jet base (the radio core) and the peculiar knot HST-1, which are currently favored as the gamma-ray emitting sites. During the VHE flaring event, the HST-1 region remains stable in terms of its structure and flux density in the optically thin regime above 2GHz, being consistent with no signs of enhanced activities reported at X-ray for this feature. The radio core shows an inverted spectrum at least up to 43GHz during this event. Astrometry of the core position, which is specified as ~20Rs from the central engine in our previous study, shows that the core position is stable on a level of 4Rs. The core at 43 and 22GHz tends to show slightly (~10%) higher flux level near the date of the VHE flux peak compared with the epochs before/after the event. The size of the 43-GHz core is estimated to be ~17Rs, which is close to the size of the emitting region suggested from the observed time scale of rapid variability at VHE. These results tend to favor the scenario that the VHE gamma-ray flare in 2010 April is associated with the radio core.

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

We introduce the GENJI program (Gamma-ray Emitting Notable AGN Monitoring by Japanese VLBI), which is a monitoring program of gamma-ray bright AGNs with the VERA array (VLBI Exploration of Radio Astrometry). The GENJI programme aims a dense monitoring at 22 GHz towards the $\gamma$-ray emitting active galactic nuclei (AGNs) to investigate the radio time variation of the core and possible ejection of new radio component, motion of jets, and their relation with the emission at other wavelengths especially in $\gamma$-rays. Currently we are monitoring 8 $\gamma$-ray-emitting notable AGNs (DA 55, 3C 84, M 87, PKS 1510-089, DA 406, NRAO 530, BL Lac, 3C 454.3) about once every two weeks. This programme is promising to trace the trend of radio time variation on shorter timescale than conventional VLBI monitoring programme and to provide complimentary data with them (e.g., MOJAVE, Boston University Blazar Project). In particular, we successfully coordinated quick follow-up observations after the GeV $\gamma$-ray flare in NRAO 530 and 3C 454.3 reported by the Fermi Gamma-ray Space Telescope. Here we present the initial results of morphology and light curves for the first 7-month operation.

We report on the VLBI follow-up observations using the Japanese VLBI Network (JVN) array at 22 GHz for the largest X-ray flare of TeV blazar Mrk 421 that occurred in mid-February, 2010. The total of five epochs of observations were performed at intervals of about 20 days between March 7 and May 31, 2010. No new-born component associated with the flare was seen directly in the total intensity images obtained by our multi-epoch VLBI observations. However, one jet component located at ~1 mas north-west from the core was able to be identified, and its proper motion can be measured as -1.66+/-0.46 mas yr^-1, which corresponds to an apparent velocity of -3.48+/-0.97 c. Here, this negative velocity indicates that the jet component was apparently moving toward the core. As the most plausible explanation, we discuss that the apparent negative velocity was possibly caused by the ejection of a new component, which could not be resolved with our observations. In this case, the obtained Doppler factor of the new component is around 10 to 20, which is consistent with the ones typically estimated by model fittings of spectral energy distribution for this source.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

Aims: We aim to constrain the structural variations within the HST-1 region downstream of the radio jet of M87, in general as well as in connection to the episodes of activity at very high energy (VHE). Methods: We analyzed and compared 26 VLBI observations of the M87 jet, obtained between 2006 and 2011 with the Very Long Baseline Array (VLBA) at 1.7 GHz and the European VLBI Network (EVN) at 5 GHz. Results: HST-1 is detected at all epochs; we model-fitted its complex structure with two or more components, the two outermost of which display a significant proper motion with a superluminal velocity around ~4c. The motion of a third feature that is detected upstream is more difficult to characterize. The overall position angle of HST-1 has changed during the time of our observations from -65deg to -90deg, while the structure has moved by over 80 mas downstream. Our results on the component evolution suggest that structural changes at the upstream edge of HST-1 can be related to the VHE events.

A Suzaku observation of a giant radio galaxy, 3C 326, which has a physical size of about 2 Mpc, was conducted on 2008 January 19 -- 21. In addition to several X-ray sources, diffuse emission was significantly detected associated with its west lobe, but the east lobe was contaminated by an unidentified X-ray source WARP J1552.4+2007. After careful evaluation of the X-ray and Non X-ray background, the 0.4 -- 7 keV X-ray spectrum of the west lobe is described by a power-law model. The photon index and 1 keV flux density was derived as $1.82_{-0.24}^{+0.26}\pm0.04$ and $19.4_{-3.2}^{+3.3}\pm 3.0$ nJy, respectively, where the first and second errors represent the statistical and systematic ones. The diffuse X-rays were attributed to be inverse Compton radiation by the synchrotron radio electrons scattering off the cosmic microwave background photons. This radio galaxy is the largest among those with lobes detected through inverse Compton X-ray emission. A comparison of the radio to X-ray fluxes yields the energy densities of electron and magnetic field as $u_e = (2.3 \pm 0.3 \pm 0.3) \times 10^{-13}$ ergs/cm3 and $u_m = (1.2_{-0.1}^{+0.2}\pm 0.2) \times 10^{-14}$ ergs/cm3, respectively. The galaxy is suggested to host a low luminosity nucleus with an absorption-corrected 2 -- 10 keV luminosity of $<2 \times 10^{42}$ ergs/s, together with a relatively weak radio core. The energetics in the west lobe of 3C 326 were compared with those of moderate radio galaxies with a size of $\sim 100$ kpc. The west lobe of 3C 326 is confirmed to agree with the correlations for the moderate radio galaxies, $u_e \propto D^{-2.2\pm0.4}$ and $u_m \propto D^{-2.4\pm0.4}$, where $D$ is their total physical size. This implies that the lobes of 3C 326 are still being energized by the jet, despite the current weakness of the nucleus.