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.

Ring-like distributions of the 6.7 GHz methanol maser spots at milliarcsecond scales represent a family of molecular structures of unknown origin associated with high-mass young stellar objects (HMYSOs). We aim to study G23.657-00.127, which has a nearly circular ring of the 6.7 GHz methanol masers, and is the most suitable target to test hypotheses on the origin of the maser rings. The European Very Long Baseline Interferometry Network (EVN) was used at three epochs spanning 10.3 yr to derive the spatio-kinematical structure of the 6.7 GHz methanol maser emission in the target. The maser cloudlets, lying in a nearly symmetric ring, expand mainly in the radial direction with a mean velocity of 3.2 km s$^{-1}$. There is an indication that the radial component of the velocity increases with cloudlet's distance from the ring centre. The tangential component does not show any clear evidence for rotation of the cloudlets or any relationship with distance from the ring centre. The blue-shifted masers may hint at an anticlockwise rotation of cloudlets in the southern part of the ring. The nearly circular structure of the ring clearly persisted for more than 10 yr. Interferometric data demonstrated that about one quarter of cloudlets show significant variability in their brightness, although the overall spectrum was non-variable in single-dish studies. Taking into account the three-dimensional motion of the maser cloudlets and their spatial distribution along a small ring, we speculate about two possible scenarios where the methanol masers trace either a spherical outflow arising from an (almost) edge-on disc, or a wide angle wind at the base of a protostellar jet. The latter is associated with near- and mid-infrared emission detected towards the ring. High angular resolution images of complementary (thermal) tracers are needed to interpret the environment of methanol masers.

Periodic flares of maser emission are thought to be induced either by variations of the seed photon flux in young binary systems or the pump rate regulated by stellar and accretion luminosities. We seek to study the variability of four maser transitions of three different species in G107.298+5.639 to constrain the dominant mechanism of periodic flares. Light curves of the 6.7 GHz methanol and 22.2 GHz water vapour maser were obtained with the Torun 32m radio telescope over 39 and 34 cycles, respectively. The target was also monitored at the 1.6 GHz hydroxyl transitions with the Nan\c{c}ay radio telescope over 13 cycles. All these maser lines were imaged using VLBI arrays. The study confirms alternating flares of the methanol and water masers with a period of 34.4d and reveals the synchronised behaviour of the methanol and hydroxyl masers in this source. The observed spatial distribution of the methanol maser cloudlets and the measured time delays of the flares of individual features imply a ring-like structure of radius 240au and thickness 30au. Internal proper motions indicate that the velocity of methanol cloudlets is dominated by a disc-wind component of about 5km/s. The methanol emission detected during only one VLBI observation is located in a region about 550au from a central star, which also exhibits OH maser flares. The erratic appearance of methanol features can be related to a powering object of relatively low luminosity which, during some variability cycles, can excite molecules only in the nearest part of the disc. A careful analysis of the maser and infrared light curves reveal a strong correlation between the 6.7 GHz line and the infrared flux densities supporting a radiative pumping of the maser. The synchronised behaviour of the hydroxyl and methanol transitions indicates a common pumping mechanism for the periodic flares of G107.298+5.639.

We compile and analyze ~200 trigonometric parallaxes and proper motions of molecular masers associated with very young high-mass stars. These measurements strongly suggest that the Milky Way is a four-arm spiral. Fitting log-periodic spirals to the locations of the masers, allows us to significantly expand our view of the structure of the Milky Way. We present an updated model for its spiral structure and incorporate it into our previously published parallax-based distance-estimation program for sources associated with spiral arms. Modeling the three-dimensional space motions yields estimates of the distance to the Galactic center, Ro = 8.15 +/- 0.15 kpc, the circular rotation speed at the Sun's position, To = 236 +/- 7 km/s, and the nature of the rotation curve. Our data strongly constrain the full circular velocity of the Sun, To + Vsun = 247 +/- 4 km/s, and its angular velocity, (To + Vsun)/Ro = 30.32 +/- 0.27 km/s/kpc. Transforming the measured space motions to a Galactocentric frame which rotates with the Galaxy, we find non-circular velocity components typically about 10 km/s. However, near the Galactic bar and in a portion of the Perseus arm, we find significantly larger non-circular motions. Young high-mass stars within 7 kpc of the Galactic center have a scale height of only 19 pc and, thus, are well suited to define the Galactic plane. We find that the orientation of the plane is consistent with the IAU-defined plane to within +/-0.1 deg., and that the Sun is offset toward the north Galactic pole by Zsun = 5.5 +/- 5.8 pc. Accounting for this offset places the central supermassive black hole, Sgr A*, in the midplane of the Galaxy. Using our improved Galactic parameters, we predict the Hulse-Taylor binary pulsar to be at a distance of 6.54 +/- 0.24 kpc, assuming its orbital decay from gravitational radiation follows general relativity.

