Waseda University Nasu telescope array is a spatial fast Fourier transform (FFT) interferometer consisting of eight linearly aligned antennas with 20 m spherical dishes. This type of interferometer was developed to survey transient radio sources with an angular resolution as high as that of a 160 m dish with a field of view as wide as that of a 20 m dish. We have been performing drift-scan-mode observations, in which the telescope scans the sky around a selected declination as the earth rotates. The black hole X-ray binary V404 Cygni underwent a new outburst in 2015 June after a quiescent period of 26 years. Because of the interest in black hole binaries, a considerable amount of data on this outburst at all wavelengths was accumulated. Using the above telescope, we had been monitoring V404 Cygni daily from one month before the X-ray outburst, and two radio flares at 1.4 GHz were detected on June 21.73 and June 26.71. The flux density and time-scale of each flare were 313+/-30 mJy and 1.50+/-0.49 days, 364+/-30 mJy and 1.70+/-0.16 days, respectively. We have also confirmed the extreme variation of radio spectra within a short period by collecting other radio data observed with several radio telescopes. Such spectral behaviors are considered to reflect the change in the opacity of the ejected blobs associated with these extreme activities in radio and X-ray. Our 1.4 GHz radio data are expected to be helpful for studying the physics of the accretion and ejection phenomena around black holes.

We present the results of CO(1-0) observations of the host galaxy of a Type I superluminous supernova (SLSN-I), SN2017egm, one of the closest SLSNe-I at z = 0.03063, by using the Atacama Large Millimeter/submillimeter Array. The molecular gas mass of the host galaxy is $M_{\rm gas} = (4.8 \pm 0.3) \times 10^9$ $M_{\odot}$, placing it on the sequence of normal star-forming galaxies in an $M_{\rm gas}$-star-formation rate (SFR) plane. The molecular hydrogen column density at the location of SN2017egm is higher than that of the Type II SN PTF10bgl, which is also located in the same host galaxy, and those of other Type II and Ia SNe located in different galaxies, suggesting that SLSNe-I have a preference for a dense molecular gas environment. On the other hand, the column density at the location of SN2017egm is comparable to those of Type Ibc SNe. The surface densities of molecular gas and the SFR at the location of SN2017egm are consistent with those of spatially resolved local star-forming galaxies and follow the Schmidt-Kennicutt relation. These facts suggest that SLSNe-I can occur in environments with the same star-formation mechanism as in normal star-forming galaxies.

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

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 investigate the molecular gas in, and star-formation properties of, the host galaxy (CGCG 137-068) of a mysterious transient, AT2018cow, at kpc and larger scales, using archival band-3 data from the Atacama Large Millimeter/submillimeter Array (ALMA). AT2018cow is the nearest Fast-Evolving Luminous Transient (FELT), and this is the very first study unveiling molecular-gas properties of FELTs. The achieved rms and beam size are 0.21 mJy beam$^{-1}$ at a velocity resolution of $40$ km s$^{-1}$ and $3".66\times2".71$ ($1.1~{\rm kpc} \times 0.8~{\rm kpc}$), respectively. CO($J$=1-0) emission is successfully detected. The total molecular gas mass inferred from the CO data is $(1.85\pm0.04)\times10^8$ M$_\odot$ with the Milky Way CO-to-H$_2$ conversion factor. The H$_2$ column density at the AT2018cow site is estimated to be $8.6\times10^{20}$ cm$^{-2}$. The ALMA data reveal that (1) CGCG 137-068 is a normal star-forming (SF) dwarf galaxy in terms of its molecular gas and star-formation properties and (2) AT2018cow is located between a CO peak and a blue star cluster. These properties suggest on-going star formation and favor the explosion of a massive star as the progenitor of AT2018cow. We also find that CGCG 137-068 has a solar or super-solar metallicity. If the metallicity of the other FELT hosts is not higher than average, then some property of SF dwarf galaxies other than metallicity may be related to FELTs.

