A far-infrared counterpart to the west hot spot of the radio galaxy Pictor A is discovered with the Spectral and Photometric Imaging REceiver (SPIRE) onboard Herschel. The color-corrected flux density of the source is measured as $70.0 \pm 9.9$ mJy at the wavelength of 350 $\mu$m. A close investigation into its radio-to-optical spectrum indicates that the mid-infrared excess over the radio synchrotron component, detected with WISE and Spitzer, significantly contributes to the far-infrared band. Thanks to the SPIRE data, it is revealed that the spectrum of the excess is described by a broken power-law model subjected to a high-energy cutoff. By applying the radiative cooling break under continuous energy injection ($\Delta \alpha = 0.5$), the broken power-law model supports an idea that the excess originates in 10-pc scale substructures within the hot spot. From the break frequency, $\nu_{\rm b} = 1.6_{-1.0}^{+3.0} \times 10^{12}$ Hz, the magnetic field was estimated as $B\simeq1$-$4$ mG. This is higher than the minimum-energy magnetic field of the substructures by a factor of $3$--$10$. Even if the origin of the excess is larger than $\sim 100$ pc, the magnetic field stronger than the minimum-energy field is confirmed. It is proposed that regions with a magnetic field locally boosted via plasma turbulence are observed as the substructures. The derived energy index below the break, $\alpha \sim 0.22$ (conservatively $<0.42$), is difficult to be attributed to the strong-shock acceleration ($\alpha = 0.5$). Stochastic acceleration and magnetic reconnection are considered as a plausible alternative mechanism.

We characterize the accuracy of linear-polarization mosaics made using the Atacama Large Millimeter/submillimeter Array (ALMA). First, we observed the bright, highly linearly polarized blazar 3C 279 at Bands 3, 5, 6, and 7 (3 mm, 1.6 mm, 1.3 mm, and 0.87 mm, respectively). At each band, we measured the blazar's polarization on an 11$\times$11 grid of evenly-spaced offset pointings covering the full-width at half-maximum (FWHM) area of the primary beam. After applying calibration solutions derived from the on-axis pointing of 3C 279 to all of the on- and off-axis data, we find that the residual polarization errors across the primary beam are similar at all frequencies: the residual errors in linear polarization fraction $P_\textrm{frac}$ and polarization position angle $\chi$ are $\lesssim$0.001 ($\lesssim$0.1% of Stokes $I$) and $\lesssim$1$^\circ$ near the center of the primary beam; the errors increase to $\sim$0.003-0.005 ($\sim$0.3-0.5% of Stokes $I$) and $\sim$1-5$^\circ$ near the FWHM as a result of the asymmetric beam patterns in the (linearly polarized) $Q$ and $U$ maps. We see the expected double-lobed "beam squint" pattern in the circular polarization (Stokes $V$) maps. Second, to test the polarization accuracy in a typical ALMA project, we performed observations of continuum linear polarization toward the Kleinmann-Low nebula in Orion (Orion-KL) using several mosaic patterns at Bands 3 and 6. We show that after mosaicking, the residual off-axis errors decrease as a result of overlapping multiple pointings. Finally, we compare the ALMA mosaics with an archival 1.3 mm CARMA polarization mosaic of Orion-KL and find good consistency in the polarization patterns.

We report results on deep Jansky Very Large Array A-configuration observations at 22 GHz of the hotspots of the radio galaxies 3C227 and 3C445. Synchrotron emission in the optical on scales up to a few kpc was reported for the four hotspots. Our VLA observations point out the presence of unresolved regions with upper limit to their linear size of about 100 pc. This is the first time that such compact components in hotspots have been detected in a mini-sample, indicating that they are not a peculiar characteristic of a few individual hotspots. The polarization may reach values up to 70 per cent in compact (about 0.1 kpc scale) regions within the hotspot, indicating a highly ordered magnetic field with size up to a hundred parsecs. On larger scales, the average polarization of the hotspot component is about 30-45 per cent, suggesting the presence of a significant random field component, rather than an ordered magnetic field. This is further supported by the displacement between the peaks in polarized intensity and in total intensity images that is observed in all the four hotspots. The electric vector position angle is not constant, but changes arbitrarily in the central part of the hotspot regions, whereas it is usually perpendicular to the total intensity contours of the outermost edge of the hotspot structure, likely marking the large-scale shock front. The misalignment between X-ray and radio-to-optical emission suggests that the former is tracing the current particle acceleration, whereas the latter marks older shock fronts.

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 report the first detection of radio polarization of a GRB afterglow with the first intensive combined use of telescopes in the millimeter and submillimeter ranges for GRB171205A. The linear polarization degree in the millimeter band at the sub-percent level ($0.27 \pm 0.04\%$) is lower than those observed in late-time optical afterglows (weighted average of $\sim 1\%$). The Faraday depolarization by non-accelerated, cool electrons in the shocked region is one of possible mechanisms for the low value. In this scenario, larger total energy by a factor of $\sim 10$ than ordinary estimate without considering non-accelerated electrons is required. The polarization position angle varies by at least 20 degrees across the millimeter band, which is not inconsistent with this scenario. This result indicates that polarimetry in the millimeter and submillimeter ranges is a unique tool for investigating GRB energetics, and coincident observations with multiple frequencies or bands would provide more accurate measurements of the non-accelerated electron fraction.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

We present ALMA Band 6 observations of the CO(2-1), HCN(3-2), and HCO$^{+}$(3-2) lines in the nearby radio galaxy / brightest cluster galaxy (BCG) of NGC 1275 with the spatial resolution of $\sim20$ pc. In the previous observations, CO(2-1) emission was detected as radial filaments lying in the east-west direction. We resolved the inner filament and found that the filament cannot be represented by a simple infalling stream both morphologically and kinematically. The observed complex nature of the filament resembles the cold gas structure predicted by recent numerical simulations of cold chaotic accretion. A crude estimate suggests that the accretion rate of the cold gas can be higher than that of hot gas. Within the central 100 pc, we detected a rotational disk of the molecular gas whose mass is $\sim10^{8} M_{\sun}$. This is the first evidence of the presence of massive cold gas disk on this spatial scale for BCGs. The disk rotation axis is approximately consistent with the axis of the radio jet on subpc scales. This probably suggests that the cold gas disk is physically connected to the innermost accretion disk which is responsible for jet launching. We also detected absorption features in the HCN(3-2) and HCO$^{+}$(3-2) spectra against the radio continuum emission mostly radiated by $\sim1.2$-pc size jet. The absorption features are blue-shifted from the systemic velocity by $\sim$300-600~km~s$^{-1}$, which suggests the presence of outflowing gas from the active galactic nucleus (AGN). We discuss the relation of the AGN feeding with cold accretion, the origin of blue-shifted absorption, and estimate of black hole mass using the molecular gas dynamics.

