We present the result of a search for orphan Gamma-Ray Burst (GRB) afterglows in the Pan-STARRS1 (PS1) COSMOS survey. There is extensive theoretical and observational evidence suggesting that GRBs are collimated jets; the direct observation of orphan GRB afterglows would further support this model. An optimal survey strategy is designed by coupling the PS1 survey with the Subaru/Hyper-Suprime-Cam (HSC) survey. The PS1 COSMOS survey, one of the survey fields in the PS1 Medium Deep Survey (PS1/MDS), searches a field of 7 deg$^2$ from December 2011 to January 2014, reaching a limiting magnitude R $\sim$ 23. The dense cadence of PS1/MDS is crucial for identifying transients, and the deep magnitude reached by the HSC survey (R $\sim$ 26) is important for evaluating potential GRB hosts. A transient classification method is employed to select potential orphan GRB afterglow candidates. After a thorough analysis of the transient and host galaxy properties, we conclude that there are no candidates in this survey field. The null result implies that the consideration of jet structures is essential for further orphan GRB afterglow surveys.

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.

Two components of jets associated with the afterglow of the gamma-ray burst GRB 160623A were observed with multi-frequency observations including long-term monitoring in a sub-millimetre range (230 GHz) using the SMA. The observed light curves with temporal breaks suggests on the basis of the standard forward-shock synchrotron radiation model that the X-ray radiation is narrowly collimated with an opening angle $\theta_{n,j}<\sim6^{\circ}$ whereas the radio radiation originated from wider jets ($\sim27^{\circ}$). The temporal and spectral evolutions of the radio afterglow agree with those expected from a synchrotron radiation modelling with typical physical parameters except for the fact that the observed wide jet opening angle for the radio emission is significantly larger than the theoretical maximum opening angle. By contrast, the opening angle of the X-ray afterglow is consistent with the typical value of GRB jets. Since the theory of the relativistic cocoon afterglow emission is similar to that of a regular afterglow with an opening angle of $\sim30^{\circ}$, the observed radio emission can be interpreted as the shocked jet cocoon emission. This result therefore indicates that the two components of the jets observed in the GRB 160623A afterglow is caused by the jet and the shocked jet cocoon afterglows.

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

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

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

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

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

We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $\alpha \approx -0.3$, where the flux density $S \propto \nu^\alpha$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales.

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.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud active galactic nuclei (AGNi) represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with a roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral distribution properties and it has been recently found that, at least in some nearby radio galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNi. Data combined from the sources considered here are not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNi needs to be re-examined by millimeter/sub-millimeter VLBI observations.

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.

The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided.

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

Centaurus A (Cen A) is a bright radio source associated with the nearby galaxy NGC 5128 where high-resolution radio observations can probe the jet at scales of less than a light-day. The South Pole Telescope (SPT) and the Atacama Pathfinder Experiment (APEX) performed a single-baseline very-long-baseline interferometry (VLBI) observation of Cen A in January 2015 as part of VLBI receiver deployment for the SPT. We measure the correlated flux density of Cen A at a wavelength of 1.4 mm on a $\sim$7000 km (5 G$\lambda$) baseline. Ascribing this correlated flux density to the core, and with the use of a contemporaneous short-baseline flux density from a Submillimeter Array observation, we infer a core brightness temperature of $1.4 \times 10^{11}$ K. This is close to the equipartition brightness temperature, where the magnetic and relativistic particle energy densities are equal. Under the assumption of a circular Gaussian core component, we derive an upper limit to the core size $\phi = 34.0 \pm 1.8~\mu\textrm{as}$, corresponding to 120 Schwarzschild radii for a black hole mass of $5.5 \times 10^7 M_{\odot}$.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of $\sim 30$% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that $\ell_1$+TSV regularization outperforms $\ell_1$+TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is $\lesssim 10$% of the traditional CLEAN beam size for $\ell_1$+TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the $\ell_1$+TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of $\sim 10-20$% in present simulations for Sgr A*.

