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.

We wish to perform a statistical study of the location and motion of individual 22 GHz water maser cloudlets, characterized by sizes that are within a few au, with respect to the radio thermal emission from young stellar objects (YSO). For this purpose, we have been carrying out the Protostellar Outflows at the EarliesT Stages (POETS) survey of a sample (38) of high-mass YSOs. The water maser positions and three-dimensional (3D) velocities were determined through Very Long Baseline Array observations with accuracies of a few milliarcsec (mas) and a few km/s, respectively. The position of the ionized core of the protostellar wind, marking the YSO, was determined through sensitive continuum Jansky Very Large Array observations with a typical error of 20 mas. The statistic of the separation of the water masers from the radio continuum shows that 84% of the masers are found within 1000 au from the YSO and 45% of them are within 200 au. Therefore, we can conclude that the 22 GHz water masers are a reliable proxy for the YSO position. The distribution of maser luminosity is strongly peaked towards low values, indicating that about half of the maser population is still undetected with the current Very Long Baseline Interferometry detection thresholds of 50-100 mJy/beam. Next-generation, sensitive radio interferometers will exploit these weak masers for an improved sampling of the velocity and magnetic fields around the YSOs. The average direction of the water maser proper motions provides a statistically-significant estimate for the orientation of the jet emitted by the YSO: 55% of the maser proper motions are directed on the sky within an angle of 30 deg from the jet axis. Finally, we show that our measurements of 3D maser velocities statistically support models in which water maser emission arises from planar shocks with propagation direction close to the plane of the sky.

We present trigonometric parallax and proper motion measurements toward 22 GHz water and 6.7 GHz methanol masers in 16 high-mass star-forming regions. These sources are all located in the Scutum spiral arm of the Milky Way. The observations were conducted as part of the Bar and Spiral Structure Legacy (BeSSeL) survey. A combination of 14 sources from a forthcoming study and 14 sources from the literature, we now have a sample of 44 sources in the Scutum spiral arm, covering a Galactic longitude range from 0$^\circ$ to 33$^\circ$. A group of 16 sources shows large peculiar motions of which 13 are oriented toward the inner Galaxy. A likely explanation for these high peculiar motions is the combined gravitational potential of the spiral arm and the Galactic bar.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

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.

The goal of the Protostellar Outflows at the EarliesT Stages (POETS) survey is to image the disk-outflow interface on scales of 10-100 au in a statistically significant sample (36) of luminous young stellar objects (YSO), targeting both the molecular and ionized components of the outflows. The outflow kinematics is studied at milliarcsecond scales through VLBI observations of the 22 GHz water masers. We employed the JVLA at 6, 13, and 22 GHz in the A- and B-Array configurations to determine the spatial structure and the spectral index of the radio continuum emission. In about half of the targets, the water masers observed at separation <= 1000 au from the YSOs trace either or both of these kinematic structures: 1) a spatially elongated distribution oriented at close angle with the direction of collimation of the maser proper motions (PM), and 2) a linear LSR velocity (Vlsr) gradient across the YSO position. The kinematic structure (1) is readily interpreted in terms of a protostellar jet, as confirmed in some targets via the comparison with independent observations of the YSO jets, in thermal (continuum and line) emissions, reported in the literature. The kinematic structure (2) is interpreted in terms of a disk-wind (DW) seen almost edge-on on the basis of several pieces of evidence: first, it is invariably directed perpendicular to the YSO jet; second, it agrees in orientation and polarity with the Vlsr gradient in thermal emissions (when reported in the literature) identifying the YSO disk at scales of <= 1000~au; third, the PMs of the masers delineating the Vlsr gradients hint at flow motions at a speed of 10-20 km/s directed at large angles with the disk midplane. In the remaining targets, the maser PMs are not collimated but rather tend to align along two almost perpendicular directions, and could originate in DW-jet systems slightly inclined (<= 30 deg) with respect to edge-on.

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

Centimeter continuum observations of protostellar jets have revealed the presence of knots of shocked gas where the flux density decreases with frequency. This spectrum is characteristic of nonthermal synchrotron radiation and implies the presence of both magnetic fields and relativistic electrons in protostellar jets. Here, we report on one of the few detections of nonthermal jet driven by a young massive star in the star-forming region G035.02$+$0.35. We made use of the NSF's Karl G. Jansky Very Large Array (VLA) to observe this region at C, Ku, and K bands with the A- and B-array configurations, and obtained sensitive radio continuum maps down to a rms of 10 $\mu$Jy beam$^{-1}$. These observations allow for a detailed spectral index analysis of the radio continuum emission in the region, which we interpret as a protostellar jet with a number of knots aligned with extended 4.5 $\mu$m emission. Two knots clearly emit nonthermal radiation and are found at similar distances, of approximately 10,000 au, each side of the central young star, from which they expand at velocities of hundreds km s$^{-1}$. We estimate both the mechanical force and the magnetic field associated with the radio jet, and infer a lower limit of $0.4\times10^{-4} $M$_{\odot}$ yr$^{-1}$ km s$^{-1}$ and values in the range $0.7-1.3 $mG, respectively.

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.

