We observed polarization of the SiO rotational transitions from Orion Source I (SrcI) to probe the magnetic field in bipolar outflows from this high mass protostar. Both 43 GHz $J$=1-0 and 86 GHz $J$=2-1 lines were mapped with $\sim$20 AU resolution, using the Very Large Array (VLA) and Atacama Large Millimeter/Submillimeter Array (ALMA), respectively. The $^{28}$SiO transitions in the ground vibrational state are a mixture of thermal and maser emission. Comparison of the polarization position angles in the $J$=1-0 and $J$=2-1 transitions allows us to set an upper limit on possible Faraday rotation of $10^{4}$ radians m$^{-2}$, which would twist the $J$=2-1 position angles typically by less than 10 degrees. The smooth, systematic polarization structure in the outflow lobes suggests a well ordered magnetic field on scales of a few hundred AU. The uniformity of the polarization suggests a field strength of $\sim$30 milli-Gauss. It is strong enough to shape the bipolar outflow and possibly lead to sub-Keplerian rotation of gas at the base of the outflow. The strikingly high fractional linear polarizations of 80-90% in the $^{28}$SiO $v$=0 masers require anisotropic pumping. We measured circular polarizations of 60% toward the strongest maser feature in the $v$=0 $J$=1-0 peak. Anisotropic resonant scattering (ARS) is likely to be responsible for this circular polarization. We also present maps of the $^{29}$SiO $v$=0 $J$=2-1 maser and several other SiO transitions at higher vibrational levels and isotopologues.

We present results of continuum and spectral line observations with ALMA and 22 GHz water (H$_2$O) maser observations using KaVA and VERA toward a high-mass star-forming region, G25.82-0.17. Multiple 1.3 mm continuum sources are revealed, indicating the presence of young stellar objects (YSOs) at different evolutionary stages, namely an ultra-compact HII region, G25.82-E, a high-mass young stellar object (HM-YSO), G25.82-W1, and starless cores, G25.82-W2 and G25.82-W3. Two SiO outflows, at N-S and SE-NW orientations, are identified. The CH$_3$OH 8$_{-1}$-7$_{0}$ E line, known to be a class I CH$_3$OH maser at 229 GHz is also detected showing a mixture of thermal and maser emission. Moreover, the H$_2$O masers are distributed in a region ~0.25" shifted from G25.82-W1. The CH$_3$OH 22$_{4}$-21$_{5}$ E line shows a compact ring-like structure at the position of G25.82-W1 with a velocity gradient, indicating a rotating disk or envelope. Assuming Keplerian rotation, the dynamical mass of G25.82-W1 is estimated to be $>$25 M$_{\odot}$ and the total mass of 20 M$_\odot$-84 M$_\odot$ is derived from the 1.3 mm continuum emission. The driving source of the N-S SiO outflow is G25.82-W1 while that of the SE-NW SiO outflow is uncertain. Detection of multiple high-mass starless$/$protostellar cores and candidates without low-mass cores implies that HM-YSOs could form in individual high-mass cores as predicted by the turbulent core accretion model. If this is the case, the high-mass star formation process in G25.82 would be consistent with a scaled-up version of low-mass star formation.

With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/submillimeter Array (ALMA). The two data sets used in this study were taken independently using different array configurations at Band 7 (330 GHz), targeting the protoplanetary disk around HD 142527; one in the shorter-baseline array configuration (~ 430 m), and the other in the longer-baseline array configuration (~ 1570 m). The image resolutions reconstructed from the two data sets are different by a factor of ~ 3. We confirm that the previously known disk structures appear on the images produced by both SpM and CLEAN at the standard beam size. The image reconstructed from the shorter-baseline data using the SpM matches that obtained with the longer-baseline data using CLEAN, achieving a super-resolution image from which a structure finer than the beam size can be reproduced. Our results demonstrate that on-going intensive development in the SpM imaging technique is beneficial to imaging with ALMA.

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

We report results of astrometric VLBI observations toward a Mira variable star BX~Cam using the VLBI array "VERA". The observations were performed from February 2012 to November 2014. Obtained parallax is 1.73$\pm$0.03 mas corresponding to a distance of 0.58$\pm$0.01 kpc. Parallax of this source was also reported in Gaia DR2 as 4.13$\pm$0.25 mas, and there is a 240 \% difference between these two measurements. Astrometric results from our VLBI observations show that we exactly traced angular motions of the seven maser spots in BX~Cam. We calculated stellar luminosities using both parallaxes, and obtained luminosities of $L_{\ast}^{\mathrm{VERA}} = 4950\pm170 L_{\odot}$ and $L_{\ast}^{\mathrm{Gaia}} = 870\pm110 L_{\odot}$. Deduced luminosities also support a validity of the parallax that we determined with VERA. Evaluating the two parallaxes, we concluded that the parallax of 1.73$\pm$0.03 mas from the VERA observations is correct for BX~Cam. We obtained a systemic motion of BX~Cam as (${\mu_{\alpha}}\cos{\delta}^{\mathrm{sys}}$, ${\mu}_{\delta}^{\mathrm{sys}}$) $=$ (13.48$\pm$0.14, $-$34.30$\pm$0.18) mas\,yr$^{-1}$. A total of 73 H$_2$O maser spots detected from our VLBI observations show a spatial distribution of 30 au $\times$ 80 au with a strong elongation along north-south direction. They show outflows with a three-dimensional velocity of 14.79$\pm$1.40 km\,s$^{-1}$. From a comparison between time variations of $V$-band magnitudes and H$_2$O maser, we found that variation of the H$_2$O maser is relevant to that seen in $V$-band even though the H$_2$O maser does not recover its maximum flux in each cycle.

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.

The black hole spacetime is described by general relativity and characterized by two quantities: the black hole mass and spin. Black hole spin measurement requires information from the vicinity of the event horizon, which is spatially resolved for the Galactic center SagittariusA* (SgrA*) and nearby radio galaxy M87 by means of very long baseline interferometry (VLBI) observations with the Event Horizon Telescope (EHT). In this paper, we simulate EHT observations for a gas cloud intermittently falling onto a black hole, and construct a method for spin measurement based on its relativistic flux variation. The light curve of the infalling gas cloud is composed of peaks formed by photons which directly reach a distant observer and by secondary ones reaching the observer after more than one rotation around the black hole. The time interval between the peaks is determined by a period of photon rotation near the photon circular orbit which uniquely depends on the spin. We perform synthetic EHT observations for SgrA* under a more realistic situation that a number of gas clouds intermittently fall towards the black hole with various initial parameters. Even for this case, the black hole spin dependence is detectable in correlated flux densities which are accurately calibrated by baselines between sites with redundant stations. The synthetic observations indicate that our methodology can be applied to EHT observations of Sgr A* since April 2017.

We report measurements of parallax and proper motion for four 22 GHz water maser sources as part of VERA Outer Rotation Curve project. All sources show Galactic latitudes of $>$ 2$^{\circ}$ and Galactocentric distances of $>$ 11 kpc at the Galactic longitude range of 95$^{\circ}$ $< l <$ 126$^{\circ}$. The sources trace the Galactic warp reaching to 200$\sim$400 pc, and indicate the signature of the warp to 600 pc toward the north Galactic pole. The new results along with previous results in the literature show the maximum height of the Galactic warp is increased with Galactocentric distance. Also, we examined velocities perpendicular to the disk for the sample, and found an oscillatory behavior between the vertical velocities and Galactic heights. This behavior suggests the existence of the bending (vertical density) waves, possibly induced by a perturbing satellite (e.g. passage of the Sagittarius dwarf galaxy).

As part our investigation into the Galactic rotation curve, we carried out Very Long Baseline Interferometry (VLBI) observations towards the star-forming region IRAS 01123+6430 using VLBI Exploration of Radio Astrometry (VERA) to measure its annual parallax and proper motion. The annual parallax was measured to be 0.151+/-0.042 mas, which corresponds to a distance of D=6.61^{+2.55}_{-1.44} kpc, and the obtained proper motion components were ({\mu}_{\alpha}cos{\delta}, {\mu}_{\delta})=(-1.44+/-0.15, -0.27+/-0.16) mas yr^{-1} in equatorial coordinates. Assuming Galactic constants of (R_0, {\Theta}_0)=(8.05+/-0.45 kpc, 238+/-14 km s^{-1}), the Galactocentric distance and rotation velocity were measured to be (R, {\Theta})=(13.04+/-2.24 kpc, 239+/-22 km s^{-1}), which are consistent with a flat Galactic rotation curve. The newly estimated distance provides a more accurate bolometric luminosity of the central young stellar object, L_Bol=(3.11+/-2.86)\times 10^3 L_solar, which corresponds to a spectral type of B1-B2. The analysis of 12CO(J=1-0) survey data obtained with the Five College Radio Astronomical Observatory (FCRAO) 14 m telescope shows that the molecular cloud associated with IRAS 01123+6430 consists of arc-like and linear components, which well matches a structure predicted by numerical simulation of the cloud-cloud collision (CCC) phenomenon. The coexistence of arc-like and linear components implies that the relative velocity of initial two clouds was as slow as 3-5 km s^{-1}, which meets the expected criteria of massive star formation where the core mass is effectively increased in the presence of low relative velocity (~3-5 km s^{-1}), as suggested by Takahira et al.(2014).

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

Using the cross-matched data of Gaia DR2 and 2MASS Point Source Catalog, we investigated the surface density distribution of stars aged ~1 Gyr in the thin disk in the range of 90{\deg} <= l <= 270{\deg}. We selected 4,654 stars above the turnoff corresponding to the age ~1 Gyr, that fall within a small box region in the color-magnitude diagram, (J-Ks)0 versus M(Ks), for which the distance and reddening are corrected. The selected sample shows an arm-like overdensity at 90{\deg} <= l <= 190{\deg}.This overdensity is located close to the Local arm traced by high-mass star forming regions (HMSFRs), but its pitch angle is slightly larger than that of the HMSFR-defined arm. Although the significance of the overdensity we report is marginal, its structure poses questions concerning both of the competing scenarios of spiral arms, the density-wave theory and the dynamic spiral arm model. The offset between the arms traced by stars and HMSFRs, i.e., gas, is difficult to be explained by the dynamic arm scenario. On the other hand, the pitch angle of the stellar Local arm, if confirmed, larger than that of the Perseus arm is difficult to be explained with the classical density-wave scenario. The dynamic arm scenario can explain it if the Local arm is in a growing up phase, while the Perseus arm is in a disrupting phase. Our result provides a new and complex picture of the Galactic spiral arms, and encourages further studies.

This paper presents a survey of microwave front-end receivers installed at radio telescopes throughout the World. This unprecedented analysis was conducted as part of a review of front-end developments for Italian radio telescopes, initiated by the Italian National Institute for Astrophysics in 2016. Fifteen international radio telescopes have been selected to be representative of the instrumentation used for radio astronomical observations in the frequency domain from 300 MHz to 116 GHz. A comprehensive description of the existing receivers is presented and their characteristics are compared and discussed. The observing performances of the complete receiving chains are also presented. An overview of on-going developments illustrates and anticipates future trends in front-end projects to meet the most ambitious scientific research goals.

We report on the first bird's-eye view of the innermost accretion disk around the high-mass protostellar object G353.273+0.641, taken by Atacama Large Millimter/submillimeter Array long-baselines. The disk traced by dust continuum emission has a radius of 250 au, surrounded by the infalling rotating envelope traced by thermal CH$_3$OH lines. This disk radius is consistent with the centrifugal radius estimated from the specific angular momentum in the envelope. The lower-limit envelope mass is $\sim$5-7 M$_{\odot}$ and accretion rate onto the stellar surface is 3 $\times$ 10$^{-3}$ M$_{\odot}$ yr$^{-1}$ or higher. The expected stellar age is well younger than 10$^{4}$ yr, indicating that the host object is one of the youngest high-mass objects at present. The disk mass is 2-7 M$_{\odot}$, depending on the dust opacity index. The estimated Toomre's $Q$ parameter is typically 1-2 and can reach 0.4 at the minimum. These $Q$ values clearly satisfy the classical criteria for the gravitational instability, and are consistent with the recent numerical studies. Observed asymmetric and clumpy structures could trace a spiral arm and/or disk fragmentation. We found that 70$\%$ of the angular momentum in the accretion flow could be removed via the gravitational torque in the disk. Our study has indicated that the dynamical nature of a self-gravitating disk could dominate the early phase of high-mass star formation. This is remarkably consistent with the early evolutionary scenario of a low-mass protostar.

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.

With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI astrometry provides unique data for studying stellar populations and Galactic structure. The BeSSel project not only provides parallax and proper motions at much larger distances than Gaia can reach, but it also uniquely samples the spiral arms of the Galaxy. The evolved stars in the BAaDE sample can potentially constrain the dynamics and stellar content of the inner bulge and bar of the Milky Way, not reachable in the optical.