Discovery of periodic maser emission was an unexpected result from monitoring observations of methanol transitions in high-mass young stellar objects. We report on the detection of five new periodic sources from a monitoring program with the Torun 32 m telescope. Variability with a period of 149 to 540 d and different patterns from sinusoidal-like to intermittent was displayed. Three dimensional structure of G59.633$-$0.192 determined from the time delays of burst peaks of the spectral features and high angular resolution map implies that the emission traces a disk. For this source the 6.7 GHz light curve followed the infrared variability supporting a radiative scheme of pumping. An unusual time delay of $\sim$80 d occurred in G30.400$-$0.296 could not be explained by the light travel time and may suggest a strong differentiation of physical conditions and excitation in this deeply embedded source. Our observations suggest the intermittent variability may present a simple response of maser medium to the underlying variability induced by the accretion luminosity while other variability patterns may reflect more complex changes in the physical conditions.

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.

MHD simulations show that the magnetic field can drive molecular outflows during the formation of massive protostars. The best probe to observationally measure both the morphology and the strength of this magnetic field at scales of 10-100 au is maser polarization. We measure the direction of magnetic fields at mas resolution around a sample of massive star-forming regions to determine whether there is a relation between the orientation of the magnetic field and of the outflows. In addition, by estimating the magnetic field strength via the Zeeman splitting measurements, the role of magnetic field in the dynamics of the massive star-forming region is investigated. We selected a flux-limited sample of 31 massive star-forming regions to perform a statistical analysis of the magnetic field properties with respect to the molecular outflows characteristics. We report the linearly and circularly polarized emission of 6.7 GHz CH3OH masers towards seven massive star-forming regions of the total sample with the EVN. The sources are: G23.44-0.18, G25.83-0.18, G25.71-0.04, G28.31-0.39, G28.83-0.25, G29.96-0.02, and G43.80-0.13. We identified a total of 219 CH3OH maser features, 47 and 2 of which showed linearly and circularly polarized emission, respectively. We measured well-ordered linear polarization vectors around all the massive YSOs and Zeeman splitting towards G25.71-0.04 and G28.83-0.25. Thanks to recent theoretical results, we were able to provide lower limits to the magnetic field strength from our Zeeman splitting measurements. We further confirm (based on ~80% of the total flux-limited sample) that the magnetic field on scales of 10-100 au is preferentially oriented along the outflow axes. The estimated magnetic field strength of |B_{||}|>61 mG and >21 mG towards G25.71-0.04 and G28.83-0.2, respectively, indicates that it dominates the dynamics of the gas in both regions.

We present new evidence for superradiance in the methanol 6.7 GHz spectral line for three different star-forming regions: S255IR-NIRS3, G24.329+0.144, and Cepheus A. Our analysis shows that some of the intensity flares exhibiting fast flux rise times and asymmetric light curves reported in these sources can naturally be explained within the context of superradiance. When a threshold for the inverted population column density is exceeded in a maser-hosting region, the radiation mode switches from one regulated by stimulated emission (maser) to superradiance. Superradiance, as a more efficient energy release mechanism, manifests itself through strong bursts of radiation emanating from spatially compact regions. Elevated inverted population densities and the triggering of superradiance can be due to a change in radiative pumping. Here, we show that an increase in the pump rate and the inverted population density of only a factor of a few results in a significant increase in radiation. While the changes in the pump rate can take place over a few hundred days, the rise in radiation flux density when superradiance is initiated is drastic and happens over a much shorter time-scale.

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.

