This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

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

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

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 the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded.

Star formation is thought to be driven by two groups of mechanisms; spontaneous collapse and triggered collapse. Triggered star formation mechanisms further diverge into cloud-cloud collision (CCC), "collect and collapse" (C\&C) and shock induced collapse of pre-existing, gravitationally stable cores, or 'radiation driven implosion' (RDI). To evaluate the contributions of these mechanisms and establish whether these processes can occur together within the same star forming region we performed mapping observations of radio frequency ammonia, and water maser emission lines in the S235 massive star forming region. Via spectral analyses of main, hyperfine and multi-transitional ammonia lines we explored the distribution of temperature and column density in the dense gas in the S235 and S235AB star forming region. The most remarkable result of the mapping observations is the discovery of high density gas in inter-core bridges which physically link dense molecular cores that house young proto-stellar clusters. The presence of dense gas implies the potential for future star formation within the system of cores and gas bridges. Cluster formation implies collapse and the continuous physical links, also seen in re-imaged archival CS and $^{13}$CO maps, suggests a common origin to the molecular cores housing these clusters, i.e the structure condensed from a single, larger parent cloud, brought about by the influence of a local expanding H${\rm II}$ region. An ammonia absorption feature co-locating with the center of the extended H${\rm II}$ region may be attributed to an older gas component left over from the period prior to formation of the H${\rm II}$ region. Our observations also detail known and new sites of water maser emission, highlighting regions of active ongoing star formation.

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.

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

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.

As a product of the maser monitoring program with the 26m telescope of the Hartebeesthoek Radio Astronomy Observatory (HartRAO), we present an unprecedented, contemporaneous flaring event of 10 maser transitions in hydroxyl, methanol, and water that began in 2015 January in the massive star-forming region NGC6334I in the velocity range -10 to -2 km/s. The 6.7 GHz methanol and 22.2 GHz water masers began flaring within 22 days of each other, while the 12.2 GHz methanol and 1665 MHz hydroxyl masers flared 80 and 113 days later respectively. The 1665 MHz, 6.7 GHz, and 22.2 GHz masers have all remained in their flared state for nearly 3 years. The brightest flaring components increased by factors of 66, 21, 26, and 20 in the 12.2 and 6.7 GHz methanol, 1665 MHz hydroxyl and 22.2 GHz water maser transitions respectively; some weaker components increased by up to a factor of 145. We also report new maser emission in the 1720, 6031, and 6035 MHz OH lines and the 23.1 GHz methanol line, along with the detection of only the fifth 4660 MHz OH maser. We note the correlation of this event with the extraordinary (sub)millimeter continuum outburst from the massive protostellar system NGC6334I-MM1 and discuss the implications of the observed time lags between different maser velocity components on the nature of the outburst. Finally, we identify two earlier epoch maser flaring events likely associated with this object, which suggest a recurring accretive phenomenon that generates powerful radiative outbursts.

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.

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

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

Based on sub-arcsecond Atacama Large Millimeter/submillimeter Array (ALMA) and Submillimeter Array (SMA) 1.3 mm continuum images of the massive protocluster NGC 6334I obtained in 2015 and 2008, we find that the dust emission from MM1 has increased by a factor of 4.0$\pm$0.3 during the intervening years, and undergone a significant change in morphology. The continuum emission from the other cluster members (MM2, MM4 and the UCHII region MM3 = NGC 6334F) has remained constant. Long term single-dish maser monitoring at HartRAO finds that multiple maser species toward NGC 6334I flared beginning in early 2015, a few months before our ALMA observation, and some persist in that state. New ALMA images obtained in 2016 July-August at 1.1 and 0.87 mm confirm the changes with respect to SMA 0.87 mm images from 2008, and indicate that the (sub)millimeter flaring has continued for at least a year. The excess continuum emission, centered on the hypercompact HII region MM1B, is extended and elongated ($1.6" \times 1.0" \approx 2100 \times 1300$ au) with multiple peaks, suggestive of general heating of the surrounding subcomponents of MM1, some of which may trace clumps in a fragmented disk rather than separate protostars. In either case, these remarkable increases in maser and dust emission provide direct observational evidence of a sudden accretion event in the growth of a massive protostar yielding a sustained luminosity surge by a factor of $70\pm20$, analogous to the largest events in simulations by Meyer et al. (2017). This target provides an excellent opportunity to assess the impact of such a rare event on a protocluster over many years.