Characterizing the variability of the extragalactic sources used for calibration in the Atacama Large Millimeter/Sub-millimeter Array (ALMA) is key to assess the flux scale uncertainty of science observations. To this end, we model the variability of 39 quasars which have been used by ALMA as secondary flux calibrators using continuous time stochastic processes. This formalism is specially adapted to the multi-frequency, quasi-periodic sampling which characterizes the calibration monitoring of ALMA. We find that simple mixtures of Ornstein-Uhlenbeck processes can describe well the flux and spectral index variability of these sources for Bands 3 and 7 (91.5 and 103.5, and 343.5 GHz, respectively). The spectral shape of the calibrators are characterized by negative spectral indices, mostly between $-0.35$ and $-0.80$, and with additional concavity. The model provides forecasts, interpolations, and uncertainty estimations for the observed fluxes that depend on the intrinsic variability of the source. These can be of practical use for the ALMA data calibrator survey and data quality assurance.

We present 870 $\mu$m ALMA polarization observation toward the Class I\hspace{-.1em}I protoplanetary disk around AS 209, which has concentric, multiple gaps and rings. We successfully detect the polarized emission and find that the polarization orientations and fractions have distinct characteristics between the inner and outer regions. In the inner region, the polarization orientations are parallel to the minor axis of the disk, which is consistent with the self-scattering model. The mean polarization fraction in the region is $\sim$0.2\%, which is lower than the expected value when the grains have the maximum polarization efficiency, which corresponds to $\lambda$/2$\pi$ $\sim$ 140 $\mu$m in grain radius. In the outer region, we detect $\sim$1.0\% polarization and find that the polarization orientations are almost in the azimuthal directions. Moreover, the polarization orientations have systematic angular deviations from the azimuthal directions with $\Delta$$\theta$ $\sim$ 4\fdg5 $\pm$ 1\fdg6. The pattern is consistent with a model that radially drifting dust grains are aligned by the gas flow against the dust grains. We consider possible scenarios of the grain dynamics at the AS 209 ring which can reproduce the polarization pattern. However, the directions of the observed angular deviations are opposite to what is predicted under the fact that the disk rotates clockwise. This poses a question in our understandings of the alignment processes and/or grain dynamics in protoplanetary disks.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary system containing the pulsar PSR B1259-63 orbiting around a Be star LS 2883 after the 2017 periastron passage. We detected radio continuum emission from the binary system in the millimeter/submillimeter wavelengths for the first time. At Band 3 (97 GHz), the flux 84 days after the periastron is almost the same as that 71 days after the periastron. Although the binary system showed intense GeV gamma-ray flares during our observations, the Band 3 flux did not indicate any time correlation with them. The Band 3 fluxes are consistent with an extrapolation of the radio spectrum at lower frequencies. Assuming that it is synchrotron emission, we constrain magnetic fields ($\lesssim 0.6$ G) and the high-energy cutoff of the electrons ($\gamma \gtrsim 360$). The flux at Band 7 (343 GHz) 69 days after the periastron shows a significant excess from the extrapolation of the radio spectrum at lower frequencies. The flux may be associated with the circumstellar disk around the Be star. We also present the results of Australian Telescope Compact Array (ATCA) observations at 94 GHz for the 2014 periastron passage, which show that the radio spectrum was relatively soft when the pulsar passed the disk.

We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model reported by Ho et al. (2019) and indicate that the spectral peaks ($\sim110$ GHz at the first epoch and $\sim67$ GHz at the second epoch) represent the synchrotron self-absorption frequency. The non-detection of linear polarization with $<$0.15% in the 230 GHz band at the phase when the effect of synchrotron self-absorption was quite small in the band may be explained by internal Faraday depolarization with high circumburst density and strong magnetic field. This result supports the stellar explosion scenario rather than the tidal disruption model. The maximum energy of accelerating particles at the shocks of AT2018cow-like objects is also discussed.

The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved structure and steep spectrum suggest a terminal shock in a radio lobe energized by the jet from the PMN J0948+0022 nucleus. This large-scale radio component is found at almost the same position angle as that of the pc-scale jet, indicating a stable jet direction up to ~1 Mpc. Its apparent one-sidedness implies a moderate advancing speed ($\beta>0.27$). The kinematic age of 1 x 10^7 year of the jet activity is consistent with the expected NLS1 phase of ~10^7 year in the AGN lifetime. The VLBA image reveals the jet structure at distances ranging from r=0.82 milliarcsec to 3.5 milliarcsec, corresponding to approximately 100 pc, where superluminal motions were found. The jet width profile ($\propto r^{1.12}$) and flux-density profile ($\propto r^{-1.44}$) depending on the distance $r$ along the jet suggest that the jet kinetic energy is converted into internal energy in this region. The jet is causally connected in a nearly conical streamline, which is supported by ambient pressure at 100 pc scales in the host galaxy of PMN J0948+0022.

NGC 1275 is one of the most conspicuous active galactic nuclei (AGN) in the local Universe. The radio jet currently emits a flux density of $\sim 10$ Jy at $\sim 1$ mm wavelengths, down from the historic high of $\sim 65$ Jy in 1980. Yet, the nature of the AGN in NGC 1275 is still controversial. It has been debated whether this is a broad emission line (BEL) Seyfert galaxy, an obscured Seyfert galaxy, a narrow line radio galaxy or a BL-Lac object. We clearly demonstrate a persistent H$\beta$ BEL over the last 35 years with a full width half maximum (FWHM) of 4150 - 6000 km/s. We also find a prominent P$\alpha$ BEL (FWHM $\approx 4770 $ km/s) and a weak CIV BEL (FWHM $\approx 4000 $ km/s), H$\beta$/CIV $\approx 2$. A far UV HST observation during suppressed jet activity reveals a low luminosity continuum. The H$\beta$ BEL luminosity is typical of broad line Seyfert galaxies with similar far UV luminosity. X-ray observations indicate a softer ionizing continuum than expected for a broad line Seyfert galaxy with similar far UV luminosity. This is opposite of the expectation of advection dominated accretion. The AGN continuum appears to be thermal emission from a low luminosity, optically thick, accretion flow with a low Eddington ratio, $\sim 0.0001$. The soft, weak ionizing continuum is consistent with the relatively weak CIV BEL. Evidence that the BEL luminosity is correlated with the jet mm wave luminosity is presented. Apparently, the accretion rate regulates jet power.

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.

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.