Context: High-mass young stellar objects (HMYSOs) can undergo accretion episodes that strongly affect the star evolution, the dynamics of the disk, and its chemical evolution. Recently reported extraordinary bursts in the methanol maser emission may be the observational signature of accretion events in deeply embedded HMYSOs. Aims: We analyze the light curve of 6.7 GHz methanol masers in S255IR-NIRS3 during the 2015-2016 burst. Methods: 8.5-year monitoring data with an average sampling interval of 5 days were obtained with the Torun 32 m radio telescope. Archival data were added, extending the time series to ~27 years. Results: The maser emission showed moderate (25-30%) variability on timescales of months to years over ~23 years since its discovery. The main burst was preceded by a one-year increase of the total flux density by a factor of 2.5, then it grew by a factor of 10 over ~0.4 years and declined by a factor of 8 during the consecutive 2.4 years. The peak maser luminosity was a factor of 24.5 higher than the pre-burst quiescent value. The light curves of individual features showed considerable diversity but indicated a general trend of suppression of the maser emission at blueshifted (<4.7 km s$^{-1}$) velocities when the redshifted emission rapidly grew and new emission features appeared at velocities >5.8 km s$^{-1}$. This new emission provided a contribution of about 80% to the maser luminosity around the peak of the burst. Conclusions: The onset of the maser burst exactly coincides with that of the infrared burst estimated from the motion of the light echo. This strongly supports the radiative pumping scheme of the maser transition. The growth of the maser luminosity is the result of an increasing volume of gas where the maser inversion is achieved.

We report results of 6.7 GHz methanol maser monitoring of 139 star-forming sites with the Torun 32 m radio telescope from June 2009 to February 2013. The targets were observed at least once a month, with higher cadences of 2-4 measurements per week for circumpolar objects. Nearly 80 percent of the sources display variability greater than 10 per cent on a time-scale between a week and a few years but about three quarters of the sample have only 1-3 spectral features which vary significantly. Irregular intensity fluctuation is the dominant type of variability and only nine objects show evidence for cyclic variations with periods of 120 to 416 d. Synchronised and anti-correlated variations of maser features are detected in four sources with a disc-like morphology. Rapid and high amplitude bursts of individual features are seen on 3-5 occasions in five sources. Long (>50 d to 20 months) lasting bursts are observed mostly for individual or groups of features in 19 sources and only one source experienced a remarkable global flare. A few flaring features display a strong anti-correlation between intensity and line-width that is expected for unsaturated amplification. There is a weak anti-correlation between the maser feature luminosity and variability measure, i.e. maser features with low luminosity tend to be more variable than those with high luminosity. The analysis of the spectral energy distribution and continuum radio emission reveals that the variability of the maser features increases when the bolometric luminosity and Lyman flux of the exciting object decreases. Our results support the concept of a major role for infrared pumping photons in triggering outburst activity of maser emission.

The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.

We conducted VLA C-configuration observations to measure positions and luminosities of Galactic Class II 6.7 GHz methanol masers and their associated ultra-compact H II regions. The spectral resolution was 3.90625 kHz and the continuum sensitivity reached 45 \uJypb. We mapped 372 methanol masers with peak flux densities of more than 2 Jy selected from the literature, 367 of them were detected. Absolute positions have nominal uncertainties of 0.3 arcsec. In this first paper on the data analysis, we present three catalogs, the first gives information on the strongest feature of 367 methanol maser sources, and the second on all detected maser spots. The third catalog present derived data of the 279 radio continuum sources found in the vicinity of maser sources. Among them, 140 show evidence of physical association with maser sources. Our catalogs list properties including distance, flux density, radial velocity and the distribution of masers on the Galactic plane is then provided as well. We found no significant relationship between luminosities of masers and their associated radio continuum counterparts.