We investigated the star formation activities in the AFGL333 region, which is in the vicinity of the W4 expanding bubble, by conducting NH3 (1,1), (2,2), and (3,3) mapping observations with the 45 m Nobeyama Radio Telescope at an angular resolution of 75". The morphology of the NH3 (1,1) map shows a bow-shape structure with the size of 2.0 x 0.6 pc as seen in the dust continuum. At the interface between the W4 bubble and the dense NH3 cloud, the compact HII region G134.2+0.8, associated with IRAS02245+6115, is located. Interestingly, just north and south of G134.2+0.8 we found NH3 emission exhibiting large velocity widths of ~ 2.8 km/s, compared to 1.8 km/s at the other positions. As the possibility of mechanical energy injection through the activity of YSO(s) is low, we considered the origin of the large turbulent gas motion as indication of interaction between the compact HII region and the periphery of the dense molecular cloud. We also found expanding motion of the CO emission associated with G134.2+0.8. The overall structure of the AFGL333-Ridge might have been formed by the expanding bubble of W4. However, the small velocity widths observed west of IRAS02245+6115, around the center of the dense molecular cloud, suggest that interaction with the compact HII region is limited. Therefore the YSOs (dominantly Class 0/I) in the core of the AFGL333-Ridge dense molecular cloud most likely formed in quiescent mode. As has been previously suggested for the large scale star formation in the W3 giant molecular cloud, our results show an apparent coexistence of induced and quiescent star formation in this region. It appears that star formation in the AFGL333 region has proceeded without significant external triggers, but accompanying stellar feedback environment.

The scaling relations and the star formation laws for molecular cloud complexes in the Milky Way is investigated. We compare their masses $M_{\rm gas}$, mass surface densities $\Sigma_{M_{\rm gas}}$, radii $R$, velocity dispersions $\sigma$, star formation rates $SFR$, and SFR densities $\Sigma_{\rm SFR}$ with those of structures ranging from cores, clumps, Giant Molecular Clouds (GMCs), to Molecular Cloud Complexes (MCCs), and to Galaxies, spanning 8 orders of magnitudes in size and 13 orders of magnitudes in mass. MCC are mostly large ($R>50$ pc), massive ($\sim 10^{6}$\,\msun) gravitationally unbound cloud structures. This results in the following universal relations: $\sigma\sim R^{0.5}$, $M_{\rm gas}\sim R^{2}$, $\Sigma_{\rm SFR}\sim \Sigma_{M_{\rm gas}}^{1.5}$, ${SFR}\sim {M_{\rm gas}}^{0.9}$, and ${SFR}\sim {\sigma}^{2.7}$. Variations in the slopes and the coefficients of these relations are found at individual scales signifying different physics acting at different scales. Additionally, there are breaks at the MCC scale in the $\sigma-R$ relation and between the starburst and the normal star-forming objects in the $SFR-M_{\rm gas}$ and $\Sigma_{\rm SFR}$-$\Sigma_{\rm M_{\rm gas}}$ relations. We propose to use the Schmidt-Kennicutt diagram to distinguish the starburst from the normal star-forming structures by applying a $\Sigma_{M_{\rm gas}}$ threshold of $\sim100$\,\msun pc$^{-2}$ and a $\Sigma_{\rm SFR}$ threshold of 1\,\msun yr$^{-1}$ kpc$^{-2}$. Mini-starburst complexes have enhanced $\Sigma_{\rm SFR}$ ($>$1\,\msun yr$^{-1}$ kpc$^{-2}$), probably caused by dynamic events such as radiation pressure, colliding flows, or spiral arm gravitational instability. Because of the dynamical evolution, gravitational boundedness does not play a significant role in characterizing the star formation activity of MCCs, especially the mini-starburst complexes.