In this work we interpret and discuss the time variable rotation measure (RM) found, for the first time over a 1-yr period, in the core region of a blazar. These results are based on a one-year, multi-frequency (15, 24, and 43 GHz) Very Long Baseline Array (VLBA) monitoring of the TeV blazar Markarian 421 (Mrk 421). We investigate the Faraday screen properties and its location with respect to the jet emitting region. Given that the 43 GHz radio core flux density and the RM time evolution suggest a similar trend, we explore the possible connection between the RM and the accretion rate. Among the various scenarios that we explore, the jet sheath is the most promising candidate for being the main source of Faraday rotation. During the one-year observing period the RM trend shows two sign reversals, which may be qualitatively interpreted within the context of the magnetic tower models. We invoke the presence of two nested helical magnetic fields in the relativistic jet with opposite helicities, whose relative contribution produce the observed RM values. The inner helical field has the poloidal component ($B_{\rm p}$) oriented in the observer's direction and produces a positive RM, while the outer helical field, with $B_{\rm p}$ in the opposite direction, produces a negative RM. We assume that the external helical field dominates the contribution to the observed RM, while the internal helical field dominates when a jet perturbation arises during the second observing epoch. Being the intrinsic polarization angle parallel to the jet axis, a pitch angle of the helical magnetic field $\phi\gtrsim 70^\circ$ is required. Additional scenarios are also considered to explain the observed RM sign reversals.

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 general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon originate. We derive a conservative formulation for the transport of energetic electrons which are initially injected at the stagnation surface and subsequently transported along flow streamlines. With this formulation the energy spectra evolution of the electrons along the flow in the presence of radiative and adiabatic cooling is determined. For flows regulated by synchrotron radiative losses and adiabatic cooling, the effective radio emission region is found to be finite, and geometrically it is more extended along the jet central axis. Moreover, the emission from regions adjacent to the stagnation surface is expected to be the most luminous as this is where the freshly injected energetic electrons concentrate. An observable stagnation surface is thus a strong prediction of the GRMHD jet model with the prescribed non-thermal electron injection. Future millimeter/sub-millimeter (mm/sub-mm) very-long-baseline interferometric (VLBI) observations of supermassive black hole candidates, such as the one at the center of M87, can verify this GRMHD jet model and its associated non-thermal electron injection mechanism.

Around a rapidly rotating black hole (BH), when the plasma accretion rate is much less than the Eddington rate, the radiatively inefficient accretion flow (RIAF) cannot supply enough MeV photons that are capable of materializing as pairs. In such a charge-starved BH magnetosphere, the force-free condition breaks down in the polar funnels. Applying the pulsar outer-magnetospheric lepton accelerator theory to super-massive BHs, we demonstrate that a strong electric field arises along the magnetic field lines in the direct vicinity of the event horizon in the funnels, that the electrons and positrons are accelerated up to 100~TeV in this vacuum gap, and that these leptons emit copious photons via inverse-Compton (IC) process between 0.1~TeV and 30~TeV for a distant observer. It is found that these IC fluxes will be detectable with Imaging Atmospheric Cherenkov Telescopes, provided that a low-luminosity active galactic nucleus is located within 1~Mpc for a million-solar-mass central BH or within 30~Mpc for a billion-solar-mass central BH. These very-high-energy fluxes are beamed in a relatively small solid angle around the rotation axis because of the inhomogeneous and anisotropic distribution of the RIAF photon field, and show an anti-correlation with the RIAF submillimeter fluxes. The gap luminosity little depends on the three-dimensional magnetic-field configuration, because the Goldreich-Julian charge density, and hence the exerted electric field is essentially governed by the frame-dragging effect, not by the magnetic field configuration.

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.