Methanol and water vapour masers are signposts of early stages of high-mass star formation but it is generally thought that due to different excitation processes they probe distinct parts of stellar environments. Here we present observations of the intermediate-mass young stellar object G107.298+5.639, revealing for the first time that 34.4 d flares of the 6.7 GHz methanol maser emission alternate with flares of individual features of the 22 GHz water maser. High angular resolution data reveal that a few components of both maser species showing periodic behaviour coincide in position and velocity and all the periodic water maser components appear in the methanol maser region of size of 360 au. The maser flares could be caused by variations in the infrared radiation field induced by cyclic accretion instabilities in a circumstellar or protobinary disc. The observations do not support either the stellar pulsations or the seed photon flux variations as the underlying mechanisms of the periodicity in the source.

Methanol masers at 6.7 GHz are well known tracers of high-mass star-forming regions. However, their origin is still not clearly understood. We aimed to determine the morphology and velocity structure for a large sample of the maser emission with generally lower peak flux densities than those in previous surveys. Using the European VLBI Network we imaged the remaining sources (17) from a sample of sources that were selected from the unbiased survey using the Torun 32 m dish. Together they form a database of a total of 63 source images with high sensitivity, milliarcsecond angular resolution and very good spectral resolution for detailed studies. We studied in detail the properties of the maser clouds and calculated the mean and median values of the projected size (17.4 au and 5.5 au, respectively) as well as the FWHM of the line (0.373 km s$^{-1}$ and 0.315 km s$^{-1}$ for the mean and median values, respectively), testing whether it was consistent with Gaussian profile. We also found maser clouds with velocity gradients (71 per cent) that ranged from 0.005 km s$^{-1}$ au$^{-1}$ to 0.210 km s$^{-1}$ au$^{-1}$. We tested the kinematic models to explain the observed structures of the 6.7 GHz emission. There were targets where the morphology supported the scenario of a rotating and expanding disk or a bipolar outflow. Comparing the interferometric and single-dish spectra we found that, typically, 50-70 per cent of the flux was missing. This phenomena is not strongly related to the distance of the source. The EVN imaging reveals that in the complete sample of 63 sources the ring-like morphology appeared in 17 per cent of sources, arcs were seen in a further 8 per cent, and the structures were complex in 46 per cent cases. The UC HII regions coincide in position in the sky for 13 per cent of the sources. They are related both to extremely high and low luminosity masers from the sample.

Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented along the outflow axis or randomly. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the orientation of the outflow axis or not. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly polarized emission of 6.7 GHz methanol masers around a sample of massive star-forming regions. This paper focuses on the first seven observed sources, G24.78+0.08, G25.65+1.05, G29.86-0.04, G35.03+0.35, G37.43+1.51, G174.20-0.08, and G213.70-12.6. For all these sources, molecular outflows have been detected in the past. We detected a total of 176 methanol masing cloudlets toward the seven massive star-forming regions, 19% of which show linearly polarized emission. The methanol masers around the massive young stellar object MM1 in G174.20-0.08 show neither linearly nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We measured significant Zeeman splitting toward both A1 and A2 in G24.78+0.08, and toward G29.86-0.04 and G213.70-12.6. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field (on scales 10-100 au) is preferentially oriented along the outflow axes.

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields can play a crucial role not only in the formation and dynamics of molecular outflows, but also in the evolution of circumstellar disks. Therefore, new measurements of magnetic fields at milliarcsecond resolution close to massive young stellar objects (YSOs) are fundamental for providing new input for numerical simulations and for understanding the formation process of massive stars. The polarized emission of 6.7 GHz CH3OH masers allows us to investigate the magnetic field close to the massive YSO where the outflows and disks are formed. Recently, we have detected with the EVN CH3OH maser polarized emission towards 10 massive YSOs. From a first statistical analysis we have found evidence that magnetic fields are primarily oriented along the molecular outflows. To improve our statistics we are carrying on a large observational EVN campaign for a total of 19 sources, the preliminary results of the first seven sources are presented in this contribution. Furthermore, we also describe our efforts to estimate the Lande' g-factors of the CH3OH maser transition to determine the magnetic field strength from our Zeeman-splitting measurements.

