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.

The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region.

Outside the Milky Way, the most luminous H2O masers at 22 GHz, called 'megamasers' because of their extreme luminosity with respect to the Galactic and extragalactic H2O masers associated with star formation, are mainly detected in active galactic nuclei. In the case of the H2O maser detected in the nuclear region of the galaxy TXS2226-184 for the first time the term 'gigamaser' was used. However, the origin of this very luminous H2O maser emission has never been investigated into details. We study the nature of the H2O gigamaser in TXS2226-184 by measuring for the first time its absolute position at mas resolution, by comparing the morphology and characteristics of the maser emission at the VLBI scales after about 20 years, and by trying to detect its polarized emission. We observed the maser emission towards TXS2226-184 three times: the very first one with the VLBA (epoch 2017.45) and the next two times with the EVN (epochs 2017.83 and 2018.44). The first two epochs were observed in phase-reference mode, while the last epoch was observed in full-polarization mode but not in phase-reference mode to increase the on-source integration time. We also retrieved and analyzed the VLBA archival data at 22 GHz of TXS2226-184 observed in 1998.40. We detected 6 H2O maser features in epoch 2017.45 (VLBA), one in epoch 2017.83 (EVN), and two in epoch 2018.44 (EVN). All of them but one are red-shifted with respect to the systemic velocity of TXS2226-184, we detected only one blue-shifted maser feature and it is the weakest one. For the first time, we were able to measure the absolute position of the H2O maser features with errors below 1 mas. No linear and circular polarization was detected. We were able to associate the H2O maser features in TXS2226-184 with the most luminous radio continuum clump reported in the literature.

Almost all dwarf galaxies in the Local Group that are not satellites of the Milky Way or M31, belong to either one of two highly-symmetric planes. It is still a matter of debate, whether these planar structures are dynamically stable or whether they only represent a transient alignment. Proper motions, if they could be measured, could help to discriminate between these scenarios. Such motions could be determined with multi-epoch Very Long Baseline Interferometry (VLBI) of sources that show emission from water and methanol at frequencies of 22 and 6.7 GHz, respectively. We report searches for such masers. We have mapped three Local Group galaxies, NGC6822, IC1613 and WLM in the bands covering the water vapor and methanol lines. These systems are members of the two above mentioned planes of galaxies. We have produced deep radio continuum (RC) maps and spectral line cubes. The former have been used to identify star forming regions and to derive global galactic star formation rates (SFRs). These SFRs turn out to be lower than those determined at other wavelengths in two of our sources. This indicates that dwarf galaxies may follow predictions on the RC-SFR relation only in individual regions of enhanced radio continuum emission, but not when considering the entire optical body of the sources. No methanol or water maser emission has been confidently detected, down to line luminosity limits of ~0.004 and 0.01 solar luminosities, respectively. This finding is consistent with the small sizes, low SFRs and metallicities of these galaxies.

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

Investigations of H$_2$O maser galaxies at X-ray energies reveal that most AGN associated with water masers are characterized by high levels of absorption. With the aim of finding new maser sources for possible interferometric follow-ups, we have searched for water maser emission in a well-defined sample of heavily absorbed AGN ($N_{\rm H} > 10^{23}$ cm$^{-2}$), including Compton-thick (CT) sources. All the galaxies in the sample were already searched for 22 GHz water maser emission in previous surveys. With the goal of providing a detection or a stringent upper limit on the H$_2$O luminosity, we re-observed some of the non-detected sources with the Green Bank Telescope. A new luminous H$_2$O maser ($L_{\rm H2O} \sim 200\,$L$_\odot$) was detected in the mid-IR-bright Seyfert 2 galaxy IRAS 15480-0344 and then followed-up with the Very Long Baseline Array. In order to shed light on the origin of the maser (jet/outflow vs. disc), we recently observed the radio continuum emission in IRAS 15480-0344 with the European VLBI network (EVN) at 1.7 and 5.0 GHz. With the newly discovered megamaser in IRAS 15480-0344 revealing a narrow ($\sim$0.6 km s$^{-1}$) and a broad ($\sim$90 km s$^{-1}$) component, the maser detection rate of the CT AGN sample is 50% (18/36), which is one of the highest ever found in maser surveys. The EVN maps show two bright sources (labeled SW and NE) in the nuclear region of IRAS 15480-0344, which we interpret as jet knots tracing regions where the radio plasma impacts dense molecular clouds. [abridged]

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