An observational result of a radio continuum and H92$\alpha$ radio recombination line of the Galactic Center Lobe (GCL), using the Yamaguchi 32 m radio telescope, is reported. The obtained spatial intensity distribution of the radio recombination line shows two distinctive ridge-like structures extending from the galactic plane vertically to the north at the eastern and western sides of the galactic center, which are connected to each other at a latitude of $1.2^{\circ}$ to form a loop-like structure as a whole. This suggests that most of the radio continuum emission of the GCL is free-free emission, and that the GCL is filled with thermal plasma. The east ridge of the GCL observed with the radio recombination line separates 30 pc from the radio arc, which has been considered as a part of the GCL, but coincides with a ridge of the radio continuum at a galactic longitude of $0^{\circ}$. The radial velocity of the radio recombination line is found to be between $-4$ and $+10$ km s$^{-1}$ across the GCL. This velocity is much smaller than the one expected from the galactic rotation, and hence indicates that the GCL exists apart from the galactic center. These characteristics of the GCL suggest that the long-standing hypothesis that the GCL was created by an explosive activity in the galactic center is unlikely, but favor that the GCL is a giant HII region.

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

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

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

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

We present the results of 3 GHz radio continuum observations of the 8 host galaxies of super-luminous supernovae (SLSNe) at $0.1 < z < 0.3$ by using the Karl G. Jansky Very Large Array. Four host galaxies are detected significantly, and two of them are found to have high star-formation rates (SFRs $>$ 20 $M_{\odot}$ yr$^{-1}$) derived from radio emission, making them the most intensely star-forming host galaxies among SLSN host galaxies. We compare radio SFRs and optical SFRs, and find that three host galaxies have an excess in radio SFRs by a factor of $>$2, suggesting the existence of dust-obscured star formation, which cannot be traced by optical studies. Two of the three host galaxies, which are located in the galaxy main sequence based on optical SFRs, are found to be above the main sequence based on their radio SFRs. This suggests a higher fraction of starburst galaxies in SLSN hosts than estimated in previous studies. We calculate extinction from the ratio between radio SFRs and dust-uncorrected optical SFRs and find that the hosts are on the trend of increasing extinction with metallicity, which is consistent with the relation in local star-forming galaxies. We also place a constraint on a pulsar-driven SN model, which predicts quasi-steady synchrotron radio emission.

We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN.

Mid-infrared properties are reported of the west hot spot of the radio galaxy Pictor A with the Wide-field Infrared Survey Explorer (WISE). The mid-infrared counterpart to the hot spot, WISE J051926.26-454554.1, is listed in the AllWISE source catalog. The source was detected in all the four WISE photometric bands. A comparison between the WISE and radio images reinforces the physical association of the \wise\ source to the hot spot. The WISE flux density of the source was carefully evaluated. A close investigation of the multi-wavelength synchrotron spectral energy distribution from the object reveals a mid-infrared excess at the wavelength of $\lambda=22$ $\mu$m with a statistical significance of $4.8 \sigma$ over the simple power-law extrapolation from the synchrotron radio spectrum. The excess is reinforced by single and double cutoff power-law modeling of the radio-to-optical spectral energy distribution. The synchrotron cutoff frequency of the main and excess components was evaluated as $7.1 \times 10^{14}$ Hz and $5.5 \times 10^{13}$ Hz, respectively. From the cutoff frequency, the magnetic field of the emission region was constrained as a function of the region size. In order to interpret the excess component, an electron population different from the main one dominating the observed radio spectrum is necessary. The excess emission is proposed to originate in a sub structure within the hot spot, in which the magnetic field is by a factor of a few stronger than that in the minimum energy condition. The relation of the mid-infrared excess to the X-ray emission is briefly discussed.

We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass of the dusty system is $\sim$ 0.2 $M_{\sun}$ at minimum and up to a few $M_{\sun}$ depending on the dust parameters. 6.7 GHz CH$_{3}$OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. This is clearly smaller than that of typical accretion disks in High-Mass star formation, indicating that the initial angular momentum was very small, or the CH$_{3}$OH masers selectively trace accreting material that has small angular momentum. In the former case, the initial specific angular momentum is estimated to be 8 $\times$ 10$^{20}$ ($M_{*}$$/$10 $M_{\sun}$)$^{0.5}$ cm$^{2}$ s$^{-1}$, or a significant fraction of the initial angular momentum was removed outside of 100 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution ($\sim$ 10 mas) thermal continuum and line observations with ALMA long baselines.

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations.

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

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