The structure of the NGC 6251 jet on the milliarcsecond scale is investigated using images taken with the European VLBI Network and the Very Long Baseline Array. We detect a structural transition of the jet from a parabolic to a conical shape at a distance of (1-2) x 10^5 times the Schwarzschild radius from the central engine, which is close to the sphere of gravitational influence of the supermassive black hole (SMBH). We also examine the jet pressure profiles with the synchrotron minimum energy assumption to discuss the physical origin of the structural transition. The NGC 6251 jet, together with the M87 jet, suggests a fundamental process of structural transition in the jets of active galactic nuclei (AGNs). Collimated AGN jets are characterized by their external galactic medium, showing that AGN jets interplay with the SMBH and its host galaxy.

Since the ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO), SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA) are working jointly to relocate the antenna to Greenland. This paper shows the status of the antenna retrofit and the work carried out after the recommissioning and subsequent disassembly of the antenna at the VLA has taken place. The next coming months will see the start of the antenna reassembly at Thule Air Base. These activities are expected to last until the fall of 2017 when commissioning should take place. In parallel, design, fabrication and testing of the last components are taking place in Taiwan.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud AGNs represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral energy distribution properties and it has been recently found that, at least in some nearby radio-galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNs. Data combined from the sources considered here is not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNs needs to be re-examined by mm/sub-mm VLBI observations.

We present the 3.5-yr monitoring results of 225 GHz opacity at the summit of the Greenland ice sheet (Greenland Summit Camp) at an altitude of 3200 m using a tipping radiometer. We chose this site as our submillimeter telescope (Greenland Telescope; GLT) site, because its location offers favorable baselines to existing submillimeter telescopes for global-scale VLBI. The site shows a clear seasonal variation with the average opacity lower by a factor of two during winter. For the winter quartiles of 25% and 50%, the Greenland site is about 10%-30% worse than the ALMA or the South Pole sites. Estimated atmospheric transmission spectra in winter season are similar to the ALMA site at lower frequencies (<450 GHz), which are transparent enough to perform astronomical observations almost all of the winter time with opacities <0.5, but 10%-25% higher opacities at higher frequencies (>450 GHz) than those at the ALMA site. This is due to the lower altitude of the Greenland site. Nevertheless, half of the winter time at the Greenland site can be used for astronomical observations at frequencies between 450 GHz and 1000 GHz with opacities <1.2, and 10% of the time show >10% transmittance in the THz (1035 GHz, 1350 GHz, and 1500 GHz) windows. One major advantage of the Greenland site in winter is that there is no diurnal variation due to the polar night condition, and therefore the durations of low-opacity conditions are significantly longer than at the ALMA site. Opacities lower than 0.05 or 0.04 can continue for more than 100 hours. Such long stable opacity conditions do not occur as often even at the South Pole; it happens only for the opacity lower than 0.05. Since the opacity variation is directly related to the sky temperature (background) variation, the Greenland site is suitable for astronomical observations that need unusually stable sky background.

A radiatively inefficient accretion flow (RIAF), which is commonly characterized by its sub-Keplerian nature, is a favored accretion model for the supermassive black hole at Galactic center, Sagittarius A$^{*}$. To investigate the observable features of a RIAF, we compare the modeled shadow images, visibilities, and spectra of three flow models with dynamics characterized by (i) a Keplerian shell which is rigidly-rotating outside the innermost stable circular orbit (ISCO) and infalling with a constant angular momentum inside ISCO, (ii) a sub-Keplerian motion, and (iii) a free-falling motion with zero angular momentum at infinity. At near-mm wavelengths the emission is dominated by the flow within several Schwarzschild radii. The energy shift due to the flow dynamics becomes important and distinguishable, suggesting that the flow dynamics are an important model parameter for interpreting the mm/submillimeter very long baseline interferometric observations with the forthcoming, fully assembled Event Horizon Telescope (EHT). As an example, we demonstrate that structural variations of Sagittarius A$^{*}$ on event horizon-scales detected in previous EHT observations can be explained by the non-stationary dynamics of a RIAF.

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.