We present high-resolution images of the submillimeter SiO line emissions of a massive young stellar object Orion Source I using the Atacama Large Millimeter/ Submillimeter Array (ALMA) at band 8. We detected the 464 GHz SiO v=4 J=11-10 line in Source I, which is the first detection of the SiO v=4 line in star-forming regions, together with the 465 GHz 29SiO v=2 J=11-10 and the 428 GHz SiO v=2 J=10-9 lines with a resolution of 50 AU. The 29SiO v=2 J=11-10 and SiO v=4 J=11-10 lines have compact structures with the diameter of <80 AU. The spatial and velocity distribution suggest that the line emissions are associated with the base of the outflow and the surface of the edge-on disk. In contrast, SiO v=2 J=10-9 emission shows a bipolar structure in the direction of northeast-southwest low-velocity outflow with ~200 AU scale. The emission line exhibits a velocity gradient along the direction of the disk elongation. With the assumption of the ring structure with Keplerian rotation, we estimated the lower limit of the central mass to be 7 solar mass and the radius of 12 AU< r <26 AU.

The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20 GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (`uv-coverage') by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging.

Many studies have shown that there are clear sequences in the period-luminosity relationship (PLR) for Mira variables and semiregular variables (SRVs) in the Large Magellanic Cloud (LMC). To investigate the PLR for SRVs in our galaxy, we examined the annual parallax measurement and conducted K'-band photometric monitoring of an SRV star SV Pegasus (SV Peg). We measured the position change of the associating H$_2$O maser spots by phase-referencing VLBI observations with VERA at 22 GHz, spanning approximately 3 yr, and detected an annual parallax of $\pi = 3.00 \pm 0.06$ mas, corresponding to a distance of $D=333 \pm 7$ pc. This result is in good agreement with the Hipparcos parallax and improves the accuracy of the distance from 35 $\%$ to 2 $\%$. However, the GAIA DR2 catalog gave a parallax of $\pi=1.12\pm0.28$ mas for SV Peg. This indicates that the GAIA result might be blurred by the effect of the stellar size because the estimated stellar radius was $\sim 5$ mas, which is comparable to the parallax. We obtained a K'-band mean magnitude of $m_{K'} = -0.48$ mag and a period of $P=177$ days from our photometric monitoring with a 1-m telescope. Using the trigonometric distance, we derived an absolute magnitude of $M_{K'}=-8.09 \pm 0.05$ mag. This result shows that the position of SV Peg in the PLR falls on the C' sequence found in the PLR in the LMC, which is similar to other SRVs in our galaxy.

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.

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

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

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

One of the outstanding problems in star-formation theory concerns the transfer of angular momentum such that mass can accrete onto a newly born young stellar object (YSO). From a theoretical standpoint, outflows and jets are predicted to play an essential role in angular momentum transfer and their rotation motions have been reported for both low- and high-mass YSOs. However, little quantitative discussion on outflow launching mechanisms have been presented for high-mass YSOs due to a lack of observational data. Here we present a clear signature of rotation in the bipolar outflow driven by Orion Source I, a high-mass YSO candidate, using the Atacama Large Millimeter/Submillimeter Array (ALMA). A rotational transition of silicon monoxide (Si18O) reveals a velocity gradient perpendicular to the outflow axis which is consistent with that of the circumstellar disk traced by a high-excitation water (H2O) line. The launching radii and outward velocity of the outflow are estimated to be >10 au and 10 km s-1, respectively. These parameters rule out a possibility that the observed outflow is produced by entrainment of a high-velocity jet, and that contribution from stellar-wind or X-wind which have smaller launching radii are significant in the case of Source I. Thus, present results provide a convincing evidence of a rotating outflow directly driven by the magneto-centrifugal disk wind launched by a high-mass YSO candidate.

We present observational results of the submillimeter H2O and SiO lines toward a candidate high-mass young stellar object Orion Source I using ALMA. The spatial structures of the high excitation lines at lower-state energies of >2500 K show compact structures consistent with the circumstellar disk and/or base of the northeast-southwest bipolar outflow with a 100 au scale. The highest excitation transition, the SiO (v=4) line at band 8, has the most compact structure. In contrast, lower-excitation transitions are more extended than 200 au tracing the outflow. Almost all the line show velocity gradients perpendicular to the outflow axis suggesting rotation motions of the circumstellar disk and outflow. While some of the detected lines show broad line profiles and spatially extended emission components indicative of thermal excitation, the strong H2O lines at 321 GHz, 474 GHz, and 658 GHz with brightness temperatures of >1000 K show clear signatures of maser action.

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

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.

The FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project is one of the legacy projects using the new multi-beam FOREST receiver installed on the Nobeyama 45-m telescope. This project aims to investigate the distribution, kinematics, and physical properties of both diffuse and dense molecular gas in the Galaxy at once by observing 12CO, 13CO, and C18O J=1-0 lines simultaneously. The mapping regions are a part of the 1st quadrant (10d < l < 50d, |b| < 1d) and the 3rd quadrant (198d < l <236d, |b| < 1d) of the Galaxy, where spiral arms, bar structure, and the molecular gas ring are included. This survey achieves the highest angular resolution to date (~20") for the Galactic plane survey in the CO J=1-0 lines, which makes it possible to find dense clumps located farther away than the previous surveys. FUGIN will provide us with an invaluable dataset for investigating the physics of the galactic interstellar medium (ISM), particularly the evolution of interstellar gas covering galactic scale structures to the internal structures of giant molecular clouds, such as small filament/clump/core. We present an overview of the FUGIN project, observation plan, and initial results, which reveal wide-field and detailed structures of molecular clouds, such as entangled filaments that have not been obvious in previous surveys, and large-scale kinematics of molecular gas such as spiral arms.

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.

From measuring the annual parallax of water masers over one and a half years with VERA, we present the trigonometric parallax and corresponding distance of another newly identified water maser source in the region of IRAS 20231+3440 as $\pi=0.611\pm0.022$ mas and $D=1.64\pm0.06$ kpc respectively. We measured the absolute proper motions of all the newly detected maser spots (30 spots) and presented two pictures describing the possible spatial distribution of the water maser as the morphology marks out an arc of masers whose average proper motion velocity in the jet direction was 14.26 km s$^{-1}$. As revealed by the ALLWISE composite image, and by applying the colour-colour method of YSO identification and classification on photometric archived data, we identified the driving source of the north maser group to be a class I, young stellar object. To further probe the nature of the progenitor, we used the momentum rate maximum value (1.2$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$ km s$^{-1}$) of the outflow to satisfy that the progenitor under investigation is a low mass young stellar object concurrently forming alongside an intermediate-mass YSO $\sim 60,000$ au ($\sim 37$ arcsecs) away from it.

We detected flaring flux variability that regularly occurred with the period of 23.9 days on a 6.7 GHz methanol maser emission at Vlsr = 25.30 km/s in G 014.23-00.50 through highly frequent monitoring using the Hitachi 32-m radio telescope. By analyzing data from 05 January 2013 to 21 January 2016, the periodic variability has persisted in at least 47 cycles, corresponding to approximately 1,100 days. The period of 23.9 days is the shortest one observed in masers at around high-mass young stellar objects so far. The flaring component normally falls below the detection limit (3 sigma) of 0.9 Jy. In the flaring periods, the component rises above the detection limit with the ratio of the peak flux density more than 180 in comparison with a quiescent phase, showing intermittent periodic variability. The time-scale of the flux rise was typically two days or shorter, and both symmetric and asymmetric profiles of flux variability were observed through intraday monitoring. These characteristics might be explained by a change in the flux of seed photons by a colliding-wind binary (CWB) system or a variation of the dust temperature by an extra heating source of a shock formed by the CWB system within a gap region in a circumbinary disk, in which the orbital semi-major axes of the binary are 0.26-0.34 au.

We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.

We have measured continuum flux densities of a high-mass protostar candidate, a radio source I in the Orion KL region (Orion Source I) using the Atacama Large Millimeter/Submillimeter Array (ALMA) at band 8 with an angular resolution of 0.1". The continuum emission at 430, 460, and 490 GHz associated with Source I shows an elongated structure along the northwest-southeast direction perpendicular to the so-called low-velocity bipolar outflow. The deconvolved size of the continuum source, 90 au times 20 au, is consistent with those reported previously at other millimeter/submillimeter wavelength. The flux density can be well fitted to the optically thick black-body spectral energy distribution (SED), and the brightness temperature is evaluated to be 700-800 K. It is much lower than that in the case of proton-electron or H- free-free radiations. Our data are consistent with the latest ALMA results by Plambeck & Wright (2016), in which the continuum emission have been proposed to arise from the edge-on circumstellar disk via thermal dust emission, unless the continuum source consists of an unresolved structure with the smaller beam filling factor.

We present observations made with the Australia Telescope Compact Array (ATCA), the Jansky Very Large Array (JVLA) and the Giant Metre-Wave Telescope of the radio source within the galaxy WISE~J071634.59-190039.2, claimed to be host of FRB~150418 by Keane et al. (2016). We have established a common flux density scale between the ATCA and JVLA observations, the main result of which is to increase the flux densities obtained by Keane et al. At a frequency of 5.5 GHz, the source has a mean flux density of 140uJy and is variable on short timescales with a modulation index of 0.36. Statistical analysis of the flux densities shows that the variations seen are consistent with refractive interstellar scintillation of the weak active galactic nucleus at the centre of the galaxy. It may therefore be the case that the FRB and the galaxy are not associated. However, taking into account the rarity of highly variable sources in the radio sky, and our lack of knowledge of the progenitors of FRBs as a class, the association between WISE~J071634.59-190039.2 and FRB~150418 remains a possibility.

Recent experimental results have confirmed a possible reduction in the GT$_+$ strengths of pf-shell nuclei. These proton-rich nuclei are of relevance in the deflagration and explosive burning phases of Type Ia supernovae. While prior GT strengths result in nucleosynthesis predictions with a lower-than-expected electron fraction, a reduction in the GT$_+$ strength can result in an slightly increased electron fraction compared to previous shell model predictions, though the enhancement is not as large as previous enhancements in going from rates computed by Fuller, Fowler, and Newman based on an independent particle model. A shell model parametrization has been developed which more closely matches experimental GT strengths. The resultant electron-capture rates are used in nucleosynthesis calculations for carbon deflagration and explosion phases of Type Ia supernovae, and the final mass fractions are compared to those obtained using more commonly-used rates.

We have measured an annual parallax of the Mira variable R~Ursae~Majoris (R~UMa) with the VLBI exploration for Radio Astronomy (VERA). From the monitoring VLBI observations spanning about two years, we detected H$_2$O maser spots in the LSR velocities ranges from 37 to 42 km\,s$^{-1}$. We derived an annual parallax of 1.97$\pm$0.05\,mas, and it gives a corresponding distance of 508$\pm$13\,pc. The VLBI maps revealed 72 maser spots distributed in $\sim$110 au area around an expected stellar position. Circumstellar kinematics of the maser spots were also revealed by subtracting a systemic motion in the Hipparcos catalog from proper motions of each maser spots derived from our VLBI observations. Infrared photometry is also conducted to measure a $K$ band apparent magnitude, and we obtained a mean magnitude of $m_K$ = 1.19$\pm$0.02\,mag. Using the trigonometric distance, the $m_K$ is converted to a $K$ band absolute magnitude of $M_K = -$7.34$\pm$0.06\,mag. This result gives a much more accurate absolute magnitude of R~UMa than previously provided. We solved a zero-point of $M_K - \log P$ relation for the Galactic Mira variables and obtained a relation of $M_K = -$3.52 $\log P$ + (1.09 $\pm$ 0.14). Other long period variables including red supergiants, whose distances were determined from astrometric VLBI, were also compiled to explore the different sequences of $M_K - \log P$ relation.

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

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

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

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.

We present an observational study of the vibrationally excited H$_{2}$O line at 658 GHz ($\nu_{2}$=1, 1$_{1, 0}$-1$_{0, 1}$) toward Orion KL using the Atacama Large Millimeter/Submillimeter Array (ALMA). This line is clearly detected at the position of the massive protostar candidate, the Source I. The spatial structure is compact with a size of about 100 AU and is elongated along the northeast-southwest low-velocity (18 km s$^{-1}$) bipolar outflow traced by 22 GHz H$_{2}$O masers, SiO masers, and thermal SiO lines. A velocity gradient can be seen perpendicular to the bipolar outflow. Overall spatial and velocity structure seems analogous to that of the 321 GHz H$_{2}$O maser line previously detected with ALMA and vibrationally excited SiO maser emission. The brightness temperature of the 658 GHz H$_{2}$O line is estimated to be higher than 2$\times$10$^{4}$ K, implying that it is emitted via maser action. Our results suggest that the 658 GHz H$_{2}$O maser line is emitted from the base of the outflow from a rotating and expanding accretion disk as observed for the SiO masers and the 321 GHz H$_{2}$O maser. We also search for two other H$_{2}$O lines at 646 GHz (9$_{7, 3}$-8$_{8, 0}$ and 9$_{7, 2}$-8$_{8, 1}$), but they are not detected in Orion KL.

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

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.

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 the trigonometric parallax of IRAS 07427-2400 with VERA to be 0.185 $\pm$ 0.027 mas, corresponding to a distance of 5.41$^{+0.92}_{-0.69}$ kpc. The result is consistent with the previous result of 5.32$^{+0.49}_{-0.42}$ kpc obtained by Choi et al. (2014) within error. To remove the effect of internal maser motions (e.g., random motions), we observed six maser features associated with IRAS 07427-2400 and determined systematic proper motions of the source by averaging proper motions of the six maser features. The obtained proper motions are ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$1.79 $\pm$ 0.32, 2.60 $\pm$ 0.17) mas yr$^{-1}$ in equatorial coordinates, while Choi et al. (2014) showed ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$2.43 $\pm$ 0.02, 2.49 $\pm$ 0.09) mas yr$^{-1}$ with one maser feature. Our astrometry results place the source in the Perseus arm, the nearest main arm in the Milky Way. Using our result with previous astrometry results obtained from observations of the Perseus arm, we conducted direct (quantitative) comparisons between 27 astrometry results and an analytic gas dynamics model based on the density-wave theory and obtained two results. First is the pitch angle of the Perseus arm determined by VLBI astrometry, 11.1 $\pm$ 1.4 deg, differing from what is determined by the spiral potential model (probably traced by stars), $\sim$ 20 deg. The second is an offset between a dense gas region and the bottom of the spiral potential model. The dense gas region traced by VLBI astrometry is located downstream of the spiral potential model, which was previously confirmed in the nearby grand-design spiral galaxy M51 in Egusa et al. (2011).

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

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

We conducted astrometric VLBI observations of water-vapor maser emission in the massive star forming region IRAS 21379+5106 to measure the annual parallax and proper motion, using VERA. The annual parallax was measured to be $0.262 \pm 0.031$ mas corresponding to a trigonometric distance of $3.82^{+0.51}_{-0.41}$ kpc. The proper motion was $(\mu_\alpha\cos{\delta}, \mu_\delta)=(-2.74 \pm 0.08, -2.87 \pm 0.18)$ mas yr$^{-1}$. Using this result, the Galactic rotational velocity was estimated to be $V_\theta=218\pm 19$ km s$^{-1}$ at the Galactocentric distance $R=9.22\pm0.43$ kpc, when we adopted the Galactic constants $R_0=8.05\pm 0.45$ kpc and $V_0=238\pm 14$ km s$^{-1}$. With newly determined distance, {the bolometric luminosity of the central young stellar object was re-evaluated to $(2.15\pm 0.54)\times 10^3 L_\odot$, which corresponds to spectral type of} B2--B3. Maser features were found to be distributed along a straight line from south-west to north-east. In addition, a vector map of the internal motions constructed from the residual proper motions implies that maser features trace a bipolar flow and that it cannot be explained by simple ballistic motion.

We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.

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 have carried out high resolution observations with Atacama Large Millimeter/Submillimeter Array (ALMA) of continuum emission from Orion KL region. We identify 11 compact sources at ALMA band 6 (245 GHz) and band 7 (339 GHz), including Hot Core, Compact Ridge, SMA1, IRc4, IRc7, and a radio source I (Source I). Spectral energy distribution (SED) of each source is determined by using previous 3 mm continuum emission data. Physical properties such as size, mass, hydrogen number density and column density are discussed based on the dust graybody SED. Among 11 identified sources, Source I, a massive protostar candidate, is a dominant energy source in Orion KL. We extensively investigate its SED from centimeter to submillimeter wavelengths. The SED of Source I can be fitted with a single power-law index of 1.97 suggesting an optically thick emission. We employ the H$^{-}$ free-free emission as an opacity source of this optically thick emission. The temperature, density, and mass of the circumstellar disk associated with Source I are constrained by the SED of H$^{-}$ free-free emission. Still the fitting result shows a significant deviation from the observed flux densities. Combined with the thermal dust graybody SED to explain excess emission at higher frequency, a smaller power-law index of 1.60 for the H$^{-}$ free-free emission is obtained in the SED fitting. The power-law index smaller than 2 would suggest a compact source size or a clumpy structure unresolved with the present study. Future higher resolution observations with ALMA are essential to reveal more detailed spatial structure and physical properties of Source I.

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.

The 22 GHz H2O maser in Orion KL has shown extraordinary burst events in 1979-1985 and 1998-1999, sometimes called supermaser. We have conducted monitoring observations of the supermaser in Orion KL using VERA, VLBI Exploration of Radio Astrometry, in the current third burst since 2011 March. Three flux maxima are detected in 2011 and 2012 with rising and falling timescales of 2-7 months. Time variations of the supermaser seem symmetric for all of the active phases. The maximum total flux density of 135000 Jy is observed in 2012 June while it is still one order of magnitude lower than those in previous bursts. The supermaser consists of two spatially different components at different velocities. They are elongated along a northwest-southeast direction perpendicular to the low-velocity outflow driven by Source I. Proper motions of the supermaser features with respect to Source I are measured toward west and southwest directions, almost parallel to the low-velocity outflow. The flux density and linewidth show an anti-correlation as expected for an unsaturated maser emission. The supermaser is located close to the methylformate (HCOOCH3) line and continuum emission peaks in the Orion Compact Ridge detected by ALMA. The broader velocity range of the weak HCOOCH3 emission at the supermaser position would be an evidence of a shock front. On the other hand, the 321 GHz H2O line is not detected at the position of the supermaser. It can be explained qualitatively by one of the theoretical H2O excitation models without extraordinary conditions. Our results support a scenario that the supermaser is excited in the dense gas interacting with the low-velocity outflow in the Compact Ridge. The extremely high flux density and its symmetric time variation for rising and falling phases could be explained by a beaming effect during the amplification process rather than changes in physical conditions.

We propose a new technique to obtain super-resolution images with radio interferometer using sparse modeling. In standard radio interferometry, sampling of ($u$, $v$) is quite often incomplete and thus obtaining an image from observed visibilities becomes an underdetermined problem, and a technique so-called "zero-padding" is often used to fill up unsampled grids in ($u$, $v$) plane, resulting in image degradation by finite beam size as well as numerous side-lobes. In this paper we show that directly solving such an underdetermined problem based on sparse modeling (in this paper LASSO) avoids the above problems introduced by zero-padding, leading to super-resolution images in which structure finer than the standard beam size (diffraction limit) can be reproduced. We present results of one-dimensional and two-dimensional simulations of interferometric imaging, and discuss its implications to super-resolution imaging, particularly focusing on imaging of black hole shadows with millimeter VLBI.

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 carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We report the results of the measurement of the trigonometric parallax of an H$_2$O maser source in \r22 with the VLBI Exploration of Radio Astrometry (VERA). The annual parallax was determined to be 0.314$\pm$0.070~mas, corresponding to a distance of 3.18$^{+0.90}_{-0.66}$~kpc. Our results confirm \r22 to be located in the Perseus arm. We computed the peculiar motion of \r22 to be $(U_{src}, V_{src}, W_{src}) = (4\pm4, -32\pm6, 8\pm6)$~\kms, where $U_{src}$, $V_{src}$, and $W_{src}$ are directed toward the Galactic center, in the direction of Galactic rotation and toward the Galactic north pole, respectively. \r22, NGC\,7538 and Cepheus A lie along the same line of sight, and are within 2$^{\circ}$ on the sky. Their parallax distances with which we derived their absolute position in the Milky Way show that \r22 and \ngc~are associated with the Perseus Arm, while \cepa~is located in the Local Arm. We compared the kinematic distances of \iras, \ga, \gb, \gc, \gd, \cepa, \ngc, \ge, \gf~and \r22~derived with flat and non-flat rotation curve with its parallax distance and found the kinematic distance derived from the non-flat rotation assumption ($-$5 to $-$39 \kms~lag) to be consistent with the parallax distance.

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 conducted phase referencing VLBI observations of the Mira variable T~Lepus (T~Lep) using VERA, from 2003 to 2006. The distance to the source was determined from its annual parallax which was measured to be 3.06$\pm$0.04 mas, corresponding to a distance of 327$\pm$4\,pc. Our observations revealed the distribution and internal kinematics of H$_2$O masers in T~Lep, and we derived a source systemic motions of 14.60$\pm$0.50 mas yr$^{-1}$ and $-$35.43$\pm$0.79 mas yr$^{-1}$ in right ascension and declination, respectively. We also determined a LSR velocity of $V_\mathrm{LSR}^{\ast} = -$27.63 km s$^{-1}$. Comparison of our result with an image recently obtained from the VLTI infrared interferometer reveals a linear scale picture of the circumstellar structure of T~Lep. Analysis of the source systemic motion in the Galacto-centric coordinate frame indicates a large peculiar motion, which is consistent with the general characteristics of AGB stars. This source makes a contribution to the calibration of the period-luminosity relation of Galactic Mira variables. From the compilation of data for nearby Mira variables found in the literature, whose distances were derived from astrometric VLBI observations, we have calibrated the Galactic Mira period-luminosity relation to a high degree of accuracy.

Multi-epoch phase-referencing VLBI (Very Long Baseline Interferometry) observations with VERA (VLBI Exploration of Radio Astrometry) were performed for the symbiotic star R Aquarii (R Aqr) from September 2005 to Oct 2006. Tracing one of the $v=2$, $J=1-0$ SiO maser spots, we measured an annual parallax of $\pi = 4.59\pm0.24$ mas, corresponding to a distance of $218_{-11}^{+12}$ pc. Our result is consistent with earlier distance measurements, but yields the highest accuracy of about $5\%$ level. Applying our distance, we derived an absolute K-band magnitude of $M_{\mathrm{K}} = -7.71 \pm 0.11$, which is consistent with the recent Period-Luminosity relation by VLBI parallax measurements for 5 OH-Mira variables. In addition, the expansion age of an inner nebulae around R Aqr is found to be about 240 years, corresponds to about the year 1773.

We report new Atacama Large Millimeter/Submillimeter Array (ALMA) observations of a circumstellar disk around Source I in Orion KL, an archetype of massive protostar candidate. We detected two ortho-H$_{2}$O lines at 321 GHz ($10_{2,9}$-$9_{3,6}$) and 336 GHz ($\nu_{2}=1, 5_{2,3}$-$6_{1,6}$) for the first time in Source I. The latter one is in vibrationally excited state at the lower state energy of 2939 K, suggesting an evidence of hot molecular gas close to Source I. The integrated intensity map of the 321 GHz line is elongated along the bipolar outflow while the 336 GHz line map is unresolved with a beam size of 0.4\arcsec. Both of these maps show velocity gradient perpendicular to the bipolar outflow. The velocity centroid map of the 321 GHz line implies spatial and velocity structure similar to that of vibrationally-excited SiO masers tracing the root of the outflow emanating from the disk surface. In contrast, the 336 GHz line is most likely emitted from the disk midplane with a diameter of 0.2\arcsec (84 AU) as traced by a radio continuum emission and a dark lane devoid of the vibrationally-excited SiO maser emission. The observed velocity gradient and the spectral profile of the 336 GHz H$_{2}$O line can be reconciled with a model of an edge-on ring-like structure with an enclosed mass of $>$7$M_{\odot}$ and an excitation temperature of $>$3000 K. The present results provide a further evidence of hot and neutral circumstellar disk rotating around Source I with a diameter of $\sim$100 AU scale.

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.

Astrometry provides the foundation for astrophysics. Accurate positions are required for the association of sources detected at different times or wavelengths, and distances are essential to estimate the size, luminosity, mass, and ages of most objects. Very Long Baseline Interferometry at radio wavelengths, with diffraction-limited imaging at sub-milliarcsec resolution, has long held the promise of micro-arcsecond astrometry. However, only in the past decade has this been routinely achieved. Currently, parallaxes for sources across the Milky Way are being measured with ~10 uas accuracy and proper motions of galaxies are being determined with accuracies of ~1 uas/y. The astrophysical applications of these measurements cover many fields, including star formation, evolved stars, stellar and super-massive black holes, Galactic structure, the history and fate of the Local Group, the Hubble constant, and tests of general relativity. This review summarizes the methods used and the astrophysical applications of micro-arcsecond radio astrometry.

We present the results of radio monitoring observations of Sgr A* at 7 mm (i.e. 43 GHz) with VLBI Exploration of Radio Astrometry (VERA), which is a VLBI array in Japan. VERA provides angular resolutions on millisecond scales, resolving structure within ~100 Schwarzschild radii of Sgr A* similar to Very Large Baseline Array (VLBA). We performed multi-epoch observations of Sgr A* in 2005 - 2008, and started monitoring it again with VERA from January 2013 for tracing the current G2 encounter event. Our preliminary results in 2013 show that Sgr A* on mas scales has been in ordinary state as of August 2013, although some fraction of the G2 cloud already passed pericenter of Sgr A* in April 2013. We will continue on monitoring Sgr A* with VERA and newly developed KaVA (KVN and VERA Array).

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.

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

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

We report the results of multi-epoch observations of Sgr A* with VLBI Exploration of Radio Astrometry (VERA) at 43 GHz, carried out from 2004 to 2008. We detected a time variation of flux at 11 % level and intrinsic size at 19 %. In addition, comparisons with previous Very Long Baseline Array (VLBA) results shows that Sgr A* underwent the flaring event at least longer than 10 days in May 2007. The intrinsic size of Sgr A* remained unchanged within 1 $\sigma$ level from the size before/after the flaring event, indicating that the brightness temperature of Sgr A* was increased. The flaring event occurred within 31 d, which is shorter than the refractive time scale. Moreover it is difficult to explain the increase in the spectral index at the flaring event by the simple interstellar scattering model. Hence, the flaring event is most likely associated with the changes in intrinsic properties of Sgr A*. We considered the origin of the brightness temperature variation, and concluded that the flaring event of Sgr A* could be explained by the continuous heating of electrons, such as a standing shock in accretion flow.

We present the first measurement of the absolute proper motions of IRAS 00259+5625 (CB3, LBN594) associated with the HI loop called the "NGC281 superbubble" that extends from the Galactic plane over ~300 pc toward decreasing galactic latitude. The proper motion components measured with VERA are (mu_alpha cos(delta), mu_delta) = (-2.48 +/- 0.32, -2.85 +/- 0.65) mas yr^{-1}, converted into (mu_l cos(b), mu_b) = (-2.72 +/- 0.32, -2.62 +/- 0.65) mas yr^{-1} in the Galactic coordinates. The measured proper motion perpendicular to the Galactic plane (mu_b) shows vertical motion away from the Galactic plane with a significance of about ~4-sigma. As for the source distance, the distance measured with VERA is marginal, 2.4^{+1.0}_{-0.6} kpc. Using the distance, an absolute vertical motion (v_{b}) of -17.9 +/- 12.2 km s^{-1} is determined with ~1.5-sigma significance. The tendency of the large vertical motion is consistent with previous VLBI results for NGC 281 associated with the same superbubble. Thus, our VLBI results indicate the superbubble expansion motion whose origin is believed to be sequential supernova explosions.

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

We present the first internal motion measurement of the 6.7-GHz methanol maser within S269, a small HII region in the outer Galaxy, which was carried out in 2006 and 2011 using the Japanese VLBI Network (JVN). Several maser groups and weak isolated spots were detected in an area spanning by ~200 mas (1000 AU). Three remarkable maser groups are aligned at a position angle of 80 degree. Two of three maser groups were also detected by a previous observation in 1998, which allowed us to study a long-term position variation of maser spots from 1998 to 2011. The angular separation between the two groups increased ~10 mas, which corresponds to an expansion velocity of ~10 km s^{-1}. Some velocity gradient (~10^{-2} km s^{-1} mas^{-1}) in the overall distribution was found. The internal motion between the maser groups support the hypothesis that the methanol masers in S269 could trace a bipolar outflow.

We present a radio observation of microquasar Cyg X-3 during an X-ray state transition from ultrasoft to hard state in the 2007 May$-$June flare using the VLBI Exploration of Radio Astrometry (VERA) at 22 GHz. During the transition, a short-lived mini-flare of ${\lesssim}3$ hr was detected prior to the major flare. In such a transition, a jet ejection is believed to occur, but there have been no direct observations to support it. An analysis of Gaussian fits to the observed visibility amplitudes shows a time variation of the source axis, or a structural change, during the mini-flare. Our model fits, together with other multiwavelength observations in the radio, soft, and hard X-rays, and the shock-in-jet models for other flaring activities at GHz wavebands, suggest a high possibility of synchrotron flares during the mini-flare, indicative of a predominant contribution from jet activity. Therefore, the mini-flare with an associated structural change is indicative of a jet ejection event in the state transition from ultrasoft to hard state.

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 report results from 5-day VLBI observations of the well-known quasar 3C 279 at 1.3 mm (230 GHz) in 2011. The measured nonzero closure phases on triangles including stations in Arizona, California and Hawaii indicate that the source structure is spatially resolved. We find an unusual inner jet direction at scales of $\sim$1 parsec extending along the northwest-southeast direction (PA = $127^{\circ}\pm3^{\circ}$), as opposed to other (previously) reported measurements on scales of a few parsecs showing inner jet direction extending to the southwest. The 1.3 mm structure corresponds closely with that observed in the central region of quasi-simultaneous super-resolution VLBA images at 7 mm. The closure phase changed significantly on the last day when compared with the rest of observations, indicating that the inner jet structure may be variable on daily timescales. The observed new direction of the inner jet shows inconsistency with the prediction of a class of jet precession models. Our observations indicate a brightness temperature of $\sim 8\times10^{10}$ K in the 1.3 mm core, much lower than that at centimeter wavelengths. Observations with better uv coverage and sensitivity in the coming years will allow the discrimination between different structure models and will provide direct images of the inner regions of the jet with 20--30 $\mu$as (5--7 light months) resolution.

We present Very Long Baseline Interferometry (VLBI) observations of 22 GHz H$_2$O masers in the high-mass star-forming region of \objectname{W75N}, carried out with VLBI Exploration of Radio Astrometry (VERA) for three-epochs in 2007 with an angular resolution of $\sim$ 1 mas. We detected H$_2$O maser emission toward the radio jet in VLA 1 and the expanding shell-like structure in VLA 2. We have made elliptical fits to the VLA 2 H$_2$O maser shell-like structure observed in the different epochs (1999, 2005, and 2007), and found that the shell is still expanding eight years after its discovery. From the difference in the size of the semi-major axes of the fitted ellipses in the epochs 1999 ($\simeq$ 71$\pm$1 mas), 2005 ($\simeq$ 97$\pm$3 mas), and 2007 ($\simeq$ 111$\pm$1 mas), we estimate an average expanding velocity of $\sim$ 5 mas yr$^{-1}$, similar to the proper motions measured in the individual H$_2$O maser features. A kinematic age of $\sim$ 20 yr is derived for this structure. In addition, our VERA observations indicate an increase in the ellipticity of the expanding shell around VLA 2 from epochs 1999 to 2007. In fact, the elliptical fit of the VERA data shows a ratio between the minor and major axes of $\sim$ 0.6, in contrast with a most circular shape for the shell detected in 1999 and 2005 (b/a $\sim$ 0.9). This suggests that we are probably observing the formation of a jet-driven H$_2$O maser structure in VLA2, evolving from a non-collimated pulsed-outflow event during the first stages of evolution of a massive young stellar object (YSO). This may support predictions made earlier by other authors on this issue, consistent with recent magnetohydrodynamical simulations. We discuss possible implications of our results in the study of the first stages of evolution of massive YSOs.

We report on results of astrometric observations of water vapor masers in the "water fountain" source IRAS 18286-0959 (I18286) with the VLBI Exploration of Radio Astrometry (VERA). These observations yielded an annual parallax of IRAS 18286-0959, pi=0.277+/-0.041 mas, corresponding to a heliocentric distance of D=3.61(+0.63)(-0.47) kpc. The maser feature, whose annual parallax was measured, showed the absolute proper motion of (mu_alpha, mu_delta)=(-3.2 +/- 0.3, -7.2 +/- 0.2) [mas/yr]. The intrinsic motion of the maser feature in the internal motions of the cluster of features in I18286 does not seem to trace the motion of the bipolar jet of I18286. Taking into account this intrinsic motion, the derived motion of the maser feature is roughly equal to that of the maser source I18286 itself. The proximity of I18286 to the Galactic midplane (z~10 pc) suggests that the parental star of the water fountain source in I18286 should be intermediate-mass AGB/post-AGB star, but the origin of a large deviation of the systemic source motion from that expected from the Galactic rotation curve is still unclear.

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 present a measurement of the trigonometric parallax of IRAS 05168+3634 with VERA. The parallax is 0.532 +/- 0.053 mas, corresponding to a distance of 1.88 +0.21/-0.17 kpc. This is significantly closer than the previous distance estimate of 6 kpc based on a kinematic distance measurement. This drastic change in the source distance implies the need for revised values of not only the physical parameters of IRAS 05168+3634, but it also implies a different location in the Galaxy, placing it in the Perseus arm rather than the Outer arm. We also measured the proper motion of the source. A combination of the distance and proper motion with the systemic velocity yields a rotation velocity {\Theta} = 227 +9/-11 km s^-1 at the source position, assuming {\Theta}_0 = 240 km s^-1. Our result, combined with previous VLBI results for six sources in the Perseus arm, indicates that the sources rotate systematically more slowly than the Galactic rotation velocity at the local standard of rest. In fact, we derive peculiar motions in the disk averaged over the seven sources in the Perseus arm of (U_mean, V_mean) = (11 +/- 3, -17 +/- 3) km s^-1, which indicates that these seven sources are moving systematically toward the Galactic Center and lag behind the overall Galactic rotation.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

Approximately 10% of active galactic nuclei exhibit relativistic jets, which are powered by accretion of matter onto super massive black holes. While the measured width profiles of such jets on large scales agree with theories of magnetic collimation, predicted structure on accretion disk scales at the jet launch point has not been detected. We report radio interferometry observations at 1.3mm wavelength of the elliptical galaxy M87 that spatially resolve the base of the jet in this source. The derived size of 5.5 +/- 0.4 Schwarzschild radii is significantly smaller than the innermost edge of a retrograde accretion disk, suggesting that the M87 jet is powered by an accretion disk in a prograde orbit around a spinning black hole.

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

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

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

We report the 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 present the astrometric VLBI observations of water masers associated with IRAS 22480+6002 (IRC+60370, hereafter I22480) with the VLBI Exploration of Radio Astrometry (VERA). The stellar type of I22480 looks unusual as a stellar maser source and has been debated since the 1970s. We successfully determined the annual parallax of a group of the water maser spots, pi=0.400+/-0.025 mas, corresponding to a distance to I22480 of D=2.50 /+0.17/-0.15 kpc. This suggests that the estimated bolometric luminosity of I22480 should be revised to 35 000 L_sun, favoring a K-type supergiant rather than an RV Tau-type variable star previously suggested. Although the spectral type is unusual as a stellar maser source, the internal motions of the water maser features suggest that the water masers are associated with the circumstellar envelope of this star. Taking into account a possible stellar motion with respect to the maser feature motions, we derived a secular proper motion of I22480, (mu_alpha, mu_delta)=(-2.58+/-0.33, -1.91+/-0.17) [mas yr^-1]. The derived motion of I22480 in the Milky Way has a deviation by ~ -30 km/s in the galactic azimuthal direction from a circular motion estimated from the galactocentric distance to I22480 and assumption of a flat Galactic rotation curve. This peculiar motion is still comparable to those typically seen in the water maser sources located in the Perseus spiral arm. Taking into account the peculiar motion and the proximity to the Galactic midplane (z ~60 pc), I22480 may be a member of the Galactic thin disk.

We investigated the ALMA science verification data of Orion KL and found a spectral signature of the vibrationally excited H2O maser line at 232.68670 GHz (nu2=1, 5,5,0-6,4,3). This line has been detected in circumstellar envelopes of late-type stars so far but not in young stellar objects including Orion KL. Thus, this is the first detection of the 232 GHz vibrationally excited H2O maser in star-forming regions. The distribution of the 232 GHz maser is concentrated at the position of the radio Source I, which is remarkably different from other molecular lines. The spectrum shows a double-peak structure at the peak velocities of -2.1 and 13.3 km s-1. It appears to be consistent with the 22 GHz H2O masers and 43 GHz SiO masers observed around Source I. Thus, the 232 GHz H2O maser around Source I would be excited by the internal heating by an embedded protostar, being associated with either the root of the outflows/jets or the circumstellar disk around Source I, as traced by the 22 GHz H2O masers or 43 GHz SiO masers, respectively.

We present the first 1.3 mm (230 GHz) very long baseline interferometry model image of an AGN jet using closure phase techniques with a four-element array. The model image of the quasar 1924-292 was obtained with four telescopes at three observatories: the James Clerk Maxwell Telescope (JCMT) on Mauna Kea in Hawaii, the Arizona Radio Observatory's Submillimeter Telescope (SMT) in Arizona, and two telescopes of the Combined Array for Research in Millimeterwave Astronomy (CARMA) in California in April 2009. With the greatly improved resolution compared with previous observations and robust closure phase measurement, the inner jet structure of 1924-292 was spatially resolved. The inner jet extends to the northwest along a position angle of $-53^\circ$ at a distance of 0.38\,mas from the tentatively identified core, in agreement with the inner jet structure inferred from lower frequencies, and making a position angle difference of $\sim 80^{\circ}$ with respect to the cm-jet. The size of the compact core is 0.15\,pc with a brightness temperature of $1.2\times10^{11}$\,K. Compared with those measured at lower frequencies, the low brightness temperature may argue in favor of the decelerating jet model or particle-cascade models. The successful measurement of closure phase paves the way for imaging and time resolving Sgr A* and nearby AGN with the Event Horizon Telescope.

We report measurement of trigonometric parallax of IRAS 05168+3634 with VERA. The parallax is 0.532 +/- 0.053 mas, corresponding to a distance of 1.88+0.21/-0.17 kpc. This result is significantly smaller than the previous distance estimate of 6 kpc based on kinematic distance. This drastic change in the source distance revises not only physical parameters of IRAS 05168+3634, but also its location of the source, placing it in the Perseus arm rather than the Outer arm. We also measure proper motions of the source. A combination of the distance and the proper motions with systemic velocity yields rotation velocity ({\Theta}) of 227+9/-11 km s-1 at the source, assuming {\Theta}0 = 240 km s-1. Our result combined with previous VLBI results for six sources in the Perseus arm indicates that the sources rotate systematically slower than the Galactic rotation velocity at the LSR. In fact, we show observed disk peculiar motions averaged over the seven sources in the Perseus arm as (Umean, Vmean) = (11 +/- 3, -17 +/- 3) km s-1, indicating that these seven sources are systematically moving toward the Galactic center, and lag behind the Galactic rotation.

We present distance measurement of the semiregular variable RX Bootis (RX Boo) with its annual parallax. Using the unique dual-beam system of the VLBI Exploration of Radio Astrometry (VERA) telescope, we conducted astrometric observations of a water maser spot accompanying RX Boo referred to the quasar J1419+2706 separated by 1.69 degrees from RX Boo. We have measured the annual parallax of RX Boo to be 7.31 +/- 0.50 mas, corresponding to a distance of 136 +10/-9 pc, from the one-year monitoring observation data of one maser spot at VLSR = 3.2 km/s. The distance itself is consistent with the one obtained with Hipparcos. The distance uncertainty is reduced by a factor of two, allowing us to determine the stellar properties more accurately. Using our distance, we discuss the location of RX Boo in various sequences of Period-Luminosity (PL) relations. We found RX Boo is located in the Mira sequence of PL relation. In addition, we calculated the radius of photosphere and the mass limitation of RX Boo and discussed its evolutionary status.

Whether the Cygnus X complex consists of one physically connected region of star formation or of multiple independent regions projected close together on the sky has been debated for decades. The main reason for this puzzling scenario is the lack of trustworthy distance measurements. We aim to understand the structure and dynamics of the star-forming regions toward Cygnus X by accurate distance and proper motion measurements. To measure trigonometric parallaxes, we observed 6.7 GHz methanol and 22 GHz water masers with the European VLBI Network and the Very Long Baseline Array. We measured the trigonometric parallaxes and proper motions of five massive star-forming regions toward the Cygnus X complex and report the following distances within a 10% accuracy: 1.30+-0.07 kpc for W 75N, 1.46^{+0.09}_{-0.08} kpc for DR 20, 1.50^{+0.08}_{-0.07} kpc for DR 21, 1.36^{+0.12}_{-0.11} kpc for IRAS20290+4052, and 3.33+-0.11kpc for AFGL 2591. While the distances of W 75N, DR 20, DR 21, and IRAS 20290+4052 are consistent with a single distance of 1.40+-0.08 kpc for the Cygnus X complex, AFGL 2591 is located at a much greater distance than previously assumed. The space velocities of the four star-forming regions in the Cygnus X complex do not suggest an expanding Stroemgren sphere.

We observed at 22 GHz with the VLBI array VERA a sample of 1536 sources with correlated flux densities brighter than 200 mJy at 8 GHz. One half of target sources has been detected. The detection limit was around 200 mJy. We derived the correlated flux densities of 877 detected sources in three ranges of projected baseline lengths. The objective of these observations was to determine the suitability of given sources as phase calibrators for dual-beam and phase-referencing observations at high frequencies. Preliminary results indicate that the number of compact extragalactic sources at 22 GHz brighter than a given correlated flux density level is twice less than at 8 GHz.

Charged-current neutrino-nucleus cross sections for 54,56Fe and 58,60Ni are calculated and compared using frameworks based on relativistic and Skyrme energy density functionals, and the shell model. The current theoretical uncertainties in modeling neutrino-nucleus cross sections are assessed in relation to the predicted Gamow-Teller transition strength and available data, multipole decomposition of the cross sections, and cross sections averaged over the Michel flux and Fermi-Dirac distribution. Employing different microscopic approaches and models, the DAR neutrino-56Fe cross section and its theoretical uncertainty are estimated: <sigma>_th=(258+-57) 10^{-42} cm^2, in very good agreement with the experimental value: <sigma>_exp=(256+-108+-43) 10^{-42} cm^2.

We present the results of astrometic observations of H2O masers associated with the star forming region G192.16-3.84 with the VLBI Exploration of Radio Astrometry (VERA). The H2O masers seem to be associated with two young stellar objects (YSOs) separated by \sim1200 AU as reported in previous observations. In the present observations, we successfully detected an annual parallax of 0.66 \pm 0.04 mas for the H2 O masers, which corresponds to a distance to G192.16-3.84 of D = 1.52 \pm 0.08 kpc from the Sun. The determined distance is shorter than the estimated kinematic distance. Using the annual parallax distance and the estimated parameters of the millimeter continuum emission, we estimate the mass of the disk plus circumstellar cloud in the southern young stellar object to be 10.0+4.3M\cdot. We also estimate the galactocentric distance and the peculiar motion -3.6 of G192.16-3.84, relative to a circular Galactic rotation: R\star = 9.99 \pm 0.08 kpc, Z\star = -0.10 \pm 0.01 kpc, and (U\star,V\star,W\star)=(-2.8\pm1.0,-10.5\pm0.3,4.9\pm2.7)[kms-1]respectively. The peculiar motion of G192.16-3.84 is within that typically found in recent VLBI astrometric results. The angular distribution and three-dimensional velocity field of H2O maser features associated with the northern YSO indicate the existence of a bipolar outflow with a major axis along the northeast-southwest direction.

We report the detection of bursts of 6.7 GHz methanol maser emission in a high-mass star-forming region, G33.64-0.21. One of the spectral components of the maser in this source changed its flux density by 7 times that of the previous day, and it decayed with a timescale of 5 days. The burst occurred repeatedly in the spectral component, and no other components showed such variability. A VLBI observation with the Japanese VLBI Network (JVN) showed that the burst location was at the southwest edge of a spatial distribution, and the bursting phenomenon occurred in a region much smaller than 70 AU. We suggest an impulsive energy release like a stellar flare as a possible mechanism for the burst. These results imply that 6.7 GHz methanol masers could be a useful new probe for studying bursting activity in the process of star formation of high-mass YSOs with a high-resolution of AU scale.

In February 2011, a burst event of the H$_{2}$O maser in Orion KL (Kleinmann-Low object) has started after 13-year silence. This is the third time to detect such phenomena in Orion KL, followed by those in 1979-1985 and 1998. We have carried out astrometric observations of the bursting H$_{2}$O maser features in Orion KL with VERA (VLBI Exploration of Radio Astrometry), a Japanese VLBI network dedicated for astrometry. The total flux of the bursting feature at the LSR velocity of 7.58 km s$^{-1}$ reaches 4.4$\times10^{4}$ Jy in March 2011. The intensity of the bursting feature is three orders of magnitudes larger than that of the same velocity feature in the quiescent phase in 2006. Two months later, another new feature appears at the LSR velocity of 6.95 km s$^{-1}$ in May 2011, separated by 12 mas north of the 7.58 km s$^{-1}$ feature. Thus, the current burst occurs at two spatially different features. The bursting masers are elongated along the northwest-southeast direction as reported in the previous burst in 1998. We determine the absolute positions of the bursting features for the first time ever with a submilli-arcsecond (mas) accuracy. Their positions are coincident with the shocked molecular gas called the Orion Compact Ridge. We tentatively detect the absolute proper motions of the bursting features toward southwest direction. It is most likely that the outflow from the radio source I or another young stellar object interacting with the Compact Ridge is a possible origin of the H$_{2}$O maser burst.

We have performed VLBI and single-dish monitoring of 22 GHz H$_{2}$O maser emission from the high mass young stellar object G353.273+0.641 with VERA (VLBI Exploration of Radio Astrometry) and Tomakamai 11-m radio telescope. Two maser flares have been detected, separated almost two years. Frequent VLBI monitoring has revealed that these flare activities have been accompanied by structural change of the prominent shock front traced by H2O maser alignments. We have detected only blue-shifted emissions and all maser features have been distributed within very small area of 200 $\times$ 200 au$^{2}$ in spite of wide velocity range (> 100 km s$^{-1}$). The light curve shows notably intermittent variation and suggests that the H$_{2}$O masers in G353.273+0.641 are excited by episodic radio jet. The time-scale of \sim2 yr and characteristic velocity of \sim500 km s$^{-1}$ also support this interpretation. Two isolated velocity components of C50 (-53 \pm 7 km s$^{-1}$) and C70 (-73 \pm 7 km s$^{-1}$) have shown synchronised linear acceleration of the flux weighted $V_{\rmn{LSR}}$ values (\sim-5 km s$^{-1}$ yr$^{-1}$) during the flare phase. This can be converted to the lower-limit momentum rate of 1.1 \times 10$^{-3}$ M$_{\sun}$ km s$^{-1}$ yr$^{-1}$. Maser properties are quite similar to that of IRAS 20126+4104 especially. This corroborates the previous suggestion that G353.273+0.641 is a candidate of high mass protostellar object. The possible pole-on geometry of disc-jet system can be suitable for direct imaging of the accretion disc in this case.

We performed the astrometry of H2O masers in the Galactic star-forming region G48.61+0.02 with the VLBI Exploration of Radio Astrometry (VERA). We derived a trigonometric parallax of 199+/-7 micro as, which corresponds to a distance of 5.03+/-0.19 kpc. The distance to G48.61+0.02 is about a half of its far kinematic distance, which was often assumed previously. This distance places G48.61+0.02 in the Sagittarius-Carina arm and near the active star forming region and the supernova remnant W51. We also obtained the three dimensional motion of G48.61+0.02, and found that it has a large peculiar motion of 40+/-5 km/s. This peculiar motion would be originated with the multiple supernovae explosions in W51, or the streaming motion across the Sagittarius-Carina arm.

We have measured the annual parallax of the H2O maser source associated with an infrared dark cloud MSXDC G034.43+00.24 from the observations with VERA (VLBI Exploration of Radio Astrometry). The parallax is 0.643 +/- 0.049 mas, corresponding to the distance of 1.56 +0.12/-0.11 kpc. This value is less than the half of the previous kinematic distance of 3.7 kpc. We revise the core mass estimates of MSXDC G034.43+00.24, based on virial masses, LTE masses and dust masses and show that the core masses decrease from the previous estimations of ~1000 Mo to hundreds of Mo. The spectral type derived from the luminosity also changes from O9.5 to B1 in the case of MM1. This spectral type is still consistent with that of the massive star. The radial velocity derived from the flat rotation model is smaller than the observed velocity, which corresponds to the peculiar motion of ~40 km/s in the line-of-sight direction.

We present the results of astrometric observations with VERA toward the H2O maser sources in IRAS 05137+3919, which is thought to be located in the far outer Galaxy. We have derived the parallax of \pi = 0.086 +/- 0.027 mas, which corresponds to the source distance of D=11.6+5.3-2.8 kpc. Although the parallax measurement is only 3-sigma level and thus the distance uncertainty is considerably large, we can strongly constrain the minimum distance to this source, locating the source at the distance from the Sun greater than 8.3 kpc (or 16.7 kpc from the Galaxy's center) at 90% confidence level. Our results provide an astrometric confirmation that this source is located in the far outer Galaxy beyond 15 kpc from the Galaxy center, indicating that IRAS 05137+3919 is one of the most distant star-forming regions from the Galaxy center.

In this paper we present the results of very long baseline interferometry (VLBI) ob- servations carried out with the VLBI Exploration of Radio Astrometry (VERA) array and the Very Long Baseline Array (VLBA) toward H2O masers in a young planetary nebula K 3-35. From the VERA observations we measured the annual parallax and proper mo- tion of a bright water maser spot in K 3-35. The resulting distance is D = 3.9+0.7 kpc. -0.5 This is the first time that the parallax of a planetary nebula is obtained by observations of its maser emission. On the other hand, the proper motion of K 3-35 as a whole was esti- mated to be {\mu}{\alpha} = -3.34+/-0.10 mas yr-1, {\mu}{\delta} = -5.93+/-0.07 mas yr-1. From these results we determined the position and velocity of K 3-35 in Galactic cylindrical coordinates: (R,{\theta},z) = (7.11+0.08-0.06 kpc, 27+/-5{\circ}, 140+25-18 pc) and (VR, V{\theta}, Vz) = (33+/-16, 233+/-11, 11+/-2) km s-1, respectively. Additionally, from our VLBA observations we measured the relative proper motions among the water maser spots located in the central region of the nebula, which have been proposed to be tracing a toroidal structure. The distribution and relative proper motions of the masers, compared with previous reported observed epochs, suggest that such structure could be totally destroyed within a few years, due to the action of high velocity winds and the expansion of the ionization front in the nebula.

We conducted the astrometry of H2O masers in the Galactic star-forming region Onsala 1 (ON1) with VLBI Exploration of Radio Astrometry (VERA). We measured a trigonometric parallax of 0.404+/-0.017 mas, corresponding to a distance of 2.47+/-0.11 kpc. ON1 is appeared to be located near the tangent point at the Galactic longitude of 69.54 deg. We estimate the angular velocity of the Galactic rotation at Sun, the ratio of the distance from Sun to the Galactic center and the Galactic rotation velocity at Sun, to be 28.7+/-1.3 km/s/kpc using the measured distance and proper motion of ON1. This value is larger than the IAU recommended value of 25.9 km/s/kpc, but consistent with other results recently obtained with the VLBI technique.

We have performed the astrometry of H2O masers in the Galactic star-forming region Onsala 2 North (ON2N) with the VLBI Exploration of Radio Astrometry.We obtained the trigonometric parallax of 0.261+/-0.009 mas, corresponding to the heliocentric distance of 3.83+/-0.13 kpc. ON2N is expected to be on the solar circle, because its radial velocity with respect to the Local Standard of Rest (LSR) is nearly zero. Using present parallax and proper motions of the masers, the Galactocentric distance of the Sun and the Galactic rotation velocity at the Sun are found to be R_0 = 7.80+/-0.26 kpc and Theta_0 = 213+/-5 km/s, respectively. The ratio of Galactic constants, namely the angular rotation velocity of the LSR can be determined more precisely, and is found to be Omega_0=Theta_0/R_0 = 27.3+/-0.8 km/s/kpc, which is consistent with the recent estimations but different from 25.9 km/s/kpc derived from the recommended values of Theta_0 and R_0 by the International Astronomical Union (1985).

We report on the astrometric observations of the 22 GHz H2O masers in the high mass star-forming region G5.89-0.39 with VERA (VLBI Exploration of Radio Astrometry). Newly derived distance of 1.28^{+0.09}_{-0.08} kpc is the most precise and significantly nearer than previous values. We revised physical parameters and reconsidered nature of G5.89-0.39 based on the new distance as follows. (1) The ionizing star of the ultra compact (UC) HII region is a late O-type (O8 - 8.5) zero age main sequence (ZAMS) star, consistent with previously established limits based on its infrared spectral line emission. (2) Crescent-like maser alignment at the position of the O type ZAMS star may trace accretion disk (or its remnant), which suggests that the star is still young and before complete evaporation of circumstellar materials. (3) Although the revised mass for the east-west outflow has been reduced, it still quite large (100 Msun) which indicates that a significant fraction of the mass is entrained material and that the dynamical age significantly underestimates the actual outflow age. Our newly-derived distance emphasizes that G5.89-0.39 is one of the nearest targets to investigate ongoing high-mass star formation and evolution in a compact cluster containing a young O-type star.

Four close radio sources in the International Celestial Reference Frame (ICRF) catalog were observed using phase referencing with the VLBA at 43, 23 and 8.6 GHz, and with VERA at 23 GHz over a one year period. The goal was to determine the stability of the radio cores, and to assess structure effects associated with positions in the ICRF. Although the four sources were compact at 8.6 GHz, the VLBA images at 43 GHz with 0.3-mas resolution showed that all were composed of several components. A component in each source was identified as the radio core using some or all of the following emission properties: compactness, spectral index, location at the end of the extended emission region, and stationary in the sky. Over the observing period, the relative positions between the four radio cores were constant to 0.02 mas---the phase referencing positional accuracy obtained at 23 and 43 GHz among the sources---suggesting that once a radio core is identified, it remains stationary in the sky to this accuracy. Other radio components in two of the four sources had detectable motion in the radio jet direction. Comparison of the 23 and 43 GHz VLBA images with the VLBA 8.6 GHz images and the ICRF positions suggests that some ICRF positions are dominated by a moving jet component; hence, they can be displaced up to 0.5 mas from the radio core, and may also reflect the motion of the jet component. Future astrometric efforts to determine a more accurate quasar reference frame at 23 and 43 GHz and from the VLBI2010 project are discussed, and supporting VLBA or EVN observations of ICRF sources at 43 GHz are recommended in order to determine the internal structure of the sources. A future collaboration between the radio (ICRF) and the optical frame of Gaia is discussed.

Cygnus X is one of the closest giant molecular cloud complexes and therefore an extensively studied region of ongoing high mass star formation. However, the distance to this region has been a long-standing issue, since sources at galactic longitude of ~80 degrees could be in the Local Arm nearby (1-2 kpc), in the Perseus Arm at ~5 kpc, or even in the outer arm (~10 kpc). We use combined observations of the EVN plus two Japanese stations to measure very accurate parallaxes of methanol masers in five star-forming regions in Cygnus X to understand if they belong to one large star-forming complex or if they are separate entities located at different distances. Here we report our preliminary result for W75N based on six epochs of VLBI observations: we find that W75N is at a distance of 1.32^{+0.11}_{-0.09} kpc, which is significantly closer than the reported values in the literature (1.5-2 kpc).

Sagittarius A*, the ~4 x 10^6 solar mass black hole candidate at the Galactic Center, can be studied on Schwarzschild radius scales with (sub)millimeter wavelength Very Long Baseline Interferometry (VLBI). We report on 1.3 mm wavelength observations of Sgr A* using a VLBI array consisting of the JCMT on Mauna Kea, the ARO/SMT on Mt. Graham in Arizona, and two telescopes of the CARMA array at Cedar Flat in California. Both Sgr A* and the quasar calibrator 1924-292 were observed over three consecutive nights, and both sources were clearly detected on all baselines. For the first time, we are able to extract 1.3 mm VLBI interferometer phase information on Sgr A* through measurement of closure phase on the triangle of baselines. On the third night of observing, the correlated flux density of Sgr A* on all VLBI baselines increased relative to the first two nights, providing strong evidence for time-variable change on scales of a few Schwarzschild radii. These results suggest that future VLBI observations with greater sensitivity and additional baselines will play a valuable role in determining the structure of emission near the event horizon of Sgr A*.

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.

We have conducted VLBA phase-referencing monitoring of H2O masers around the red supergiant, S Persei, for six years. We have fitted maser motions to a simple expanding-shell model with a common annual parallax and stellar proper motion, and obtained the annual parallax as 0.413+/-0.017 mas, and the stellar proper motion as (-0.49+/-0.23 mas/yr, -1.19+/-0.20 mas/yr) in right ascension and declination, respectively. The obtained annual parallax corresponds to the trigonometric distance of 2.42+0.11-0.09 kpc. Assuming the Galactocentric distance of the Sun of 8.5 kpc, the circular rotational velocity of the LSR at the distance of the Sun of 220 km/s, and a flat Galactic rotation curve, S Persei is suggested to have a non-circular motion deviating from the Galactic circular rotation for 15 km/s, which is mainly dominated by the anti rotation direction component of 12.9+/-2.9 km/s. This red supergiant is thought to belong to the OB association, Per OB1, so that this non-circular motion is representative of a motion of the OB association in the Milky Way. This non-circular motion is somewhat larger than that explained by the standard density-wave theory for a spiral galaxy, and is attributed to either a cluster shuffling of the OB association, or to non-linear interactions between non-stationary spiral arms and multi-phase interstellar media. The latter comes from a new view of a spiral arm formation in the Milky Way suggested by recent large N-body/smoothed particle hydrodynamics numerical simulations.

We have carried out multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22~GHz H$_{2}$O masers associated with a Class 0 protostar L1448C in the Perseus molecular cloud. The maser features trace the base of collimated bipolar jet driven by one of the infrared counter parts of L1448C named as L1448C(N) or L1448-mm A. We detected possible evidences for apparent acceleration and precession of the jet according to the three-dimensional velocity structure. Based on the phase-referencing VLBI astrometry, we have successfully detected an annual parallax of the H$_{2}$O maser in L1448C to be 4.31$\pm$0.33~milliarcseconds (mas) which corresponds to a distance of 232$\pm$18~pc from the Sun. The present result is in good agreement with that of another H$_{2}$O maser source NGC~1333 SVS13A in the Perseus molecular cloud, 235~pc. It is also consistent with the photometric distance, 220~pc. Thus, the distance to the western part of the Perseus molecular cloud complex would be constrained to be about 235~pc rather than the larger value, 300~pc, previously reported.

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 report on trigonometric parallax measurements for the Galactic star forming region G14.33-0.64 toward the Sagittarius spiral arm. We conducted multi-epoch phase-referencing observations of an H2O maser source in G14.33-0.64 with the Japanese VLBI array VERA. We successfully detected a parallax of 0.893+/-0.101 mas, corresponding to a source distance of 1.12+/-0.13 kpc, which is less than half of the kinematic distance for G14.33-0.64. Our new distance measurement demonstrates that the Sagittarius arm lies at a closer distance of ~1 kpc, instead of previously assumed ~2-3 kpc from kinematic distances. The previously suggested deviation of the Sagittarius arm toward the Galactic center from the symmetrically fitted model (Taylor & Cordes 1993) is likely due to large errors of kinematic distances at low galactic longitudes. G14.33-0.64 most likely traces the near side of the Sagittarius arm. We attempted fitting the pitch angle of the arm with other parallax measurements along the arm, which yielded two possible pitch angles of i=34.7+/-2.7 degrees and i=11.2+/-10.5 degrees. Our proper motion measurements suggest G14.33-0.64 has no significant peculiar motion relative to the differential rotation of the Galaxy (assumed to be in a circular orbit), indicating that the source motion is in good agreement with the Galactic rotation.

We report on recent VLBA/VERA/IVS observational tests of General Relativity. First, we will summarize the results from the 2005 VLBA experiment that determined gamma with an accuracy of 0.0003 by measuring the deflection of four compact radio sources by the solar gravitational field. We discuss the limits of precision that can be obtained with VLBA experiments in the future. We describe recent experiments using the three global arrays to measure the aberration of gravity when Jupiter and Saturn passed within a few arcmin of bright radio sources. These reductions are still in progress, but the anticipated positional accuracy of the VLBA experiment may be about 0.01 mas.

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

We are using the VLBA and the Japanese VERA project to measure trigonometric parallaxes and proper motions of masers found in high-mass star-forming regions across the Milky Way. Early results from 18 sources locate several spiral arms. The Perseus spiral arm has a pitch angle of 16 +/- 3 degrees, which favors four rather than two spiral arms for the Galaxy. Combining positions, distances, proper motions, and radial velocities yields complete 3-dimensional kinematic information. We find that star forming regions on average are orbiting the Galaxy ~15 km/s slower than expected for circular orbits. By fitting the measurements to a model of the Galaxy, we estimate the distance to the Galactic center R_o = 8.4 +/- 0.6 kpc and a circular rotation speed Theta_o = 254 +/- 16 km/s. The ratio Theta_o/R_o can be determined to higher accuracy than either parameter individually, and we find it to be 30.3 +/- 0.9 km/s/kpc, in good agreement with the angular rotation rate determined from the proper motion of Sgr A*. The data favor a rotation curve for the Galaxy that is nearly flat or slightly rising with Galactocentric distance. Kinematic distances are generally too large, sometimes by factors greater than two; they can be brought into better agreement with the trigonometric parallaxes by increasing Theta_o/R_o from the IAU recommended value of 25.9 km/s/kpc to a value near 30 km/s/kpc. We offer a "revised" prescription for calculating kinematic distances and their uncertainties, as well as a new approach for defining Galactic coordinates. Finally, our estimates of Theta_o and To/R_o, when coupled with direct estimates of R_o, provide evidence that the rotation curve of the Milky Way is similar to that of the Andromeda galaxy, suggesting that the dark matter halos of these two dominant Local Group galaxy are comparably massive.

We present a unified rotation curve of the Galaxy re-constructed from the existing data by re-calculating the distances and velocities for a set of galactic constants R_0=8 kpc and V_0=200 km/s. We decompose it into a bulge with de Vaucouleurs-law profile of half-mass scale radius 0.5 kpc and mass 1.8 x 10^{10}M_{sun}, an exponential disk of scale radius 3.5 kpc of 6.5 x 10^{10}M_{sun}, and an isothermal dark halo of terminal velocity 200 km/s. The r^{1/4}-law fit was obtained for the first time for the Milky Way's rotation curve. After fitting by these fundamental structures, two local minima, or the dips, of rotation velocity are prominent at radii 3 and 9 kpc. The 3-kpc dip is consistent with the observed bar. It is alternatively explained by a massive ring with the density maximum at radius 4 kpc. The 9-kpc dip is clearly exhibited as the most peculiar feature in the galactic rotation curve. We explain it by a massive ring of amplitude as large as 0.3 to 0.4 times the disk density with the density peak at radius 11 kpc. This great ring may be related to the Perseus arm, while no peculiar feature of HI-gas is associated.

We have made multi-epoch VLBI observations of H2O maser emission in the massive star forming region IRAS 06061+2151 with the Japanese VLBI network (JVN) from 2005 May to 2007 October. The detected maser features are distributed within an 1\arcsec$\times$1\arcsec (2000 au$\times$2000 au at the source position) around the ultra-compact H {\small\bf II} region seen in radio continuum emission. Their bipolar morphology and expanding motion traced through their relative proper motions indicate that they are excited by an energetic bipolar outflow. Our three-dimensional model fitting has shown that the maser kinematical structure in IRAS 06061+2151 is able to be explained by a biconical outflow with a large opening angle ($>$ 50\degr). The position angle of the flow major axis coincides very well with that of the large scale jet seen in 2.1$\:\mu\rmn{m}$ hydrogen emission. This maser geometry indicates the existence of dual structures composed of a collimated jet and a less collimated massive molecular flow. We have also detected a large velocity gradient in the southern maser group. This can be explained by a very small (on a scale of several tens of au) and clumpy (the density contrast by an order of magnitude or more) structure of the parental cloud. Such a structure may be formed by strong instability of shock front or splitting of high density core.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We present the results of daily monitoring of 6.7 GHz methanol maser in Cepheus A (Cep A) using Yamaguchi 32-m radio telescope as well as the results of imaging observations conducted with the JVN (Japanese VLBI Network). We indentified five spectral features, which are grouped into red-shifted (-1.9 and -2.7 km/s) and blue-shifted (-3.8, -4.2, and -4.9 km/s), and we detected rapid variabilities of these maser features within a monitoring period of 81 days. The red-shifted features decreased in flux density to 50% of its initial value, while the flux density of the blue-shifted features rapidly increased within a 30 days. The time variation of these maser features showed two remarkable properties; synchronization and anti-correlation between the red-shifted and the blue-shifted. The spatial distribution of the maser spots obtained by the JVN observation showed an arclike structure with a scale of $\sim$1400 AU, and separations of the five maser features were found to be larger than 100 AU. The absolute position of the methanol maser was also obtained based on the phase-referencing observations, and the arclike structure were found to be associated with the Cep A-HW2 object, with the elongation of the arclike structure nearly perpendicularly to the radio continuum jet from the Cep A-HW2 object. These properties of the masers, namely, the synchronization of flux variation, and the spectral and spatial isolation of features, suggest that the collisional excitation by shock wave from a common exciting source is unlikely. Instead, the synchronized time variation of the masers can be explained if all the maser features are excited by infrared radiation from dust which is heated by a common exciting source with a rapid variability.

We present measurements with the VLBA of the variability in the centroid position of Sgr A* relative to a background quasar at 7-mm wavelength. We find an average centroid wander of 71 +/- 45 micro-arcsec for time scales between 50 and 100 min and 113 +/- 50 micro-arcsec for timescales between 100 and 200 min, with no secular trend. These are sufficient to begin constraining the viability of the hot-spot model for the radio variability of Sgr A*. It is possible to rule out hot spots with orbital radii above 15GM_SgrA*/c^2 that contribute more than 30% of the total 7-mm flux. However, closer or less luminous hot spots remain unconstrained. Since the fractional variability of Sgr A* during our observations was ~20% on time scales of hours, the hot-spot model for Sgr A*'s radio variability remains consistent with these limits. Improved monitoring of Sgr A*'s centroid position has the potential to place significant constraints upon the existence and morphology of inhomogeneities in a supermassive black hole accretion flow.

We have observed 13 methanol maser sources associated with massive star-forming regions; W3(OH), Mon R2, S 255, W 33A, IRAS 18151-1208, G 24.78+0.08, G 29.95-0.02, IRAS 18556+0136, W 48, OH 43.8-0.1, ON 1, Cep A and NGC 7538 at 6.7 GHz using the Japanese VLBI Network (JVN). Twelve of the thirteen sources were detected at our longest baseline of $\sim$50 M$\lambda$, and their images are presented. Seven of them are the first VLBI images at 6.7 GHz. This high detection rate and the small fringe spacing of $\sim$4 milli-arcsecond suggest that most of the methanol maser sources have compact structure. Given this compactness as well as the known properties of long-life and small internal-motion, this methanol maser line is suitable for astrometry with VLBI.

We report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the maser features in SVS 13 of 4.25+/-0.32 mas, corresponding to the distance of 235+/-18 pc from the Sun. Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser features were detectable only for 6 months, the present results provide the distance to NGC 1333 with much higher accuracy than photometric methods. The absolute positions and proper motions have been derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s-1 are most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed proper motions of the maser features are currently difficult to attribute to either the jet or the rotating circumstellar disk associated with VLA4A, which should be investigated through future high-resolution astrometric observations of VLA4A and other radio sources in NGC 1333.

We have performed high-precision astrometry of H2O maser sources in Galactic star forming region Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189+/-8 micro-arcsec, corresponding to the source distance of 5.28 +0.24/-0.22 kpc. This is the smallest parallax ever measured, and the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3%. This gives the strongest constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence of large amount of dark matter in the Galaxy's outer disk.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

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

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

We present the basic features and the activities of Japanese VLBI network (JVN), a newly-established VLBI network with baselines ranging from 50 to 2560 km spreading across the Japanese islands, and capable of observing at 6.7, 8.4, and 22 GHz. We show a number of results of JVN observations: 8.4-GHz continuum images of a Giga-hertz Peaked Spectrum (GPS) source and radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s), the spatial and velocity structures of water masers in NML Cygni as well as methanol masers in Cep A, and demonstrative observations with the bigradient phase referencing.

This paper presents results of a survey search for bright compact radio sources at 22 GHz with the VERA radio-interferometer. Each source from a list of 2494 objects was observed in one scan for 2 minutes. The purpose of this survey was to find compact extragalactic sources bright enough at 22 GHz to be useful as phase calibrators. Observed sources were either a) within 6 degrees of the Galactic plane, or b) within 11 degrees from the Galactic center; or c) within 2 degrees from known water masers. Among the observed sources, 549 were detected, including 180 extragalactic objects which were not previously observed with the very long baseline interferometry technique. Estimates of the correlated flux densities of the detected sources are presented.

We propose bigradient phase referencing (BPR), a new radio-observation technique, and report its performance using the Japanese very-long-baseline-interferometry network (JVN). In this method, a weak source is detected by phase-referencing using a primary calibrator, in order to play a role as a secondary calibrator for phase-referencing to a weak target. We will be given the opportunity to select a calibrator from lots of milli-Jansky sources, one of which may be located at the position closer to the target. With such a smaller separation, high-quality phase-referencing can be achieved. Furthermore, a subsequent more-sophisticated calibration can relocate array's focus to a hypothetical point much closer to the target; a higher quality of phase referencing is available. Our demonstrative observations with strong radio sources have proved the capabilities of BPR in terms of image dynamic ranges and astrometric reproducibility. The image dynamic range on a target has been improved with a factor of about six compared to that of normal phase-referencing; the resultant position difference of target's emission between two epochs was only 62+-50 micro-arcsecond, even with less than 2300-km baselines at 8.4 GHz and fast-switching of a target-calibrator pair of a 2.1-degree separation.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We present the results of multi-epoch VERA observations of W49N - OH43.8-0.1 water maser pair. Based on the dual-beam VLBI observation with VERA, we successfully obtained the phase-referenced maps of OH43.8-0.1 with respect to the W49N reference spot for 3 epochs with a time span of 6 months. The maps were in good agreement with previous studies obtained with a single-beam VLBI, and were also consistent with each other with an accuracy of about 0.2 mas. Moreover, there are systematic, rather linear displacements of maser feature positions, which may be the relative proper motions of maser features caused by the Galactic rotation as well as internal motions of individual maser features.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

We present the results of first dual-beam observations with VERA (VLBI Exploration of Radio Astrometry). The first dual-beam observations of a pair of H2O maser sources W49N and OH43.8-0.1 have been carried out on 2002 May 29 and July 23, and fringes of H2O maser lines at 22 GHz have been successfully detected. While the residual fringe phases of both sources show rapid variations over 360 degree due to the atmospheric fluctuation, the differential phase between the two sources remains constant for 1 hour with r.m.s. of 8 degree, demonstrating that the atmospheric phase fluctuation is removed effectively by the dual-beam phase referencing. The analysis based on Allan standard deviation reveals that the differential phase is mostly dominated by white phase noise, and the coherence function calculated from the differential phase shows that after phase referencing the fringe visibility can be integrated for arbitrarily long time. These results demonstrate VERA's high capability of phase referencing, indicating that VERA is a promising tool for phase referencing VLBI astrometry at 10 microarcsec level accuracy.

We analyze the data of the gravitational microlensing survey carried out by by the MOA group during 2000 towards the Galactic Bulge (GB). Our observations are designed to detect efficiently high magnification events with faint source stars and short timescale events, by increasing the the sampling rate up to 6 times per night and using Difference Image Analysis (DIA). We detect 28 microlensing candidates in 12 GB fields corresponding to 16 deg^2. We use Monte Carlo simulations to estimate our microlensing event detection efficiency, where we construct the I-band extinction map of our GB fields in order to find dereddened magnitudes. We find a systematic bias and large uncertainty in the measured value of the timescale $t_{\rm Eout}$ in our simulations. They are associated with blending and unresolved sources, and are allowed for in our measurements. We compute an optical depth tau = 2.59_{-0.64}^{+0.84} \times 10^{-6} towards the GB for events with timescales 0.3<t_E<200 days. We consider disk-disk lensing, and obtain an optical depth tau_{bulge} = 3.36_{-0.81}^{+1.11} \times 10^{-6}[0.77/(1-f_{disk})] for the bulge component assuming a 23% stellar contribution from disk stars. These observed optical depths are consistent with previous measurements by the MACHO and OGLE groups, and still higher than those predicted by existing Galactic models. We present the timescale distribution of the observed events, and find there are no significant short events of a few days, in spite of our high detection efficiency for short timescale events down to t_E = 0.3 days. We find that half of all our detected events have high magnification (>10). These events are useful for studies of extra-solar planets.

We investigate properties of astrometric microlensing of distant sources (such as QSOs and radio galaxies) caused by stars in the Galaxy, mainly focusing on application to the VERA (VLBI Exploration of Radio Astrometry) project. Assuming typical parameters for the Galaxy disk and bulge, we show that the maximum optical depth for astrometric shift of 10 micro-arcsec level is 8.9x10^{-2} for QSO-disk lensing case and 3.8x10^{-2} for QSO-bulge lensing case. We also find that the maximum optical depth for QSO-disk lensing is larger by an order of magnitude than that for disk-disk or bulge-disk lensing case (assuming a typical source distance of 8 to 10 kpc). In addition to optical depth, we also calculate the event rate and find that the maximum event rate for QSO-disk lensing case is 1.2x10^{-2} event per year, which is about 30 times greater than that for disk-disk lensing case. This high event rate implies that if one monitors 10 QSOs behind the Galactic center region for 10 years, at least one astrometric microlensing event should be detected. Therefore, monitoring distant radio sources with VERA can be a new tool to study astrometric microlensing caused by stars in the Galaxy. We discuss the implications of astrometric microlensing for VERA by focusing on estimating the lens mass, and also present some possible candidates of radio sources toward which astrometric microlensing events should be searched for with VERA. (abridged)

One hundred and forty six long-period red variable stars in the Large Magellanic Cloud (LMC) from the three year MOA project database were analysed. A careful periodic analysis was performed on these stars and a catalogue of their magnitudes, colours, periods and amplitudes is presented. We convert our blue and red magnitudes to $K$ band values using 19 oxygen-rich stars. A group of red short-period stars separated from the Mira sequence has been found on a (log P, K) diagram. They are located at the short period side of the Mira sequence consistent with the work of Wood and Sebo (1996). There are two interpretations for such stars; a difference in pulsation mode or a difference in chemical composition. We investigated the properties of these stars together with their colour, amplitude and periodicity. We conclude that they have small amplitudes and less regular variability. They are likely to be higher mode pulsators. A large scatter has been also found on the long period side of the (log P, K) diagram. This is possibly a systematic spread given that the blue band of our photometric system covers both standard B and V bands and affects carbon-rich stars.

We have determined distances and Hubble ratios for galaxies at intermediate redshifts, cz=about 10,000 to 35,000 km/s, by applying the CO-line Tully-Fisher relation to our 12CO (J=1-0)-line observations using the Nobeyama 45-m telescope, and near-IR (NIR) photometry in the J and H bands using the 1.88-m telescope at Okayama Astrophysical Observatory. By averaging the Hubble ratios from J-band result, we obtain a Hubble constant of H_0 = 60 (+/- 10) km/s/Mpc. We argue that the CO line-NIR TFR can be a complimentary method to the other methods for measuring distances of galaxies at intermediate and high redshifts.

In this paper we investigate the properties of astrometric microlensing of distant radio sources (QSOs and radio galaxies) due to MACHOs, and discuss their implications for VERA (VLBI Exploration of Radio Astrometry). First we show that in case of astrometric microlensing of distant sources, the event duration is only a function of the lens mass and tangential velocity, but independent of the lens distance, in contrast to the well-known three-fold degeneracy for photometric microlensing. Moreover, the lens mass $M$ is scaled by the tangential velocity $v_\perp$ as $M\propto v_\perp$, rather than $M\propto v_\perp^2$ which is the case for photometric microlensing. Thus, in astrometric microlensing the dependence of the lens mass on the unknown parameter $v_\perp$ is weaker, indicating that the duration of astrometric microlensing event is a better quantity to study the mass of lensing objects than that of photometric microlensing. We also calculate the optical depth and event rate, and show that within 20$^\circ$ of the galactic center a typical event rate for 10 $\mu$as-level shift is larger than $2.5 \times 10^{-4}$ event per year, assuming that a quarter of the halo is made up with MACHOs. This indicates that if one monitors a few hundred sources for $\sim$20 years, such an astrometric microlensing event is detectable. Since a typical event duration is found to be fairly long (5 to 15 years), the frequency of the monitoring observation can be relatively low, i.e., once per six months, which is rather reasonable for practical observations. We discuss practical strategy for observing astrometric microlensing with VERA, and argue that an astrometric microlensing event due to MACHOs can be detected by VERA within a few decades.

We have performed high-resolution synthesis observations of the 12CO (J=1-0) line emission from the radio lobe edge-on spiral NGC 3079 using a 7-element mm-wave interferometer at Nobeyama, which consisted of the 45-m telescope and 6-element Array. The molecular nuclear disk (NMD) of 750 pc radius is found to be inclined by 20 degrees from the optical disk, and the NMD has spiral arms. An ultra-high-density molecular core (UHC) was found at the nucleus. The gaseous mass of the UHC within 125 pc radius is as large as about 3 x 10^8 Msun, an order of magnitude more massive than that in the same area of the Galactic Center, and the mean density is as high as approximately 3 x 10^3 H2 cm^-3. A position-velocity diagram along the major axis indicates that the rotation curve starts at a finite velocity exceeding 300 \kms already from the nucleus. The surface mass density in the central region is estimated to be as high as about 10^5 Msun pc^-2, producing a very deep gravitational potential. We argue that the very large differential rotation in such a deep potential will keep the UHC gravitationally stable from current star formation.

We present new measurement of $^{12}$CO($J=1-0$) emission from 16 luminous infrared galaxies (LIGs) at intermediate redshift ($cz \sim 10,000 - 50,000 {\rm km s^{-1}}$). These new data were selected by isolated and normal morphology. The CO observations were performed using the NRO 45-m telescope. Comparison of the CO and dust properties of the new result with those from other CO measurements revealed characteristics of this sample: (1) It is the deepest CO observations of IRAS galaxies at intermediate redshift without strong interaction features. (2) It has typical properties of normal IRAS galaxies in terms of star-formation efficiency, color-color diagrams and galactic nuclear activity. (3) It has smaller gas-to-dust ratio than normal IRAS galaxies. This can be explained by two-component dust model, and our sample consists of most of warm dust.

In order to search for new VLBI sources in the Galactic plane that can be used as phase reference sources in differential VLBI, we have conducted 22 GHz observations of radio sources in the Galactic plane using the Japanese VLBI Network (J-Net). We have observed 267 VLBI source candidates selected from existing radio surveys and have detected 93 sources at the signal-to-noise ratio larger than 5. While 42 of the 93 detected sources had already been detected with VLBI at relatively lower frequency (typically 2 to 8 GHz), the remaining 51 are found to be new VLBI sources detected for the first time. These VLBI sources are located within $|b|\le 5^\circ$, and have a large number of Galactic maser sources around them. Thus, they are potential candidates for phase reference sources for VLBI Exploration of Radio Astrometry (VERA), which is the first VLBI array dedicated to the phase referencing VLBI astrometry aiming at measuring the parallax and proper motion of maser sources in the whole Galaxy.

We propose a new method to constrain the location of microlensing objects using EAGLE (Extremely Amplified Gravitational LEnsing) events. We have estimated the rate of EAGLE events by taking the finite-source effect in to account. We found that the EAGLE event rate for using a 1-m class telescope w hose limiting magnitude is $V \sim 21$ is the same as or higher than that of the ordinary microlensing events which have been found to date. We have also found that the fraction of transit EAGLE events is large enough to detect: between $4 \sim 80 % $ depending on the lens location. Since the lens proper motion can be measured for a transit event, one can distinguish whether the lens is a MACHO (MAssive Compact Halo Object) in our hal o or one of the known stars in the Large Magellanic Cloud (LMC) from the proper motion measurement for each transit EAGLE event. Moreover, we show that the fraction of transit EAGLEs in all EAGLE events signif icantly depends on the lensing locations: the transit EAGLE fraction for the sel f-lensing case is $2 \sim 15$ times larger than that for halo MACHOs. Thus, one can constrain the location of lens objects by the statistics of the tr ansit events fraction. We show that we can reasonably expect $0 \sim 6$ transit events out of 21 EAGLE events in 3 years. We can also constrain the lens population properties at a gre ater than 99% confidence level depending on the number of transit events de tected. We also present the duration of EAGLE events, and show how an hourly ob servational mode is more suitable for an EAGLE event search program.

We have obtained interferometer observations of the central region of the Sb galaxy NGC 4527 in the 12CO (J=1-0) line emission using the Nobeyama Millimeter Array. We also obtained optical (Halpha, [NII]) spectroscopy using the Okayama 188-cm reflector along the major axis. The kinematical structure and the distribution of HII regions show symmetry around the nucleus, while the distribution of molecular gas is asymmetric. The molecular-gas mass shares only 10% of the dynamical mass in the central 1 kpc radius region. Using position-velocity diagrams, we have derived a center-to-disk rotation curve. It rises steeply in the central region, attaining a maximum of about 250 km/s at 400 pc radius, and then decreases to a minimum at 2 kpc radius, followed by a disk and outer flat part. The rotation curve may provide us with the most similar case to that of the Milky Way Galaxy.

We present high-resolution central-to-outer rotation curves for Sb, SBb, Sc, and SBc galaxies. We discuss their general characteristics, particularly their central behavior, as well as dependencies on morphological types, activity, and peculiarity. The rotation curves generally show a steep nuclear rise, and high-velocity central rotation, followed by a broad maximum in the disk and then by a flat rotation due to the massive halo. Since the central high velocity and steep rise are common to all massive galaxies, it cannot be due to non-circular motions. Disk rotation curves of barred galaxies show larger dispersion than those of normal galaxies which is probably due to non-circular motions. Interacting galaxies show often perturbed outer rotation curves, while their central rotation shows no particular peculiarity. Also, central activities, such as starbursts and AGN, appear to show no particular correlation with the property of rotation curves. This would suggest that the central activities are triggered by a more local effect than the global dynamical property.

The cluster Zwicky 1615.8+3505 is considered to be a dynamically young poor cluster. To investigate the morphology and star-formation activity of galaxies under the environment of a dynamically young poor cluster, we have performed V, R, and I surface photometry, optical low-resolution spectroscopy, and 12CO (J=1-0) line observations for member galaxies. Our data show that more than 90% of the observed galaxies show regular morphologies and no star-formation activities, indicating that the environment does not affect these galaxy properties. Among sixteen galaxies observed, only NGC 6104 shows a significant star-formation activity, and shows a distorted morphology, indicating a tidal interaction. This galaxy contains double knots, and only one knot possesses Seyfert activity, though the sizes and luminosities are similar to each other; we also discuss this feature.

We report on the mass-to-light ratio determination based on a newly selected binary galaxy sample, which includes a large number of pairs whose separations exceed a few hundred kpc. The probability distributions of the projected separation and the velocity difference have been calculated considering the contamination of optical pairs, and the mass-to-light ratio has been determined based on the maximum likelihood method. The best estimate of $M/L$ in the B band for 57 pairs is found to be 28 $\sim$ 36 depending on the orbital parameters and the distribution of optical pairs (solar unit, $H_0=50$ km s$^{-1}$ Mpc$^{-1}$). The best estimate of $M/L$ for 30 pure spiral pairs is found to be 12 $\sim$ 16. These results are relatively smaller than those obtained in previous studies, but consistent with each other within the errors. Although the number of pairs with large separation is significantly increased compared to previous samples, $M/L$ does not show any tendency of increase, but found to be almost independent of the separation of pairs beyond 100 kpc. The constancy of $M/L$ beyond 100 kpc may indicate that the typical halo size of spiral galaxies is less than $\sim 100$ kpc.

We investigate the possibility of lens mass determination for a caustic crossing microlensing event based on a space telescope observation. We demonstrate that the parallax due to the orbital motion of a space telescope causes a periodic fluctuation of the light curve, from which the lens distance can be derived. Since the proper motion of the lens relative to the source is also measurable for a caustic crossing event, one can find a full solution for microlensing properties of the event, including the lens mass. To determine the lens mass with sufficient accuracy, the light curve near the caustic crossing should be observed within uncertainty of $\sim$ 1%. We argue that the Hubble Space Telescope observation of the caustic crossing supplied with ground-based observations of the full light curve will enable us to determine the mass of MACHOs, which is crucial for understanding the nature of MACHOs.

Caustic crossing microlensing events provide us a unique opportunity to measure the relative proper motion of the lens to the source, and so those caused by binary MACHOs are of great importance for understanding the structure of the Galactic halo and the nature of MACHOs. The microlensing event 98-SMC-01, occurred in June 1998, is the first event for which the proper motion is ever measured through the caustic crossing, and this event may be caused by binary MACHOs as we argue in this Letter. Motivated by the possible existence of binary MACHOs, we have performed the Monte Carlo simulations of caustic crossing events by binary MACHOs and investigated the properties and detectability of the events. Our calculation shows that typical caustic crossing events have the interval between two caustic crossings ($t_{\rm cc}$) of about 5 days. We argue that with the current strategy of binary event search the proper motions of these typical events are not measurable because of the short time scale. Therefore the proper motion distribution measured from caustic crossing events suffers significantly from {`}time scale bias{'}, which is a bias toward finding long time scale events and hence slowly moving lenses. We predict there are two times more short time scale events ($t_{\rm cc}\le 10$ days) than long time scale events ($t_{\rm cc}\ge 10$ days), and propose an hourly monitoring observation instead of the nightly monitoring currently undertaken to detect caustic crossing events by binary MACHOs more efficiently.

We have obtained optical CCD spectroscopy along the major axes of 22 nearby spiral galaxies of Sb and Sc types in order to analyze their rotation curves. By subtracting the stellar continuum emission, we have obtained position velocity (PV) diagrams of the H alpha and [NII] lines. We point out that the H alphaline is often superposed by a broad stellar absorption feature (Balmer wind) in the nuclear regions, and, therefore, the [NII] line is a better tracer of kinematics in the central a few hundred pc regions. By applying the envelope-tracing technique to the H alpha and [NII] PV diagrams, we have derived nucleus-to-disk rotation curves of the observed galaxies. The rotation curves rise steeply within the central a few hundred parsecs, indicating rapidly rotating nuclear disk and mass concentration near the nucleus. [For more rotation curves : http://www.ioa.s.u-tokyo.ac.jp/]

If the Galactic rotation speed at the Solar circle is $\sim 200$ km s$^{-1}$ or smaller, which is supported by several recent studies, the rotation curve of the Galaxy could be declining in the outermost region. Motivated by this, we investigate the effect of such declining rotation curve on the estimate of the MACHO mass and the fractional contribution of the MACHOs to the Galactic dark halo. Using Hernquist and Plummer halo models instead of the standard halo model, we find that the MACHO mass could be significantly smaller than that for the standard halo case. In particular, there exists a certain set of halo parameters for which the MACHO mass is 0.1$M_\odot$ or less and at the same time the MACHO contribution to the total mass of the halo is almost 100 %. This result indicates that a halo which consists solely of brown dwarfs can be consistent with both of the observed microlensing properties and the constraints from the rotation curve, provided the outer rotation curve is indeed declining.

We propose a criterion for examining whether or not the uncertainty of the outer rotation curves is sufficiently small to distinguish a Keplerian rotation curve from flat ones. We have applied this criterion to both Keplerian and non-Keplerian rotation curves so far obtained, and investigated their relative fraction. We also studied the minimum extent of the dark halos indicated by non-Keplerian rotation curves using the criterion. We have found that one cannot rule out the possibility that a significant fraction of rotation curves become Keplerian within 10-times the disk scale length. If the Keplerian rotation curves so far observed trace the mass truncation, several galaxies may have rather small halos, the extent of which is not larger than twice that of the optical disk.

We have obtained high-resolution position-velocity (PV) diagrams along the major axes of the central regions of nearby galaxies in the CO-line emission using the Nobeyama 45-m telescope and the Millimeter Array. Nuclear rotation curves for 14 galaxies have been derived based on the PV diagrams using the envelope-tracing method. The nuclear rotation curves for most of the galaxies show a steep rise within a few hundred pc, which indicates a high-density concentration of mass. Keywords: Galaxies: general - Galaxies: structure - ISM: molecular line

Using an outer rotation curve of the Galaxy, we explore the galactic constants and the mass of the Galaxy. We show that \Theta_0 of 200 km/s is more favorable than the IAU standard value of 220 km/s, and also show that if \Theta_0 is smaller than 207 km/s the rotation curve beyond 2R_0 is declining in Keplerian fashion. In the case of \Theta_0= 200 km/s and R_0= 7.6 kpc, the total mass and the extent of the Galaxy inferred from the rotation curve are 2.0+/-0.3x 10^{11} M_\odot and 15 kpc, respectively. These results may significantly change the previous view of the Galaxy, that its outer region is dominated by a massive dark halo extending out to several tens of kpc.

We report on the first results of a long-term project to derive distances of galaxies at cosmological distances by applying the CO-line width-luminosity relation. We have obtained deep CO-line observations of galaxies at redshifts up to 29,000 km/s using the Nobeyama 45-m mm-wave Telescope, and some supplementary data were obtained by using the IRAM 30-m telescope. We have detected the CO line emission for several galaxies, and used their CO line widths to estimate the absolute luminosities using the line-width-luminosity relation. In order to obtain photometric data and inclination correction, we also performed optical imaging observations of the CO-detected galaxies using the CFHT 3.6-m telescope at high resolution. The radio and optical data have been combined to derive the distance moduli and distances of the galaxies, and Hubble ratios were estimated for these galaxies. We propose that the CO line width-luminosity relation can be a powerful method to derive distances of galaxies to redfhift of z = 0.1 and to derive the Hubble ratio in a significant volume of the universe. Key words: Cosmology - Galaxies: general - Distance scale - CO line

We have examined the distribution of HI and H$_2$ gases in four face-on galaxies by using the observed dat a of CO and HI line emissions from the literatures. We demonstrate that the gas phase transition of HI and H$_2$ occurs suddenly within a narrow range of radi us, which we call the molecular front. We have tried to explain such phase transition in galactic scale with a help of the phase transition theory proposed by Elmegreen. The crucial parameters for determinating the molecular fraction $f_{\rm mol}$ are interstellar pressure $P $, UV radiation field $U$, and metallicity $Z$, and we have constructed a model galaxy in which $P$, $U$ a nd $Z$ obey an exponential function of the galacto-centric radius. The model shows that the molecular front must be a fundamental feature of galaxies which has an exponentia l disk, and that the metallicity gradient is most crucial for the formation of the molecular front. We have also tried to reproduce the observed molecular fraction $f_{\rm mol}$ by giving the set of ($P$, $ U$, $Z$) observationally, and show that the model can describe the variation of the molecular fraction $f_ {\rm mol}$ in galaxies quite well. We discuss the implication of the molecular front for the chemical evolution of galaxies.