Context. The role of magnetic fields during the formation of high-mass stars is not yet fully understood, and the processes related to the early fragmentation and collapse are largely unexplored today. The high-mass star forming region G9.62+0.19 is a well known source, presenting several cores at different evolutionary stages. Aims. We determine the magnetic field morphology and strength in the high-mass star forming region G9.62+0.19, to investigate its relation to the evolutionary sequence of the cores. Methods. We use Band 7 ALMA observations in full polarisation mode and we analyse the polarised dust emission. We estimate the magnetic field strength via the Davis-Chandrasekhar-Fermi and the Structure Function methods. Results. We resolve several protostellar cores embedded in a bright and dusty filamentary structure. The polarised emission is clearly detected in six regions. Moreover the magnetic field is oriented along the filament and appears perpendicular to the direction of the outflows. We suggest an evolutionary sequence of the magnetic field, and the less evolved hot core exhibits a magnetic field stronger than the more evolved one. We detect linear polarisation from thermal line emission and we tentatively compared linear polarisation vectors from our observations with previous linearly polarised OH masers observations. We also compute the spectral index, the column density and the mass for some of the cores. Conclusions. The high magnetic field strength and the smooth polarised emission indicate that the magnetic field could play an important role for the fragmentation and the collapse process in the star forming region G9.62+019 and that the evolution of the cores can be magnetically regulated. On average, the magnetic field derived by the linear polarised emission from dust, thermal lines and masers is pointing in the same direction and has consistent strength.

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

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

Magnetic fields play an important role during star formation. Direct magnetic field strength observations have proven specifically challenging in the extremely dynamic protostellar phase. Because of their occurrence in the densest parts of star forming regions, masers, through polarization observations, are the main source of magnetic field strength and morphology measurements around protostars. Of all maser species, methanol is one of the strongest and most abundant tracers of gas around high-mass protostellar disks and in outflows. However, as experimental determination of the magnetic characteristics of methanol has remained largely unsuccessful, a robust magnetic field strength analysis of these regions could hitherto not be performed. Here we report a quantitative theoretical model of the magnetic properties of methanol, including the complicated hyperfine structure that results from its internal rotation. We show that the large range in values of the Land\'{e} g-factors of the hyperfine components of each maser line lead to conclusions which differ substantially from the current interpretation based on a single effective g-factor. These conclusions are more consistent with other observations and confirm the presence of dynamically important magnetic fields around protostars. Additionally, our calculations show that (non-linear) Zeeman effects must be taken into account to further enhance the accuracy of cosmological electron-to-proton mass ratio determinations using methanol.

Context. The importance of the magnetic field in high-mass-star formation is not yet fully clear and there are still many open questions concerning its role in the accretion processes and generation of jets and outflows. In the past few years, masers have been successfully used to probe the magnetic field morphology and strength at scales of a few au around massive protostars, by measuring linear polarisation angles and Zeeman splitting. The massive protostar IRAS 18089-1732 is a well studied high-mass-star forming region, showing a hot core chemistry and a disc-outflow system. Previous SMA observations of polarised dust revealed an ordered magnetic field oriented around the disc of IRAS 18089-1732. Aims. We want to determine the magnetic field in the dense region probed by 6.7 GHz methanol maser observations and compare it with observations in dust continuum polarisation, to investigate how the magnetic field in the compact maser region relates to the large-scale field around massive protostars. Methods. We reduced MERLIN observations at 6.7 GHz of IRAS 18089-1732 and we analysed the polarised emission by methanol masers. Results. Our MERLIN observations show that the magnetic field in the 6.7 GHz methanol maser region is consistent with the magnetic field constrained by the SMA dust polarisation observations. A tentative detection of circularly polarised line emission is also presented. Conclusions. We found that the magnetic field in the maser region has the same orientation as in the disk. Thus the large-scale field component, even at the au scale of the masers, dominates over any small-scale field fluctuations. We obtained, from the circular polarisation tentative detection, a field strength along the line of sight of 5.5 mG which appeared to be consistent with the previous estimates.

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.

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.

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 current paradigm of star formation through accretion disks, and magnetohydrodynamically driven gas ejections, predicts the development of collimated outflows, rather than expansion without any preferential direction. We present radio continuum observations of the massive protostar W75N(B)-VLA 2, showing that it is a thermal, collimated ionized wind and that it has evolved in 18 years from a compact source into an elongated one. This is consistent with the evolution of the associated expanding water-vapor maser shell, which changed from a nearly circular morphology, tracing an almost isotropic outflow, to an elliptical one outlining collimated motions. We model this behavior in terms of an episodic, short-lived, originally isotropic, ionized wind whose morphology evolves as it moves within a toroidal density stratification.

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

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

The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

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