We investigate the magnetic field configuration of the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. By stacking images over three years in order to enhance the sensitivity, we reveal, for the first time, systematic transverse gradients of the Faraday rotation measure in several knots along the jet. Combining this result with polarization properties and the dynamics of the jet, we suggest the magnetic structure in several knots at kiloparsec scales consists of a systematically wrapped, tightly wound helical configuration. Our analysis brings us a new paradigm where the M87 jet is a fundamentally current carrying system produced in the vicinity of the supermassive black hole, transferring a huge amount of the electromagnetic energy over the host galaxy scale.

We investigate the distribution of Faraday rotation measure (RM) in the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. We resolve the structure of the RM in several knots along the jet for the first time. We derive the power spectrum in the arcsecond scale jet and find indications that the RM cannot be associated with a turbulent magnetic field with 3D Kolmogorov spectrum. Our analysis indicates that the RM probed on jet scales has a significant contribution of a Faraday screen associated with the vicinity of the jet, in contrast with that on kiloparsec scales, typically assumed to be disconnected from the jet. Comparison with previous RM analyses suggests that the magnetic fields giving rise to the RMs observed in jet scales have different properties and are well less turbulent than these observed in the lobes.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

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 report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

A planned rapid submillimeter (submm) Gamma Ray Burst (GRBs) follow-up observations conducted using the Greenland Telescope (GLT) is presented. The GLT is a 12-m submm telescope to be located at the top of the Greenland ice sheet, where the high-altitude and dry weather porvides excellent conditions for observations at submm wavelengths. With its combination of wavelength window and rapid responding system, the GLT will explore new insights on GRBs. Summarizing the current achievements of submm GRB follow-ups, we identify the following three scientific goals regarding GRBs: (1) systematic detection of bright submm emissions originating from reverse shock (RS) in the early afterglow phase, (2) characterization of forward shock and RS emissions by capturing their peak flux and frequencies and performing continuous monitoring, and (3) detections of GRBs as a result of the explosion of first-generation stars result of GRBs at a high redshift through systematic rapid follow ups. The light curves and spectra calculated by available theoretical models clearly show that the GLT could play a crucial role in these studies.

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 ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011. SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA), SAO's main partner for this project, are working jointly to relocate the antenna to Greenland to carry out millimeter and submillimeter VLBI observations. This paper presents the work carried out on upgrading the antenna to enable operation in the Arctic climate by the GLT Team to make this challenging project possible, with an emphasis on the unexpected telescope components that had to be either redesigned or changed. Five-years of inactivity, with the antenna laying idle in the desert of New Mexico, coupled with the extreme weather conditions of the selected site in Greenland have it necessary to significantly refurbish the antenna. We found that many components did need to be replaced, such as the antenna support cone, the azimuth bearing, the carbon fiber quadrupod, the hexapod, the HVAC, the tiltmeters, the antenna electronic enclosures housing servo and other drive components, and the cables. We selected Vertex, the original antenna manufacturer, for the main design work, which is in progress. The next coming months will see the major antenna components and subsystems shipped to a site of the US East Coast for test-fitting the major antenna components, which have been retrofitted. The following step will be to ship the components to Greenland to carry out VLBI and single dish observations. Antenna reassembly at Summit Station should take place during the summer of 2018.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

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.

We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87.

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.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

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.