We report the discovery of 6.7 GHz methanol maser periodic flares in four massive star forming regions and the updated light curve for the known periodic source G22.357+0.066. The observations were carried out with the Torun 32 m radio telescope between June 2009 and April 2014. Flux density variations with period of 120 to 245 d were detected for some or all spectral features. A variability pattern with a fast rise and relatively slow fall on time-scale of 30-60 d dominated. A reverse pattern was observed for some features of G22.357+0.066, while sinusoidal-like variations were detected in G25.411+0.105. A weak burst lasting ~520 d with the velocity drift of 0.24 km/s/yr occurred in G22.357+0.066. For three sources for which high resolution maps are available, we found that the features with periodic behaviour are separated by more than 500 au from those without any periodicity. This suggests that the maser flares are not triggered by large-scale homogeneous variations in either the background seed photon flux or the luminosity of the exciting source and a mechanism which is able to produce local changes in the pumping conditions is required.

Methanol masers at 6.7 GHz are associated with high-mass star-forming regions (HMSFRs) and often have mid-infrared (MIR) counterparts characterized by extended emission at 4.5 $\mu$m, which likely traces outflows from massive young stellar objects (MYSOs). Our objectives are to determine the milliarcsecond (mas) morphology of the maser emission and to examine if it comes from one or several candidate MIR counterparts in the clusters of MYSOs. The European VLBI Network (EVN) was used to image the 6.7 GHz maser line with ~2.'1 field of view toward 14 maser sites from the Torun catalog. Quasi-simultaneous observations were carried out with the Torun 32 m telescope. We obtained maps with mas angular resolution that showed diversity of methanol emission morphology: a linear distribution (e.g., G37.753-00.189), a ring-like (G40.425+00.700), and a complex one (e.g., G45.467+00.053). The maser emission is usually associated with the strongest MIR counterpart in the clusters; no maser emission was detected from other MIR sources in the fields of view of 2.'1 in diameter. The maser source luminosity seems to correlate with the total luminosity of the central MYSO. Although the Very Long Baseline Interferometry (VLBI) technique resolves a significant part of the maser emission, the morphology is still well determined. This indicates that the majority of maser components have compact cores.

6.7 GHz methanol maser emission from the well-studied star-forming region Cepheus A was monitored with the Torun 32 m radio telescope. We found synchronized and anticorrelated changes of the flux density of the two blueshifted and one redshifted maser features for ~30 per cent of 1340 d of our observations. Two of those features exhibited high amplitude flux density variations with periods of 84-87 d over the last 290 d interval of the monitoring. We also report on two flares of emission at two different redshifted velocities completely covered during the whole outburst. These flare events lasted 510-670 d and showed a very rapid linear rise and slow exponential decline, which may be caused by variability of the seed flux density. The flux density of the two strongest features dropped by a factor of 2-5 on a time-scale ~22 yr, while other features have not changed significantly during this period, but showed strong variability on time-scales <5 yr.

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.

Milliarcsecond VLBI maps of regions containing 6.7 GHz methanol maser emission have lead to the recent discovery of ring-like distributions of maser spots and the plausible hypothesis that they may be tracing circumstellar disks around forming high mass stars. We aimed to test this hypothesis by imaging these regions in the near and mid-infrared at high spatial resolution and compare the observed emission to the expected infrared morphologies as inferred from the geometries of the maser rings. In the near infrared we used the Gemini North adaptive optics system of Altair/NIRI, while in the mid-infrared we used the combination of the Gemini South instrument T-ReCS and super-resolution techniques. Resultant images had a resolution of approximately 150 mas in both the near-infrared and mid-infrared. We discuss the expected distribution of circumstellar material around young and massive accreting (proto)stars and what infrared emission geometries would be expected for the different maser ring orientations under the assumption that the masers are coming from within circumstellar disks. Based upon the observed infrared emission geometries for the four targets in our sample and the results of SED modeling of the massive young stellar objects associated with the maser rings, we do not find compelling evidence in support of the hypothesis that methanol masers rings reside in circumstellar disks.