We previously observed the star-forming region G16.59$-$0.05 through interferometric observations of both thermal and maser lines, and identified a high-mass young stellar object (YSO) which is surrounded by an accretion disk and drives a nonthermal radio jet. We performed high-angular-resolution (beam FWHM ~0.15") 1.2-mm continuum and line observations towards G16.59$-$0.05 with the Atacama Large Millimeter Array (ALMA). The main dust clump, with size ~10$^4$ au, is resolved into four relatively compact (diameter ~2000 au) millimeter (mm) sources. The source harboring the high-mass YSO is the most prominent in molecular emission. By fitting the emission profiles of several unblended and optically thin transitions of CH$_3$OCH$_3$ and CH$_3$OH, we derived gas temperatures inside the mm-sources in the range 42--131 K, and calculated masses of 1--5 $M_{\odot}$. A well-defined Local Standard of Rest velocity (Vlsr) gradient is detected in most of the high-density molecular tracers at the position of the high-mass YSO, pinpointed by compact 22-GHz free-free emission. This gradient is oriented along a direction forming a large (~70 degree) angle with the radio jet, traced by elongated 13-GHz continuum emission. The butterfly-like shapes of the P-V plots and the linear pattern of the emission peaks of the molecular lines at high velocity confirm that this Vlsr gradient is due to rotation of the gas in the disk surrounding the high-mass YSO. The disk radius is ~500 au, and the Vlsr distribution along the major axis of the disk is well reproduced by a Keplerian profile around a central mass of 10$\pm$2 $M_{\odot}$. The position of the YSO is offset by >~ 0.1" from the axis of the radio jet and the dust emission peak. To explain this displacement we argue that the high-mass YSO could have moved from the center of the parental mm source owing to dynamical interaction with one or more companions.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

Abridged. Methods. We made use of the Karl G. Jansky Very Large Array (VLA) in the B configuration at K band, and in the A configuration at both Ku and C bands, in order to image the radio continuum emission towards 25 H$_2$O maser sites with an angular resolution and thermal rms of the order of $0.1''$ and 10 $\mu$Jy beam$^{-1}$, respectively. These targets add to our pilot study of 11 maser sites presented in Moscadelli et al. (2016). The sample of H$_2$O maser sites was selected among those regions having an accurate distance measurement, obtained through maser trigonometric parallaxes, and H$_2$O maser luminosities in excess of 10$^{-6}$ L$_{\odot}$. Results. We present high-resolution radio continuum images of 33 sources belonging to 25 star-forming regions. In each region, we detect radio continuum emission within a few 1000 au of the H$_2$O masers' position; 50% of the radio continuum sources are associated with bolometric luminosities exceeding 5 $\times$ 10$^{3}$ L$_{\odot}$, including W33A and G240.32$+$0.07. We provide a detailed spectral index analysis for each radio continuum source, based on the integrated fluxes at each frequency, and produce spectral index maps with the multi-frequency-synthesis deconvolution algorithm of CASA. The radio continuum emission traces thermal bremsstrahlung in (proto)stellar winds and jets, with flux densities at 22 GHz below 3 mJy, and spectral index values between $-0.1$ and 1.3. We prove a strong correlation ($r>$0.8) between the radio continuum luminosity (L$_{\rm rad}$) and the H$_2$O maser luminosity (L$_{\rm H_2O}$) of ($\rm L_{8GHz}/mJy\times kpc^2)=10^{3.8}\times(L_{H_2O}/L_{\odot})^{0.74}$. Since H$_2$O masers are excited through shocks driven by (proto)stellar winds and jets, these results provide support to the idea that the radio continuum emission around young stars is dominated by shock-ionization.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Context. Recent observations of the massive young stellar object S255 NIRS 3 have revealed a large increase in both methanol maser flux density and IR emission, which have been interpreted as the result of an accretion outburst, possibly due to instabilities in a circumstellar disk. This indicates that this type of accretion event could be common in young/forming early-type stars and in their lower mass siblings, and supports the idea that accretion onto the star may occur in a non-continuous way. Aims. As accretion and ejection are believed to be tightly associated phenomena, we wanted to confirm the accretion interpretation of the outburst in S255 NIRS 3 by detecting the corresponding burst of the associated thermal jet. Methods. We monitored the radio continuum emission from S255 NIRS 3 at four bands using the Karl G. Jansky Very Large Array. The millimetre continuum emission was also observed with both the Northern Extended Millimeter Array of IRAM and the Atacama Large Millimeter/submillimeter Array. Results. We have detected an exponential increase in the radio flux density from 6 to 45 GHz starting right after July 10, 2016, namely about 13 months after the estimated onset of the IR outburst. This is the first ever detection of a radio burst associated with an IR accretion outburst from a young stellar object. The flux density at all observed centimetre bands can be reproduced with a simple expanding jet model. At millimetre wavelengths we infer a marginal flux increase with respect to the literature values and we show this is due to free-free emission from the radio jet. Abridged.

Methanol and water masers indicate young stellar objects. They often exhibit flares, and a fraction shows periodic activity. Several mechanisms might explain this behavior but the lack of concurrent infrared (IR) data complicates to identify the cause. Recently, 6.7 GHz methanol maser flares were observed, triggered by accretion bursts of high-mass YSOs which confirmed the IR-pumping of these masers. This suggests that regular IR changes might lead to maser periodicity. Hence, we scrutinized space-based IR imaging of YSOs associated with periodic methanol masers. We succeeded to extract the IR light curve from NEOWISE data for the intermediate mass YSO G107.298+5.639. Thus, for the first time a relationship between the maser and IR variability could be established. While the IR light curve shows the same period of ~34.6 days as the masers, its shape is distinct from that of the maser flares. Possible reasons for the IR periodicity are discussed.

This review arose from the European Radio Astronomy Technical Forum (ERATec) meeting held in Firenze, October 2015, and aims to highlight the breadth and depth of the high-impact science that will be aided and assisted by the use of simultaneous mm-wavelength receivers. Recent results and opportunities are presented and discussed from the fields of: continuum VLBI (observations of weak sources, astrometry, observations of AGN cores in spectral index and Faraday rotation), spectral line VLBI (observations of evolved stars and massive star-forming regions) and time domain observations of the flux variations arising in the compact jets of X-ray binaries. Our survey brings together a large range of important science applications, which will greatly benefit from simultaneous observing at mm-wavelengths. Such facilities are essential to allow these applications to become more efficient, more sensitive and more scientifically robust. In some cases without simultaneous receivers the science goals are simply unachievable. Similar benefits would exist in many other high frequency astronomical fields of research.

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.

Aims: In order to test the nature of an (accretion) disk in the vicinity of Cepheus A HW2, we measured the three-dimensional velocity field of the CH3OH maser spots, which are projected within 1000au of the HW2 object, with an accuracy of the order of 0.1km/s. Methods: We made use of the European VLBI Network (EVN) to image the 6.7GHz CH3OH maser emission towards Cepheus A HW2 with 4.5 milli-arcsecond resolution (3au). We observed at three epochs spaced by one year between 2013 and 2015. During the last epoch, on mid-march 2015, we benefited from the new deployed Sardinia Radio Telescope. Results: We show that the CH3OH velocity vectors lie on a preferential plane for the gas motion with only small deviations of 12+/-9 degrees away from the plane. This plane is oriented at a position angle of 134 degrees east of north, and inclined by 26 degrees with the line-of-sight, closely matching the orientation of the disk-like structure previously reported by Patel et al.(2005). Knowing the orientation of the equatorial plane, we can reconstruct a face-on view of the CH3OH gas kinematics onto the plane. CH3OH maser emission is detected within a radius of 900au from HW2, and down to a radius of about 300au, the latter coincident with the extent of the dust emission at 0.9mm. The velocity field is dominated by an infall component of about 2km/s down to a radius of 300au, where a rotational component of 4km/s becomes dominant. We discuss the nature of this velocity field and the implications for the enclosed mass. Conclusions: These findings bring direct support to the interpretation that the high-density gas and dust emission, surrounding Cepheus A HW2, trace an accretion disk.

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.

Aims: We want to study the physical properties of the ionized jet emission in the vicinity of an O-type young stellar object (YSO), and estimate how efficient is the transfer of energy and momentum from small- to large-scale outflows. Methods: We conducted Karl G. Jansky Very Large Array (VLA) observations, at both 22 and 45 GHz, of the compact and faint radio continuum emission in the high-mass star-forming region G023.01-00.41, with an angular resolution between 0.3" and 0.1", and a thermal rms of the order of 10 uJy/beam. Results: We discovered a collimated thermal (bremsstrahlung) jet emission, with a radio luminosity (L_rad) of 24 mJy kpc^2 at 45 GHz, in the inner 1000 AU from an O-type YSO. The radio thermal jet has an opening angle of 44 degrees and brings a momentum rate of 8 10^-3 M_sun yr^-1 km/s. By combining the new data with previous observations of the molecular outflow and water maser shocks, we can trace the outflow emission from its driving source through the molecular clump, across more than two order of magnitude in length (500 AU-0.2 pc). We find that the momentum-transfer efficiency, between the inner jet emission and the extended outflow of entrained ambient gas, is near unity. This result suggests that the large-scale flow is swept-up by the mechanical force of the radio jet emission, which originates in the inner 1000 AU from the high-mass YSO.

The nature of the spiral structure of the Milky Way has long been debated. Only in the last decade have astronomers been able to accurately measure distances to a substantial number of high-mass star-forming regions, the classic tracers of spiral structure in galaxies. We report distance measurements at radio wavelengths using the Very Long Baseline Array for eight regions of massive star formation near the Local spiral arm of the Milky Way. Combined with previous measurements, these observations reveal that the Local Arm is larger than previously thought, and both its pitch angle and star formation rate are comparable to those of the Galaxy's major spiral arms, such as Sagittarius and Perseus. Toward the constellation Cygnus, sources in the Local Arm extend for a great distance along our line of sight and roughly along the solar orbit. Because of this orientation, these sources cluster both on the sky and in velocity to form the complex and long enigmatic Cygnus X region. We also identify a spur that branches between the Local and Sagittarius spiral arms.

We report full polarimetric VLBA observations of water masers towards the Turner-Welch Object in the W3(OH) high-mass star forming complex. This object drives a synchrotron jet, which is quite exceptional for a high-mass protostar, and is associated with a strongly polarized water maser source, W3(H$_2$O), making it an optimal target to investigate the role of magnetic fields on the innermost scales of protostellar disk-jet systems. The linearly polarized emission from water masers provides clues on the orientation of the local magnetic field, while the measurement of the Zeeman splitting from circular polarization provides its strength. The water masers trace a bipolar, biconical outflow at the center of the synchrotron jet. Although on scales of a few thousand AU the magnetic field inferred from the masers is on average orientated along the flow axis, on smaller scales (10s to 100s of AU), we have revealed a misalignment between the magnetic field and the velocity vectors, which arises from the compression of the field component along the shock front. Our measurements support a scenario where the magnetic field would evolve from having a dominant component parallel to the outflow velocity in the pre-shock gas, with field strengths of the order of a few tens of mG (at densities of $10^7$ cm$^{-3}$), to being mainly dominated by the perpendicular component of order of a few hundred of mG in the post-shock gas where the water masers are excited (at densities of $10^9$ cm$^{-3}$). The general implication is that in the undisturbed (i.e. not-shocked) circumstellar gas, the flow velocities would follow closely the magnetic field lines, while in the gas shocked by the prostostellar jet the magnetic field would be re-configured to be parallel to the shock front.

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

We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object (YSO), at very high spatial resolution (10-100 AU). We performed milli-arcsecond polarimetric observations of the strong CH3OH maser emission observed in the vicinity of an O-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH3OH masing gas previously measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2 degree/AU). Overall, the velocity field vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30 degrees. We interpret this finding as the contribution of a turbulent velocity field of about 3.5 km/s, responsible for braking up the alignment between the velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk plane, with an average speed of 7 km/s. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 500-600 AU, and inward for shorter radii.

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.

