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 report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

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.

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

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 imaged the continuum and molecular line emission from Orion Source I (SrcI) with up to 30 mas (12 AU) resolution at 43, 99, 223, and 340 GHz in an attempt to probe the structure and chemistry of the circumstellar disk and bipolar outflow associated with this high mass protostar. The continuum spectral index ranges from $\sim$2 along the midplane of the disk to $\sim$3 along the edges, consistent with dust that is optically thick in the midplane but becomes optically thin at the periphery. Salt (NaCl) emission is visible where the dust is optically thin; it provides a unique tracer of the velocity field within the disk. All other molecules that we have mapped - H$_2$O, AlO, SiO, SiS, SO, and SO$_2$ - appear to originate primarily in the bipolar outflow. The base of the outflow is corotating with the disk. SiS shows a filamentary structure that is most prominent along the edges of the outflow. The molecular distributions suggest that Si and Al released from dust grains in the disk react with oxygen derived from H$_2$O to form SiO and AlO, and with SO and SO$_2$ to form SiS.

This paper reports observations of a 22 GHz water maser `superburst' in the G25.65+1.05 massive star forming region, conducted in response to an alert from the Maser Monitoring Organisation (M2O). Very long baseline interferometry (VLBI) observations using the European VLBI Network (EVN) recorded a maser flux density of $1.2 \times 10^{4}$ Jy. The superburst was investigated in the spectral, structural and temporal domains and its cause was determined to be an increase in maser path length generated by the superposition of multiple maser emitting regions aligning in the line of sight to the observer. This conclusion was based on the location of the bursting maser in the context of the star forming region, its complex structure, and its rapid onset and decay.

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

We present the observational result of a glycine precursor, methylamine (CH$_3$NH$_2$), together with methanol (CH$_3$OH) and methanimine (CH$_2$NH) towards high-mass star-forming regions, NGC6334I, G10.47+0.03, G31.41+0.3, and W51~e1/e2 using ALMA. The molecular abundances toward these sources were derived using the rotational diagram method and compared with our state-of-the-art chemical model. We found that the observed ratio of "CH$_3$NH$_2$/CH$_3$OH" is in between 0.11 and 2.2. We also found that the observed "CH$_3$NH$_2$/CH$_3$OH" ratio agrees well with our chemical model by considering the formation of CH$_3$NH$_2$ on the grain surface via hydrogenation process to HCN. This result clearly shows the importance of hydrogenation processes to form CH$_3$NH$_2$. NGC63343I MM3, where CH$_3$NH$_2$ was not detected in this study and showed "CH$_3$NH$_2$/CH$_3$OH" ratio of less than 0.02, is clearly distinguished from the other cores.

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements.

We have observed the Class I protostellar source Elias 29 with Atacama Large Millimeter/submillimeter Array (ALMA). We have detected CS, SO, $^{34}$SO, SO$_2$, and SiO line emissions in a compact component concentrated near the protostar and a ridge component separated from the protostar by 4\arcsec\ ($\sim 500$ au). The former component is found to be abundant in SO and SO$_2$ but deficient in CS. The abundance ratio SO/CS is as high as $3^{+13}_{-2} \times 10^2$ at the protostar, which is even higher than that in the outflow-shocked region of L1157 B1. However, organic molecules (HCOOCH$_3$, CH$_3$OCH$_3$, CCH, and c-C$_3$H$_2$) are deficient in Elias 29. We attribute the deficiency in organic molecules and richness in SO and SO$_2$ to the evolved nature of the source or the relatively high dust temperature (\protect\raisebox{-0.7ex}{$\:\stackrel{\textstyle >}{\sim}\:$} 20 K) in the parent cloud of Elias 29. The SO and SO$_2$ emissions trace rotation around the protostar. Assuming a highly inclined configuration ($i \geq 65$\degr; 0\degr\ for a face-on configuration) and Keplerian motion for simplicity, the protostellar mass is estimated to be (0.8 -- 1.0) \Msun. The $^{34}$SO and SO$_2$ emissions are asymmetric in their spectra; the blue-shifted components are weaker than the red-shifted ones. Although this may be attributed to the asymmetric molecular distribution, other possibilities are also discussed.

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.

The bright maser emission produced by several molecular species at centimeter to long millimeter wavelengths provides an essential tool for understanding the process of massive star formation. Unimpeded by the high dust optical depths that affect shorter wavelength observations, the high brightness temperature of these emission lines offers a way to resolve accretion and outflow motions down to scales below $\sim$1 au in deeply embedded Galactic star-forming regions at kiloparsec distances. The recent identification of extraordinary accretion outbursts in two high-mass protostars, both of which were heralded by maser flares, has rapidly impacted the traditional view of massive protostellar evolution, leading to new hydrodynamic simulations that can produce such episodic outbursts. In order to understand how these massive protostars evolve in response to such events, larger, more sensitive ground-based centimeter wavelength interferometers are needed that can simultaneously image multiple maser species in the molecular gas along with faint continuum from the central ionized gas. Fiducial observations of a large sample of massive protostars will be essential in order to pinpoint the progenitors of future accretion outbursts, and to quantify the outburst-induced changes in their protostellar photospheres and outflow and accretion structures. Knowledge gained from these studies will have broader impact on the general topic of accretion onto massive objects.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

In a rare and spectacular display, two well-known massive star forming regions, W49N and G25.65+1.05, recently underwent maser 'super burst' - their fluxes suddenly increasing above 30,000 and 18,000 Jy, respectively, reaching several orders of magnitude above their usual values. In quick-response, ToO observations with the EVN, VLBA and KaVA were obtained constituting a 4 week campaign - producing a high-cadence multi-epoch VLBI investigation of the maser emission. The combination of high-resolution, polarisation and flux monitoring during the burst provides one of the best accounts, to date, of the maser super burst phenomenon, aiding their use as astrophysical tools. These proceedings contain the preliminary results of our campaign.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

Imaging the bright maser emission produced by several molecular species at centimeter wavelengths is an essential tool for understanding the process of massive star formation because it provides a way to probe the kinematics of dense molecular gas at high angular resolution. Unimpeded by the high dust optical depths that affect shorter wavelength observations, the high brightness temperature of these emission lines offers a way to resolve accretion and outflow motions down to scales as fine as $\sim$1-10 au in deeply embedded Galactic star-forming regions, and at sub-pc scales in nearby galaxies. The Next Generation Very Large Array will provide the capabilities needed to fully exploit these powerful tracers.

We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a mean factor of 6.5 (to 4.2 kJy) with highly-blueshifted features (up to 70 km/s from LSR) becoming more prominent. The strongest flaring water masers reside 3000 au north of MM1B and form a remarkable bow shock pattern whose vertex coincides with CM2 and tail points back to MM1B. Excited OH masers trace a secondary bow shock located ~120 au downstream. ALMA images of CS (6-5) reveal a highly-collimated north-south structure encompassing the flaring masers to the north and the non-flaring masers to the south seen in projection toward the MM3-UCHII region. Proper motions of the southern water masers over 5.3 years indicate a bulk projected motion of 117 km/s southward from MM1B with a dynamical time of 170 yr. We conclude that CM2, the water masers, and many of the excited OH masers trace the interaction of the high velocity bipolar outflow from MM1B with ambient molecular gas. The previously-excavated outflow cavity has apparently allowed the radiative energy of the current outburst to propagate freely until terminating at the northern bow shock where it strengthened the masers. Additionally, water masers have been detected toward MM7 for the first time, and a highly-collimated CS (6-5) outflow has been detected toward MM4.

We have carried out observations in the 42$-$46 and 82$-$103 GHz bands with the Nobeyama 45-m radio telescope, and in the 338.2$-$339.2 and 348.45$-$349.45 GHz bands with the ASTE 10-m telescope toward three high-mass star-forming regions containing massive young stellar objects (MYSOs), G12.89+0.49, G16.86$-$2.16, and G28.28$-$0.36. We have detected HC$_{3}$N including its $^{13}$C and D isotopologues, CH$_{3}$OH, CH$_{3}$CCH, and several complex organic molecules (COMs). Combining our previous results of HC$_{5}$N in these sources, we compare the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratios in the three observed sources. The ratio in G28.28$-$0.36 is derived to be $0.091^{+0.109}_{-0.039}$, which is higher than that in G12.89+0.49 by one order of magnitude, and that in G16.86$-$2.16 by a factor of $\sim 5$. We investigate the relationship between the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratio and the $N$(CH$_{3}$CCH)/$N$(CH$_{3}$OH) ratio. The relationships of the two column density ratios in G28.28$-$0.36 and G16.86$-$2.16 are similar to each other, while HC$_{5}$N is less abundant when compared to CH$_{3}$CCH in G12.89+0.49. These results imply a chemical diversity in the lukewarm ($T \sim 20-30$ K) envelope around MYSOs. Besides, several spectral lines from complex organic molecules, including very-high-excitation energy lines, have been detected toward G12.89+0.49, while the line density is significantly low in G28.28$-$0.36. These results suggest that organic-poor MYSOs are surrounded by a carbon-chain-rich lukewarm envelope (G28.28$-$0.36), while organic-rich MYSOs, namely hot cores, are surrounded by a CH$_{3}$OH-rich lukewarm envelope (G12.89+0.49 and G16.86$-$2.16).

The study of the chemical evolution of glycine in the interstellar medium is one of challenging topics in astrochemistry. Here, we present the chemical modeling of glycine in hot cores using the state-of-the-art three-phase chemical model NAUTILUS, which is focused on the latest glycine chemistry. For the formation process of glycine on the grain surface, we obtained consistent results with previous studies that glycine would be formed via the reactions of COOH with CH$_2$NH$_2$. However, we will report three important findings regarding the chemical evolution and the detectability of interstellar glycine. First, with the experimentally obtained binding energy from the temperature programmed thermal desorption (TPD) experiment, a large proportion of glycine was destroyed through the grain surface reactions with NH or CH$_3$O radicals before it fully evaporates. As a result, the formation process in the gas phase is more important than thermal evaporation from grains. If this is the case, NH$_2$OH and CH$_3$COOH rather than CH$_3$NH$_2$ and CH$_2$NH would be the essential precursors to the gas phase glycine. Secondly, since the gas phase glycine will be quickly destroyed by positive ions or radicals, early evolutionary phase of the hot cores would be the preferable target for the future glycine surveys. Thirdly, we suggest the possibility that the suprathermal hydrogen atoms can strongly accelerate the formation of COOH radicals from CO$_2$, resulting in the dramatic increase of formation rate of glycine on grains. The efficiency of this process should be investigated in detail by theoretical and experimental studies in the future.

Formation processes of high-mass stars have been long-standing issues in astronomy and astrophysics. This is mainly because of major difficulties in observational studies such as a smaller number of high-mass young stellar objects (YSOs), larger distances, and more complex structures in young high-mass clusters compared with nearby low-mass isolated star-forming regions (SFRs), and extremely large opacity of interstellar dust except for centimeter to submillimeter wavelengths. High resolution and high sensitivity observations with Atacama Large Millimeter/Submillimeter Array (ALMA) at millimeter/submillimeter wavelengths will overcome these observational difficulties even for statistical studies with increasing number of high-mass YSO samples. This review will summarize recent progresses in high-mass star-formation studies with ALMA such as clumps and filaments in giant molecular cloud complexes and infrared dark clouds (IRDCs), protostellar disks and outflows in dense cores, chemistry, masers, and accretion bursts in high-mass SFRs.

The different spatial distributions of N-bearing and O-bearing species, as is well known towards Orion~KL, is one of the long-lasting mysteries. We conducted a survey observation and chemical modeling study to investigate if the different distributions of O- and N-bearing species are widely recognized in general star-forming regions. First, we report our observational results of complex organic molecules (COMs) with the 45~m radio telescope at the Nobeyama Radio Observatory towards eight star-forming regions. Through our spectral survey ranging from 80 to 108~GHz, we detected CH$_3$OH, HCOOCH$_3$, CH$_3$OCH$_3$, (CH$_3$)$_2$CO, CH$_3$CHO, CH$_3$CH$_2$CN, CH$_2$CHCN, and NH$_2$CHO. Their molecular abundances were derived via the rotation diagram and the least squares methods. We found that N-bearing molecules, tend to show stronger correlations with other N-bearing molecules rather than O-bearing molecules. While G10.47+0.03 showed high fractional abundances of N-bearing species, those in NGC6334F were not so rich, being less than 0.01 compared to CH$_3$OH. Then, the molecular abundances towards these sources were evaluated by chemical modeling with NAUTILUS three-phase gas-grain chemical code. Through the simulations of time evolutions for the abundances of COMs, we suggest that observed correlations of fractional abundances between COMs can be explained by the combination of the different temperature structures inside the hot cores and the different evolutionary phase. Since our modeling could not fully explain the observed excitation temperatures, it is important to investigate the efficiency of grain surface reactions and their activation barriers, and the binding energy of COMs to further promote our understanding.

In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$.

Chemical diversity of the gas in low-mass protostellar cores is widely recognized. In order to explore its origin, a survey of chemical composition toward 36 Class 0/I protostars in the Perseus molecular cloud complex, which are selected in an unbiased way under certain physical conditions, has been conducted with IRAM 30 m and NRO 45 m telescope. Multiple lines of C2H, c-C3H2 and CH3OH have been observed to characterize the chemical composition averaged over a 1000 au scale around the protostar. The derived beam-averaged column densities show significant chemical diversity among the sources, where the column density ratios of C2H/CH3OH are spread out by 2 orders of magnitude. From previous studies, the hot corino sources have abundant CH3OH but deficient C2H, their C2H/CH3OH column density ratios being relatively low. In contrast, the warm-carbon-chain chemistry (WCCC) sources are found to reveal the high C2H/CH3OH column density ratios. We find that the majority of the sources have intermediate characters between these two distinct chemistry types. A possible trend is seen between the C2H/CH3OH ratio and the distance of the source from the edge of a molecular cloud. The sources located near cloud edges or in isolated clouds tend to have a high C2H/CH3OH ratio. On the other hand, the sources having a low C2H/CH3OH ratio tend to be located in inner regions of the molecular cloud complex. This result gives an important clue to an understanding of the origin of the chemical diversity of protostellar cores in terms of environmental effects.

We have observed complex organic molecules (COMs) and deuterated species toward a hot core/corino (HC) associated with the infrared dark cloud (IRDC) clump G34.43+00.24 MM3 with the Atacama Large Millimeter/submillimeter Array (ALMA). We have detected six normal-COMs (CH3OH, CH3CHO, CH3CH2CN, CH3OCH3, HCOOCH3, and NH2CHO), one deuterated-COM (CH2DCN), and two deuterated fundamental molecules (D2CO and DNC) toward G34.43+00.24 MM3 HC. None of these lines, except for CH3OH, are detected toward the shocked regions in our data, which suggests that COMs do not originate in shocks. The abundance of the COMs relative to CH3OH in G34.43+00.24 MM3 HC is found to be similar to those in high-mass hot cores, rather than those in hot corinos in low-mass star-forming regions. This result suggests that the physical conditions of the warm-up phase of G34.43+00.24 MM3 HC are similar to those of high-mass sources. On the other hand, the D2CO abundance relative to CH$_3$OH in G34.43+00.24 MM3 HC is higher than that of other hot cores, and seems to be comparable to that of hot corinos. The relatively high D2CO/CH3OH ratio of G34.43+00.24 MM3 HC implies a long cold starless phase of G34.43+00.24 MM3 HC.

We have investigated the TUKH122 prestellar core in the Orion A cloud using ALMA 3 mm dust continuum, N$_2$H$^+$ ($J=1-0$), and CH$_3$OH ($J_K=2_K-1_K$) molecular line observations. Previous studies showed that TUKH122 is likely on the verge of star formation because the turbulence is almost dissipated and chemically evolved among other starless cores in the Orion A cloud. By combining ALMA 12-m and ACA data, we recover extended emission with a resolution of $\sim5"$ corresponding to 0.01 pc and identify 6 condensations with a mass range of $0.1-0.4$ $M_\odot$ and a radius of $\lesssim0.01$ pc. These condensations are gravitationally bound following a virial analysis and are embedded in the filament including the elongated core with a mass of $\sim29$ $M_\odot$ and a radial density profile of $r^{-1.6}$ derived by {\it Herschel}. The separation of these condensations is $\sim0.035$ pc, consistent with the thermal jeans length at a density of $4.4\times10^5$ cm$^{-3}$. This density is similar to the central part of the core. We also find a tendency that the N$_2$H$^+$ molecule seems to deplete at the dust peak condensation. This condensation may be beginning to collapse because the linewidth becomes broader. Therefore, the fragmentation still occurs in the prestellar core by thermal Jeans instability and multiple stars are formed within the TUKH122 prestellar core. The CH$_3$OH emission shows a large shell-like distribution and surrounds these condensations, suggesting that the CH$_3$OH molecule formed on dust grains is released into gas phase by non-thermal desorption such as photoevaporation caused by cosmic-ray induced UV radiation.

We report results of the observations of 15 regions in several molecular lines for a statistical study of massive cluster-forming clumps. We identified 24 clumps based on the C18O (J=1-0) data obtained by the NRO 45 m telescope, and found that 16 of them are associated with young clusters. The clumps associated with clusters have a typical mass, radius, and molecular density of ~1 X 10^3 Mo, ~0.5 pc, ~1 X 10^5 cm^-3, respectively. We categorized the clumps and clusters into four types according to the spatial coincidence of gas and star density, and discussed their evolutions: Clumps without clusters (Type 1), clumps showing good correlations with clusters (Type 2), clumps showing poor correlations with clusters (Type 3), and clusters with no associated clumps (Type 4). Analyses of the velocity structures and the chemical compositions imply that the clump + cluster systems should evolve from Type 1 to Type 4. We found that some of the Type 2 clumps are infalling on the clump-scale to form clusters at the clump center, which should commonly occur in the beginning of cluster formation. Interestingly, all of the identified Type 1 clumps are likely to be older than the Type 2 clumps in terms of chemical compositions, suggesting that they have been gravitationally stable for a long time possibly being supported by the strong magnetic field of > 1 mG.Type 1 clumps younger than the observed Type 2 clumps should be very rare to find because of their short lifetime.

We report the first sub-arcsecond VLA imaging of 6 GHz continuum, methanol maser, and excited-state hydroxyl maser emission toward the massive protostellar cluster NGC6334I following the recent 2015 outburst in (sub)millimeter continuum toward MM1, the strongest (sub)millimeter source in the protocluster. In addition to detections toward the previously known 6.7 GHz Class II methanol maser sites in the hot core MM2 and the UCHII region MM3 (NGC6334F), we find new maser features toward several components of MM1, along with weaker features $\sim1''$ north, west, and southwest of MM1, and toward the non-thermal radio continuum source CM2. None of these areas have heretofore exhibited Class II methanol maser emission in three decades of observations. The strongest MM1 masers trace a dust cavity, while no masers are seen toward the strongest dust sources MM1A, 1B and 1D. The locations of the masers are consistent with a combination of increased radiative pumping due to elevated dust grain temperature following the outburst, the presence of infrared photon propagation cavities, and the presence of high methanol column densities as indicated by ALMA images of thermal transitions. The non-thermal radio emission source CM2 ($2''$ north of MM1) also exhibits new maser emission from the excited 6.035 and 6.030 GHz OH lines. Using the Zeeman effect, we measure a line-of-sight magnetic field of +0.5 to +3.7 mG toward CM2. In agreement with previous studies, we also detect numerous methanol and excited OH maser spots toward the UCHII region MM3, with predominantly negative line-of-sight magnetic field strengths of -2 to -5 mG and an intriguing south-north field reversal.

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

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.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

We present sensitive high angular resolution ($\sim$ 0.1$''$ -- 0.3$''$) continuum ALMA (The Atacama Large Millimeter/Submillimeter Array) observations of the archetypal hot core located in Orion-KL. The observations were made in five different spectral bands (bands 3, 6, 7, 8, and 9) covering a very broad range of frequencies (149 -- 658 GHz). Apart of the well-know millimeter emitting objects located in this region (Orion Source I and BN), we report the first submillimeter detection of three compact continuum sources (ALMA 1-3) in the vicinities of the Orion-KL hot molecular core. These three continuum objects have spectral indices between 1.47 to 1.56, and brightness temperatures between 100 to 200 K at 658 GHz suggesting that we are seeing moderate optically thick dust emission with possible grain growth. However, as these objects are not associated with warm molecular gas, and some of them are farther out from the molecular core, we thus conclude that they cannot heat the molecular core. This result favours the hypothesis that the hot molecular core in Orion-KL core is heated externally.

We have conducted an interferometric line survey in the 0.8 mm band toward the young high-mass protostar candidate NGC 2264 CMM3 with ALMA. CMM3 is resolved into the two continuum peaks, CMM3A and CMM3B, at an angular separation of 0.9". Thus, CMM3 is found to be a binary system candidate. We have detected molecular outflows associated with CMM3A and CMM3B each, indicating active star formation. In addition to the two peaks, six faint continuum peaks are detected around CMM3A and CMM3B, most of which are thought to be evolved low-mass protostars. CMM3A is found to be rich in molecular line emission including complex organic molecules such as HCOOCH3 and CH3OCH3. The emission of complex organic molecules is distributed within a compact region around the continuum peak of CMM3A. Hence, CMM3A apparently harbors a hot core. On the other hand, CMM3B is deficient in molecular line emission, although its continuum flux is almost comparable to that of CMM3A. Possible origins of the striking difference between CMM3A and CMM3B are discussed.

It has been a long-standing problem to detect interstellar glycine (NH$_2$CH$_2$COOH), the simplest amino acid, in studying a possible relation between the Universe and origin of life. In the last about 40 years all surveys of glycine failed, and it would be an alternative strategy to search for precursor(s) to glycine. Such studies of precursors would be crucial prior to conducting sensitive surveys by ALMA. Laboratory studies have suggested that CH$_3$NH$_2$ is a possible precursor to glycine. Further theoretical study also suggested that the CH$_2$NH$_2$ radical that can be formed from CH$_3$NH$_2$ through photodissociation can be a good precursor to glycine. Thus we observed CH$_3$NH$_2$ towards several hot core sources by using the Nobeyama 45m radio telescope, and succeeded in finding a new CH$_3$NH$_2$ source, G10.47+0.03, with its fractional abundance of 2.4+-0.6x10^{-8}; at the time of writing, this source is the most abundant source of CH$_3$NH$_2$ ever known. We also found another new source of CH$_3$NH$_2$, NGC6334F, however, the detection may be tentative since the signal-to-noise ratios are not sufficient to claim secure detection. Detailed analysis of the detected data revealed the detection of the hyperfine components of CH$_3$NH$_2$ for the first time. We found that the observed abundance of CH$_3$NH$_2$ agrees fairly well with the theoretically predicted value by Garrod (2013). Detectability of interstellar glycine is discussed.

We used the Atacama Large (sub-)Millimeter Array (ALMA) and the IRAM Plateau de Bure Interferometer (PdBI) to image, with an angular resolution of 0.5$''$ (120 au) and 1$''$ (235 au), respectively, the emission from 11 different organic molecules in the protostellar binary NGC1333 IRAS 4A. We clearly disentangled A1 and A2, the two protostellar cores present. For the first time, we were able to derive the column densities and fractional abundances simultaneously for the two objects, allowing us to analyse the chemical differences between them. Molecular emission from organic molecules is concentrated exclusively in A2 even though A1 is the strongest continuum emitter. The protostellar core A2 displays typical hot corino abundances and its deconvolved size is 70 au. In contrast, the upper limits we placed on molecular abundances for A1 are extremely low, lying about one order of magnitude below prestellar values. The difference in the amount of organic molecules present in A1 and A2 ranges between one and two orders of magnitude. Our results suggest that the optical depth of dust emission at these wavelengths is unlikely to be sufficiently high to completely hide a hot corino in A1 similar in size to that in A2. Thus, the significant contrast in molecular richness found between the two sources is most probably real. We estimate that the size of a hypothetical hot corino in A1 should be less than 12 au. Our results favour a scenario in which the protostar in A2 is either more massive and/or subject to a higher accretion rate than A1, as a result of inhomogeneous fragmentation of the parental molecular clump. This naturally explains the smaller current envelope mass in A2 with respect to A1 along with its molecular richness.

We carried out line survey observations at the 26-30 GHz band toward the four high-mass star-forming regions containing hot cores, G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Robert C. Byrd Green Bank Telescope. We have detected HC5N from all of the sources, and HC7N from the three sources, except for G10.30-0.15. We further conducted observations of HC5N at the 42-46 GHz and 82-103 GHz bands toward the three sources, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Nobeyama 45 m radio telescope. The rotational lines of HC5N with the high excitation energies (Eu/k=63-100 K), which are hardly excited in the cold dark clouds, have been detected from the three sources. The rotational temperatures of HC5N are found to be 13-20 K in the three sources. The detection of the lines with the high excitation energies and the derived rotational temperatures indicate that HC5N exists in the warm gas within 0.07-0.1 pc radii around massive young stellar objects. The column densities of HC5N in the three sources are derived to be (2.0+-2.8)*10^(13) cm-2 . We compare the ratios between N(HC5N) the column density of HC5N and W(CH3OH) the integrated intensity of the thermal CH3OH emission line among the three high-mass star-forming regions. We found a possibility of the chemical differentiation in the three high-mass star-forming regions; G28.28-0.36 shows the largest N (HC5N)/W (CH3OH) ratio of > 8.0*10^(14) in units of (K km s-1 )-1 cm-2 , while G12.89+0.49 and G16.86-2.16 show the smaller values (~2*10^(13) ).

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

We have resolved for the first time the radial and vertical structure of the almost edge-on envelope/disk system of the low-mass Class 0 protostar L1527. For that, we have used ALMA observations with a spatial resolution of 0.25$^{\prime\prime}$$\times$0.13$^{\prime\prime}$ and 0.37$^{\prime\prime}$$\times$0.23$^{\prime\prime}$ at 0.8 mm and 1.2 mm, respectively. The L1527 dust continuum emission has a deconvolved size of 78 au $\times$ 21 au, and shows a flared disk-like structure. A thin infalling-rotating envelope is seen in the CCH emission outward of about 150 au, and its thickness is increased by a factor of 2 inward of it. This radius lies between the centrifugal radius (200 au) and the centrifugal barrier of the infalling-rotating envelope (100 au). The gas stagnates in front of the centrifugal barrier and moves toward vertical directions. SO emission is concentrated around and inside the centrifugal barrier. The rotation speed of the SO emitting gas is found to be decelerated around the centrifugal barrier. A part of the angular momentum could be extracted by the gas which moves away from the mid-plane around the centrifugal barrier. If this is the case, the centrifugal barrier would be related to the launching mechanism of low velocity outflows, such as disk winds.

The Class 0 protostar, L483, has been observed in various molecular lines in the 1.2 mm band at a sub-arcsecond resolution with ALMA. An infalling-rotating envelope is traced by the CS line, while a very compact component with a broad velocity width is observed for the CS, SO, HNCO, NH$_2$CHO, and HCOOCH$_3$ lines. Although this source is regarded as the warm carbon-chain chemistry (WCCC) candidate source at a 1000 au scale, complex organic molecules characteristic of hot corinos such as NH$_2$CHO and HCOOCH$_3$ are detected in the vicinity of the protostar. Thus, both hot corino chemistry and WCCC are seen in L483. Although such a mixed chemical character source has been recognized as an intermediate source in previous single-dish observations, we here report the first spatially-resolved detection. A kinematic structure of the infalling-rotating envelope is roughly explained by a simple ballistic model with the protostellar mass of 0.1--0.2 $M_\odot$ and the radius of the centrifugal barrier (a half of the centrifugal radius) of 30--200 au, assuming the inclination angle of 80\degr\ (0\degr\ for a face-on). The broad line emission observed in the above molecules most likely comes from the disk component inside the centrifugal barrier. Thus, a drastic chemical change is seen around the centrifugal barrier.

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.

Observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

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 report the first evidence of a hot corino in a Bok globule. This is based on the ALMA observations in the 1.2 mm band toward the low-mass Class 0 protostar IRAS 19347+0727 in B335. Saturated complex organic molecules (COMs), CH$_3$CHO, HCOOCH$_3$, and NH$_2$CHO, are detected in a compact region within a few 10 au around the protostar. Additionally, CH$_3$OCH$_3$, C$_2$H$_5$OH, C$_2$H$_5$CN, and CH$_3$COCH$_3$ are tentatively detected. Carbon-chain related molecules, CCH and c-C$_3$H$_2$, are also found in this source, whose distributions are extended over a few 100 au scale. On the other hand, sulfur-bearing molecules CS, SO, and SO$_2$, have both compact and extended components. Fractional abundances of the COMs relative to H$_2$ are found to be comparable to those in known hot-corino sources. Though the COMs lines are as broad as 5-8 km s$^{-1}$, they do not show obvious rotation motion in the present observation. Thus, the COMs mainly exist in a structure whose distribution is much smaller than the synthesized beam (0."58 x 0."52).

Context. ALMA observations of the high-mass star-forming region G35.20-0.74N have revealed the presence of a Keplerian disk in core B rotating about a massive object of 18 Msun, as computed from the velocity field. The luminosity of such a massive star would be comparable to (or higher than) the luminosity of the whole star-forming region. To solve this problem it has been proposed that core B could harbor a binary system. This could also explain the possible precession of the radio jet associated with this core, which has been suggested by its S-shaped morphology. Aims. To establish the origin of the free-free emission from core B and investigate the existence of a binary system at the center of this massive core and the possible precession of the radio jet. Methods. We carried out VLA continuum observations of G35.20-0.74N at 2 cm in the B configuration and at 1.3 cm and 7 mm in the A and B configurations. The bandwidth at 7 mm covers the CH3OH maser line at 44.069 GHz. Continuum images at 6 and 3.6 cm in the A configuration were obtained from the VLA archive. We also carried out VERA observations of the H2O maser line at 22.235 GHz. Results. The observations have revealed the presence of a binary system of UC/HC HII regions at the geometrical center of the radio jet in G35.20-0.74N. This binary system, which is associated with a Keplerian rotating disk, consists of two B-type stars of 11 and 6 Msun. The S-shaped morphology of the radio jet has been successfully explained as due to precession produced by the binary system. The analysis of the precession of the radio jet has allowed us to better interpret the IR emission in the region, which would be not tracing a wide-angle cavity open by a single outflow with a position angle of ~55 deg but two different flows: a precessing one in the NE-SW direction associated with the radio jet, and a second one in an almost E--W direction ...

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

We conducted survey observations of a glycine precursor, methanimine or methylenimine (CH$_2$NH), with the NRO 45 m telescope and the SMT telescope towards 12 high-mass and two low-mass star-forming regions in order to increase number of CH$_2$NH sources and to better understand the characteristics of CH2NH sources. As a result of our survey, CH$_2$NH was detected in eight sources, including four new sources. The estimated fractional abundances were ~10$^8$ in Orion KL and G10.47+0.03, while they were ~10$^9$ towards the other sources. Our hydrogen recombination line and past studies suggest that CH$_2$NH-rich sources have less evolved HII regions. The less destruction rates by UV flux from the central star would be contributed to the high CH$_2$NH abundances towards CH$_2$NH-rich sources. Our gas-grain chemical simulations suggest that CH$_2$NH is mostly formed in the gas-phase by neutral-neutral reactions rather than the product of thermal evaporation from the dust surfaces.

We report the detection of a wing component in NH$_3$ emission toward the starless core TUKH122 with subthermal turbulence in the Orion A cloud. This NH$_3$ core is suggested to be on the verge of star formation because the turbulence inside the NH$_3$ core is almost completely dissipated, and also because it is surrounded by CCS, which resembles the prestellar core L1544 in Taurus showing infall motions. Observations were carried out with the Nobeyama 45 m telescope at 0.05 km s$^{-1}$ velocity resolution. We find that the NH$_3$ line profile consists of two components. The quiescent main component has a small linewidth of 0.3 km s$^{-1}$ dominated by thermal motions, and the red-shifted wing component has a large linewidth of 1.36 km s$^{-1}$ representing turbulent motions. These components show kinetic temperatures of 11 K and $<$ 30 K, respectively. Furthermore, there is a clear velocity offset between the NH$_3$ quiescent gas ($VLSR=3.7$ km s$^{-1}$) and the turbulent gas ($VLSR=4.4$ km s$^{-1}$). The centroid velocity of the turbulent gas corresponds to that of the surrounding gas traced by the $^{13}$CO ($J=1-0$) and CS ($J=2-1$) lines. LVG model calculations for CS and CO show that the turbulent gas has a temperature of $8-13$ K and an H$_2$ density of ~ $10^4$ cm$^{-3}$, suggesting that the temperature of the turbulent component is also ~ 10 K. The detections of both NH$_3$ quiescent and wing components may indicate a sharp transition from the turbulent parent cloud to the quiescent dense core.

Sub-arcsecond images of the rotational line emission of CS and SO have been obtained toward the Class I protostar IRAS 04365$+$2535 in TMC-1A with ALMA. A compact component around the protostar is clearly detected in the CS and SO emission. The velocity structure of the compact component of CS reveals infalling-rotating motion conserving the angular momentum. It is well explained by a ballistic model of an infalling-rotating envelope with the radius of the centrifugal barrier (a half of the centrifugal radius) of 50 AU, although the distribution of the infalling gas is asymmetric around the protostar. The distribution of SO is mostly concentrated around the radius of the centrifugal barrier of the simple model. Thus a drastic change in chemical composition of the gas infalling onto the protostar is found to occur at a 50 AU scale probably due to accretion shocks, demonstrating that the infalling material is significantly processed before being delivered into the disk.

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 are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

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 have observed the HC3N (J=10-9) and N2H+ (J=1-0) lines toward the Vela C molecular clouds with the Mopra 22 m telescope to study chemical characteristics of dense cores. The intensity distributions of these molecules are similar to each other at an angular resolution of 53", corresponding to 0.19 pc suggesting that these molecules trace the same dense cores. We identified 25 local peaks in the velocity-integrated intensity maps of the HC3N and/or N2H+ emission. Assuming LTE conditions, we calculated the column densities of these molecules and found a tendency that N2H+/HC3N abundance ratio seems to be low in starless regions while it seems to be high in star-forming regions, similar to the tendencies in the NH3/CCS, NH3/HC3N, and N2H+/CCS abundance ratios found in previous studies of dark clouds and the Orion A GMC. We suggest that carbon chain molecules, including HC3N, may trace chemically young molecular gas and N-bearing molecules, such as N2H+, may trace later stages of chemical evolution in the Vela C molecular clouds. It may be possible that the N2H+/HC3N abundance ratio of ~ 1.4 divides the star-forming and starless peaks in the Vela C, although it is not as clear as those in NH3/CCS, NH3/HC3N, and N2H+/CCS for the Orion A GMC. This less clear separation may be caused by our lower spatial resolution or the misclassification of star-forming and starless peaks due to the larger distance of the Vela C. It might be also possible that the HC3N (J=10-9) transition is not a good chemical evolution tracer compared with CCS (J=4-3 and 7-6) transitions.

Spectral line survey observations are conducted toward the high-mass protostar candidate NGC 2264 CMM3 in the 4 mm, 3 mm, and 0.8 mm bands with the Nobeyama 45 m telescope and the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. In total, 265 emission lines are detected in the 4 mm and 3 mm bands, and 74 emission lines in the 0.8 mm band. As a result, 36 molecular species and 30 isotopologues are identified. In addition to the fundamental molecular species, many emission lines of carbon-chain molecules such as HC5N, C4H, CCS, and C3S are detected in the 4 mm and 3 mm bands. Deuterated molecular species are also detected with relatively strong intensities. On the other hand, emission lines of complex organic molecules such as HCOOCH3, and CH3OCH3 are found to be weak. For the molecules for which multiple transitions are detected, rotation temperatures are derived to be 7-33 K except for CH3OH. Emission lines with high upper-state energies (Eu > 150 K) are detected for CH3OH, indicating existence of a hot core. In comparison with the chemical composition of the Orion KL, carbon-chain molecules and deuterated molecules are found to be abundant in NGC 2264 CMM3, while sulfur-bearing species and complex organic molecules are deficient. These characteristics indicate chemical youth of NGC 2264 CMM3 in spite of its location at the center of the cluster forming core, NGC 2264 C.

We recently reported the first identification of rotational transitions of methyl formate (HCOOCH$_{3}$) in the second torsionally excited state toward Orion Kleinmann-Low (KL) observed with the Nobeyama 45 m telescope. In combination with the identified transitions of methyl formate in the ground state and the first torsional excited state, it was found that there is a difference in rotational temperature and vibrational temperature, where the latter is higher. In this study, high spatial resolution analysis by using Atacama Large Millimeter/Submillimeter Array (ALMA) science verification data was carried out to verify and understand this difference. Toward the Compact Ridge, two different velocity components at 7.3 and 9.1 km s$^{-1}$ were confirmed, while a single component at 7.3 km s$^{-1}$ was identified towards the Hot Core. The intensity maps in the ground, first, and second torsional excited states have quite similar distributions. Using extensive ALMA data, we determined the rotational and vibrational temperatures for the Compact Ridge and Hot Core by the conventional rotation diagram method. The rotational temperature and vibrational temperatures agree for the Hot Core and for one component of the Compact Ridge. At the 7.3 km s$^{-1}$ velocity component for the Compact Ridge, the rotational temperature was found to be higher than the vibrational temperature. This is different from what we obtained from the results by using the single-dish observation. The difference might be explained by the beam dilution effect of the single-dish data and/or the smaller number of observed transitions within the limited range of energy levels ($\leq$30 K) of $E_u$ in the previous study.

We have observed the clump G34.43+00.24 MM3 associated with an infrared dark cloud in DNC $J$=3--2, HN$^{13}$C $J$=3--2, and N$_2$H$^+$ $J$=3--2 with the Atacama Large Millimeter/submillimeter Array (ALMA). The N$_2$H$^+$ emission is found to be relatively weak near the hot core and the outflows, and its distribution is clearly anti-correlated with the CS emission. This result indicates that a young outflow is interacting with cold ambient gas. The HN$^{13}$C emission is compact and mostly emanates from the hot core, whereas the DNC emission is extended around the hot core. Thus, the DNC and HN$^{13}$C emission traces warm regions near the protostar differently. The DNC emission is stronger than the HN$^{13}$C emission toward most parts of this clump. The DNC/HNC abundance ratio averaged within a $15^{\prime\prime} \times 15^{\prime\prime}$ area around the phase center is higher than 0.06. This ratio is much higher than the value obtained by the previous single-dish observations of DNC and HN$^{13}$C $J$=1--0 ($\sim$0.003). It seems likely that the DNC and HNC emission observed with the single-dish telescope traces lower density envelopes, while that observed with ALMA traces higher density and highly deuterated regions. We have compared the observational results with chemical-model results in order to investigate the behavior of DNC and HNC in the dense cores. Taking these results into account, we suggest that the low DNC/HNC ratio in the high-mass sources obtained by the single-dish observations are at least partly due to the low filling factor of the high density regions.

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 have observed a molecular clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the CH3OH J_K=9_{-1}--8_0 E, 5_0-4_0 E, and 5_{-1}-4_{-1} E lines at sub-arcsecond resolution by using the Atacama Large Millimeter/submillimeter Array. By comparing the CH3OH J_K=9_{-1}-8_0 E emission with the CH3OH 5_0-4_0 E and 5_{-1}-4_{-1} E emission, we have found that the CH3OH J_K=9_{-1}-8_0 E emission is masing. We have clearly shown that the CH3OH J_K=9_{-1}-8_0 masers arise from the post shocked gas in the interacting regions between the outflows and ambient dense gas. Toward the strongest peak of the CH3OH maser emission, SiO J=6-5 emission is very weak. This indicates that the CH3OH maser emission traces relatively old shocks or weak shocks.

Sub-arcsecond ($0.^{\prime\prime}5$) images of H$_2$CO and CCH line emission have been obtained in the $0.8$ mm band toward the low-mass protostar IRAS 15398-3359 in the Lupus 1 cloud as one of the Cycle 0 projects of the Atacama Large Millimeter/Submillimeter Array. We have detected a compact component concentrated in the vicinity of the protostar and a well-collimated outflow cavity extending along the northeast-southwest axis. The inclination angle of the outflow is found to be about $20^\circ$, or almost edge-on, based on the kinematic structure of the outflow cavity. This is in contrast to previous suggestions of a more pole-on geometry. The centrally concentrated component is interpreted by use of a model of the infalling rotating envelope with the estimated inclination angle, and the mass of the protostar is estimated to be less than $0.09\ M_\odot$. Higher spatial resolution data are needed to infer the presence of a rotationally supported disk for this source, hinted at by a weak high-velocity H$_2$CO emission associated with the protostar.

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

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

The Serpens South infrared dark cloud consists of several filamentary ridges, some of which fragment into dense clumps. On the basis of CCS ($J_N=4_3-3_2$), HC$_3$N ($J=5-4$), N$_2$H$^+$ ($J=1-0$), and SiO ($J=2-1, v=0$) observations, we investigated the kinematics and chemical evolution of these filamentary ridges. We find that CCS is extremely abundant along the main filament in the protocluster clump. We emphasize that Serpens South is the first cluster-forming region where extremely-strong CCS emission is detected. The CCS-to-N$_2$H$^+$ abundance ratio is estimated to be about 0.5 toward the protocluster clump, whereas it is about 3 in the other parts of the main filament. We identify six dense ridges with different $V_{\rm LSR}$. These ridges appear to converge toward the protocluster clump, suggesting that the collisions of these ridges may have triggered cluster formation. The collisions presumably happened within a few $\times \ 10^5$ yr because CCS is abundant only in such a short time. The short lifetime agrees with the fact that the number fraction of Class I objects, whose typical lifetime is $0.4 \times \ 10^5$ yr, is extremely high as about 70 percent in the protocluster clump. In the northern part, two ridges appear to have partially collided, forming a V-shape clump. In addition, we detected strong bipolar SiO emission that is due to the molecular outflow blowing out of the protostellar clump, as well as extended weak SiO emission that may originate from the filament collisions.

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 imaged two starless molecular cloud cores, TUKH083 and TUKH122, in the Orion A giant molecular cloud in the CCS and ammonia (NH$_3$) emission with the Very Large Array. TUKH122 contains one NH$_3$ core "TUKH122-n," which is elongated and has a smooth oval boundary. Where observed, the CCS emission surrounds the NH$_3$ core. This configuration resembles that of the N$_2$H$^+$ and CCS distribution in the Taurus starless core L1544, a well-studied example of a dense prestellar core exhibiting infall motions. The linewidth of TUKH122-n is narrow (0.20 km s$^{-1}$) in the NH$_3$ emission line and therefore dominated by thermal motions. The smooth oval shape of the core boundary and narrow linewidth in NH$_3$ seem to imply that TUKH122-n is dynamically relaxed and quiescent. TUKH122-n is similar to L1544 in the kinetic temperature (10 K), linear size (0.03 pc), and virial mass ($\sim$ 2 $M_{\odot}$). Our results strongly suggest that TUKH122-n is on the verge of star formation. TUKH122-n is embedded in the 0.2 pc massive (virial mass $\sim$ 30 $M_{\odot}$) turbulent parent core, while the L1544 NH$_3$ core is embedded in the 0.2 pc less-massive (virial mass $\sim$ 10 $M_{\odot}$) thermal parent core. TUKH083 shows complicated distribution in NH$_3$, but was not detected in CCS. The CCS emission toward TUKH083 appears to be extended, and is resolved out in our interferometric observations.

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 used the Nobeyama 45 m radio telescope to reveal spatial distributions of CCS and HC$_3$N in L1147, one of carbon-chain producing regions (CCPRs) candidates, where carbon-chain molecules are dominant rather than NH$_3$. We found that three cores (two CCS cores and one HC$_3$N core) exist along the NE-SW filament traced by the 850 micron dust continuum, which are away from a Very Low Luminosity Object (VeLLO - a source that may turn into sub-stellar mass brown dwarf). The column densities of CCS are 3-7x10$^{12}$ cm$^{-2}$ and those of HC$_3$N are 2-6x10$^{12}$ cm$^{-2}$, respectively, much lower than those previously reported towards other CCPRs. We also found that two CCS peaks are displaced from that of HC$_3$N. In order to interpret such interleaved distributions, we conducted chemical reaction network simulations, and found that slightly different gas densities could lead to large variation of CCS-to-HC$_3$N ratio in the early phase of dark cloud evolution. Such a chemical "variation" may be seen in other CCPRs. Finally we were able to confirm that the L1147 filament can be regarded as a CCPR.

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.

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.

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 have mapped six molecular cloud cores in the Orion A giant molecular cloud (GMC), whose kinetic temperatures range from 10 to 30 K, in CCS and N2H+ with Nobeyama 45 m radio telescope to study their chemical characteristics. We identified 31 intensity peaks in the CCS and N2H+ emission in these molecular cloud cores. It is found for cores with temperatures lower than ~ 25 K that the column density ratio of N(N2H+)/N(CCS) is low toward starless core regions while it is high toward star-forming core regions, in case that we detected both of the CCS and N2H+ emission. This is very similar to the tendency found in dark clouds (kinetic temperature ~ 10 K). The criterion found in the Orion A GMC is N(N2H+)/N(CCS) ~ 2-3. In some cases, the CCS emission is detected toward protostars as well as the N2H+ emission. Secondary late-stage CCS peak in the chemical evolution caused by CO depletion may be a possible explanation for this. We found that the chemical variation of CCS and N2H+ can also be used as a tracer of evolution in warm (10-25 K) GMC cores. On the other hand, some protostars do not accompany N2H+ intensity peaks but are associated with dust continuum emitting regions, suggesting that the N2H+ abundance might be decreased due to CO evaporation in warmer star-forming sites.

We have observed a cluster forming clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the 1.3 mm continuum and the CH3OH, CS, 13CS, SiO, CH3CH2CN, and HCOOCH3 lines with the Atacama Large Millimeter/submillimeter Array and in K-band with the Keck telescope. We have found a young outflow toward the center of this clump in the SiO, CS, and CH3OH lines. This outflow is likely driven by a protostar embedded in a hot core, which is traced by the CH3CH2CN, HCOOCH3, 13CS, and high excitation CH3OH lines. The size of the hot core is about 800 x 300 AU in spite of its low mass (<1.1 M_sun), suggesting a high accretion rate or the presence of multiple star system harboring a few hot corinos. The outflow is highly collimated, and the dynamical timescale is estimated to be less than 740 yr. In addition, we have also detected extended emission of SiO, CS, and CH3OH, which is not associated with the hot core and the outflow. This emission may be related to past star formation activity in the clump. Although G34.43+00.24 MM3 is surrounded by a dark feature in infrared, it has already experienced active formation of low-mass stars in an early stage of clump evolution.

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.

The millimeter-wave rotational transition line ($J_K = 1_0 - 0_0$) of deuterated methane CH$_3$D has tentatively been detected toward the low-mass Class 0 protostar IRAS 04368+2557 in L1527 with the Heinrich Hertz Submillimeter Telescope. This is the first detection of CH$_3$D in interstellar clouds, if confirmed. The column density and fractional abundance of CH$_3$D are determined to be $(9.1\pm3.4) \times 10^{15}$ cm$^{-2}$ and $(3.0\pm1.1) \times 10^{-7}$, respectively, where we assume the rotational temperature of 25 K. The column density and fractional abundance of the gaseous CH$_4$ are estimated to be $(1.3-4.6) \times 10^{17}$ cm$^{-2}$ and $(4.3-15.2) \times 10^{-6}$, respectively, by adopting the molecular D/H ratios of 2--7% reported for various molecules in L1527. The fractional abundance of CH$_4$ is higher than or comparable to that found in high-mass star-forming cores by infrared observations. It is sufficiently high for triggering the efficient production of various carbon-chain molecules in a lukewarm region near the protostar, which supports the picture of the warm carbon-chain chemistry.

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 results of multi-epoch very long baseline interferometry (VLBI) water (H2O) maser observations carried out with the VLBI Exploration of Radio Astrometry (VERA) toward the Cepheus A HW3d object. We measured for the first time relative proper motions of the H2O maser features, whose spatio-kinematics traces a compact bipolar outflow. This outflow looks highly collimated and expanding through ~ 280 AU (400 mas) at a mean velocity of ~ 21 km/s (~ 6 mas/yr) without taking into account the turbulent central maser cluster. The opening angle of the outflow is estimated to be ~ 30{\circ}. The dynamical time-scale of the outflow is estimated to be ~ 100 years. Our results provide strong support that HW3d harbors an internal massive young star, and the observed outflow could be tracing a very early phase of star formation. We also have analyzed Very Large Array (VLA) archive data of 1.3 cm continuum emission obtained in 1995 and 2006 toward Cepheus A. The comparative result of the HW3d continuum emission suggests the possibility of the existence of distinct young stellar objects (YSOs) in HW3d and/or strong variability in one of their radio continuum emission components.

We have observed the HN13C J=1-0 and DNC J=1-0 lines toward 18 massive clumps, including infrared dark clouds (IRDCs) and high-mass protostellar objects (HMPOs), by using the Nobeyama Radio Observatory 45 m telescope. We have found that the HN13C emission is stronger than the DNC emission toward all the observed sources. The averaged DNC/HNC ratio is indeed lower toward the observed high-mass sources (0.009\pm0.005) than toward the low-mass starless and star-forming cores (0.06). The kinetic temperature derived from the NH3 (J, K) = (1, 1) and (2, 2) line intensities is higher toward the observed high-mass sources than toward the low-mass cores. However the DNC/HNC ratio of some IRDCs involving the Spitzer 24 {\mu}m sources is found to be lower than that of HMPOs, although the kinetic temperature of the IRDCs is lower than that of the HMPOs. This implies that the DNC/HNC ratio does not depend only on the current kinetic temperature. With the aid of chemical model simulations, we discuss how the DNC/HNC ratio decreases after the birth of protostars. We suggest that the DNC/HNC ratio in star-forming cores depends on the physical conditions and history in their starless-core phase, such as its duration time and the gas kinetic temperature.

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.

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.

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 found two dense cores, CB130-3 and L673-SMM4, where the carbon-chain molecules are extremely abundant relative to NH$_{3}$, during a survey observation of radio emission lines of CCS, HC$_{3}$N, HC$_{5}$N, and NH$_{3}$ toward dark cloud cores. Judging from the low NH$_{3}$/CCS ratios, they are possible candidates for "Carbon-Chain--Producing Regions (CCPRs)" recognized as chemically young dark cloud cores. The deuterium fractionation ratios DNC/HN$^{13}$C in CB130-3 and L673-SMM4 are found to be 1.28$^{+0.27}_{-0.05}$ and 1.96$^{+0.32}_{-0.01}$, respectively, which are comparable to or slightly higher than those in CCPRs found previously. We suggest that the dense cores of CB130-3 and L673-SMM4 are analogous to CCPRs but their chemical evolutionary phase would be slightly older than those of the dense cores in the Taurus region.

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 have observed molecular cloud cores in the Orion A giant molecular cloud (GMC) in CCS, HC3N, DNC, and HN13C to study their chemical characteristics. We have detected CCS in the Orion A GMC for the first time. CCS was detected in about a third of the observed cores. The cores detected in CCS are not localized but are widely distributed over the Orion A GMC. The CCS peak intensity of the core tends to be high in the southern region of the Orion A GMC. The HC3N peak intensity of the core also tends to be high in the southern region, while there are HC3N intense cores near Orion KL, which is not seen in CCS. The core associated with Orion KL shows broad HC3N line profile, and star formation activity near Orion KL seems to enhance the HC3N emission. The column density ratio of NH3 to CCS is lower near the middle of the filament, and is higher toward the northern and southern regions along the Orion A GMC filament. This ratio is known to trace the chemical evolution in nearby dark cloud cores, but seems to be affected by core gas temperature in the Orion A GMC: cores with low NH3 to CCS column density ratios tend to have warmer gas temperature. The value of the column density ratio of DNC to HN13C is generally similar to that in dark cloud cores, but becomes lower around Orion KL due to higher gas temperature.

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 carried out an observation of the CCS ($J_{N}$=2$_{1}-$1$_{0}$) line with the Very Large Array in its D-configuration toward a protostellar core L483 (IRAS~18140$-$0440). This is a candidate source of the newly found carbon-chain rich environment called "Warm-Carbon-Chain-Chemistry (WCCC)", according to the previous observations of carbon-chain molecules. The CCS distribution in L483 is found to consist of two clumps aligned in the northwest-southeast direction, well tracing the CCS ridge observed with the single-dish radio telescope. The most remarkable feature is that CCS is depleted at the core center. Such a CCS distribution with the central hole is consistent with those of previously observed prestellar and protostellar cores, but it is rather unexpected for L483. This is because the distribution of CS, which is usually similar to that of CCS, is centrally peaked. Our results imply that the CCS ($J_{N}$=2$_{1}-$1$_{0}$) line would selectively trace the outer cold envelope in the chemically less evolved phase that is seriously resolved out with the interferometric observation. Thus, it is most likely that the high abundance of CCS in L483 relative to the other WCCC sources is not due to the activity of the protostar, although it would be related to its younger chemical evolutionary stage, or a short timescale of the prestellar phase.

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 present results of a survey of CCS, HC$_{3}$N, and HC$_{5}$N toward 40 dark cloud cores to search for "Carbon-Chain--Producing Regions(CCPRs)", where carbon-chain molecules are extremely abundant relative to NH$_{3}$, as in L1495B, L1521B, L1521E, and the cyanopolyyne peak of TMC-1. We have mainly observed toward cores where the NH$_{3}$ lines are weak, not detected, or not observed in previous surveys, and the CCS, HC$_{3}$N, and HC$_{5}$N lines have been detected toward 17, 17, and 5 sources, respectively. Among them, we have found a CCPR, L492, and its possible candidates, L1517D, L530D, L1147, and L1172B. They all show low abundance ratios of [NH$_{3}$]/[CCS] (hereafter called the NH$_{3}$/CCS ratio) indicating the chemical youth. Combining our results with those of previous surveys, we have found a significant variation of the NH$_{3}$/CCS ratio among dark cloud cores and among molecular cloud complexes. Such a variation is also suggested by the detection rates of carbon-chain molecules. For instance, the NH$_{3}$/CCS ratios are higher and the detection rates of carbon-chain molecules are lower in the Ophiuchus cores than in the Taurus cores. An origin of these systematic abundance variation is discussed in terms of the difference in the evolutionary stage or the contraction timescale. We have also identified a carbon-chain-rich star-forming core, L483, where intense HC$_{3}$N and HC$_{5}$N lines are detected. This is a possible candidate for a core with "Warm Carbon-Chain Chemistry".

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 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 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 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 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 have surveyed the N2H+ J=1-0, HC3N J=5-4, CCS J_N=4_3-3_2, NH3 (J, K) = (1, 1), (2, 2), (3, 3), and CH3OH J=7-6 lines toward the 55 massive clumps associated with infrared dark clouds by using the Nobeyama Radio Observatory 45 m telescope and the Atacama Submillimeter Telescope Experiment 10 m telescope. The N2H+, HC3N, and NH3 lines are detected toward most of the objects. On the other hand, the CCS emission is detected toward none of the objects. The [CCS]/[N2H+] ratios are found to be mostly lower than unity even in the Spitzer 24 micron dark objects. This suggests that most of the massive clumps are chemically more evolved than the low-mass starless cores. The CH3OH emission is detected toward 18 out of 55 objects. All the CH3OH-detected objects are associated with the Spitzer 24 micron sources, suggesting that star formation has already started in all the CH3OH-detected objects. The velocity widths of the CH3OH J_K=7_0-6_0 A+ and 7_{-1}-6_{-1} E lines are broader than those of N2H+ J=1-0. The CH3OH J_K=7_0-6_0 A+ and 7_{-1}-6_{-1} E lines tend to have broader linewidth in the MSX dark objects than in the others, the former being younger or less luminous than the latter. The origin of the broad emission is discussed in terms of the interaction between an outflow and an ambient cloud.

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

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 report H2O masers associated with the massive-star forming region G192.16-3.84 observed with the new Japan VLBI network at three epochs spanned for two months, which have revealed the three-dimensional kinematical structure of the whole \h2o maser region in G192.16-3.84, containing two young stellar objects separated by ~1200 AU. The maser spatio-kinematical structure has well persisted since previous observations, in which the masers are expected to be associated with a highly-collimated bipolar jet and an infalling-rotating disk in the northern and southern clusters of H2O maser features, respectively. We estimated a jet expansion speed of ~100 km/s and re-estimated a dynamical age of the whole jet to be 5.6x10^4 yrs. We have investigated the spatial distribution of Doppler velocities during the previous and present observations and relative proper motions of H2O maser features in the southern cluster, and a relative bulk motion between the two maser clusters. They are well explained by a model of an infalling-rotating disk with a radius of ~1000 AU and a central stellar mass of 5-10 M_sun, rather than by a model of a bipolar jet perpendicular to the observed CO outflow. Based on the derived H2O maser spatio-kinematical parameters, we discuss the formation mechanism of the massive young stellar objects and the outflow development in G192.16-3.84.

We report on molecular abundances and distributions in a starless dense core L492. We have found that the abundances of carbon-chain molecules such as CCS, C$_{3}$S, HC$_{3}$N, HC$_{5}$N, and HC$_{7}$N are comparable to those in chemically young dark cloud cores called "carbon-chain--producing regions", such as L1495B, L1521B, L1521E, and TMC-1. This is the first dark cloud core with extremely rich in carbon-chain-molecules that is found outside the Taurus region. In addition, the deuterium fractionation ratios of DNC/HNC and DCO$^{+}$/HCO$^{+}$ are also comparable to those in carbon-chain--producing regions, being significantly lower than those in the evolved prestellar cores such as L1498 and L1544. On the other hand, the abundances of NH$_{3}$ and N$_{2}$H$^{+}$ are systematically higher than those in carbon-chain--producing regions. Our mapping observations reveal that the central hole of molecular distributions, which were reported for CCS and C$^{34}$S in evolved prestellar cores is not significant in L492, indicating that the depletion factor of molecules is not very high. Furthermore, L492 is dynamically more evolved than carbon-chain--producing regions, and the protostellar collapse has started like L1498 and L1544. Therefore, it is likely that the chemical and dynamical evolutionary stage of L492 is intermediate between carbon-chain--producing regions (L1495B, L1521B, L1521E, and TMC-1) and evolved prestellar cores (L1498 and L1544).

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 results of phase-referencing VLBI observations of water maser sources in Orion-Monoceros Molecular Clouds with VERA (VLBI Exploration of Radio Astrometry), which is newly constructed Japanese VLBI network. Main topics of this poster are (1) the aim of one of the first scientific projects for VERA "3-Dimensional Structure and Kinematics of Orion-Monoceros Molecular Cloud Complex"; (2) current status (sensitivity and astrometric accuracy) of phase-referencing VLBI observations with VERA; and (3) results of VLBI observations of water sources in Orion-Monoceros Molecular Clouds with VERA.

We present the first comprehensive study on physical and chemical properties of quiescent starless cores L1495B and L1521B, which are known to be rich in carbon-chain molecules like the cyanopolyyne peak of TMC-1 and L1521E. We have detected radio spectral lines of various carbon-chain molecules such as CCS, C$_{3}$S, C$_{4}$H, HC$_{3}$N, and HC$_{5}$N. On the other hand, the NH$_{3}$ lines are weak and the N$_{2}$H$^{+}$ lines are not detected. According to our mapping observations of the HC$_{3}$N, CCS, and C$_{3}$S lines, the dense cores in L1495B and L1521B are compact with the radius of 0.063 and 0.044 pc, respectively, and have a simple elliptical structure. The distributions of CCS seem to be different from those of well-studied starless cores, L1498 and L1544, where the distribution of CCS shows a shell-like structure. Since the H$^{13}$CO$^{+}$, HN$^{13}$C, and C$^{34}$S lines are detected in L1495B and L1521B, the densities of these cores are high enough to excite the NH$_{3}$ and N$_{2}$H$^{+}$ lines. Therefore, the abundances of NH$_{3}$ and N$_{2}$H$^{+}$ relative to carbon-chain molecules are apparently deficient, as observed in L1521E. We found that longer carbon-chain molecules such as HC$_{5}$N and C$_{4}$H are more abundant in TMC-1 than L1495B and L1521B, while those of sulfur-bearing molecules such as C$^{34}$S, CCS, and C$_{3}$S are comparable. Both distributions and abundances of the observed molecules of L1495B and L1521B are quite similar to those of L1521E, strongly suggesting that L1495B and L1521B is in a very early stage of physical and chemical evolution.

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 present results of mapping observations of the DNC, HN^13C, and H^13CO^+ lines (J=1-0) toward 4 nearby dark cloud cores, TMC-1, L1512, L1544, and L63, along with observations of the DNC and HN^13C lines (J=2-1) toward selected positions. By use of statistical equilibrium calculations based on the LVG model, the H_2 densities are derived to be (1.4-5.5)*10^5 cm^-3, and the [DNC]/[HN^13C] ratios are derived to be 1.25-5.44 with a typical uncertainty by a factor of 2. The observed [DNC]/[HNC] ratios range from 0.02 to 0.09, assuming the [^12C]/[^13C] ratio of 60. Distributions of DNC and HN^13C are generally similar to each other, whereas the distribution of H^13CO^+ is more extended than those of DNC and HN^13C, indicating that they reside in an inner part of the cores than HCO^+. The [DNC]/[HN^13C] ratio is rather constant within each core, although a small systematic gradients are observed in TMC-1 and L63. Particularly, no such systematic gradient is found in L1512 and L1544, where a significant effect of depletion of molecules is reported toward the central part of the cores. This suggests that the [DNC]/[HNC] ratio would not be very sensitive to depletion factor, unlike the [DCO^+]/[HCO^+] ratio. On the other hand, the core to core variation of the [DNC]/[HNC] ratio, which range an order of magnitude, is more remarkable than the variation within each core. These results are interpreted qualitatively by a combination of three competing time-dependent processes; gas-phase deuterium fractionation, depletion of molecules onto grain surface, and dynamical evolution of a core.

We have studied the physical and chemical properties of a quiescent starless core L1521E with various molecular lines. It is found that there exists a compact dense core traced by the H^13CO^+, HN^13C, CCS, and HC_3N lines; their distributions have a single peak at the same position. The core radius is as small as 0.031 pc, whereas the H_2 density at the peak position is as high as (1.3-5.6)times10^5 cm^-3. Although the density is high enough to excite the inversion transitions of NH_3, these lines are found to be very faint in L1521E. The distributions of NH_3 and CCS seem to be different from those of well-studied starless cores, L1498 and L1544, where the distribution of CCS shows a shell-like structure while that of NH_3 is concentrated at the center of the core. Abundances of carbon-chain molecules are higher in L1521E than the other dark cloud cores, and especially those of sulfur-bearing molecules C_nS are comparable to the cyanopolyyne peak of TMC-1. Our results suggest that L1521E would be in a very early stage of physical and chemical evolution.