We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247.

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

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

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 followed-up the massive young stellar object (MYSO) S255-NIRS3 (=S255-IRS1b) during its recent accretion outburst event in the Ks band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete NIR light-curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4 mag brighter than the quiescent phase and then a relatively moderate year-scale fading until the last observation, similar to that of the accretion burst events such as EXors found in lower-mass young stellar objects. The behavior of the Ks band light curve is similar to that observed in 6.7 GHz class II methanol maser emission, with a sudden increase followed by moderate year-scale fading. However, the maser emission peaks appear 30-50 days earlier than that of the Ks band emission. The similarities confirmed that the origins of the maser emission and the Ks band continuum emission is common as previously shown from another infrared and radio observations by Stecklum et al. (2016); Caratti o Garatti et al. (2017a); Moscadelli et al. (2017). However, the differences in energy transfer paths, such as the exciting/emitting/scattering structures, may cause the delay in the flux-peak dates.

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

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.

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

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

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

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.

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

We 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 present the results of multi-epoch VLBI observations of water masers in the AGFL 5142 massive star forming region. We measure an annual parallax of $\pi=0.467 \pm 0.010$ mas, corresponding to a source distance of $D=2.14^{+0.051}_{-0.049}$ kpc. Proper motion and line of sight velocities reveal the 3D kinematics of masers in this region, most of which associate with millimeter sources from the literature. In particular we find remarkable bipolar bowshocks expanding from the most massive member, AFGL 5142 MM1, which are used to investigate the physical properties of its protostellar jet. We attempt to link the known outflows in this region to possible progenitors by considering a precessing jet scenario and we discuss the episodic nature of ejections in AFGL 5142.

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.

Context. The S235AB star forming region houses a massive young stellar object which has recently been reported to exhibit possible evidence of jet rotation - an illusive yet crucial component of disk aided star formation theories. Aims. To confirm the presence of a molecular counterpart to the jet and to further study the molecular environment in in S235AB. Methods. We search for velocity wings in the line emission of thermal SiO (J=2-1, v=0), a tracer of shocked gas, which would indicate the presence of jet activity. Utilising other lines detected in our survey we use the relative intensities of intra species transitions, isotopes and hyperfine transitions to derive opacities, temperatures, column densities and abundances of various molecular species in S235AB. Results. The SiO (J=2-1, v=0) emission exhibits velocity wing of up to 75 km/s above and below the velocity of the star, indicating the presence of a jet. The molecular environment describes an evolutionary stage resemblant of a hot molecular core.

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

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

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

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

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

The 6.7 GHz methanol maser emission, a tracer of forming massive stars, sometimes shows enigmatic periodic flux variations over several 10-100 days. In this Letter, we propose that this periodic variations could be explained by the pulsation of massive protostars growing under rapid mass accretion with rates of Mdot > 10^-3 Msun/yr. Our stellar evolution calculations predict that the massive protostars have very large radius exceeding 100 Rsun at maximum, and we here study the pulsational stability of such the bloated protostars by way of the linear stability analysis. We show that the protostar becomes pulsationally unstable with various periods of several 10-100 days, depending on different accretion rates. With the fact that the stellar luminosity when the star is pulsationally unstable also depends on the accretion rate, we derive the period-luminosity relation log (L/Lsun) = 4.62 + 0.98log(P/100 day), which is testable with future observations. Our models further show that the radius and mass of the pulsating massive protostar should also depend on the period. It would be possible to infer such protostellar properties and the accretion rate with the observed period. Measuring the maser periods enables a direct diagnosis of the structure of accreting massive protostars, which are deeply embedded in dense gas and inaccessible with other observations.

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.

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

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