W75N(B) is a massive star-forming region that contains three radio continuum sources (VLA 1, VLA 2, and VLA 3), which are thought to be three massive young stellar objects at three different evolutionary stages. VLA 1 is the most evolved and VLA 2 the least evolved source. The 22 GHz H2O masers associated with VLA 1 and VLA 2 have been mapped at several epochs over eight years. While the H2O masers in VLA 1 show a persistent linear distribution along a radio jet, those in VLA 2 are distributed around an expanding shell. Furthermore, H2O maser polarimetric measurements revealed magnetic fields aligned with the two structures. Using new polarimetric observations of H2O masers, we aim to confirm the elliptical expansion of the shell-like structure around VLA 2 and, at the same time, to determine if the magnetic fields around the two sources have changed. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz H2O masers towards the massive star-forming region W75N(B). The H2O maser distribution around VLA 1 is unchanged from that previously observed. We made an elliptical fit of the H2O masers around VLA 2. We find that the shell-like structure is still expanding along the direction parallel to the thermal radio jet of VLA 1. While the magnetic field around VLA 1 has not changed in the past 7 years, the magnetic field around VLA 2 has changed its orientation according to the new direction of the major-axis of the shell-like structure and it is now aligned with the magnetic field in VLA 1.

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.

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields play an important role in the formation and dynamics of molecular outflows. Masers, in particular 6.7-GHz CH3OH masers, are the best probes of the magnetic field morphologies around massive young stellar objects on the smallest scales of 10-100 AU. This paper focuses on 4 massive young stellar objects, IRAS06058+2138-NIRS1, IRAS22272+6358A, S255-IR, and S231, which complement our previous 2012 sample (the first EVN group). From all these sources, molecular outflows have been detected in the past. Seven of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz CH3OH masers in the star-forming regions in this second EVN group. We detected a total of 128 CH3OH masing cloudlets. Fractional linear polarization (0.8%-11.3%) was detected towards 18% of the CH3OH masers in our sample. The linear polarization vectors are well ordered in all the massive young stellar objects. We measured significant Zeeman-splitting in IRAS06058+2138-NIRS1 (DVz=3.8+/-0.6 m/s) and S255-IR (DVz=3.2+/-0.7 m/s). By considering the 20 massive young stellar objects towards which the morphology of magnetic fields was determined by observing 6.7-GHz CH3OH masers in both hemispheres, we find no evident correlation between the linear distributions of CH3OH masers and the outflows or the linear polarization vectors. On the other hand, we present first statistical evidence that the magnetic field (on scales 10-100 AU) is primarily oriented along the large-scale outflow direction. Moreover, we empirically find that the linear polarization fraction of unsaturated CH3OH masers is P_l<4.5%.

We present Very Long Baseline Interferometry (VLBI) observations of 22 GHz H$_2$O masers in the high-mass star-forming region of \objectname{W75N}, carried out with VLBI Exploration of Radio Astrometry (VERA) for three-epochs in 2007 with an angular resolution of $\sim$ 1 mas. We detected H$_2$O maser emission toward the radio jet in VLA 1 and the expanding shell-like structure in VLA 2. We have made elliptical fits to the VLA 2 H$_2$O maser shell-like structure observed in the different epochs (1999, 2005, and 2007), and found that the shell is still expanding eight years after its discovery. From the difference in the size of the semi-major axes of the fitted ellipses in the epochs 1999 ($\simeq$ 71$\pm$1 mas), 2005 ($\simeq$ 97$\pm$3 mas), and 2007 ($\simeq$ 111$\pm$1 mas), we estimate an average expanding velocity of $\sim$ 5 mas yr$^{-1}$, similar to the proper motions measured in the individual H$_2$O maser features. A kinematic age of $\sim$ 20 yr is derived for this structure. In addition, our VERA observations indicate an increase in the ellipticity of the expanding shell around VLA 2 from epochs 1999 to 2007. In fact, the elliptical fit of the VERA data shows a ratio between the minor and major axes of $\sim$ 0.6, in contrast with a most circular shape for the shell detected in 1999 and 2005 (b/a $\sim$ 0.9). This suggests that we are probably observing the formation of a jet-driven H$_2$O maser structure in VLA2, evolving from a non-collimated pulsed-outflow event during the first stages of evolution of a massive young stellar object (YSO). This may support predictions made earlier by other authors on this issue, consistent with recent magnetohydrodynamical simulations. We discuss possible implications of our results in the study of the first stages of evolution of massive YSOs.