This chapter summarizes radio astrometry in relation to current very long baseline interferometry (VLBI) projects and describes its perspectives with the SKA. The scientific goals of the astrometry with the SKA have been discussed in the international and Japanese communities of researchers, whose major issues are shown here. We have demonstrated some of the issues, such as censuses of possible targets and the technical feasibility of astrometry in the SKA frequency bands. The preliminary results of our case studies on SKA astrometry are also presented. In addition, possible synergy and commensality of the SKA astrometric projects with those in the optical and infrared astrometric missions, especially JASMINE (Japan Astrometry Satellite Mission for INfrared Exploration) are discussed.

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

KaVA (KVN and VERA Array) is a new combined VLBI array with KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). First, we briefly review the imaging capabilities of KaVA array which actually achieves more than three times better dynamic range than that achieved by VERA alone. The KaVA images clearly show detailed structures of extended radio jets in AGNs. Next, we represent the key science program to be led by KaVA AGN sub working group. We will conduct the monitoring observations of Sgr A* and M87 because of the largeness of their central super-massive black hole angular sizes. The main science goals of the program are (i) testing magnetically-driven-jet paradigm by mapping velocity fields of the M87 jet, and (ii) obtaining tight constraints on physical properties of radio emitting region in Sgr A*.

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 report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individual maser features. This fact suggests that the H$_{2}$O maser traces a region where molecular gas is steadily accelerated. There were five maser flares during five-years monitoring, and maser distributions in four of them were densely monitored by the VLBI Exploration of Radio Astrometry (VERA). The overall distribution of the maser features suggests the presence of a bipolar jet, with the 3D kinematics indicating that it is almost face-on (inclination angle of $\sim$ 8$^{\fdg}$--17$^{\fdg}$ from the line-of-sight). Most of maser features were recurrently excited within a region of 100$\times$100 AU$^{2}$ around the radio continuum peak, while their spatial distributions significantly varied between each flare. This confirms that episodic propagations of outflow shocks recurrently invoke intermittent flare activities. We also measured annual parallax, deriving the source distance of 1.70 $^{+0.19}_{-0.16}$ kpc that is consistent with the commonly-used photometric distance.

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 have observed bursting variability of the 6.7 GHz methanol maser of G33.641-0.228. Five bursts were detected in the observation period of 294 days from 2009 to 2012. The typical burst is a large flux density rise in about one day followed by a slow fall. A non-typical burst observed in 2010 showed a large and rapid flux density enhancement from the stable state, but the rise and fall of the flux density were temporally symmetric and a fast fluctuation continued 12 days. On average, the bursts occurred once every 59 days, although bursting was not periodic. Since the average power required for causing the burst of order of 10^21 Js^-1 is far smaller than the luminosity of G33.641-0.228, a very small fraction of the source's power would be sufficient to cause the burst occasionally. The burst can be explained as a solar-flare like event in which the energy is accumulated in the magnetic field of the circumstellar disk, and is released for a short time. However, the mechanism of the energy release and the dust heating process are still unknown.

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 have detected periodic flares of the 6.7 GHz methanol maser from an intermediate-mass star-forming region IRAS22198+6336. The maser was monitored daily in 2011, 2012, and 2013. Six flares were observed with a period of 34.6 days. The variation pattern is intermittent, and the flux ratio of the flaring and the quiescent states exceeds 30. Such intermittent variation with the short period uniquely characterizes the variation of the IRAS22198+6336 maser. At least five spectral components were identified. The spectral components varied almost synchronously, but their peak times differed by 1.8 days. These characteristics can be explained by the colliding-wind binary model.

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 have measured the internal proper motions of the 6.7 GHz methanol masers associated with Cepheus A (Cep A) HW2 using Very Long Baseline Interferometery (VLBI) observations. We conducted three epochs of VLBI monitoring observations of the 6.7 GHz methanol masers in Cep A-HW2 with the Japanese VLBI Network (JVN) over the period between 2006-2008. In 2006, we were able to use phase-referencing to measure the absolute coordinates of the maser emission with an accuracy of a few milliarcseconds. We compared the maser distribution with other molecular line observations that trace the rotating disk. We measured the internal proper motions for 29 methanol maser spots, of which 19 were identified at all three epochs and the remaining ten at only two epochs. The magnitude of proper motions ranged from 0.2 to 7.4 km/s, with an average of 3.1 km/s. Although there are large uncertainties in the observed internal proper motions of the methanol maser spots in Cep A, they are well fitted by a disk that includes both rotation and infall velocity components. The derived rotation and infall velocities at the disk radius of 680 au are 0.5 +- 0.7 and 1.8 +- 0.7 km/s, respectively. Assuming that the modeled disk motion accurately represents the accretion disk around the Cep A-HW2 high-mass YSO, we estimated the mass infall rate to be 3 x 10^{-4} n_8 Msun/yr (n_8 is the gas volume density in units of 10^{8} cm^{-3}). The combination of the estimated mass infall rate and the magnitude of the fitted infall velocity suggests that Cep A-HW2 is at an evolutionary phase of active gas accretion from the disk onto the central high-mass YSO. The infall momentum rate is estimated to be 5 x 10^{-4} n_8 Msun/yr km/s, which is larger than the estimated stellar radiation pressure of the HW2 object, supporting the hypothesis that this object is in an active gas accretion phase.

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