Maser emission plays an important role as a tool in star formation studies. It is widely used for deriving kinematics, as well as the physical conditions of different structures, hidden in the dense environment very close to the young stars, for example associated with the onset of jets and outflows. We will summarize the recent observational and theoretical progress on this topic since the last maser symposium: the IAU Symposium 242 in Alice Springs.

The 22.2 GHz water masers are often associated with the 6.7 GHz methanol masers but owing to the different excitation conditions they likely probe independent spatial and kinematic regions around the powering young massive star. We compared the emission of these two maser species on milliarcsecond scales to determine in which structures the masers arise and to test a disc-outflow scenario where the methanol emission arises in a circumstellar disc while the water emission comes from an outflow. We obtained high-angular and spectral resolution 22.2 GHz water maser observations of the two sources G31.581+00.077 and G33.641-00.228 using the EVN. In both objects the water maser spots form complex and filamentary structures of sizes 18-160 AU. The emission towards the source G31.581+00.077 comes from two distinct regions of which one is related to the methanol maser source of ring-like shape. In both targets the main axis of methanol distribution is orthogonal to the water maser distribution. Most of water masers appear to trace shocks on a working surface between an outflow/jet and a dense envelope. Some spots are possibly related to the disc-wind interface which is as close as 100-150 AU to the regions of methanol emission.

We report the discovery of periodic flares of 6.7GHz methanol maser in the young massive stellar object G22.357+0.066. The target was monitored in the methanol maser line over 20 months with the Torun 32m telescope. The emission was also mapped at two epochs using the EVN. The 6.7GHz methanol maser shows periodic variations with a period of 179 days. The periodic behavior is stable for the last three densely sampled cycles and has even been stable over ~12 years, as the archival data suggest. The maser structure mapped with the EVN remains unchanged at two epochs just at the putative flare maxima separated by two years. The time delays of up to ~16 days seen between maser features are combined with the map of spots to construct the 3-dimensional structure of the maser region. The emission originating in a single ~100 AU layer can be modulated by periodic changes in the infrared pumping radiation or in the free-free background emission from an HII region.

We present the results of a survey of radio continuum sources near the Galactic plane using the Very Long Baseline Array (VLBA). Our observations are designed to identify compact extragalactic sources of milliarcsecond size that can be used for parallax measurements in the Bar and Spiral Structure Legacy Survey. We selected point sources from the NVSS and CORNISH catalogs with flux densities above 30 mJy and within $1.5\degr$ of known maser targets. Of the 1529 sources observed, 199 were detected. For sources detected on 3 or more baselines, we determined accurate positions and evaluated their quality as potential calibrators. Most of the 1330 sources that were not detected with the VLBA are probably of extragalactic origin.

Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity ($\Theta_0$) and curve and the distance to the Galactic center (R$_0$). Here, we present an update on our first results, implementing a recent change in the knowledge about the Solar motion. It seems unavoidable that the IAU recommended values for R$_0$ and $\Theta_0$ need a substantial revision. In particular the combination of 8.5 kpc and 220 \kms\, can be ruled out with high confidence. Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr\,A* gives best values of R$_0$ = 8.3 $\pm$ 0.23 kpc and $\Theta_0$ = 239 or 246 $\pm$ 7 \kms, for Solar motions of V$_ \odot$ = 12.23 and 5.25 \kms, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) Survey.

In recent studies of methanol masers, a substantial fraction of the objects show maser components aligned in large-scale elliptical configurations. These can be readily interpreted as rings centred on a high mass star in formation, seen in projection. Remarkably, most of these rings do not show signs of rotation, but rather the radial motions dominate. This must mean that their dynamics are governed by other than gravitational forces. In particular, we have studied the methanol masers around Cep A in detail, where it can be argued that the methanol masers show signs of infall. In this paper we discuss the dynamics of the Cep A methanol maser and sources from the Torun blind survey to argue that at least in a fraction of sources methanol masers could be associated with the shock interface between the large scale accretion, regulated by the magnetic field, and a 1000-AU scale circumstellar disk. We discuss the validity of such a model for the overall population of methanol maser sources.