The role of magnetic fields in the formation of high-mass stars is still under debate, and recent measurements of their orientation and strength by using polarized maser emissions are contributing new insights. Masers polarization, in particular of the 6.7-GHz methanol masers, are one of the best probes of the magnetic field morphologies around massive protostars. Determining the magnetic field morphology around an increasing number of massive protostars at milliarcsecond resolution by observing 6.7-GHz methanol masers is crucial to better understand the role of magnetic fields in massive star formation.The First EVN Group consists of 4 massive star-forming complexes: W51, W48, IRAS18556+0138, and W3(OH). These contain well-studied \hii ~regions from some of which molecular bipolar outflows were also detected (W51-e2, G35.20-0.74N). Nine of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz methanol masers in the star-forming regions of the First EVN Group. We detected a total of 154 CH3OH masers, one third of these towards W3(OH). Fractional linear polarization (1.2-11.5%) was detected towards 55 masers. The linear polarization vectors are well-ordered in all the massive star-forming regions. We measured significant Zeeman-splitting in 3 massive star-forming regions (W51, W48, and W3(OH)) revealing a range of separations -3.5 m/s<\Delta V_{z}<3.8 m/s with the smallest |\Delta V_{z}|=0.4m/s. We were also able to compare our magnetic field results with those obtained from submillimeter wavelength dust observation in W51 and show that the magnetic field at low and high resolutions are in perfect agreement.

NGC7538 is a complex massive star-forming region. The region is composed of several radio continuum sources, one of which is IRS1, a high-mass protostar, from which a 0.3 pc molecular bipolar outflow was detected. Several maser species have been detected around IRS1. The CH3OH masers have been suggested to trace a Keplerian-disk, while the H2O masers are almost aligned to the outflow. More recent results suggested that the region hosts a torus and potentially a disk, but with a different inclination than the Keplerian-disk that is supposed to be traced by the CH3OH masers. Tracing the magnetic field close to protostars is fundamental for determining the orientation of the disk/torus. Recent studies showed that during the protostellar phase of high-mass star formation the magnetic field is oriented along the outflows and around or on the surfaces of the disk/torus. The observations of polarized maser emissions at milliarcsecond resolution can make a crucial contribution to understanding the orientation of the magnetic field and, consequently, the orientation of the disk/torus in NGC7538-IRS1. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22GHz H2O masers toward NGC7538-IRS1. The European VLBI Network and the MERLIN telescopes were used to measure the linear polarization and the Zeeman-splitting of the 6.7GHz CH3OH masers toward the same region. We detected 17 H2O masers and 49 CH3OH masers at high angular resolution. We detected linear polarization emission toward two H2O masers and toward twenty CH3OH masers. The CH3OH masers, most of which only show a core structure, seem to trace rotating and potentially infalling gas in the inner part of a torus. Significant Zeeman-splitting was measured in three CH3OH masers. [...] We also propose a new description of the structure of the NGC7538-IRS1 maser region.

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

A debated topic in star formation theory is the role of magnetic fields during the protostellar phase of high-mass stars. It is still unclear how magnetic fields influence the formation and dynamics of massive disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from polarized maser emissions, which allows us to investigate the magnetic field on small scales by using very long-baseline interferometry. The massive star-forming region W75N contains three radio continuum sources (VLA1, VLA2, and VLA3), at three different evolutionary stages, and associated masers, while a large-scale molecular bipolar outflow is also present. Very recently, polarization observations of the 6.7 GHz methanol masers at milliarsecond resolution have been able to probe the strength and structure of the magnetic field over more than 2000 AU around VLA1. The magnetic field is parallel to the outflow, suggesting that VLA1 is its powering source. The observations of water masers at 22 GHz can give more information about the gas dynamics and the magnetic fields around VLA1 and VLA2. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz water masers in the star-forming region W75N. We detected 124 water masers, 36 around VLA1 and 88 around VLA2 of W75N, which indicate two different physical environments around the two sources, where VLA1 is in a more evolved state. The linear polarization of the masers confirms the tightly ordered magnetic field around VLA1, which is aligned with the large-scale molecular outflow, and also reveals an ordered magnetic field around VLA2, which is not parallel to the outflow. [abridged]