The origin of polarized emission from protoplanetary disks is uncertain. Three mechanisms have been proposed for such polarized emission: grain alignment with magnetic fields, grain alignment with radiation gradients, and self-scattering of thermal dust emission. Aiming to observationally identify the polarization mechanisms, we present ALMA polarization observations of the 0.87 mm dust continuum emission toward the circumstellar disk around HD 142527 with high spatial resolution. We confirm that the polarization vectors in the northern region are consistent with self-scattering. Furthermore, we show that the polarization vectors in the southern region are consistent with grain alignment by magnetic fields, although self-scattering cannot be ruled out. To understand the differences between the polarization mechanisms, we propose a simple grain size segregation model: small dust grains ($\lesssim$ 100 microns) are dominant and aligned with magnetic fields in the southern region, and middle-sized ($\sim100$ microns) grains in the upper layer emit self-scattered polarized emission in the northern region. The grain size near the middle plane in the northern region cannot be measured because the emission at 0.87 mm is optically thick. However, it can be speculated that larger dust grains ($\gtrsim$ cm) may accumulate near this plane. These results are consistent with those of a previous analysis of the disk, in which large grain accumulation and optically thick emission from the northern region were found. This model is also consistent with theories where smaller dust grains are aligned with magnetic fields. The magnetic fields are toroidal, at least in the southern region.

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.

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.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

We present the large statistics of the galaxy effective radius in the rest-frame far-infrared (FIR) wavelength Re(FIR) obtained from 1627 deep ALMA 1-mm band maps that become public by 2017 July. Our ALMA sample consists of 1034 sources at z=0-6 that typically have the star-formation rate of ~100-1000 M$_{\odot}$/yr and the stellar mass of ~10$^{10}$-10$^{11.5}$ M$_{\odot}$. We homogeneously derive Re(FIR) and FIR luminosity L_FIR of our ALMA sources via the same uv-visibility method with the exponential disk fitting, carefully evaluating the selection and measurement incompletenesses by realistic Monte-Carlo simulations. We find that there is a positive correlation between Re(FIR) and L_FIR at the >99% significance level. The best-fit power-law function, Re(FIR) $\propto$ L_FIR$^{\alpha}$, provides $\alpha=0.28+/-0.07$, and shows that the average Re(FIR) at a fixed L_FIR decreases toward high redshifts. The best-fit $\alpha$ and the redshift evolution of Re(FIR) are similar to those of the galaxy effective radius in the rest-frame UV (optical) wavelength Re(UV) (Re(Opt.)) revealed by Hubble Space Telescope (HST) studies. We identify that our ALMA sources have significant trends of Re(FIR) <~ Re(UV) and Re(Opt.), suggesting that dusty starbursts take place in a compact region. Moreover, Re(FIR) of our ALMA sources is comparable to Re(Opt.) of quiescent galaxies at z~1-3 as a function of stellar mass, supporting the evolutionary connection between these two galaxy populations. We also investigate rest-frame UV and optical morphologies of our ALMA sources with deep HST images, and find that ~30-40% of our ALMA sources are classified as major mergers. This indicates that dusty starbursts are triggered not only by the major mergers but also the other mechanism(s).

We present Very Long Baseline Array polarimetric observations of the innermost jet of 3C$\sim$84 (NGC$\sim$1275) at 43$\sim$GHz. A significant polarized emission is detected at the hotspot of the innermost re-started jet, which is located $\sim$1 pc south from the radio core. While the previous report presented a hotspot at the southern end of the western limb, the hotspot location has been moved to the southern end of the eastern limb. Faraday rotation is detected within an entire bandwidth of the 43-GHz band. The measured rotation measure (RM) is at most (6.3$\pm$1.9)$\times10^{5}$$\sim$rad$\sim$m$^{-2}$ and might be slightly time variable on the timescale of a month by a factor of a few. Our measured RM and the RM previously reported by the CARMA and SMA observations cannot be consistently explained by the spherical accretion flow with a power-law profile. We propose that a clumpy/inhomogeneous ambient medium is responsible for the observed rotation measure. Using equipartition magnetic field, we derive the electron density of $2\times10^{4}$$\sim$cm$^{-3}$. Such an electron density is consistent with the cloud of narrow line emission region around the central engine. We also discuss the magnetic field configuration from black hole scale to pc scale and the origin of low polarization.

The mechanisms causing millimeter-wave polarization in protoplanetary disks are under debate. To disentangle the polarization mechanisms, we observe the protoplanetary disk around HL Tau at 3.1 mm with the Atacama Large Millimeter/submillimeter Array (ALMA), which had polarization detected with CARMA at 1.3 mm. We successfully detect the ring-like azimuthal polarized emission at 3.1 mm. This indicates that dust grains are aligned with the major axis being in the azimuthal direction, which is consistent with the theory of radiative alignment of elongated dust grains, where the major axis of dust grains is perpendicular to the radiation flux. Furthermore, the morphology of the polarization vectors at 3.1 mm is completely different from those at 1.3 mm. We interpret that the polarization at 3.1 mm to be dominated by the grain alignment with the radiative flux producing azimuthal polarization vectors, while the self-scattering dominates at 1.3 mm and produces the polarization vectors parallel to the minor axis of the disk. By modeling the total polarization fraction with a single grain population model, the maximum grain size is constrained to be $100{\rm~\mu m}$, which is smaller than the previous predictions based on the spectral index between ALMA at 3 mm and VLA at 7 mm.

We explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching and we found that a fossil shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the followings: (i) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm$^{-3}$ radiates bright Inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (ii) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of $10^{3}$ cm$^{-3}$ and the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. A possible production of ultra high energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