We present the results of studies of a new class of 6.7 GHz methanol maser sources with a ring-like emission structure discovered recently with the EVN. We have used the VLA to search for water masers at 22 GHz and radio continuum at 8.4 GHz towards a sample of high-mass star forming regions showing a ring-like distribution of methanol maser spots. Using the Gemini telescopes we found mid-infrared (MIR) counterparts of five methanol rings with a resolution of 0."15. The centres of methanol maser rings are located within, typically, only 0."2 of the MIR emission peak, implying their physical relation with a central star. These results strongly support a scenario wherein the ring-like structures appear at the very early stage of massive star formation before either water-maser outflows or H II regions are seen.

22 GHz water and 6.7 GHz methanol masers are usually thought as signposts of early stages of high-mass star formation but little is known about their associations and the physical environments they occur in. The aim was to obtain accurate positions and morphologies of the water maser emission and relate them to the methanol maser emission recently mapped with Very Long Baseline Interferometry. A sample of 31 methanol maser sources was searched for 22 GHz water masers using the VLA and observed in the 6.7 GHz methanol maser line with the 32 m Torun dish simultaneously. Water maser clusters were detected towards 27 sites finding 15 new sources. The detection rate of water maser emission associated with methanol sources was as high as 71%. In a large number of objects (18/21) the structure of water maser is well aligned with that of the extended emission at 4.5 $\mu$m confirming the origin of water emission from outflows. The sources with methanol emission with ring-like morphologies, which likely trace a circumstellar disk/torus, either do not show associated water masers or the distribution of water maser spots is orthogonal to the major axis of the ring. The two maser species are generally powered by the same high-mass young stellar object but probe different parts of its environment. The morphology of water and methanol maser emission in a minority of sources is consistent with a scenario that 6.7 GHz methanol masers trace a disc/torus around a protostar while the associated 22 GHz water masers arise in outflows. The majority of sources in which methanol maser emission is associated with the water maser appears to trace outflows. The two types of associations might be related to different evolutionary phases.

We investigate which structures the 6.7 GHz methanol masers trace in the environment of high-mass protostar candidates by observing a homogenous sample of methanol masers selected from Torun surveys. We also probed their origins by looking for associated H II regions and IR emission. We selected 30 methanol sources with improved position accuracies achieved using MERLIN and another 3 from the literature. We imaged 31 of these using the European VLBI Network's expanded array of telescopes with 5-cm (6-GHz) receivers. We used the VLA to search for 8.4 GHz radio continuum counterparts and inspected Spitzer GLIMPSE data at 3.6-8 um from the archive. High angular resolution images allowed us to analyze the morphology and kinematics of the methanol masers in great detail and verify their association with radio continuum and mid-infrared emission. A new class of "ring-like" methanol masers in star--forming regions appeared to be suprisingly common, 29 % of the sample. The new morphology strongly suggests that methanol masers originate in the disc or torus around a proto- or a young massive star. However, the maser kinematics indicate the strong influence of outflow or infall. This suggests that they form at the interface between the disc/torus and a flow. This is also strongly supported by Spitzer results because the majority of the masers coincide with 4.5 um emission to within less than 1 arcsec. Only four masers are associated with the central parts of UC H II regions. This implies that 6.7 GHz methanol maser emission occurs before H II region observable at cm wavelengths is formed.

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.

Methanol masers are associated with young high-mass stars and are an important tool for investigating the process of massive star formation. The recently discovered methanol maser ring in G23.657-00.127 provides an excellent ``laboratory'' for a detailed study of the nature and physical origin of methanol maser emission, as well as parallax and proper motion measurements. Multi-epoch observations of the 12.2 GHz methanol maser line from the ring were conducted using the Very Long Baseline Array. Interferometric observations with milliarcsecond resolution enabled us to track single maser spots in great detail over a period of 2 years. We have determined the trigonometric parallax of G23.657-00.127 to be 0.313+/-0.039 mas, giving a distance of 3.19{+0.46}{-0.35} kpc. The proper motion of the source indicates that it is moving with the same circular velocity as the LSR, but it shows a large peculiar motion of about 35 km/s toward the Galactic center.