Very luminous extragalactic water masers, the megamasers, are associated with active galactic nuclei (AGN) in galaxies characterized by accretion disks, radio jets, and nuclear outflows. Weaker masers, the kilomasers, seem to be mostly related to star formation activity, although the possibility exists that some of these sources may belong to the weak tail of the AGN maser distribution. It is of particular importance to accurately locate the water maser emission to reveal its origin and shed light onto extragalactic star forming activity or to elucidate the highly obscured central regions of galaxies. We performed interferometric observations of three galaxies, NGC3556, Arp299, and NGC4151, where water emission was found. Statistical tools have been used to study the relation between OH and water maser emission in galaxies. The maser in NGC3556 is associated with a compact radio continuum source that is most likely a supernova remnant or radio supernova. In Arp299, the luminous water maser has been decomposed in three main emitting regions associated with the nuclear regions of the two main galaxies of the system, NGC3690 and IC694, and the region of overlap. In NGC4151, only one of the two previously observed maser components has been tentatively detected. This feature, if real, is associated with the galaxy's central region. The only galaxy, so far, where luminous maser emission from two maser species, OH and H2O has been confidently detected is Arp299. Weaker masers from these two species do instead coexist in a number of objects. A larger number of objects searched for both maser species are, however, necessary to better assess these last two results.

We present polarization observations of the 6.7-GHz methanol masers around the massive protostar Cepheus A HW2 and its associated disc. The data were taken with the Multi-Element Radio Linked Interferometer Network. The maser polarization is used to determine the full three-dimensional magnetic field structure around Cepheus A HW2. The observations suggest that the masers probe the large scale magnetic field and not isolated pockets of a compressed field. We find that the magnetic field is predominantly aligned along the protostellar outflow and perpendicular to the molecular and dust disc. From the three-dimensional magnetic field orientation and measurements of the magnetic field strength along the line of sight, we are able to determine that the high density material, in which the masers occurs, is threaded by a large scale magnetic field of ~23 mG. This indicates that the protostellar environment at ~1000 AU from Cepheus A HW2 is slightly supercritical (lambda approximately 1.7) and the relation between density and magnetic field is consistent with collapse along the magnetic field lines. Thus, the observations indicate that the magnetic field likely regulates accretion onto the disc. The magnetic field dominates the turbulent energies by approximately a factor of three and is sufficiently strong to be the crucial component stabilizing the massive accretion disc and sustaining the high accretion rates needed during massive star-formation.

The role of magnetic fields during the protostellar phase of high-mass star-formation is a debated topic. In particular, it is still unclear how magnetic fields influence the formation and dynamic of disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from H2O and OH maser observations. Recently, the first 6.7 GHz methanol maser polarization observations were made, and they reveal strong and ordered magnetic fields. The morphology of the magnetic field during high-mass star-formation needs to be investigated on small scales, which can only be done using very long baseline interferometry observations. The massive star-forming regionW75N contains three radio sources and associated masers, while a large-scale molecular bipolar outflow is also present. Polarization observations of the 6.7 GHz methanol masers at high angular resolution probe the strength and structure of the magnetic field and determine its relation to the outflow. Eight of the European VLBI network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7 GHz methanol masers in the star-forming region W75N. We detected 10 methanol maser features, 4 of which were undetected in previous work. All arise near the source VLA1 of W75N. The linear polarization of the masers reveals a tightly ordered magnetic field over more than 2000 AU around VLA1 that is exactly aligned with the large-scale molecular outflow. This is consistent with the twisted magnetic field model proposed for explaining dust polarization observations. The Zeeman-splitting measured on 3 of the maser features indicates a dynamically important magnetic field in the maser region of the order of 50mG. We suggest VLA1 is the powering sources of the bipolar outflow.

Recently, a relationship between the water maser detection rate and far infrared (FIR) flux densities has been established as a result of two 22 GHz maser surveys in a complete sample of galaxies (Dec>-30 degree) with 100 micron flux densities of > 50 Jy and > 30 Jy. This survey has been extended to the southern galaxies in order to discover new maser sources and to investigate the galaxies hosting the maser spots with particular emphasis on their nuclear regions. A sample of 12 galaxies with Dec<-30 degree and S(100 micron)>50 Jy was observed with the 70-m telescope of the Canberra Deep Space Communication Complex (CDSCC) at Tidbinbilla (Australia) in a search for water maser emission. The average 3 sigma noise level of the survey is 15 mJy for a 0.42 km/s channel, corresponding to a detection threshold of ~0.1 solar luminosities for the isotropic maser luminosity at a distance of 25 Mpc. Two new detections are reported: a kilomaser with an isotropic luminosity L_H2O ~5 solar luminosities in NGC3620 and a maser with about twice this luminosity in the merger system NGC3256. The detections have been followed-up through continuum and spectral line interferometric observations with the Australia Telescope Compact Array (ATCA). In NGC3256, a fraction (about a third) of the maser emission arises from two hot spots associated with star formation activity, which are offset from the galactic nuclei of the system. The remaining emission may arise from weaker centers of maser activity distributed over the central 50 arcsec. [abridged]