We present Atacama Large Millimeter Array (ALMA) polarization observations at 97.5 GHz of the southern hot spot of the radio galaxy 3C 445. The hot spot structure is dominated by two bright components enshrouded by diffuse emission. Both components show fractional polarization between 30 and 40 per cent, suggesting the presence of shocks. The polarized emission of the western component has a displacement of about 0.5 kpc outward with respect to the total intensity emission, and may trace the surface of a front shock. Strong polarization is observed in a thin strip marking the ridge of the hot spot structure visible from radio to optical. No significant polarization is detected in the diffuse emission between the main components, suggesting a highly disordered magnetic field likely produced by turbulence and instabilities in the downstream region that may be at the origin of the extended optical emission observed in this hot spot. The polarization properties support a scenario in which a combination of both multiple and intermittent shock fronts due to jet dithering, and spatially distributed stochastic second-order Fermi acceleration processes are present in the hot spot complex.

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 report the detection of a new feature at the centre of NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. This feature emerges ~ 2 mas (~ 0.8 pc) north of the central core in recent 15 and 43 GHz VLBA images, and seems to be the counterjet to a known radio jet expanding to the south of the core. Apparently, the two jets were born through an outburst around 2005. From the ratio of the apparent lengths of the two jets from the core, we found that the jet angle to the line of sight is \theta=65^\circ\pm 16^\circ, which is not much different from the angle of the outer jets generated by an activity around 1959 and constrains theories on gamma-ray emission from jets. The new northern jet has a strongly inverted spectrum in contrast with the southern jet. This suggests that the central black hole is surrounded by a subparsec-scale accretion disk with the density of >~ 10^5 cm^-3. The brightness of the counterjet suggests that the disk is highly inhomogeneous. The ambient gas density in the direction of the jet is ~ 8 cm^-3 if the current jet activity is similar to the past average.

We present the polarization observations toward the circumstellar disk around HD 142527 by using Atacama Large Millimeter/submillimeter Array (ALMA) at the frequency of 343 GHz. The beam size is $0.51 " \times 0.44 "$, which corresponds to the spatial resolution of $\sim$ 71 $\times$ 62 AU. The polarized intensity displays a ring-like structure with a peak located on the east side with a polarization fraction of $P= 3.26 \pm 0.02$ %, which is different from the peak of the continuum emission from the northeast region. The polarized intensity is significantly weaker at the peak of the continuum where $P= 0.220 \pm 0.010$ %. The polarization vectors are in the radial direction in the main ring of the polarized intensity, while there are two regions outside at the northwest and northeast areas where the vectors are in the azimuthal direction. If the polarization vectors represent the magnetic field morphology, the polarization vectors indicate the toroidal magnetic field configuration on the main ring and the poloidal fields outside. On the other hand, the flip of the polarization vectors is predicted by the self-scattering of thermal dust emission due to the change of the direction of thermal radiation flux. Therefore, we conclude that self-scattering of thermal dust emission plays a major role in producing polarization at millimeter wavelengths in this protoplanetary disk. Also, this puts a constraint on the maximum grain size to be approximately 150 ${\rm \mu m}$ if we assume compact spherical dust grains.

We analyzed Suzaku/XIS data of 2006--2015 observations of a gamma-ray emitting radio galaxy NGC 1275, and brightening of the nucleus in the X-ray band was found in 2013--2015, correlating with GeV Gamma-ray brightening. This is the first evidence of variability with correlation between GeV gamma-ray and X-ray for NGC 1275. We also analyzed Swift/XRT data of NGC 1275, and found that X-ray was flaring by a factor of $\sim$5 in several days in 2006, 2010, and 2013. The X-ray spectrum during the flare was featureless and somewhat steeper with a photon index of $\sim$2 against $\sim$1.7 in the normal state, indicating that a synchrotron component became brighter. A large Xray to GeV gamma-ray flux ratio in the flare could be explained by the shock-in-jet scenario. On the other hand, a long-term gradual brightening of radio, X-ray, and GeV gamma-ray with a larger gamma-ray amplitude could be origin of other than internal shocks, and then we discuss some possibilities.

We report on the 1.6 GHz (18 cm) VLBI observations of the unresolved, steady TeV source HESS J1943+213 located in the Galactic plane, performed with the European VLBI Network (EVN) in 2014. Our new observations with a nearly full EVN array provide the deepest image of HESS J1943+213 at the highest resolution ever achieved, enabling us to resolve the long-standing issues of the source identification. The milliarcsecond-scale structure of HESS J1943+213 has a clear asymmetric morphology, consisting of a compact core and a diffuse jet-like tail. This is broadly consistent with the previous e-EVN observations of the source performed in 2011, and re-analyzed in this work. The core component is characterized by the brightness temperature of $\gtrsim1.8 \times 10^9$ K, which is typical for low-luminosity blazars in general. Overall, radio properties of HESS J1943+213 are consistent with the source classification as an "extreme high-frequency-peaked BL Lac object". Remarkably, we note that since HESS J1943+213 does not reveal any optical or infrared signatures of the AGN activity, it would never be recognized and identified as a BL Lac object, if not its location close to the Galactic plane where the High Energy Stereoscopic System has surveyed, and the follow-up dedicated X-ray and radio studies triggered by the source detection in the TeV range. Our results suggest therefore a presence of an unrecognized, possibly very numerous population of particularly extreme HBLs, and simultaneously demonstrate that the low-frequency VLBI observations with high-angular resolution are indispensable for a proper identification of such objects.

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.

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 present statistics of 133 faint 1.2-mm continuum sources detected in about 120 deep ALMA pointing data that include all the archival deep data available by 2015 June. We derive number counts of 1.2 mm continuum sources down to 0.02 mJy partly with the assistance of gravitational lensing, and find that the total integrated 1.2 mm flux of the securely identified sources is 22.9^(+6.7)_(-5.6) Jy deg^(-2) that corresponds to 104^(+31)_(-25) % of the extragalactic background light (EBL) measured by COBE observations. These results suggest that the major 1.2 mm EBL contributors are sources with 0.02 mJy, and that very faint 1.2 mm sources with ~< 0.02 mJy contribute negligibly to the EBL with the possible flattening and/or truncation of number counts in this very faint flux regime. To understand the physical origin of our faint ALMA sources, we measure the galaxy bias bg by the counts-in-cells technique, and place a stringent upper limit of bg < 3.5 that is not similar to bg values of massive DRGs and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near-infrared (NIR) deep fields, we identify optical-NIR counterparts for 59% of our faint ALMA sources, majority of which have luminosities, colors, and the IRX-beta relation same as sBzKs and LBGs. We thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as sBzKs and LBGs in optical studies, and that these faint ALMA sources are not miniature (U)LIRGs simply scaled down with the infrared brightness.