The structure and dynamics of the M87 jet from sub-milli-arcsec to arcsecond scales are continuously examined. We analysed the VLBA archival data taken at 43 and 86 GHz to measure the size of VLBI cores. Millimeter/sub-mm VLBI cores are considered as innermost jet emissions, which has been originally suggested by Blandford & K\"onigl. Those components fairly follow an extrapolated parabolic streamline in our previous study so that the jet has a single power-law structure with nearly five orders of magnitude in the distance starting from the vicinity of the supermassive black hole (SMBH), less than 10 Schwarzschild radius ($r_{\rm s}$). We further inspect the jet parabolic structure as a counterpart of the magnetohydrodynamic (MHD) nozzle in order to identify the property of a bulk acceleration. We interpret that the parabolic jet consists of Poynting-flux dominated flows, powered by large amplitude, nonlinear torsional Alfv\'en waves. We examine the non-relativistic MHD nozzle equation in a parabolic shape. The nature of trans-fast magnetosonic flow is similar to the one of transonic solution of Parker's hydrodynamic solar wind; the jet becomes super-escape as well as super-fast magnetosonic at around $\sim 10^{3}\,r_{\rm s}$, while the upstream trans-Alfv\'enic flow speed increases linearly as a function of the distance at $\sim 10^{2}$ - $10^{3}\, r_{\rm s}$. We here point out that this is the first evidence to identify these features in astrophysical jets. We propose that the M87 jet is magnetically accelerated, but thermally confined by the stratified ISM inside the sphere of gravitational influence of the SMBH potential, which may be a norm in AGN jets.

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

We report the latest results of 225 GHz atmospheric opacity measurements from two arctic sites; one on high coastal terrain near the Eureka weather station, on Ellesmere Island, Canada, and the other at the Summit Station near the peak of the Greenland icecap. This is a campaign to search for a site to deploy a new telescope for submillimeter Very Long Baseline Interferometry and THz astronomy in the northern hemisphere. Since 2011, we have obtained 3 months of winter data near Eureka, and about one year of data at the Summit Station. The results indicate that these sites offer a highly transparent atmosphere for observations in submillimeter wavelengths. The Summit Station is particularly excellent, and its zenith opacity at 225 GHz is statistically similar to the Atacama Large Milllimeter/submillimeter Array in Chile. In winter, the opacity at the Summit Station is even comparable to that observed at the South Pole.

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.

The structure of the M87 jet, from milli-arcsec to arcsecond scales, is extensively investigated, utilizing the images taken with the EVN, MERLIN and VLBA. We discover that the jet maintains a parabolic streamline over a range in sizescale equal to $10^{5}$ times the Schwarzschild radius. The jet then transitions into a conical shape further downstream. This suggests that the magnetohydrodynamic jet is initially subjected to the confinement by the external gas which is dominated by the gravitational influence of the supermassive black hole. Afterwards the jet is then freely expanding with a conical shape. This geometrical transition indicates that the origin of the HST-1 complex may be a consequence of the over-collimation of the jet. Our result suggests that when even higher angular resolution is provided by a future submm VLBI experiment, we will be able to explore the origin of active galactic nuclei jets.

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

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.

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

Both in the Faraday Rotation Measure and the intrinsic polarization angle, new features are revealed to indicate a helical magnetic field operating along the jet of the bright active galactic nuclei 3C 273. The helical field has been suggested to be related to the formation and collimation of jets by magnetohydrodynamic models. The distribution of the RM shows a systematic gradient with respect to the jet axis, which is expected by a helical magnetic field. In addition, the helical field can consistently explain two types in the direction of the projected magnetic field: parallel and perpendicular to the jet axis. Further, if the helical magnetic field is generated by winding up of an initial field by rotation of the accretion disk, we can uniquely determine the direction of the disk rotation, since the jet is approaching us.

We consider the use of polarization properties as a means to discriminate between Synchrotron Self-Absorption (SSA) and Free--Free Absorption (FFA) in GHz-Peaked Spectrum (GPS) sources. The polarization position angle (PA) of synchrotron radiation at high frequencies for the optically thin regime is perpendicular to the magnetic field, whereas it is parallel to the magnetic field at low frequencies for the optically thick regime. Therefore, SSA produces a change in PA of $90^{\circ}$ across the spectral peak, while FFA does not result in such a change. We analyzed polarization data from VLA observations for six GPS sources to see if such a change in PA was present. Our results indicate that there is no significant evidence for $90^{\circ}$ change in PA across the spectral peak, suggesting that FFA is more likely than SSA for low-frequency cutoffs in these sources.