We report on the results of X-ray and radio follow-up observations of two GeV gamma-ray sources 2FGL J0923.5+1508 and 2FGL J1502.1+5548, selected as candidates for high-redshift blazars from unassociated sources in the {\it Fermi} Large Area Telescope Second Source Catalog. We utilize the Suzaku satellite and the VLBI Exploration of Radio Astrometry (VERA) telescopes for X-ray and radio observations, respectively. For 2FGL J0923.5+1508, a possible radio counterpart NVSS J092357+150518 is found at 1.4 GHz from an existing catalog, but we do not detect any X-ray emission from it and derive a flux upper limit $F_{\rm 2-8 keV} <$ 1.37 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$. Radio observations at 6.7 GHz also result in an upper limit of $S_{\rm 6.7 GHz}$ $<$ 19 mJy, implying a steep radio spectrum that is not expected for a blazar. On the other hand, we detect X-rays from NVSS J150229+555204, the potential 1.4 GHz radio counterpart of 2FGL J1502.1+5548. The X-ray spectrum can be fitted with an absorbed power-law model with a photon index $\gamma$=1.8$^{+0.3}_{-0.2}$ and the unabsorbed flux is $F_{\rm 2-8 keV}$=4.3$^{+1.1}_{-1.0}$ $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$. Moreover, we detect unresolved radio emission at 6.7 GHz with flux $S_{\rm 6.7 GHz}$=30.1 mJy, indicating a compact, flat-spectrum radio source. If NVSS J150229+555204 is indeed associated with 2FGL J1502.1+5548, we find that its multiwavelength spectrum is consistent with a blazar at redshift $z \sim 3-4$.

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 first detection of a radio-continuum and molecular jet associated with a dominant blue-shifted maser source, G353.273+0.641. A radio jet is extended 3000 au along NW-SE direction. H$_{2}$O masers are found to be clustered in the root of a bipolar radio jet. A molecular jet is detected by thermal SiO ($\upsilon$ = 0, $J$ = 2-1) emission. The SiO spectrum is extremely wide (-120 -- +87 km s$^{-1}$) and significantly blue-shift dominated, similar to the maser emission. The observed geometry and remarkable spectral similarity between H$_{2}$O maser and SiO strongly suggests the existence of a maser-scale ($\sim$ 340 au) molecular jet that is enclosed by the extended radio jet. We propose a "disc-masking" scenario as the origin of the strong blue-shift dominance, where an optically thick disc obscures a red-shifted lobe of a compact jet.

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.

We measured the trigonometric annual parallax of H$_2$O maser source associated with the massive star-forming regions IRAS 06061+2151 with VERA. The annual parallax of $0.496\pm0.031$ mas corresponding to a distance of $2.02^{+0.13}_{-0.12}$ kpc was obtained by 10 epochs of observations from 2007 October to 2009 September. This distance was obtained with a higher accuracy than the photometric distance previously measured, and places IRAS 06061+2151 in the Perseus spiral arm. We found that IRAS 06061+2151 also has a peculiar motion of larger than 15 km s$^{-1}$ counter to the Galactic rotation. That is similar to five sources in the Perseus spiral arm, whose parallaxes and proper motions have already been measured with higher accuracy. Moreover, these sources move at on average 27 km s$^{-1}$ toward the Galactic center and counter to the Galactic rotation.