The high-mass star-forming site G009.62+00.20E hosts the 6.7 GHz methanol maser source with the greatest flux density in the Galaxy which has been flaring periodically over the last ten years. We performed high-resolution astrometric measurements of the CH3OH, H2O, and OH maser emission and 7 mm continuum in the region. The radio continuum emission was resolved in two sources separated by 1300 AU. The CH3OH maser cloudlets are distributed along two north-south ridges of emission to the east and west of the strongest radio continuum component. This component likely pinpoints a massive young stellar object which heats up its dusty envelope, providing a constant IR pumping for the Class II CH3OH maser transitions. We suggest that the periodic maser activity may be accounted for by an independent, pulsating, IR radiation field provided by a bloated protostar in the vicinity of the brightest masers. We also report about the discovery of an elliptical distribution of CH3OH maser emission in the region of periodic variability.

To constrain theoretical models of high-mass star formation, observational signatures of mass accretion in O-type forming stars are desirable. Using the JVLA, we have mapped the hot and dense molecular gas in the hot core NGC7538 IRS1, with 0.2'' angular resolution, in seven metastable (J=K) inversion transitions of ammonia: (J,K)=(6,6), (7,7), (9,9), (10,10), (12,12), (13,13), and (14,14). These lines arise from energy levels between ~400 K and ~1950 K above the ground state, and are observed in absorption against the HC-HII region associated with NGC7538 IRS1. With a 500 AU linear resolution, we resolve the elongated North-South ammonia structure into two compact components: the main core and a southernmost component. Previous observations of the radio continuum with a 0.08'' (or 200 AU) resolution, resolved in turn the compact core in two (northern and southern) components. These features correspond to a triple system of high-mass YSOs IRS1a, IRS1b, and IRS1c identified with VLBI measurements of methanol masers. The velocity maps of the compact core show a clear velocity gradient in all lines, which is indicative of rotation in a (circumbinary) envelope, containing ~40 solar masses (dynamical mass). In addition, we derived physical conditions of the molecular gas: rotational temperatures ~280 K, ammonia column densities ~1.4-2.5 x 10^19 cm-2, H_2 volume densities ~3.5-6.2 x 10^10 cm-3, and a total gas mass in the range of 19-34 solar masses, for the main core. We conclude that NGC7538 IRS1 is the densest hot molecular core known, containing a rotating envelope which hosts a multiple system of high-mass YSOs, possibly surrounded by accretion disks. Future JVLA observations in the A-configuration are needed to resolve the binary system in the core and may allow to study the gas kinematics in the accretion disks associated with individual binary members.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

We report measurements of parallaxes and proper motions of ten high-mass star-forming regions in the Sagittarius spiral arm of the Milky Way as part of the BeSSeL Survey with the VLBA. Combining these results with eight others from the literature, we investigated the structure and kinematics of the arm between Galactocentric azimuth around -2 and 65 deg. We found that the spiral pitch angle is 7.3 +- 1.5 deg; the arm's half-width, defined as the rms deviation from the fitted spiral, is around 0.2 kpc; and the nearest portion of the Sagittarius arm is 1.4 +- 0.2 kpc from the Sun. Unlike for adjacent spiral arms, we found no evidence for significant peculiar motions of sources in the Sagittarius arm opposite to Galactic rotation.

NGC7538 IRS1 is considered the best high-mass accretion disk candidate around an O-type young star in the northern hemisphere. We investigated the 3D kinematics and dynamics of circumstellar gas with very high linear resolution, from tens to 1500 AU, with the ultimate goal of building a comprehensive dynamical model for this YSO. We employed four different observing epochs of EVN data at 6.7 GHz, spanning almost eight years, which enabled us to measure, besides line-of-sight (l.o.s.) velocities and positions, also l.o.s. accelerations and proper motions of methanol masers. In addition, we imaged with the JVLA-B array highly-excited ammonia inversion lines, from (6,6) to (13,13), which enabled us to probe the hottest molecular gas very close to the exciting source(s). We found five 6.7 GHz maser clusters which are distributed over a region extended N-S across ~1500 AU and are associated with three peaks of the radio continuum. We proposed that these maser clusters identify three individual high-mass YSOs, named IRS1a, IRS1b, and IRS1c. We modeled the maser clusters in IRS1a and IRS1b in terms of edge-on disks in centrifugal equilibrium. In the first case, masers may trace a quasi-Keplerian thin disk, orbiting around a high-mass YSO, IRS1a, of up to 25 solar masses. This YSO dominates the bolometric luminosity of the region. The second disk is both massive (<16 Msun within ~500 AU) and thick, and the mass of the central YSO, IRS1b, is constrained to be at most a few solar masses. In summary, we present compelling evidence that NGC7538 IRS1 is not forming just one single high-mass YSO, but consists of a multiple system of high-mass YSOs, which are surrounded by accretion disks, and are probably driving individual outflows. This new model naturally explains all the different orientations and disk/outflow structures proposed for the region in previous models.

(Abridged) METHODS: We performed SMA observations at 1.3 mm with both the most extended and compact array configurations, providing sub-arcsecond and high sensitivity maps of various molecular lines, including both hot-core and outflow tracers. We also reconstruct the spectral energy distribution of the region from millimeter to near infrared wavelengths, using the Herschel/Hi-GAL maps, as well as archival data. RESULTS: From the spectral energy distribution, we derive a bolometric luminosity of about 4x10^4 Lsun. Our interferometric observations reveal that the distribution of dense gas and dust in the HMC is significantly flattened and extends up to a radius of 8000 AU from the center of radio continuum and maser emission in the region. The equatorial plane of this HMC is strictly perpendicular to the elongation of the collimated bipolar outflow, as imaged on scales of about 0.1-0.5 pc in the main CO isotopomers as well as in the SiO(5-4) line. In the innermost HMC regions (ca. 1000 AU), the velocity field traced by the CH3CN(12_K-11_K) line emission shows that molecular gas is both expanding along the outflow direction following a Hubble-law, and rotating about the outflow axis, in agreement with the (3-D) velocity field traced by methanol masers. The velocity field associated with rotation indicates a dynamical mass of 19 Msun at the center of the core. The latter is likely to be concentrated in a single O9.5 ZAMS star, consistent with the estimated bolometric luminosity of G023.01-00.41. The physical properties of the CO(2-1) outflow emission, such as its momentum rate 6x10^-3 Msun km/s /yr and its outflow rate 2x10^-4 Msun/yr, support our estimates of the luminosity (and mass) of the embedded young stellar object.

IRAS20126+4104 is a well studied B0.5 protostar that is surrounded by a ~1000 au Keplerian disk and is where a large-scale outflow originates. Both 6.7-GHz CH3OH masers and 22-GHz H2O masers have been detected toward this young stellar object. The CH3OH masers trace the Keplerian disk, while the H2O masers are associated with the surface of the conical jet. Recently, observations of dust polarized emission (350 um) at an angular resolution of 9 arcseconds (~15000 au) have revealed an S-shaped morphology of the magnetic field around IRAS20126+4104. The observations of polarized maser emissions at milliarcsecond resolution (~20 au) can make a crucial contribution to understanding the orientation of the magnetic field close to IRAS20126+4104. This will allow us to determine whether the magnetic field morphology changes from arcsecond resolution to milliarcsecond resolution. The European VLBI Network was used to measure the linear polarization and the Zeeman splitting of the 6.7-GHz CH3OH masers toward IRAS20126+4104. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman splitting of the 22-GHz H2O masers toward the same region. We detected 26 CH3OH masers and 5 H2O masers at high angular resolution. Linear polarization emission was observed toward three CH3OH masers and toward one H2O maser. Significant Zeeman splitting was measured in one CH3OH maser (\Delta V_{Z}=-9.2 +/- 1.4 m/s). No significant (5 sigma) magnetic field strength was measured using the H2O masers. We found that in IRAS20126+4104 the rotational energy is less than the magnetic energy.

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.

We report trigonometric parallax measurements of masers in the massive star forming complex W43 from VLBA observations as part of the BeSSeL Survey. Based on measurements of three 12 GHz methanol maser sources (G029.86-00.04, G029.95-00.01 and G031.28+00.06) and one 22 GHz water maser source (G031.58+00.07) toward W43, we derived a distance of $5.49^{+0.39}_{-0.34}$ kpc to W43. By associating the masers with CO molecular clouds, and associating the clouds kinematically with CO longitude-velocity spiral features, we assign W43 to the Scutum spiral arm, close to the near end of the Galactic bar. The peculiar motion of W43 is about 20 km/s toward the Galactic Center and is very likely induced by the gravitational attraction of the bar.