We investigate shell emission associated with dying radio loud AGNs. First, based on our recent work by Ito et al. (2015), we describe the dynamical and spectral evolutions of shells after stopping the jet energy injection. We find that the shell emission overwhelms that of the radio lobes soon after stopping the jet energy injection because fresh electrons are continuously supplied into the shell via the forward shock while the radio lobes rapidly fade out without jet energy injection. We find that such fossil shells can be a new class of target sources for SKA telescope. Next, we apply the model to the nearby radio source 3C84. Then, we find that the fossil shell emission in 3C84 is less luminous in radio band while it is bright in TeV gamma-ray band and it can be detectable by CTA.

We report on Atacama Large Millimeter/submillimeter Array (ALMA) detections of molecular absorption lines in Bands 3, 6 and 7 toward four radio-loud quasars, which were observed as the bandpass and complex gain calibrators. The absorption systems, three of which are newly detected, are found to be Galactic origin. Moreover, HCO absorption lines toward two objects are detected, which almost doubles the number of HCO absorption samples in the Galactic diffuse medium. In addition, high HCO to H13CO+ column density ratios are found, suggesting that the interstellar media (ISM) observed toward the two calibrators are in photodissociation regions, which observationally illustrates the chemistry of diffuse ISM driven by ultraviolet (UV) radiation. These results demonstrate that calibrators in the ALMA Archive are potential sources for the quest for new absorption systems and for detailed investigation of the nature of the ISM.

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $\mu$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

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.

We report on VLBI observations of NGC 1275 with VLBI Exploration of Radio Astrometry (VERA) and Very Long Baseline Array (VLBA) during 2008 and 2013 at 15, 22 and 43 GHz. Our observations provide long-term variations in the radio flux and structure of the inner jet in NGC 1275 on sub-parsec scales, particularly at C3 which is a new radio component accountable for the recent active state in radio regime since early-2000s. We found the apparent velocity of C3 was $\beta_{\rm{app}} = 0.267 \pm 0.007 $, which was consistent with typical velocities of the hot spot and/or radio lobes in the young radio sources. Furthermore, the radio flux and size of C3 had simultaneously increased during September 2007 and May 2012 with keeping optically-thin spectra, suggesting that the particle acceleration occurred inside C3 as generally seen in the hot spots of the radio lobes. Our new observations suggest that C3 is a new-born hot spot and/or radio lobe at a very early stage of its evolution. Our results imply the radio lobes in the radio galaxies might be already formed on sub-parsec scales from the central SMBH and might play an important role to produce $\rm \gamma$-ray emission in the young radio sources.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

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 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 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.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

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 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.

The velocity field of the M87 jet from milli-arcsecond (mas) to arcsecond scales is extensively investigated together with new radio images taken by EVN observations. We detected proper motions of components located at between 160 mas from the core and the HST-1 complex for the first time. Newly derived velocity fields exhibits a systematic increase from sub-to-superluminal speed in the upstream of HST-1. If we assume that the observed velocities reflect the bulk flow, we here suggest that the M87 jet may be gradually accelerated through a distance of 10^6 times of the Schwarzschild radius of the supermassive black hole. The acceleration zone is co-spatial with the jet parabolic region, which is interpreted as the collimation zone of the jet (Asada & Nakamura 2012). The acceleration and collimation take place simultaneously, which we suggest a characteristic of magnetohydrodynamic flows. Distribution of the velocity field has a peak at HST-1, which is considered as the site of over-collimation, and shows a deceleration downstream of HST-1 where the jet is conical. Our interpretation of the velocity map in the M87 jet gives a hypothesis in AGNs that the acceleration and collimation zone of relativistic jets extends over the whole scale within the sphere of influence of the supermassive black hole.

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.

We conducted the first multi-frequency polarimetric imaging of four broad absorption line (BAL) quasars using Very Long Baseline Array at milli-arcsecond resolutions to investigate the inclination of the non-thermal jet and test the hypothesis that radio sources in BAL quasars are still young. Among these four sources, J0928+446, J1018+0530, and J1405+4056 show one-sided structures in parsec scales, and polarized emission detected in the core. These characteristics are consistent with those of blazars. We set constraints on viewing angles to $<$66 deg for these jets, in the framework of a Doppler beaming effect. J1159+0112 exhibits an unpolarized gigahertz peaked spectrum component and several discrete blobs with steep spectra on both sides of the central component across $\sim$1 kpc. These properties are consistent with those of young radio sources. We discuss the structures of jets and AGN wind.

We conducted Very Long Baseline Array (VLBA) observations of seven nearby narrow-line Seyfert 1 (NLS1) galaxies at 1.7 GHz (18cm) with milli-arcsecond resolution. This is the first systematic very long baseline interferometry (VLBI) study focusing on the central parsec-scale regions of radio-quiet NLS1s. Five of the seven were detected at a brightness temperature of >~5x10^6 K and contain radio cores with high brightness temperatures of >6x10^7 K, indicating a nonthermal process driven by jet-producing central engines as is observed in radio-loud NLS1s and other active galactic nucleus (AGN) classes. VLBA images of MRK 1239, MRK 705, and MRK 766 exhibit parsec-scale jets with clear linear structures. A large portion of the radio power comes from diffuse emission components that are distributed within the nuclear regions (<~300 pc), which is a common characteristic throughout the observed NLS1s. Jet kinetic powers limited by the Eddington limit may be insufficient to allow the jets escape to kiloparsec scales for these radio-quiet NLS1s with low-mass black holes of <~10^7 M_sun.

This paper reports very-long-baseline interferometry observations of the radio-loud broad absorption line (BAL) quasar J1020+4320 at 1.7, 2.3, 6.7, and 8.4 GHz using the Japanese VLBI network (JVN) and European VLBI network (EVN). The radio morphology is compact with a size of ~10 pc. The convex radio spectrum is stable over the last decade; an observed peak frequency of 3.2 GHz is equivalent to 9.5 GHz in the rest frame, suggesting an age of the order of ~100 years as a radio source, according to an observed correlation between linear size and peak frequency of compact steep spectrum (CSS) and giga-hertz peaked spectrum (GPS) radio sources. A low-frequency radio excess suggests relic of past jet activity. J1020+4320 may be one of the quasars with recurrent and short-lived jet activity during a BAL-outflowing phase.

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.