We present the results of SiO (2-1) and SO2 (12-13) line observations of Sgr B2(N) made with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~2arcsec. Our analysis of the SiO and SO2 line emission reveals a bipolar molecular outflow in an east-west direction whose driving source is located at K2. In addition, SO2 line core shows a north-south velocity gradient most probably indicating a hot core of molecular gas rotating around K2. Fractional abundances of SO2 and SiO (X(SO2) and X(SiO), respectively) in the outflowing molecular gas are derived from comparisons with the C18O emission. Assuming an excitation temperature of 100+-50 K, we calculate X(SO2) = 2.3X10^-8 and X(SiO) = 1.2X10^-9. The outflow from SgrB2(N) K2 is characterized as a young (5X10^3yr) and massive (~2000Msun), but moderately collimated (~60deg) outflow. We also report a possible detection of the SiO (v=2, J=2-1) maser emission from the position of K2. If confirmed, it would make Sgr B2(N) the 4th star forming region associated with SiO masers.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present the results of continuum and 12CO(3-2) and CH3OH(7-6) line observations of IRAS16547-4247 made with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~0.5". The 12CO(3-2) emission shows two high-velocity outflows whose driving sources are located within the dust continuum peak. The alignment of these outflows do not coincide with that of the wide-angle, large scale, bipolar outflow detected with APEX in previous studies. The CH3OH(7-6) line emission traces an hourglass structure associated with the cavity walls created by the outflow lobes. Taking into account our results together with the position of the H2O and class I CH3OH maser clusters, we discuss two possible scenarios that can explain the hourglass structure observed in IRAS16547-4247: (1) precession of a biconical jet, (2) multiple, or at least two, driving sources powering intersecting outflows. Combining the available evidence, namely, the presence of two cross-aligned bipolar outflows and two different H2O maser groups, we suggest that IRAS16547-4247 represents an early formation phase of a protocluster.

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

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

We report the results of VLBI observations of H$_{2}$O masers in the IRAS 20143+3634 star forming region using VERA (VLBI Exploration of Radio Astronomy). By tracking masers for a period of over two years we measured a trigonometric parallax of $\pi = 0.367 \pm 0.037$ mas, corresponding to a source distance of $D = 2.72 ^{+0.31}_{-0.25}$ kpc and placing it in the Local spiral arm. Our trigonometric distance is just 60% of the previous estimation based on radial velocity, significantly impacting the astrophysics of the source. We measured proper motions of $-2.99 \pm 0.16$ mas yr$^{-1}$ and $-4.37 \pm 0.43$ mas yr$^{-1}$ in R.A. and Decl. respectively, which were used to estimate the peculiar motion of the source as $(U_{s},V_{s},W_{s}) = (-0.9 \pm 2.9, -8.5 \pm 1.6, +8.0 \pm 4.3)$ km s$^{-1}$ for $R_0=8$ kpc and $\Theta_0=221$ km s$^{-1}$, and $(U_{s},V_{s},W_{s}) = (-1.0 \pm 2.9, -9.3 \pm 1.5, +8.0 \pm 4.3)$ km s$^{-1}$ for $R_0=8.5$ kpc and $\Theta_0=235$ km s$^{-1}$. IRAS 20143+3634 was found to be located near the tangent point in the Cygnus direction. Using our observations we derived the angular velocity of Galactic rotation of the local standard of rest (LSR), $\Omega_{0} = 27.3 \pm 1.6$ km s$^{-1}$ kpc$^{-1}$, which is consistent with previous values derived using VLBI astrometry of SFRs at the tangent points and Solar circle. It is higher than the value recommended by the IAU of $\Omega_{0} = 25.9$ km s$^{-1}$ kpc$^{-1}$ which was calculated using the Galactocentric distance of the Sun and circular velocity of the LSR.

We present the results of sulfur monoxide, SO line emission observations of G0.253+0.016 with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of 1.7 arcsec. The dense and massive molecular cloud of G0.253+0.016 is highly sub-structured, yet shows no obvious signs of cluster formation. We found three outstanding features of the cloud from the SO emission, namely, shell structure of radius 1.3 pc, large velocity gradients of 20 km/s/pc with the cloud, and cores with large velocity dispersions (30-40 km/s) around the shell structure. We suggest that these large-velocity dispersion cores will form high-mass stars in the future. In attempt to explore the formation scenario of the dense cloud, we compared our results with numerical simulations, thus, we propose that G0.253+0.016 may have formed due to a cloud-cloud collision process.

Planetary nebulae (PNe) have diverse morphological shapes, including point-symmetric and multipolar structures. Many PNe also have complicated internal structures such as torus, lobes, knots, and ansae. A complete accounting of all the morphological structures through physical models is difficult. A first step toward such an understanding is to derive the true three-dimensional structure of the nebulae. In this paper, we show that a multipolar nebula with three pairs of lobes can explain many of such features, if orientation and sensitivity effects are taken into account. Using only six parameters - the inclination and position angles of each pair - we are able to simulate the observed images of 20 PNe with complex structures. We suggest that the multipolar structure is an intrinsic structure of PNe and the statistics of multipolar PNe has been severely underestimated in the past.

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.