We report four trigonometric parallaxes for high-mass star-forming regions within 4 kpc of the Galactic center. These measurements were made with the VLBA as part of the BeSSeL Survey. By associating these sources kinematically with large-scale features in CO and HI longitude-velocity diagrams, we begin to outline some major features of the inner Milky Way: the Connecting arm, the near and far 3 kpc arms, and the Norma arm. The Connecting arm in the first Galactic quadrant lies closer to the Galactic center than the far 3 kpc arm and is offset by the long-bar's major axis near its leading edge, supporting the presence of an inner Lindblad resonance. Assuming the 3 kpc arms are a continuous physical structure, the relative Galactocentric distance of its near and far sides suggests highly elliptical streamlines of gas around the bar(s) and a bar corotation radius, r_CR > 3.6 kpc. At a Galactic longitude near 10{\deg} and a heliocentric distance of about 5 kpc, the near 3 kpc arm and the Norma arm intersect on a face-on view of our Galaxy, while passing at different Galactic latitudes. We provide an accurate distance measurement to the W31 star-forming complex of 4.95 (+0.5;-0.43) kpc from the Sun, which associates it with a bright CO feature belonging to the near 3 kpc arm.

We wish to investigate the structure, velocity field, and stellar content of the G16.59-0.05 high-mass star-forming region, where previous studies have established the presence of two almost perpendicular (NE-SW and SE-NW), massive outflows, and a rotating disk traced by methanol maser emission. We performed Very Large Array observations of the radio continuum and ammonia line emission, complemented by COMICS/Subaru and Hi-GAL/Herschel images in the mid- and far-infrared (IR). Our centimeter continuum maps reveal a collimated radio jet that is oriented E-W and centered on the methanol maser disk, placed at the SE border of a compact molecular core. The spectral index of the jet is negative, indicating non-thermal emission over most of the jet, except the peak close to the maser disk, where thermal free-free emission is observed. We find that the ammonia emission presents a bipolar structure consistent (on a smaller scale) in direction and velocity with that of the NE-SW bipolar outflow detected in previous CO observations. After analyzing our previous N2H+(1-0) observations again, we conclude that two scenarios are possible. In one case both the radio jet and the ammonia emission would trace the root of the large-scale CO bipolar outflow. The different orientation of the jet and the ammonia flow could be explained by precession and/or a non-isotropic density distribution around the star. In the other case, the N2H+(1-0) and ammonia bipolarity is interpreted as two overlapping clumps moving with different velocities along the line of sight. The ammonia gas also seems to undergo rotation consistent with the maser disk. Our IR images complemented by archival data allow us to derive a bolometric luminosity of about 10^4 L_sun and to conclude that most of the luminosity is due to the young stellar object associated with the maser disk.

Maser lines from different molecular species, including water, hydroxyl, and methanol, are common observational phenomena associated with massive star forming regions. In particular, the methanol maser appears as an ideal tool to study the early phases of massive star formation. However, it is difficult to establish the exact start of the methanol maser phase, and it would then be interesting to detect and study low-flux density methanol masers (i.e., < 0.1 Jy or even << 0.1 Jy), in order to determine if they can effectively be used to mark a specific evolutionary phase in high-mass star formation. Past surveys have been unable to systematically detect many low-flux density methanol masers, and thus we do not yet know how many such masers exist in the Galaxy and what is their physical nature. Out of a sample of 107 observed Hi-GAL sources we detected a total of 32 methanol masers, with 22 sources being new and weak (median peak flux density 0.07 Jy) detections, in the Galactic longitude range [32.0, 59.8]deg. We also detected 12 6.035-GHz OH maser, with 9 objects being new detections. Our survey covers a similar range of source distances as the "Arecibo Methanol Maser Galactic Plane Survey", but the methanol masers detected by us are clearly shifted towards lower integrated flux densities. The newly detected methanol masers are mostly of low-luminosity and, except for some sources, their weakness is not due to distance effects or positional offsets. No specific correlation is found with the physical parameters of the Hi-GAL clumps, except for sources with both CH3OH and OH masers which tend to have higher mass and luminosity. The intensity of the methanol masers correlates well with the velocity range of the maser emission, which suggests that the low brightness of these masers is related to the number of maser spots in the emitting region and their evolution with time.

Trigonometric parallax measurements of nine water masers associated with the Local arm of the Milky Way were carried out as part of the BeSSeL Survey using the VLBA. When combined with 21 other parallax measurements from the literature, the data allow us to study the distribution and 3-dimensional motions of star forming regions in the spiral arm over the entire northern sky. Our results suggest that the Local arm does not have the large pitch angle characteristic of a short spur. Instead its active star formation, overall length (>5 kpc), and shallow pitch angle (~10 degrees) suggest that it is more like the adjacent Perseus and Sagittarius arms; perhaps it is a branch of one of these arms. Contrary to previous results, we find the Local arm to be closer to the Perseus than to the Sagittarius arm, suggesting that a branching from the former may be more likely. An average peculiar motion of near-zero toward both the Galactic center and north Galactic pole, and counter rotation of ~ 5 km/s were observed, indicating that the Local arm has similar kinematic properties as found for other major spiral arms.

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

We imaged the molecular outflows towards the cluster of high-mass young stellar objects G24.78+0.08 at high-angular resolution using SiO emission, which is considered the classical tracer of protostellar jets. We performed SiO observations with the VLA interferometer in the J = 1-0 v=0 transition and with the SMA array in the 5-4 transition. A complementary IRAM 30-m single-dish survey in the (2-1), (3-2), (5-4), and (6-5) SiO lines was also carried out. Two collimated SiO high-velocity outflows driven by the A2 and C millimeter continuum massive cores have been imaged. On the other hand, we detected no SiO outflow driven by the young stellar objects in more evolved evolutionary phases that are associated with ultracompact (B) or hypercompact (A1) HII regions. The LVG analysis reveals high-density gas (10^3-10^4 cm-3), with well constrained SiO column densities (0.5-1 10^15 cm-2). The driving source of the A2 outflow is associated with typical hot core tracers such as methyl formate, vinyl cyanide, cyanoacetilene, and acetone. The driving source of the main SiO outflow in G24 has an estimated luminosity of a few 10^4 Lsun (typical of a late O-type star) and is embedded in the 1.3 mm continuum core A2, which in turn is located at the centre of a hot core that rotates on a plane perpendicular to the outflow main axis. The present SiO images support a scenario similar to the low-mass case for massive star formation, where jets that are clearly traced by SiO emission, create outflows of swept-up ambient gas usually traced by CO.