The properties of tidally induced arms provide a means to study molecular cloud formation and the subsequent star formation under environmental conditions which in principle are different from quasi stationary spiral arms. We report the properties of a newly discovered molecular gas arm of likely tidal origin at the south of NGC 4039 and the overlap region in the Antennae galaxies, with a resolution of 1"68 x 0"85, using the Atacama Large Millimeter/submillimeter Array science verification CO(2-1) data. The arm extends 3.4 kpc (34") and is characterized by widths of ~ 200 pc (2") and velocity widths of typically \DeltaV ~ 10-20 km/s . About 10 clumps are strung out along this structure, most of them unresolved, with average surface densities of \Sigma_gas ~ 10-100 Msun pc^{-2}, and masses of (1-8) x 10^6 Msun. These structures resemble the morphology of beads on a string, with an almost equidistant separation between the beads of about 350 pc, which may represent a characteristic separation scale for giant molecular associations. We find that the star formation efficiency at a resolution of 6" (600 pc) is in general a factor of 10 higher than in disk galaxies and other tidal arms and bridges. This arm is linked, based on the distribution and kinematics, to the base of the western spiral arm of NGC 4039, but its morphology is different to that predicted by high-resolution simulations of the Antennae galaxies.

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 the finding of kiloparsec (kpc)-scale radio structures in three radio-loud narrow-line Seyfert 1 (NLS1) galaxies from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) of the Very Large Array (VLA), which increases the number of known radio-loud NLS1s with kpc-scale structures to six, including two gamma-ray emitting NLS1s (PMN J0948+0022 and 1H 0323+342) detected by the Fermi Gamma-ray Space Telescope. The detection rate of extended radio emissions in NLS1s is lower than that in broad-line active galactic nuclei (AGNs) with a statistical significance. We found both core-dominated (blazar-like) and lobe-dominated (radio-galaxy-like) radio structures in these six NLS1s, which can be understood in the framework of the unified scheme of radio-loud AGNs that considers radio galaxies as non-beamed parent populations of blazars. Five of the six NLS1s have (i) extended radio luminosities suggesting jet kinetic powers of >~10^44 erg/s, which is sufficient to make jets escape from hosts' dense environments, (ii) black holes of >~10^7 solar mass, which can generate the necessary jet powers from near-Eddington mass accretion, and (iii) two-sided radio structures at kpc scales, requiring expansion rates of ~0.01c--0.3c and kinematic ages of >~10^7 years. On the other hand, most typical NLS1s would be driven by black holes of <~10^7 solar mass in a limited lifetime of ~10^7 years. Hence the kpc-scale radio structures may originate in a small window of opportunity during the final stage of the NLS1 phase just before growing into broad-line AGNs.

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.

The radio galaxy 3C 84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGNs) and one of the best laboratories to study the radio properties of the sub-pc jet in connection with the gamma-ray emission. In order to identify possible radio counterparts of the gamma-ray emissions in 3C 84, we study the change in structure within the central 1 pc and the light curve of sub-pc-size components C1, C2, and C3. We search for any correlation between changes in the radio components and the gamma-ray flares by making use of VLBI and single dish data. Throughout the radio monitoring spanning over two GeV gamma-ray flares detected by the {\it Fermi}-LAT and the MAGIC Cherenkov Telescope in the periods of 2009 April to May and 2010 June to August, total flux density in radio band increases on average. This flux increase mostly originates in C3. Although the gamma-ray flares span on the timescale of days to weeks, no clear correlation with the radio light curve on this timescale is found. Any new prominent components and change in morphology associated with the gamma-ray flares are not found on the VLBI images.

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.

Following the discovery of a new radio component right before the GeV \gamma-ray detection since 2008 August by Fermi Gamma-ray Space Telescope, we present a detailed study of the kinematics and lightcurve on the central sub-pc scale of 3C 84 using the archival VLBA 43-GHz data covering the period between 2002 January to 2008 November. We find that the new component "C3", previously reported by the observations with the VLBI Exploration of Radio Astrometry (VERA), was already formed in 2003. The flux density of C3 increases moderately until 2008, and then it becomes brighter rapidly after 2008. The radio core, C1, also shows a similar trend. The apparent speed of C3 with reference to the core C1 shows moderate acceleration from 0.10c to 0.47c between 2003 November to 2008 November, but is still sub-relativistic. We further try to fit the observed broadband spectrum by the one-zone synchrotron self-Compton (SSC) model using the measured apparent speed of C3. The fit can reproduce the observed \gamma-ray emission, but does not agree with the observed radio spectral index between 22 and 43 GHz.

We conducted millimeter continuum observations for samples of nearby early-type galaxies (21 sources) and nearby low-luminosity active galactic nuclei (LLAGN; 16 sources) at 100 GHz (3 mm) using the Nobeyama Millimeter Array (NMA). In addition, we performed quasi-simultaneous observations at 150 GHz (2 mm) and 100 GHz for five LLAGNs. Compact nuclear emissions showing flat or inverted spectra at centimeter-to-millimeter wavelengths were found in many LLAGNs and several early-type galaxies. Moreover, significant flux variability was detected in three LLAGNs. These radio properties are similar to Sgr A*. The observed radio luminosities are consistent with the fundamental plane of black hole activity that has suggested on the basis of samples with black hole masses ranging from 10 to 10^10 M_sun. This implies nuclear jets powered by radiatively inefficient accretion flows onto black holes.

We conducted very long baseline interferometry (VLBI) observations of five radio-loud narrow-line Seyfert 1 (NLS1) galaxies in milliarcsecond resolutions at 1.7 GHz (18 cm) using the Very Long Baseline Array (VLBA). Significant parsec-scale structures were revealed for three out of the five sources with high brightness temperature by direct imaging; this is convincing evidence for nonthermal jets. FBQS J1644+2619 with an inverted spectrum showed a prominent one-sided linear structure, indicating Doppler beaming with an intrinsic jet speed of >0.74c. FBQS J1629+4007, also with an inverted spectrum, showed rapid flux variability, indicating Doppler beaming with an intrinsic jet speed of >0.88c. Thus, we found convincing evidence that these two NLS1s can generate at least mildly or highly relativistic jets, which may make them apparently radio loud even if they are intrinsically radio quiet. On the other hand, the other three NLS1s had steep spectra and two of them showed significantly diffuse pc-scale structures, which were unlikely to be strongly beamed. Thus, some NLS1s have the ability to generate jets strong enough to make them intrinsically radio loud without Doppler beaming. NLS1s as a class show a number of extreme properties and radio-loud ones are very rare. We build on these radio results to understand that the central engines of radio-loud NLS1s are essentially the same as that of other radio-loud AGNs in terms of the formation of nonthermal jets.