We have observed the OH 1612-MHz maser emission towards the proto-planetary nebula candidate OH17.7-2.0 that underwent a very strong and unusual outburst in 2003. Phase-referencing data were obtained with the EVN in order to localize the outburst and to examine its possible causes. The majority of the emission comes from an incomplete spherical shell with inner and outer radii of 220 and 850 mas, respectively. There is a strong evidence for maser components that arise due to the interaction of a jet-like post-AGB outflow with the remnant outer AGB shell. The most prominent signature of such an interaction is the strongly bursting polarized emission near 73.3km/s coming from two unresolved components of brightness temperature up to 10^11K located at the edge of the biconal region 2500 AU from the central star. It is remarkable that this OH biconal region is well-aligned with the polar outflow inferred from the near-infrared image.

We characterize the molecular environment of candidate massive young stellar objects (MYSOs) signposted by methanol masers. Single pixel observations of 10 transitions of HCO^+, CO and CS isotopomers were carried out, using the IRAM 30m telescope. We studied a sample of 28 targets for which the 6.7GHz maser emission positions are known with a sub-arcsecond accuracy. The systemic velocity inferred from the optically thin lines agrees within 3km/s with the central velocity of the maser emission for most of the sources. About 64% of the sources show line wings in one or more transitions of CO, HCO^+ and CS species, indicating the presence of molecular outflows. Comparison of the widths of line wings and methanol maser emission suggests that the 6.7GHz maser line traces the environment of MYSO of various kinematic regimes. Therefore conditions conducive for the methanol maser can exist in the inner parts of molecular clouds or circumstellar discs as well as in the outer parts associated with molecular outflows. Calculations of the physical conditions based on the CO and HCO^+ lines and the CS line intensity ratios refine the input parameters for maser models.

We report on EVN imaging of the 6.7 GHz methanol maser emission from the candidate high-mass protostar G23.657-0.127. The masers originate in a nearly circular ring of 127 mas radius and 12 mas width. The ring structure points at a central exciting object which characteristics are typical for a young massive star; its bolometric luminosity is estimated to be <3.2*10^4 L_sun and <1.2*10^5 L_sun for near (5.1 kpc) and far (10.5 kpc) kinematic distances, respectively. However, the spatial geometry of the underlying maser region remains ambiguous. We consider scenarios in which the methanol masers originate in a spherical bubble or in a rotating disc seen nearly face-on.

We present the results of MERLIN polarization mapping of OH masers at 1665 and 1667 MHz towards the Cepheus A star-forming region. The maser emission is spread over a region of 6 arcsec by 10 arcsec, twice the extent previously detected. In contrast to the 22 GHz water masers, the OH masers associated with H II regions show neither clear velocity gradients nor regular structures. We identified ten Zeeman pairs which imply a magnetic field strength along the line-of-sight from -17.3 to +12.7 mG. The magnetic field is organised on the arcsecond scale, pointing towards us in the west and away from us in the east side. The linearly polarized components, detected for the first time, show regularities in the polarization position angles depending on their position. The electric vectors of OH masers observed towards the outer parts of H II regions are consistent with the interstellar magnetic field orientation, while those seen towards the centres of H II regions are parallel to the radio-jets. A Zeeman quartet inside a southern H II region has now been monitored for 25 years; we confirm that the magnetic field decays monotonically over that period.

We present the first EVN maps of 6.7 GHz methanol masers of three high-mass protostar candidates selected from the Torun unbiased survey of the Galactic plane. A variety of linear and arc like structures was detected. A number of maser clusters with projected sizes of 20-100 AU show monotonic velocity gradients. Some of them are roughly perpendicular to the major axes of these structures and can arise behind shock fronts.