This work presents a detailed study of the gas kinematics towards the "Hot Molecular Core" (HMC) G31.41+0.31 via multi-epoch VLBI observations of the H2O 22 GHz and CH3OH 6.7 GHz masers, and single-epoch VLBI of the OH 1.6 GHz masers. Water masers present a symmetric spatial distribution with respect to the HMC center, where two nearby (0.2" apart), compact, VLA sources (labeled "A" and "B") are previously detected. The spatial distribution of a first group of water masers, named "J1", is well fit with an elliptical profile, and the maser proper motions mainly diverge from the ellipse center, with average speed of 36 km s-1. These findings strongly suggest that the "J1" water maser group traces the heads of a young (dynamical time of 1.3 10^3 yr), powerful (momentum rate of ~0.2 M_sun yr-1 km s-1), collimated (semi-opening angle ~10 deg) jet emerging from a MYSO located close (within 0.15") to the VLA source "B". Most of the water features not belonging to "J1" present an elongated (about 2" in size), NE--SW oriented (PA = 70 deg), S-shape distribution, which we denote with the label "J2". The elongated distribution of the "J2" group and the direction of motion, approximately parallel to the direction of elongation, of most "J2" water masers suggests the presence of another collimated outflow, emitted from a MYSO near the VLA source "A". The orientation of the "J2" jet agrees well with that (PA = 68 deg) of the well-defined V_LSR gradient across the HMC revealed by previous interferometric, thermal line observations. Furthermore, the "J2" jet is powerful enough to sustain the large momentum rate, 0.3 M_sun yr-1 km s-1, estimated assuming that the V_LSR gradient represents a collimated outflow. These two facts lead us to favour the interpretation of the V_LSR gradient across the G31.41+0.31 HMC in terms of a compact and collimated outflow.

We used the Very Long Baseline Array (VLBA) and the European VLBI Network (EVN) to perform phase-referenced VLBI observations of the three most powerful maser transitions associated with the high-mass star-forming region G28.87+0.07: the 22.2 GHz H$_{2}$O, 6.7 GHz CH$_{3}$OH, and 1.665 GHz OH lines. We also performed VLA observations of the radio continuum emission at 1.3 and 3.6 cm and Subaru observations of the continuum emission at 24.5 $\mu$m. Two centimeter continuum sources are detected and one of them (named "HMC") is compact and placed at the center of the observed distribution of H$_{2}$O, CH$_{3}$OH and OH masers. The bipolar distribution of line-of-sight (l.o.s) velocities and the pattern of the proper motions suggest that the water masers are driven by a (proto)stellar jet interacting with the dense circumstellar gas. The same jet could both excite the centimeter continuum source named "HMC" (interpreted as free-free emission from shocked gas) and power the molecular outflow observed at larger scales -- although one cannot exclude that the free-free continuum is rather originating from a hypercompact \ion{H}{2} region. At 24.5 $\mu$m, we identify two objects separated along the north-south direction, whose absolute positions agree with those of the two VLA continuum sources. We establish that $\sim$90% of the luminosity of the region ($\sim$\times10^{5} L_\sun$) is coming from the radio source "HMC", which confirms the existence of an embedded massive young stellar object (MYSO) exciting the masers and possibly still undergoing heavy accretion from the surrounding envelope.

We report on a detailed study of the water maser kinematics and radio continuum emission toward the most massive and young object in the star-forming region AFGL 2591. Our analysis shows at least two spatial scales of multiple star formation, one projected across 0.1 pc on the sky and another one at about 2000 AU from a ZAMS star of about 38 Solar masses. This young stellar object drives a powerful jet- and wind-driven outflow system with the water masers associated to the outflow walls, previously detected as a limb-brightened cavity in the NIR band. At about 1300 AU to the north of this object a younger protostar drives two bow shocks, outlined by arc-like water maser emission, at 200 AU either side of the source. We have traced the velocity profile of the gas that expands along these arc-like maser structures and compared it with the jet-driven outflow model. This analysis suggests that the ambient medium around the northern protostar is swept up by a jet-driven shock (>66 km/s) and perhaps a lower-velocity (~10 km/s) wind with an opening angle of about 20 degrees from the jet axis.

We report a trigonometric parallax measurement with the Very Long Baseline Array for the water maser in the distant high-mass star-forming region G75.30+1.32. This source has a heliocentric distance of 9.25+-0.45 kpc, which places it in the Outer arm in the first Galactic quadrant. It lies 200 pc above the Galactic plane and is associated with a substantial HI enhancement at the border of a large molecular cloud. At a Galactocentric radius of 10.7 kpc, G75.30+1.32 is in a region of the Galaxy where the disk is significantly warped toward the North Galactic Pole. While the star-forming region has an instantaneous Galactic orbit that is nearly circular, it displays a significant motion of 18 km/s toward the Galactic plane. The present results, when combined with two previous maser studies in the Outer arm, yield a pitch angle of about 12 degrees for a large section of the arm extending from the first quadrant to the third.

Observational signatures of infalling envelopes and outflowing material in early stages of protostellar evolution, and at small radii from the protostar, are essential to progress in the understanding of the mass-accretion process in star formation. In this letter, we report a detailed study of the accretion and outflow structure around a protostar in the well-known high-mass star-forming region AFGL 5142. We focus on the mm source MM-1, which exhibits hot-core chemistry, radio continuum emission, and strong water and methanol masers. Remarkably, our Very Long Baseline Interferometry (VLBI) observations of molecular masers over six years provided us with the 3-D velocity field of circumstellar molecular gas with a resolution of 0.001-0.005 arcseconds and at radii <0.23 arcseconds (or 400 AU) from the protostar. In particular, our measurements of methanol maser emission enabled, for the first time, a direct measurement of infall of a molecular envelope (radius of 300 AU and velocity of 5 km/s) onto an intermediate- to high-mass protostar. We estimate an infall rate of 0.0006 n_8 Msun/year, where n_8 is the ambient volume density in units of 10^8 cm-3 (required for maser excitation). In addition, our measurements of water maser (and radio continuum) emission identify a collimated bipolar molecular outflow (and ionized jet) from MM-1. The evidence of simultaneous accretion and outflow at small spatial scales, makes AFGL 5142 an extremely compelling target for high-angular resolution studies of high-mass star formation.

This work presents a study of the milliarcsecond structure of the 6.7 GHz methanol masers at high-velocity resolution (0.09 km s^(-1)) in four high-mass star-forming regions: G16.59-0.05, G23.01-0.41, IRAS20126+4104, and AFGL5142. We studied these sources by means of multi-epoch VLBI observations in the 22 GHz water and 6.7 GHz methanol masers, to determine the 3-D gas kinematics within a few thousand AU from the (proto)star. The present work exploits the 6.7 GHz maser data collected so far to investigate the milliarcsecond structure of this maser emission at high-velocity resolution. Most of the detected 6.7 GHz maser features present an ordered (linear, or arc-like) distribution of maser spots on the plane of the sky, together with a regular variation in the spot LSR velocity (V_LSR) with position. Typical values for the amplitude of the V_LSR gradients (defined in terms of the derivative of the spot V_LSR with position) are found to be 0.1-0.2 km s^(-1) mas^(-1). In each of the four target sources, the orientation and the amplitude of most of the feature V_LSR gradients remain remarkably stable in time, on timescales of (at least) several years. We also find that the data are consistent with having the V_LSR gradients and proper motion vectors in the same direction on the sky, considered the measurement uncertainties. The time persistency, the ordered angular and spatial distribution, and the orientation generally similar to the proper motions, altogether suggest a kinematical interpretation for the origin of the 6.7 GHz maser V_LSR gradients. This work shows that the organized motions (outflow, infall, and rotation) revealed by the (22 GHz water and 6.7 GHz methanol) masers on large scales (~100-1000 AU) also persist to very small (~10 AU) scales.

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.

We report trigonometric parallaxes for three massive star forming regions, corresponding to distances of $2.34^{+0.13}_{-0.11}$ kpc for G12.89+0.49 (also known as IRAS 18089-1732),$1.98^{+0.14}_{-0.12}$ kpc for G15.03-0.68 (in the M17 region), and $8.0^{+4.0}_{-2.0}$ kpc for G27.36-0.16. Both G12.89+0.49 and G15.03-0.68 are located in the Carina-Sagittarius spiral arm.

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.

IRAS 20126+4104 is one of the best candidates for a high-mass (proto)star surrounded by an accretion disk. Such a fact may be used to set constraints on theories of high-mass star formation, but requires confirmation that the mass and luminosity of IRAS 20126+4104 are indeed typical of a B0.5 star, which in turn requires an accurate estimate of the distance. We used the Very Long Baseline Array and the European VLBI Network to observe the 22.2 GHz water and 6.7 GHz methanol masers in IRAS 20126+4104 at a number of epochs suitably distributed in time. The absolute positions of the maser features were established with respect to reference quasars, which allowed us to derive absolute proper motions. From the parallax of the water masers we obtain a distance of 1.64 \pm 0.05 kpc, which is very similar to the value adopted so far in the literature (1.7 kpc) and confirms that IRAS 20126+4104 is a high-mass (proto)star. From the methanol masers we derive the component in the plane of the sky of the systemic velocity of the disk+star system (-16 km/s in right-ascension and +7.6 km/s in declination). Accurate knowledge of the distance and systemic velocity allows us to improve on the model fit to the water maser jet presented in a previous study. Finally, we identify two groups of methanol maser features, one undergoing rotation in the disk and possibly distributed along a narrow ring centered on the star, the other characterised by relative proper motions indicating that the features are moving away from the disk, perpendicular to it. We speculate that the latter group might be tracing the disk material marginally entrained by the jet.

The present paper focuses on the high-mass star-forming region G23.01-0.41. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G23.01-0.41: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the HMC at 1.3 and 3.6 cm. Results: We have detected H2O, CH3OH, and OH maser emission clustered within 2000 AU from the center of a flattened HMC, oriented SE-NW, from which emerges a massive 12CO outflow, elongated NE-SW, extended up to the pc-scale. Although the three maser species show a clearly different spatial and velocity distribution and sample distinct environments around the massive YSO, the spatial symmetry and velocity field of each maser specie can be explained in terms of expansion from a common center, which possibly denotes the position of the YSO driving the maser motion. Water masers trace both a fast shock (up to 50 km/s) closer to the YSO, powered by a wide-angle wind, and a slower (20 km/s) bipolar jet, at the base of the large-scale outflow. Since the compact free-free emission is found offset from the putative location of the YSO along a direction consistent with that of the maser jet axis, we interpret the radio continuum in terms of a thermal jet. The velocity field of methanol masers can be explained in terms of a composition of slow (4 km/s in amplitude) motions of radial expansion and rotation about an axis approximately parallel to the maser jet. Finally, the distribution of line of sight velocities of the hydroxyl masers suggests that they can trace gas less dense (n(H2) < 10^6 cm^-3) and more distant from the YSO than that traced by the water and methanol masers, which is expanding toward the observer. (Abridged)

We derive accurate proper motions of the CH3OH 12 GHz masers towards the W3(OH) UCHII region, employing seven epochs of VLBA observations spanning a time interval of about 10 yr. The achieved velocity accuracy is of the order of 0.1 km/s, adequate to precisely measure the relative velocities of most of the 12 GHz masers in W3(OH), with amplitude varying in the range 0.3 - 3 km/s. Towards W3(OH), the most intense 12 GHz masers concentrate in a small area towards the north (the northern clump) of the UCHII region. We have compared the proper motions of the CH3OH 12 GHz masers with those (derived from literature data) of the OH 6035 MHz masers, emitting from the same region of the methanol masers. In the northern clump, the two maser emissions emerge from nearby (but likely distinct) cloudlets of masing gas with, in general, a rather smooth variation of line-of-sight and sky-projected velocities, which suggests some connection of the environments and kinematics traced by both maser types. The conical outflow model, previously proposed to account for the 12 GHz maser kinematics in the northern clump, does not reproduce the new, accurate measurements of 12 GHz maser proper motions and has to be rejected. We focus on the subset of 12 GHz masers of the northern clump belonging to the "linear structure at P.A. = 130-140 degree", whose regular variation of LSR velocities with position presents evidence for some ordered motion. We show that the 3-dimensional velocities of this "linear distribution" of 12GHz masers can be well fitted considering a flat, rotating disk, seen almost edge-on.

The present paper focuses on the high-mass star-forming region G16.59-0.05. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G16.59-0.05: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the star-forming region at 1.3 and 3.6 cm. Results: This is the first work to report accurate measurements of the "relative" proper motions of the 6.7 GHz CH3OH masers. The different spatial and 3-D velocity distribution clearly indicate that the 22 GHz water and 6.7 GHz methanol masers are tracing different kinematic environments. The bipolar distribution of 6.7 GHz maser l.o.s. velocities and the regular pattern of observed proper motions suggest that these masers are tracing rotation around a central mass of about 35 solar masses. The flattened spatial distribution of the 6.7 GHz masers, oriented NW-SE, suggests that they can originate in a disk/toroid rotating around the massive YSO which drives the 12CO(2-1) outflow, oriented NE-SW, observed on arcsec scale. The extended, radio continuum source observed close to the 6.7 GHz masers could be excited by a wide-angle wind emitted from the YSO associated with the methanol masers, and such a wind is proven to be sufficiently energetic to drive the NE-SW 12CO(2-1) outflow. The H2O masers distribute across a region offset about 0".5 to the NW of the CH3OH masers, in the same area where emission of high-density molecular tracers, typical of HMCs, was detected. We postulate that a distinct YSO, possibly in an earlier evolutionary phase than that exciting the methanol masers, is responsible for the excitation of the water masers and the HMC molecular lines. (Abridged)

We report a trigonometric parallax of 12 GHz methanol masers associated with the massive star forming region G9.62+0.20, corresponding to a distance of 5.2+-0.6 kpc. With a local standard of rest velocity of about 2 km/s, the region's kinematic distances of 0.5 and 16 kpc differ greatly from the distance derived here. Our measurement of the peculiar motion of the star forming region shows a very large deviation from a circular Galactic orbit: 41 km/s radially outward from the Galactic center and 60 km/s counter to Galactic rotation. The combination of its radial velocity and distance places G9.62+0.20 in the inner region of the Galaxy close to the Expanding Near 3 kpc-Arm, where the bulge/bar potential has strong gravitational influence. We also map the distribution of 12 GHz methanol masers, locate them with respect to a hypercompact H II region, and compare our data with the periodic flare phenomenon reported previously for this source.

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 report trigonometric parallaxes for the high-mass star forming regions G35.20-0.74 and G35.20-1.74, corresponding to distances of 2.19 (+0.24 -0.20) kpc and 3.27 (+0.56 -0.42) kpc, respectively. The distances to both sources are close to their near kinematic distances and place them in the Carina-Sagittarius spiral arm. Combining the distances and proper motions with observed radial velocities gives the locations and full space motions of the star forming regions. Assuming a standard model of the Galaxy, G35.20-0.74 and G35.20-1.74 have peculiar motions of ~13 km/s and ~16 km/s counter to Galactic rotation and ~9 km/s toward the North Galactic Pole.

We report trigonometric parallaxes for the massive star-forming regions G23.01-0.41 and G23.44-0.18, corresponding to distances of 4.59 [+0.38/-0.33] kpc and 5.88 [+1.37/-0.93] kpc, respectively. The distance to G23.01-0.41 is smaller than its near kinematic distance assuming a standard model of the Milky Way and less than half of its far kinematic distance, which has usually been assumed. This places it in the Crux-Scutum spiral arm. The distance to G23.44-0.18 is close to its near kinematic distance and most likely places it in the Norma spiral arm near the end of the Galactic bar. Combining the distance and proper motions with observed radial velocities gives the location and full space motion of the star forming regions. We find large deviations from circular Galactic orbits for these sources: both sources show peculiar motions of 20 to 30 km/s counter to Galactic rotation and toward the Galactic center. These large peculiar motions might be the result of gravitational perturbations from the Galactic bar.

We report trigonometric parallaxes for the sources NGC 7538 and Cep A, corresponding to distances of 2.65 [+0.12/-0.11] kpc and 0.70 [+0.04/-0.04] kpc, respectively. The distance to NGC 7538 is considerably smaller than its kinematic distance and places it in the Perseus spiral arm. The distance to Cep A is also smaller than its kinematic distance and places it in the Local arm or spur. Combining the distance and proper motions with observed radial velocities gives the location and full space motion of the star forming regions. We find significant deviations from circular Galactic orbits for these sources: both sources show large peculiar motions (> 10 km/s) counter to Galactic rotation and NGC 7538 has a comparable peculiar motion toward the Galactic center.

We report trigonometric parallaxes for G59.7+0.1 and W 51 IRS2, corresponding to distances of 2.16^{+0.10}_{-0.09} kpc and 5.1^{+2.9}_{-1.4} kpc, respectively. The distance to G59.7+0.1 is smaller than its near kinematic distance and places it between the Carina-Sagittarius and Perseus spiral arms, probably in the Local (Orion) spur. The distance to W 51 IRS2, while subject to significant uncertainty, is close to its kinematic distance and places it near the tangent point of the Carina-Sagittarius arm. It also agrees well with a recent estimate based on O-type star spectro/photometry. Combining the distances and proper motions with observed radial velocities gives the full space motions of the star forming regions. We find modest deviations of 5 to 10 km/s from circular Galactic orbits for these sources, both counter to Galactic rotation and toward the Galactic center.

We are conducting a large program with the NRAO Very Long Baseline Array (VLBA) to measure trigonometric parallaxes of massive star-forming regions across the Milky Way. Here we report measurement of the parallax and proper motion of methanol masers in S 252 and G232.6+1.0. The parallax of S 252 is 0.476 +/- 0.006 mas (2.10 [+0.027/-0.026] kpc), placing it in the Perseus spiral arm. The parallax of G232.6+1.0 is 0.596 +/- 0.035 mas (1.68 [+0.11/-0.09] kpc), placing it between the Carina-Sagittarius and Perseus arms, possibly in a Local (Orion) spur of the Carina-Sagittarius arm. For both sources, kinematic distances are significantly greater than their parallax distances. Our parallaxes and proper motions yield full space motions accurate to about 1 km/s. Both sources orbit the Galaxy about 13 km/s slower than circular rotation.

We present new radio continuum observations, taken at 1.4 GHz with the Very Large Array, and new observations in the CO line, taken with the Plateau de Bure interferometer, of two regions of the spiral galaxy NGC 3627. Comparing these observations with archival Spitzer and H_{alpha} data we study the correlation of the radio continuum (RC), infrared-8 um and CO emissions at a spatial resolution of 100 pc. We compare the point-by-point variations of the RC, CO, and 8um brightnesses in two distinct regions of 2 kpc \times 2 kpc in size of NGC 3627. We examined scale much lower than the electron diffusion scale, where a break down of the correlations would be expected. However no evidence for such correlation break down has been found. The RC emission follows well the distribution of CO and the widths of several structures, measured along slices across them, are comparable. Furthermore, we found that down to a spatial scale of 100 pc the radio continuum emission is correlated with dust emissions at 4.5, 5.8, and 8 um, that traces different dust temperatures. We present a new perspective with which visualizing and studying the RC-CO-24 um and RC-CO-8 um correlations, by using a three-dimensional representation. We fit the observed data with a three dimensional line, obtaining a rms of 0.25 dex. The observed correlation enhances the complexity of the electrons diffusion, losses and injection mechanisms, and of their connection to star formation processes, described by molecular and dust emissions. We plan to further investigate this connection using spatially resolved spectral index studies at low radio frequencies where the thermal emission component is seemingly negligible.

Context. Water vapor emission at 22 GHz from masers associated with star-forming regions is highly variable. Aims. We present a database of up to 20 years of monitoring of a sample of 43 masers within star-forming regions. The sample covers a large range of luminosities of the associated IRAS source and is representative of the entire population of H2O masers of this type. The database forms a good starting point for any further study of H2O maser variability. Methods. The observations were obtained with the Medicina 32-m radiotelescope, at a rate of 4-5 observations per year. Results. To provide a database that can be easily accessed through the web, we give for each source: plots of the calibrated spectra, the velocity-time-flux density plot, the light curve of the integrated flux, the lower and upper envelopes of the maser emission, the mean spectrum, and the rate of the maser occurrence as a function of velocity. Figures for just one source are given in the text for representative purposes. Figures for all the sources are given in electronic form in the on-line appendix. A discussion of the main properties of the H2O variability in our sample will be presented in a forthcoming paper.

Context. G24.78+0.08 A1 is a 20 Msun star surrounded by a hypercompact (HC) HII region, driving a CO bipolar outflow, and located at the center of a massive rotating toroid undergoing infall towards the HC region. Recent water maser observations suggest that the HC region is expanding and accretion onto the star is halted. Aims. This study aims to confirm the expansion scenario proposed for the HC region on the basis of recent H2O maser observations. Methods. We carried out continuum VLA observations at 1.3cm and 7mm with the A array plus Pie Town configuration to map the HC region towards G24 A1. Results. The emission of the HC region has been resolved and shows a ring shape structure. The profiles of the emission obtained by taking slices at different angles passing through the barycenter of the HC region confirm the shell structure of the emission. The ratio between the inner and the outer radius of the shell, Ri/Ro, derived fitting the normalized brightness temperature profile passing through the peak of the 7mm emission, is 0.9, which indicates that the shell is thin. The deconvolved outer radius estimated from the fit is 590 AU. These results imply that the HC region in G24 A1 cannot be described in terms of a classical, homogeneous HII region but is instead an ionized shell. This gives support to the model of an expanding wind-driven, ionized shell suggested by the kinematics and distribution of the H2O masers associated with the HC region. According to this model, the HC region is expanding on very short times scales, 21-66 yr.

Most previous high-angular (<0.1 arcsec) resolution studies of molecular masers in high-mass star forming regions (SFRs) have concentrated mainly on either water or methanol masers. While high-angular resolution observations have clarified that water masers originate from shocks associated with protostellar jets, different environments have been proposed in several sources to explain the origin of methanol masers. Tha aim of the paper is to investigate the nature of the methanol maser birthplace in SFRs and the association between the water and methanol maser emission in the same young stellar object. We have conducted phase-reference Very Long Baseline Interferometry (VLBI) observations of water and methanol masers toward two high-mass SFRs, Sh 2-255 IR and AFGL 5142. In Sh 2-255 IR water masers are aligned along a direction close to the orientation of the molecular outflow observed on angular scales of 1-10 arcsec, tracing possibly the disk-wind emerging from the disk atmosphere. In AFGL 5142 water masers trace expansion at the base of a protostellar jet, whilst methanol masers are more probably tracing infalling than outflowing gas. The results for AFGL 5142 suggest that water and methanol masers trace different kinematic structures in the circumstellar gas.

Aims: We present new radio continuum observations of ten BIMA SONG galaxies, taken at 1.4 GHz with the Very Large Array. These observations allow us to extend the study of the relationships between the radio continuum (RC) and CO emission to 22 CO luminous galaxies for which single dish CO images have been added to interferometric data. New Spitzer infrared (IR) images of six of these galaxies have been released. The analysis of these high resolution images allowed us to probe the RC-IR-CO correlations down to linear scales of a few hundred pc. Results: for the 22 galaxies analysed, the RC-CO correlation on scales from $\sim10$ kpc down to $\sim100$ pc is nearly linear and has a scatter of a factor of two, i.e. comparable to that of the global correlations. There is no evidence for any severe degradation of the scatter below the kpc scale. This also applies to the six galaxies for which high-resolution mid-IR data are available. In the case of NGC 5194, we find that the non-thermal radio spectral index is correlated with the RC/FIR ratio. Conclusions: The scatter of the point-by-point correlations does not increase significantly with spatial resolution. We thus conclude that we have not yet probed the physical scales at which the correlations break down. However, we observe local deviations from the correlations in regions with a high star formation rate, such as the spiral arms, where we observe a flat radio spectrum and a low RC/FIR ratio. In the intra-arm regions and in the peripheral regions of the disk, the RC/FIR is generally higher and it is characterized by a steepening of the radio spectrum.

We have conducted phase-reference multi-epoch observations of the 22.2 GHz water masers using the VLBA and multi-frequency study of the continuum emission using the VLA towards the high-mass SFR AFGL 5142. The water maser emission comes from two elongated structures (indicated as Group I and Group II), with the measured proper motions aligned along the structures' elongation axes. Each group consists of two (blue- and red-shifted) clusters of features separated by a few hundreds and thousands of AU respectively for Group I and Group II. The maser features of Group II have both positions and velocities aligned along a direction close to the axis of the outflow traced by HCO+ and SiO emission on angular scales of tens of arcsec. We predict that the maser emission arises from dense, shocked molecular clumps displaced along the axis of the molecular outflow. The two maser clusters of Group I are oriented on the sky along a direction forming a large angle (> 60 degrees) with the axis of the jet/outflow traced by Group II maser features. We have detected a compact (8.4 and 22 GHz) continuum source that falls close to the centroid of Group I masers, indicating that the source ionizing the gas is also responsible for the excitation of the water masers. The kinematic analysis indicates that the Group I masers trace outflowing rather than rotating gas, discarding the Keplerian disk scenario proposed in a previous paper for Group I. Since the axis joining the two maser clusters of Group II does not cross the position of the continuum source, Group II masers might be excited by an (undetected) massive YSO, distinct from the one (pinpointed by the VLA continuum emission) responsible for the excitation of the Group I masers.

This article reports first results of a long-term observational program aimed to study the earliest evolution of jet/disk systems in low-mass YSOs by means of VLBI observations of the 22.2 GHz water masers. We report here data for the cluster of low-mass YSOs in the Serpens molecular core and for the single object RNO~15-FIR. Towards Serpens SMM1, the most luminous sub-mm source of the Serpens cluster, the water maser emission comes from two small (< 5 AU in size) clusters of features separated by ~25 AU, having line of sight velocities strongly red-shifted (by more than 10 km/s) with respect to the LSR velocity of the molecular cloud. The two maser clusters are oriented on the sky along a direction that is approximately perpendicular to the axis of the radio continuum jet observed with the VLA towards SMM1. The spatial and velocity distribution of the maser features lead us to favor the interpretation that the maser emission is excited by interaction of the receding lobe of the jet with dense gas in the accretion disk surrounding the YSO in SMM1. Towards RNO~15-FIR, the few detected maser features have both positions and (absolute) velocities aligned along a direction that is parallel to the axis of the molecular outflow observed on much larger angular scales. In this case the maser emission likely emerges from dense, shocked molecular clumps displaced along the axis of the jet emerging from the YSO. The protostar in Serpens SMM1 is more massive than the one in RNO~15-FIR. We discuss the case where a high mass ejection rate can generate jets sufficiently powerful to sweep away from their course the densest portions of circumstellar gas. In this case, the excitation conditions for water masers might preferably occur at the interface between the jet and the accretion disk, rather than along the jet axis.

We report on the first multi-epoch, phase referenced VLBI observations of the water maser emission in a high-mass protostar associated with a disk-jet system. The source under study, IRAS 20126+4104, has been extensively investigated in a large variety of tracers, including water maser VLBA data acquired by us three years before the present observations. The new findings fully confirm the interpretation proposed in our previous study, namely that the maser spots are expanding from a common origin coincident with the protostar. We also demonstrate that the observed 3-D velocities of the maser spots can be fitted with a model assuming that the spots are moving along the surface of a conical jet, with speed increasing for increasing distance from the cone vertex. We also present the results of single-dish monitoring of the water maser spectra in IRAS 20126+4104. These reveal that the peak velocity of some maser lines decreases linearly with time. We speculate that such a deceleration could be due to braking of the shocks from which the maser emission originates, due to mass loading at the shock front or dissipation of the shock energy.

We have studied the relationships between the radio continuum (RC) and CO emission for a set of galaxies selected from the BIMA Survey of Nearby Galaxies. We find that the global CO-RC correlation is as tight as the global FIR-RC correlation for the 24 galaxies studied. Within 9 galaxies with ~6 arcsec CO and RC data available, the CO and RC emission is as tightly correlated as its global value; the radially averaged correlation is nearly linear, extends over four order of magnitude and holds down to the smallest linear resolution of the observations, which is ~100 pc. We define qco as the log of the ratio of the CO to RC flux as a way to characterize the CO-RC correlation. Combining 6 arcsec pixel-by-pixel comparisons across all sources yields an average small-scale correlation of qco = 1.1 +/- 0.28; that is, the spatially resolved correlation has a dispersion that is less than a factor of 2. There are however systematic variations in the CO/RC ratio; the strongest organized structures in qco tend to be found along spiral arms and on size scales much larger than the resolution of the observations. We do not measure any systematic trend in CO/RC ratio as a function of radius in galaxies. The constancy of the CO/RC ratio stands in contrast to the previously measured decrease in the FIR/RC ratio as a function of radius in galaxies. We suggest that the excellent correlation between the CO, RC and FIR emission in galaxies is a consequence of regulation by hydrostatic pressure; this model links all three emissions without invoking an explicit dependence on a star formation scenario.

Extragalactic water vapor masers with 50, 1000, 1, and 230 solar (isotropic) luminosities were detected toward Mrk1066 (UGC2456), Mrk34, NGC3556 and Arp299, respectively. The interacting system Arp299 appears to show two maser hotspots separated by 20 arcsec. A statistical analysis of 53 extragalactic H2O sources indicates (1) that the correlation between IRAS Point Source and H2O luminosities, established for individual star forming regions in the galactic disk, also holds for AGN dominated megamaser galaxies, (2) that maser luminosities are not correlated with 60/100 micron color temperatures and (3) that only a small fraction of the luminous megamasers detectable with 100-m sized telescopes have so far been identified. The slope of the H2O luminosity function, -1.5, indicates that the number of detectable masers is almost independent of their luminosity. If the LF is not steepening at very high maser luminosities, H2O megamasers at significant redshifts should be detectable with present day state-of-the-art facilities.

We present VLBI observations of the H$_2$O maser emission towards a selection of massive young stellar objects (YSOs). In one of these sources, IRAS20126+4104, the measured proper motions confirm that the water masers spots are tracing the very root (around 100 AU) of a conical bipolar outflow emerging near the position of the embedded YSO, and are an excellent tool to investigate the structure and kinematics of the outflow/jet. We also present preliminary results of a first epoch EVN observations of a selection of 5 high-mass YSO to assess the precise association of maser spots and molecular outflows in these sources.

We have conducted multi-epoch EVN observations of the 22.2 GHz water masers towards four high-mass star forming regions (Sh 2-255 IR, IRAS 23139+5939, WB89-234, and OMC2). The (three) observing epochs span a time range of 6 months. In each region, the water maser emission likely originates close (within a few hundreds of AU) to a forming high-mass YSO. Several maser features (~10) have been detected for each source and, for those features persistent over the three epochs, proper motions have been derived. The amplitudes of the proper motions are found to be larger than the range of variation of the line-of-sight velocities and in each of the observed sources the proper motion orientation seems to indicate an expansion motion. The gas kinematics traced by the 22.2 GHz water masers is compatible with the shock-excited nature of water maser emission. Three different kinematic models (a spherical expanding shell, a Keplerian rotating disk, and a conical outflow) were fitted to the 3-dimensional velocity field of the detected maser features. The results of these fits, together with the comparison of the VLBI maps with the highest-resolution images of the sources in several thermal tracers, suggest that the water maser features are most likely tracing the inner portion of the molecular outflows detected at much larger-scales.

We report the first detection of a water megamaser in a radio-loud galaxy, 3C403, and present a follow-up study using the VLA. 3C403 has been observed as a part of a small sample of FRII galaxies with evidence of nuclear obscuration. The isotropic luminosity of the maser is ~1200 solar luminosities. With a recessional velocity of cz~17680 km/s it is the most distant water maser so far reported. The line arises from the densest (> 10^(8) cm^(-3)) interstellar gas component ever observed in a radio-loud galaxy. Two spectral features are identified, likely bracketing the systemic velocity of the galaxy. Our interferometric data clearly indicate that these arise from a location within 0.1" (~110 pc) from the active galactic nucleus. We conclude that the maser spots are most likely associated with the tangentially seen parts of a nuclear accretion disk, while an association with dense warm gas interacting with the radio jets cannot yet be ruled out entirely.

We have conducted multi-epoch EVN observations of the 22.2 GHz water masers towards the high-mass young stellar object in AFGL 5142. With four observing epochs, spanning a time of ~1 year, 12 distinct maser features have been detected, 7 out of these detected in more than one epoch. The positions and the velocities of the VLBI features agree well with those of the emission centers previously identified by means of VLA observations. For a few features, persistent over three or four epochs, accurate values of the proper motions are derived. On the basis of their spatial distribution, the observed maser features can be divided into two groups. A model fit to the positions and velocities of the maser features of Group I, detected in the same region (within ~500 mas) where the massive YSO should be located, demonstrates that these might arise on the surface of a nearly edge-on Keplerian disk, rotating around a massive young stellar object. The maser features of Group II, found at large distances from the YSO (> 1 arcsec), have positions and line-of-sight velocities in agreement with the blue-shifted lobe of a large scale molecular outflow (traced by the HCO+ and SiO emission), and might result from the interaction between the gas flowing away from the young stellar object and the ambient gas of the progenitor molecular core.

We report the results of a search for 22 GHz water maser emission in IRAS-bright galaxies, using the 100-m Effelsberg telescope. In particular, we present the details of four new maser detections (IC342, NGC2146, NGC3556, and Arp299) and follow-up interferometric studies. A comparison between water maser detection rates derived in the present study and those in previous similar surveys is also presented.

We report the first results of an ongoing survey at 22 GHz with the 100-m Effelsberg telescope to search for water maser emission in bright IRAS sources. We have detected water vapor emission in IC 342. The maser, associated with a star forming region ~10-15 arcsec west of the nucleus, consists of a single 0.5 km/s wide feature and reaches an isotropic luminosity of 0.01 L(sun) (D = 1.8 Mpc). Our detection raises the detection rate among northern galaxies with IRAS point source fluxes S(100micron) > 50 Jy to 16%.

We combine the radio continuum images from the NRAO VLA Sky Survey with the CO-line observations from the extragalactic CO survey of the Five College Radio Astronomy Observatory to study the relationship between molecular gas and the star formation rate within the disks of 180 spiral galaxies at 45" resolution. We find a tight correlation between these quantities. On average, the ratio between the radio continuum and the CO emission is constant, within a factor of 3, both inside the same galaxy and from galaxy to galaxy. The mean star formation efficiency deduced from the radio continuum corresponds to convert 3.5% of the available molecular gas into stars on a time scale of 10^8 yr and depends weakly on general galaxy properties, such as Hubble type or nuclear activity. A comparison is made with another similar analysis performed using the H_alpha luminosity as star formation indicator. The overall agreement we find between the two studies reinforces the use of the radio luminosity as star formation rate indicator not only on global but also on local scales.