There is growing evidence of relativistic jets in radio-loud narrow-line Seyfert 1 (RL-NLS1) galaxies. We constrain the observational properties of the radio emission in the first RL-NLS1 galaxy ever detected in gamma-rays, PMN J0948+0022, i.e., its flux density and structure in total intensity and in polarization, its compactness, and variability. We performed three real-time e-VLBI observations of PMN J0948+0022 at 22 GHz, using a global array including telescopes in Europe, East Asia, and Australia. These are the first e-VLBI science observations ever carried out with a global array, reaching a maximum baseline length of 12458 km. The observations were part of a large multiwavelength campaign in 2009. The source is detected at all three epochs. The structure is dominated by a bright component, more compact than 55 microarcsec, with a fainter component at a position angle theta~ 35deg. Relativistic beaming is required by the observed brightness temperature of 3.4x10^11 K. Polarization is detected at a level of about 1%. The parameters derived by the VLBI observations, in addition to the broad-band properties, confirm that PMN J0948+0022 is similar to flat spectrum radio quasars. Global e-VLBI is a reliable and promising technique for future studies.

The detection of gamma-ray emission by Fermi-LAT from the radio loud Narrow Line Seyfert 1 PMN J0948+0022 (Abdo et al. 2009, ApJ 699, 976) triggered a multi-wavelength campaign between March and July 2009. Given its high compactness (Doi et al. 2006, PASJ 58, 829), inverted spectrum, and 0deg declination, the source was an ideal target to observe at 22 GHz with a Global VLBI array extending from Europe to East Asia and Australia. In order to deliver prompt results to be analysed in combination with the other instruments participating in the campaign, the observations were carried out with real time VLBI, for the first time on a Global scale. Indeed, the main results have been published just a few months after the campaign (Abdo et al. 2009, ApJ 707, 727). Here we present additional details about the e-VLBI observations.

We report multi-epoch Very Long Baseline Interferometry (VLBI) study of the innermost jet of 3C 84. We carried out 14-epoch VLBI observations during 2006-2009 with the Japanese VLBI Network (JVN) and the VLBI Exploration of Radio Astrometry (VERA), immediately following the radio outburst that began in 2005. Comparison between VLBI lightcurve and single-dish lightcurve indicates that the outburst was associated with the central ~1 pc core. We found that this outburst accompanied the emergence of a new component, and the projected speed of this new component was 0:23c from 2007/297 (2007 October 24) to 2009/114 (2009 April 24). We argue the site of gamma-ray emission detected by Fermi/LAT and jet kinematics in connection with gamma-ray emission mechanism.

Using deep Chandra ACIS observation data for Cygnus A, we report evidence of non-thermal X-ray emission from radio lobes surrounded by a rich intra-cluster medium (ICM). The diffuse X-ray emission, which are associated with the eastern and western radio lobes, were observed in a 0.7--7 keV Chandra$ ACIS image. The lobe spectra are reproduced with not only a single-temperature Mekal model, such as that of the surrounding ICM component, but also an additional power-law (PL) model. The X-ray flux densities of PL components for the eastern and western lobes at 1 keV are derived as 77.7^{+28.9}_{-31.9} nJy and 52.4^{+42.9}_{-42.4} nJy, respectively, and the photon indices are 1.69^{+0.07}_{-0.13} and 1.84^{+2.90}_{-0.12}, respectively. The non-thermal component is considered to be produced via the inverse Compton (IC) process, as is often seen in the X-ray emission from radio lobes. From a re-analysis of radio observation data, the multiwavelength spectra strongly suggest that the seed photon source of the IC X-rays includes both cosmic microwave background radiation and synchrotron radiation from the lobes. The derived parameters indicate significant dominance of the electron energy density over the magnetic field energy density in the Cygnus A lobes under the rich ICM environment.

Multi-epoch Very Long Baseline Interferometry (VLBI) study of the sub-pc scale jet of 3C 84 is presented. We carried out 14-epoch VLBI observations during 2006-2009 with the Japanese VLBI Network (JVN) and the VLBI Exploration of Radio Astrometry (VERA), immediately following the radio outburst that began in 2005. We confirmed that the outburst was associated with the central ~1 pc core, accompanying the emergence of a new component. This is striking evidence of the recurrence of jet activity. The new component became brighter during 2008, in contrast to the constant gamma-ray emission that was observed with the Fermi Gamma-ray Space Telescope during the same time. We found that the projected speed of the new component is 0.23c from 2007/297 (2007 October 24) to 2009/114 (2009 April 24). The direction of movement of this component differs from that of the pre-existing component by ~40 degree. This is the first measurement of kinematics of a sub-pc jet in a gamma-ray active phase. Possible detection of jet deceleration and the jet kinematics in connection with the gamma-ray emission is discussed.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

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.

We argue that the origin of "FRI/FRI{-.1em}I dichotomy" -- the division between Fanaroff-Riley class I (FRI) with subsonic lobes and class I{-.1em}I (FRI{-.1em}I) radio sources with supersonic lobes is sharp in the radio-optical luminosity plane (Owen-White diagram) -- can be explained by the deceleration of advancing radio lobes. The deceleration is caused by the growth of the effective cross-sectional area of radio lobes. We derive the condition in which an initially supersonic lobe turns into a subsonic lobe, combining the ram-pressure equilibrium between the hot spots and the ambient medium with the relation between "the hot spot radius" and "the linear size of radio sources" obtained from the radio observations. We find that the dividing line between the supersonic lobes and subsonic ones is determined by the ratio of the jet power $L_{\rm j}$ to the number density of the ambient matter at the core radius of the host galaxy $\bar{n}_{\rm a}$. It is also found that there exists the maximal ratio of $(L_{\rm j}/\bar{n}_{\rm a})$ and its value resides in $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}\approx 10^{44-47} {\rm erg} {\rm s}^{-1} {\rm cm}^{3}$, taking account of considerable uncertainties. This suggests that the maximal value $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}$ separates between FRIs and FRI{-.1em}Is.

We present a new prediction of GeV $\gamma$-ray emission from radio lobes of young AGN jets. In the previous work of Kino et al. (2007), MeV $\gamma$-ray bremsstrahlung emission was predicted from young cocoons/radio-lobes in the regime of no coolings. In this study, we include cooling effects of bremsstrahlung emission and adiabatic loss. With the initial conditions determined by observed young radio lobes, we solve a set of equations describing the expanding lobe evolution. Then we find that the lobes initially have electron temperature of $\sim$GeV, and they cool down to $\sim$MeV by the adiabatic loss. Correspondingly, the lobes initially yield bright bremsstrahlung luminosity in $\sim$GeV range and they fade out. We estimate these $\gamma$-ray emissions and show that nearby young radio lobes could be detected with Fermi Gamma-ray Space Telescope.

Gamma-ray emission from cocoons of young radio galaxies is predicted. Considering the process of adiabatic injection of the shock dissipation energy and mass of the relativistic jet into the cocoon, we find that the thermal electron temperature of the cocoon is typically predicted to be of the order of $\sim$ MeV, and is determined only by the bulk Lorentz factor of the jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal Bremsstrahlung emissions at energies $\sim$MeV. Hotter cocoons (i.e., GeV) for younger sources are also discussed.

In order to reveal the long-term evolution of relativistic jets in active galactic nuclei (AGNs), we examine the dynamical evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)]. By comparing the observed relation between the hot spot size and the linear size of radio source with a coevolution model of hot spot and cocoon, we find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). The deceleration is caused by the growth of the cross-sectional area of the cocoon head. Moreover, by comparing the hot spot speed with the sound speed of the ambient medium, we predict that only CSOs whose initial advance speed is higher than 0.3-0.5c can evolve into FRIIs.

We examine the evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)], by comparing the relation between the hot spot size and the projected linear size with a coevolution model of hot spots and a cocoon. We take account of the deceleration effect by the cocoon head growth. We find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). This is ascribed to the change of the power-law index of ambient density profile in the MSO phase ($\sim $1 kpc). It is also found that the cocoon shape becomes nearly spherical or disrupted for MSOs, while an elongated morphology is predicted for CSOs and FRIIs. This seems to be consistent with the higher fraction of distorted morphology of MSOs than that of CSOs and FRI. Finally, we predict that only CSOs whose initial advance speed is higher than about 0.1c can evolve into FRIIs, comparing the hot spot speed with the sound speed of the ambient medium.

We present results from multifrequency polarimetry of NRAO 140 using the Very Long Baseline Array. These observations allow us to reveal the distributions of both the polarization position angle and the Faraday rotation measure (RM). These distributions are powerful tools to discern the projected and line-of-sight components of the magnetic field, respectively. We find a systematic gradient in the RM distribution, with its sign being opposite at either side of the jet with respect to the jet axis. The sign of the RM changes only with the direction of the magnetic field component along the line of sight, so this can be explained by the existence of helical magnetic components associated with the jet itself. We derive two constraints for the pitch angle of the helical magnetic field from the distributions of the RM and the projected magnetic field; the RM distribution indicates that the helical fields are tightly wound, while that of the projected magnetic field suggests they are loosely wound around the jet axis. This inconsistency may be explained if the Faraday rotator is not cospatial with theemitting region. Our results may point toward a physical picture in which an ultra-relativistic jet (spine) with a loosely wound helical magnetic field is surrounded by a sub-relativistic wind layer (sheath) with a tightly wound helical magnetic field.

The existence of a gradient in the Faraday rotation measure (RM) of the quasar 3C 273 jet is confirmed by follow-up observations. A gradient transverse to the jet axis is seen for more than 20 mas in projected distance. Taking account of the viewing angle, we estimate it to be more than 100 pc. Comparing to the distribution of the RM in 1995, we detect a time variation of it at the same distance from the core over 7 yr. We discuss the origin of the Faraday rotation based on this rapid time variation. We rule out foreground media such as a narrow-line region, and suggest a helical magnetic field in the sheath region as the origin of this gradient of the RM.

We performed phase-reference very long baseline interferometry (VLBI) observations on five radio-loud narrow-line Seyfert 1 galaxies (NLS1s) at 8.4 GHz with the Japanese VLBI Network (JVN). Each of the five targets (RXS J08066+7248, RXS J16290+4007, RXS J16333+4718, RXS J16446+2619, and B3 1702+457) in milli-Jansky levels were detected and unresolved in milli-arcsecond resolutions, i.e., with brightness temperatures higher than 10^7 K. The nonthermal processes of active galactic nuclei (AGN) activity, rather than starbursts, are predominantly responsible for the radio emissions from these NLS1s. Out of the nine known radio-loud NLS1s, including the ones chosen for this study, we found that the four most radio-loud objects exclusively have inverted spectra. This suggests a possibility that these NLS1s are radio-loud due to Doppler beaming, which can apparently enhance both the radio power and the spectral frequency.

We observed the narrow-line quasar SDSS J094857.3+002225, which has the highest known radio loudness for a narrow-line Seyfert~1 galaxy (NLS1), at 1.7--15.4 GHz with the Very Long Baseline Array (VLBA). This is the first very-long-baseline interferometry (VLBI) investigation for a radio-loud NLS1. We independently found very high brightness temperatures from (1) its compactness in a VLBA image and (2) flux variation among the VLBA observation, our other observations with the VLBA, and the Very Large Array (VLA). A Doppler factor larger than 2.7--5.5 was required to meet an intrinsic limit of brightness temperature in the rest frame. This is evidence for highly relativistic nonthermal jets in an NLS1. We suggest that the Doppler factor is one of the most crucial parameters determining the radio loudness of NLS1s. The accretion disk of SDSS J094857.3+002225 is probably in the very high state, rather than the high/soft state, by analogy with X-ray binaries with strong radio outbursts and superluminal jets such as GRS 1915+105.

We present a study of the kinematic and spectral ages of the Gigahertz-Peaked Spectrum (GPS) source CTD~93. Measurements of the hot spot separation over 8.5 yr show evidence of an increase. The separation rate along the source axis is 0.34$\pm0.11c$ (H$_{0}$=72 km s$^{-1}$ Mpc$^{-1}$), which results in a kinematic age of 2200$\pm$700 yr. Assuming that two hot spots are moving apart at equal speeds, we derive an advance speed of 0.17$\pm0.06c$. The radio lobe spectra show a high frequency steepening, as expected if energetic electrons lose energy by synchrotron radiation. The spectral break decreases with the distance from the hot spot in the northern component of CTD~93. This tendency is expected from the basic scenario of radio lobe evolution involving particle acceleration at the hot spots, with the radio lobes populated by high energy electrons which have leaked from the hot spots. Although a core-jet morphology for CTD~93 has previously been proposed, these results indicate that the morphology is similar to that of Compact Symmetric Object (CSO). From the spectral fits in the northern component we found a break frequency of 3.7 GHz at the edge of the lobe. The resultant spectral age is $\sim$300 yr assuming the equipartition magnetic field. This requires the advance speed of 0.26$c$, which shows a good agreement of the hot spot motion of 0.17$\pm0.06c$. Our results strongly support the hypothesis that CSOs are young radio sources.